EXTERNAL  PLTPO,VISSP
	DIMENSION EF(181,71),O(3,3)
	DATA	O/0.408,-0.707,0.577,0.816,0.,-0.577,0.408,0.707,0.577/

	AS(J,I)=0.1*FLOAT(MOD(I+J,2))
	NT=71
	NP=181
	S=1.0
	DT=3.14159/FLOAT(NT-1)
	DP=6.28318/FLOAT(NP)
	T=3.14159
	DO 20 I=1,NT
	P=0.35
	DO 10 J=1,NP
	TF=0.67*(1.0+0.5*COS(3.0*T)+S*0.1667*COS(9.0*T))
	PF=0.67*(1.0+0.5*SIN(3.0*P)-S*0.1667*SIN(9.0*P))
	EF(J,I)=TF*PF
 10	P=P+DP
 20	T=T-DT

	CALL PFNAME ('DEM19.BLD')
	CALL PLT00
	CALL PLTFR
	CALL PLTLA('DEM19')
	CALL PLTSV(EF,NP,NT,O,VISSP,PLTPO)
	CALL PLTEND

	STOP
	END