Skip to content

Latest commit

 

History

History
33 lines (26 loc) · 1.34 KB

README.md

File metadata and controls

33 lines (26 loc) · 1.34 KB

insol

Currently there is one public function available in the insol module: calc_insol_day. This function returns the average insolation for a given day of year, latitude, and time between 5 Ma years before present and 1 Ma into the future. The function can calculate the insolation for a point, a vector of latitudes, or a 2D array of latitudes, returning a double with the same size as the input.

double precision function calc_insol_day(day,lats,time_bp,[S0],[day_year],[fldr])
day        Day of the year
lats       point, vector or 2D array of latitude values of interest 
           (-90:90 degrees)
time_bp    "time before present", ie, year relative to 1950
S0         Solar constant (optional, default 1365.0 W/m2)
day_year   Total number of days in the year, currently only
           handles day_year=360 (optional)
fldr       Folder containing table of precomputed orbital 
           parameter values (optional, default=input/)

Additional interfaces to be developed may include:

double precision function calc_insol_daily(days,lat,time_bp,[S0],[day_year],[fldr])

... to calculate daily insolation over a year for a given latitude.

double precision function calc_insol_hour(hour,day,lats,time_bp,[S0],[day_year],[fldr])

... to calculate the hourly insolation for a given day and latitude(s).

More?