source: LMDZ5/trunk/libf/phymar/Mod_PHY_S0_dat.f90 @ 5440

Last change on this file since 5440 was 2089, checked in by Laurent Fairhead, 10 years ago

Inclusion de la physique de MAR


Integration of MAR physics

File size: 4.3 KB
RevLine 
[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
Note: See TracBrowser for help on using the repository browser.