DIMENSION EYE(3),ST1(20,20,2),IS2(4,75),IOBJ(20,20)
	CALL OPENPL(2,5,'3DTEST2.PLT',-1)
	CALL LIMITS(2.5,17.5,2.5,17.5)
	IU=1
	OPEN(UNIT=IU, ACCESS='SEQUENTIAL', FORM='UNFORMATTED',
     *	     NAME='SY:3D.TMP', TYPE='SCRATCH')
	EYE(1)=250.
	EYE(2)=150.
	EYE(3)=100.
	CALL INIT3D(0,EYE,20,20,20,ST1,4,75,IS2,IU)
	DO 120 I=1,20
	  U=(10.5-FLOAT(I))/20.
	  UU=U*U
	  DO 110 J=1,20
	    V=(FLOAT(J)-10.5)/20.
	    VV=V*V
	    A=1./SQRT(UU+VV)
	    DO 100 K=1,20
	      IF(I*J.EQ.1 .OR. I*K.EQ.1 .OR. J*K.EQ.1) GOTO 80
	      W=(FLOAT(K)-10.5)/20.
	      IF(1./SQRT(UU+VV+W*W)+(0.5-V)**2*A .LE. 4.) GOTO 90
 80	      IOBJ(J,K)=1
	      GOTO 100
 90	      IOBJ(J,K)=0
 100	  CONTINUE
 110	CONTINUE
	CALL DANDR(20,20,ST1,4,75,75,IS2,IU,IOBJ)
 120	CONTINUE
	CALL GREND
 999	STOP
	END