C	[DEM17.FOR]
C	STRAWBERRY
C	DEMONSTRATIONFOR THE PROGRAM PLTOV, WHICH CALCULATES
C	THE OUTER BOUND OF TWO FUNCTIONS DEFINED OVER A SPHERICAL
C	SURFACE. THE DEMONSTRATION SHOWS A "STRAWBERRY" SURROUNDED
C	BY A SPHERE.
C	[22-MAY-74]

	EXTERNAL  PLTPO,VISSP
	DIMENSION O1(3,3),O2(3,3)
	DIMENSION EFF(16,7),GEE(16,7)
	DIMENSION EF(151,61)
	DATA  O1/0.577,0.577,0.577,
     *       -0.707,0.000,0.707,
     *         0.408,-0.816,0.408/
	DATA  O2/0.577,-0.577,0.577,
     *         0.707,0.000,-0.707,
     *          0.408,0.816,0.408/

	AS(I,J)=0.1*FLOAT(MOD(I+J,2))
	L=2
	F=1.0
	NT=61
	NP=151
	NTT=7
	NPP=16
	DT=3.14159/FLOAT(NT-1)
	DP=6.28318/FLOAT(NP)
	T=0.0
	DO 20 I=1,NP
	P=0.5
	DO 10 J=1,NP
	TF=0.2+0.8*SIN(0.5*T)
	PF=0.2*(4.0+SIN(P)+0.1*SIN(3.0*P))
	EF(J,I)=TF*PF+AS(J,I)
   10	P=P+DP
   20	T=T+DT

	CALL PFNAME('DEM17.BLD')
	CALL PLT00
	CALL PLTFR
	CALL PLTLA ('DEM17')

	DO 30 I=1,NTT
	DO 30 J=1,NPP
	EFF(J,I)=EF(10*(J-1)+1,10*(I-1)+1)
   30	GEE(J,I)=F
	IF (L.EQ.1) CALL PLTOV (EFF,GEE,NPP,NTT,1.0,O1,PLTPO)
	IF (L.EQ.2) CALL PLTOV (EFF,GEE,NPP,NTT,1.0,O2,PLTPO)

	IF (L.EQ.1) CALL PLTSV (EF,NP,NT,O1,VISSP,PLTPO)
	IF (L.EQ.2) CALL PLTSV (EF,NP,NT,O2,VISSP,PLTPO)
	CALL PLTEND
	CALL EXIT
	END