Solar Insolation Computer Codes for Windows

There are two FORTRAN codes for calculating Solar Insolation:

Berger (http://www.astr.ucl.ac.be/research/astro.html) code Insol.f: This solution is suitable to compute for a 1,000,000 period centered on the present day. The code was meant to be compiled on a Unix machine and calculates only for one year. I revised the code to calculate a range of years and compiled it to run on a Windows machine (InsoMY.f=Inso for Many Years). The files in the zip file InsolMY.zip are: InsoMY.f, InsoMY.exe, Insol.in (contains the input parameters), Inso.out, InsoMYSession.txt, READMEBerger.txt (this file). I use the solar constant 1367 watts/m^2. (See http://www.ras.org.uk/pdfs/Solanki.pdf.)
I also have developed an Excel workbook for doing the Berger calculation. It is 7.5 Mbytes zipped and 17.5 Mbytes unzipped. If you want a copy of it, let me know (roperld@vt.edu) and I will put a copy of it on a disk and send it to you.

Laskar (http://xml.gsfc.nasa.gov/archive/catalogs/6/6063/index_long.html) code Insola.f: This solution is suitable to compute for -20,000,000 in the past to 10,000,000 in the future. The code was meant to be compiled on a Unix machine and required that another code, prepinsol.f, be compiled and run first. I revised the code and compiled it to run on a Windows machine (InsolaSA.f=Insola Stand Alone) without the need for running the prepinsol.f code first. The files in the zip file InsolaSA.zip are: InsolaSA.f, InsolSA.exe, InsolaSA.par, InsolaSA.out, InsolaSession.txt. I also created a version that writes the output file with only the insolation values in a column, to make it easy to import the data into a spreadsheet. Those files are also in the zip file: InsolaWOnly.f and InsolaWOnly.exe. The input data files must be downloaded from ftp://ftp.imcce.fr/pub/ephem/sun/la93: They are ORBELN.ASC.Z, ORBELP.ASC.Z, PREC0N.ASC.Z and PREC0P.ASC.Z, which must be unzipped and placed in the same directory as the InsolaSA.exe and the InsolaSA.par files. The parameter file, InsolaSA.par, has the solar constant 1367 watts/m^2, which can be changed by an editor. (See http://www.ras.org.uk/pdfs/Solanki.pdf.)