EXTERNAL  PLTPO,VISSP
	DIMENSION EF(181,71),O(3,3)

C	RX=RZ=RY=0	DEM000
	DATA	O/0.408,-0.707,0.577,
     *		0.816,0.,-0.577,
     *		0.408,0.707,0.577/

C RX=90,RY=RZ=0		DEMX00
C	DATA	O/0.408,-0.577,-0.707,
C     *		0.816,0.577,0.,
C     *		0.408,-0.577,0.707/

C RY=90,RX=RZ=0		DEMY00
C	DATA	O/0.577,-0.707,-0.408,
C     *		-0.577,0.,-0.816,
C     *		0.577,0.707,-0.408/

C RZ=90,RX=RY=0		DEMZ00
C	DATA	O/0.707,0.408,0.577,
C     *		0.,0.816,-0.577,
C     *		-0.707,0.408,0.577/

C RX=RZ=90,RY=0		DEMXZ0
C	DATA	O/0.577,0.408,-0.707,
C     *		-0.577,0.816,0.,
C     *		0.577,0.408,0.707/

C RZ=RX=90,RY=0		DEMZX0
C	DATA	O/0.707,-0.577,0.408,
C     *		0.,0.577,0.816,
C     *		-0.707,-0.577,0.408/

	AS(J,I)=0.1*FLOAT(MOD(I+J,2))
	NT=71
	NP=181
	DT=3.14159/FLOAT(NT-1)
	DP=6.28318/FLOAT(NP)
	T=0.0
	DO 20 I=1,NT
	P=0.0
	DO 10 J=1,NP
	R1=0.9+SIN(3.0*T+P)
	R2=1.0+SIN(5.0*T-P)
	E1=0.75*SIN(T)*EXP(-R1*R1)
	E2=0.25*SIN(T)*EXP(-R2*R2)
	EF(J,I)=0.33+E1-E2+0.25*AS(J,I)
 10	P=P+DP
 20	T=T+DT

	DO 30 I=1,NT,10
	DO 30 J=1,NP,10
 30	EF(J,I)=1.0

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

	STOP
	END