[2089] | 1 | module Mod_PHY_S0_dat |
---|
| 2 | |
---|
| 3 | !--------------------------------------------------------------------------+ |
---|
| 4 | ! Sat 27-Apr-2013 MAR | |
---|
| 5 | ! module Mod_PHY_S0_dat contains the main data of | |
---|
| 6 | ! Insolation Scheme | |
---|
| 7 | ! | |
---|
| 8 | ! version 3.p.4.1 created by H. Gallee, Thu 25-Apr-2013 | |
---|
| 9 | ! Last Modification by H. Gallee, Sat 27-Apr-2013 | |
---|
| 10 | ! | |
---|
| 11 | !--------------------------------------------------------------------------+ |
---|
| 12 | |
---|
| 13 | |
---|
| 14 | use Mod_Real |
---|
| 15 | |
---|
| 16 | |
---|
| 17 | IMPLICIT NONE |
---|
| 18 | |
---|
| 19 | |
---|
| 20 | |
---|
| 21 | integer, SAVE :: kBP = 0 ! Time [kyr B.P.] |
---|
| 22 | real(kind=real8), SAVE :: rsunS0 ! Insolation, Top of the Atmosphere [W/m2] |
---|
| 23 | real(kind=real8), SAVE :: cszEPS = 1.e-3 ! Minimum accepted cos(Sun zenith.Dist.) [-] |
---|
| 24 | real(kind=real8), SAVE :: ecc ! Earth Orbit Eccentricity [-] |
---|
| 25 | real(kind=real8), SAVE :: perh ! Earth Orbit Longitude of the Perihelion [degree] |
---|
| 26 | real(kind=real8), SAVE :: xob ! Earth Orbit Obliquity [degree] |
---|
| 27 | real(kind=real8), SAVE :: ecc_EO(-10:0) ! Earth Orbit Eccentricity [-] |
---|
| 28 | real(kind=real8), SAVE :: perhEO(-10:0) ! Earth Orbit Longitude of the Perihelion [degree] |
---|
| 29 | real(kind=real8), SAVE :: xob_EO(-10:0) ! Earth Orbit Obliquity [degree] |
---|
| 30 | |
---|
| 31 | data ecc_EO( 0) / 0.01673 / ! Eccentricity |
---|
| 32 | data perhEO( 0) /102.4 / ! Longitude of the Perihelion [degree] |
---|
| 33 | data xob_EO( 0) / 23.445 / ! Obliquity [degree] |
---|
| 34 | |
---|
| 35 | data ecc_EO( -6) / 0.018682/ ! Eccentricity |
---|
| 36 | data perhEO( -6) / 0.87 / ! Longitude of the Perihelion [degree] |
---|
| 37 | data xob_EO( -6) / 24.105 / ! Obliquity [degree] |
---|
| 38 | |
---|
| 39 | data ecc_EO(-10) / 0.019419/ ! Eccentricity |
---|
| 40 | data perhEO(-10) /294.81 / ! Longitude of the Perihelion [degree] |
---|
| 41 | data xob_EO(-10) / 24.226 / ! Obliquity [degree] |
---|
| 42 | |
---|
| 43 | |
---|
| 44 | real(kind=real8), SAVE :: pirr ! 180 / pi / 3600 [degree/radian/sec] |
---|
| 45 | |
---|
| 46 | real(kind=real8), SAVE :: om = 0.0172142 ! Radian advance in one day [radian] |
---|
| 47 | real(kind=real8), SAVE :: Tyear = 365.25 ! Length of one year [day] |
---|
| 48 | real(kind=real8), SAVE :: Step ! Advance on Earth Orbit in 1 day [degree/day] |
---|
| 49 | |
---|
| 50 | real(kind=real8), SAVE :: xee ! Square of Eccentricity (Ecc) [-] |
---|
| 51 | real(kind=real8), SAVE :: xse ! Square Root (1-Ecc**2) [-] |
---|
| 52 | real(kind=real8), SAVE :: xe3 ! Square Root (1-Ecc**2) X Ecc [-] |
---|
| 53 | real(kind=real8), SAVE :: xl ! Longitude of Aphelion [degree] |
---|
| 54 | real(kind=real8), SAVE :: xllp ! Longitude of Aphelion [radian] |
---|
| 55 | real(kind=real8), SAVE :: so ! sinus of Obliquity [-] |
---|
| 56 | real(kind=real8), SAVE :: xlam ! true long. sun for mean long. = 0 [degree] |
---|
| 57 | |
---|
| 58 | |
---|
| 59 | |
---|
| 60 | end module Mod_PHY_S0_dat |
---|