 This programme computes the total daily irradiation received at the top 
of the atmosphere for a given latitude and time in the year (in Watts/Meter^2)
 Enter successively:
     PHI : Latitude
     TIMESTART (Integer in 1000 Years, 0 for the present
       Negative for the past >-1000; Positive for future <1000)
     TIMEEND (Greater than TIMESTART)
     IOPT: = 1 True Long. Sun TLS; = 2  Month MA - Day J
90
1
10
1

 Enter true longitude of Sun (TLS)
TLS=  0, 30, 60, 90,120,150,180,210,240,270,300,330
for Mar,Apr,May,Jun,Jul,Aug,Sep,Oct,Nov,Dec,Jan,Feb
90
 Time =    1000.0   Eccentricity = 0.01631   Precession Param. = 0.01423
   Longitude Perihelion =  119.2   Obliquity = 23.317
Latitude =  90.0 True Long. Sun =   90.0
 Day Insol. =   520.68 Watts/m^2     Length day = 24.00 Hours
Time(x10^3 year)   Insolation(Watts/m^2)
      1.000         520.68
      2.000         521.47
      3.000         523.25
      4.000         525.62
      5.000         528.17
      6.000         530.49
      7.000         532.23
      8.000         533.17
      9.000         533.20
     10.000         532.33
