#OBSOLETE EXAMPLE SCRIPT

from math import *
from solar_flux import *

degrad = pi/180.
obliquity = degrad*22.5
L = 1360. #Solar constant for Earth
#
# Make a file of the flux factor vs. time for
# a few different latitudes
outfile = open("flux.out","w")
for season in [2*i for i in range(180)]:
    outfile.write("%f %f %f %f\n"%(season,
                                  L*solar(0.,degrad*season,obliquity),
                                  L*solar(degrad*45,degrad*season,obliquity),
                                  L*solar(degrad*60,degrad*season,obliquity)))

outfile.close()

# Now make a file of the annual mean flux vs latitude
# (meanflux function assumes a circular orbit)

outfile = open("meanflux.out","w")
for lat in range(-90,91):
    outfile.write("%f %f\n"%(lat,L*meanflux(degrad*lat,obliquity)))
outfile.close()