PROGRAM TEST

	real day,solstice,month,date,tilt,pi,lat
        

	pi=3.141592654
	day=1.
	lat=38.
	do 10, month=1.,12.,1.
           IF(MONTH.EQ.1.)DATE=DAY
	   IF(MONTH.EQ.2.)DATE=DAY+31.
	   IF(MONTH.EQ.3.)DATE=DAY+59.
	   IF(MONTH.EQ.4.)DATE=DAY+90.
	   IF(MONTH.EQ.5.)DATE=DAY+120.
	   IF(MONTH.EQ.6.)DATE=DAY+151.
	   IF(MONTH.EQ.7.)DATE=DAY+181.
	   IF(MONTH.EQ.8.)DATE=DAY+212.
	   IF(MONTH.EQ.9.)DATE=DAY+243.
	   IF(MONTH.EQ.10.)DATE=DAY+273.
	   IF(MONTH.EQ.11.)DATE=DAY+304.
	   IF(MONTH.EQ.12.)DATE=DAY+334.
	   SOLSTICE=172.
	   TILT=23.45
	   POLE=TILT*COS((DATE-SOLSTICE)*2.*PI/365.)
 10	   write(6,*)MONTH,date,tilt,pole,COS((LAT-POLE)*PI/180.)
       end