Insolation Examples

Diurnal cycle of insolation

using Dates

include("plot_diurnal_cycle.jl")

# Pasadena in January
lat, lon = [34.15, -118.14]
date = DateTime(2020, 01, 10)
timezone = +8 # Pacific Standard Time
od = Insolation.OrbitalData()
diurnal_cycle(lat, lon, date, od, timezone, "Pasadena_January.png")

# Finland in June
lat, lon = [66.50, 25.73]
date = DateTime(2020, 06, 10)
timezone = -3 # Eastern European Summer Time
diurnal_cycle(lat, lon, date, od, timezone, "Finland_June.png")

"/home/runner/work/Insolation.jl/Insolation.jl/docs/build/Finland_June.png"

Daily-mean insolation

Insolation in J2000

import Insolation
import Insolation.Parameters as IP
import ClimaParams as CP

FT = Float64
param_set = IP.InsolationParameters(FT)

include("plot_insolation.jl")

γ0 = IP.obliq_epoch(param_set)
ϖ0 = IP.lon_perihelion_epoch(param_set)
e0 = IP.eccentricity_epoch(param_set)
od = Insolation.OrbitalData()

days, lats, F0 = calc_day_lat_insolation(od, 365, 180, param_set)
plot_day_lat_insolation(days, lats, F0, "YlOrRd", title, "insol_example1.png")

"/home/runner/work/Insolation.jl/Insolation.jl/docs/build/insol_example1.png"

Insolation with smaller obliquity

import Insolation
import Insolation.Parameters as IP
import ClimaParams as CP

FT = Float64
param_set = IP.InsolationParameters(FT)

od = Insolation.OrbitalData()

# decrease γ to 20.0°
param_set = IP.InsolationParameters(FT, (; obliq_epoch = deg2rad(20.0)))
γ1 = IP.obliq_epoch(param_set)
days, lats, F2 = calc_day_lat_insolation(od, 365, 180, param_set)

plot_day_lat_insolation(days,lats,F2,"YlOrRd",  title, "insol_example2a.png")
plot_day_lat_insolation(days, lats, F2-F0, "PRGn", title, "insol_example2b.png")

"/home/runner/work/Insolation.jl/Insolation.jl/docs/build/insol_example2b.png"

Insolation with very large obliquity (like Uranus)

import Insolation
import Insolation.Parameters as IP
import ClimaParams as CP

FT = Float64
param_set = IP.InsolationParameters(FT)

od = Insolation.OrbitalData()

# now change obliquity to 97.86°
param_set = IP.InsolationParameters(FT, (;obliq_epoch = deg2rad(97.86)))
γ4 = IP.obliq_epoch(param_set)
days, lats, F5 = calc_day_lat_insolation(od, 365, 180, param_set)

plot_day_lat_insolation(days,lats,F5,"YlOrRd", title, "insol_example3a.png")
plot_day_lat_insolation(days, lats, F5-F0, "PRGn", title, "insol_example3b.png")

"/home/runner/work/Insolation.jl/Insolation.jl/docs/build/insol_example3b.png"