GEUS-Glaciology-and-Climate / AWS_v3

PROMICE AWS v3 processor - IDL/GDL
0 stars 0 forks source link

Fractional year equation #3

Open BaptisteVandecrux opened 3 years ago

BaptisteVandecrux commented 3 years ago

https://github.com/GEUS-PROMICE/AWS_v3/blob/dc958542137403b480601565726b50d62075b103/AWSdataprocessing_v3.pro#L495

shouldn't it be:

     d0_rad = 2.*!pi*(day_year+(hour+minute/60. -12)/24.-1.)/365.

Source: NOAA Solar Calculator, 2020, https://www.esrl.noaa.gov/gmd/grad/solcalc/solareqns.PDF

Dirk-van-As commented 3 years ago

Half day shift? Perhaps your d0 is defined differently. Perhaps compare/plot resulting declinations ("yours" vs "mine") to see which one makes more sense?