source: LMDZ5/branches/testing/libf/phymar/Mod_PHY____dat.f90 @ 3946

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

Merged trunk changes -r2070:2158 into testing branch. Compilation problems introduced by revision r2155 have been corrected by hand

File size: 6.7 KB
Line 
1     module Mod_PHY____dat
2
3!------------------------------------------------------------------------------+
4!                                                         Sat  8-Jun-2013  MAR |
5!    module Mod_PHY____dat  contains specific constants for                    |
6!           MAR PHYsical    Parameterizations                                  |
7!                                                                              |
8!                     These constants may be slightly different in HOST Model  |
9!                                 and may be modified in  PHY________INI       |
10!                                                                              |
11!     version 3.p.4.1 created by H. Gallee,               Tue 12-Mar-2013      |
12!           Last Modification by H. Gallee,               Sat  8-Jun-2013      |
13!                                                                              |
14!------------------------------------------------------------------------------+
15
16
17
18! General Variables
19! =================
20
21      use Mod_Real
22
23
24      IMPLICIT NONE
25
26
27! Physical Constants (may be changed in the interface with the Host Model)
28! ==================
29
30      real(kind=real8), SAVE  ::  zer0   =    0.00         !                                !    0.00
31      real(kind=real8), SAVE  ::  half   =    0.50         !                                !    0.50
32      real(kind=real8), SAVE  ::  un_1   =    1.00         !                                !    1.00
33      real(kind=real8), SAVE  ::  piNmbr =    3.1416       ! pi                             !    3.1416
34      real(kind=real8), SAVE  ::  Dg2Rad = 1745.33d-5      ! pi / 180                       !           
35      real(kind=real8), SAVE  ::  epsq   =    1.e-6        ! Arbirary Small Value, H2O      !    1.e-6
36      real(kind=real8), SAVE  ::  eps1   =    1.e-1        ! Arbirary Small Value           !    1.e-1
37      real(kind=real8), SAVE  ::  eps6   =    1.e-6        ! Arbirary Small Value           !    1.e-6
38      real(kind=real8), SAVE  ::  epsn   =    1.e-9        ! Arbirary Small Value           !    1.e-9
39      real(kind=real8), SAVE  ::  epsp   =    1.e-12       ! Arbirary Small Value           !    1.e-12
40      real(kind=real8), SAVE  ::  R_1by3 =    0.333333     ! 1 / 3                          !
41      real(kind=real8), SAVE  ::  R_5by3 =    1.666666     ! 5 / 3                          !
42      real(kind=real8), SAVE  ::  R_1000 =    1.e+3        !                                !    1.e+3
43      real(kind=real8), SAVE  ::  ea_MAX =   50.           ! MAX allowed exponential Argum. ! computed by HOST
44      real(kind=real8), SAVE  ::  ea_MIN =  -50.           ! MIN allowed exponential Argum. ! computed by HOST
45      real(kind=real8), SAVE  ::  A_MolV =    1.35d-5      ! Air Molecular Viscosity        !    1.35d-5 m2/s
46      real(kind=real8), SAVE  ::  rhoIce =  920.00d+0      ! Ice         Specific Mass      !  920.00d+0 kg/m3
47      real(kind=real8), SAVE  ::  BSnoRo =  255.00d+0      ! Blown Snow  Specific Mass      !  255.00d+0 kg/m3
48      real(kind=real8), SAVE  ::  LhvH2O = 2508.00d+3      ! Latent Heat of Vapor. of Rain  ! 2500.00d+3 J/kg
49      real(kind=real8), SAVE  ::  LhfH2O =    3.34d+5      ! Latent Heat of Fusion of Snow  !    3.34d+5 J/kg
50      real(kind=real8), SAVE  ::  LhsH2O = 2833.60d+3      ! Latent Heat of Sublim.of Snow  ! 2833.60d+3 J/kg
51      real(kind=real8), SAVE  ::  CpdAir = 1004.708845     ! dry air specific heat at cst p ! 1004.00    J/kg/K
52      real(kind=real8), SAVE  ::  R_DAir =  287.05967      ! dry air perfect gas law  cst   !  287.      J/kg/K
53      real(kind=real8), SAVE  ::  RCp    =    0.285857     ! R / Cp                         ! 287./1004. -
54      real(kind=real8), SAVE  ::  p0_kap =    3.730037     !                                ! 100 kPa ** (R/Cp)
55      real(kind=real8), SAVE  ::  Lv_CPd                   ! LhvH2O      /  CpdAir          !
56      real(kind=real8), SAVE  ::  Ls_CPd                   ! LhsH2O      /  CpdAir          !
57      real(kind=real8), SAVE  ::  Lc_CPd                   ! LhfH2O      /  CpdAir          !
58      real(kind=real8), SAVE  ::  hC_Wat = 4186.00d+0      ! Water Heat Capacity            ! 4186.00d+0 J/kg/K
59      real(kind=real8), SAVE  ::  rhoWat = 1000.00d+0      ! Water Specific Mass            ! 1000.00d+0 kg/m3
60      real(kind=real8), SAVE  ::  Tf_Sno =  273.16         ! Snow  Melting  Point           !  273.16    K
61      real(kind=real8), SAVE  ::  Tf_Sea =  271.2          ! Sea   Melting  Point           !  271.2     K
62      real(kind=real8), SAVE  ::  StefBo =    5.67d-8      ! Stefan-Boltzman Constant       !    5.67d-8 W/m2/K4
63      real(kind=real8), SAVE  ::  Grav_F =    9.81         !     Gravitational  Force       !    9.81    m/s2
64!     real(kind=real8), SAVE  ::  Grav_I =    0.101937     ! 1 /(Gravitational  Force)      ! 1 /9.81    s2/m
65      real(kind=real8), SAVE  ::  Grav_I                   ! 1 /(Gravitational  Force)      ! 1 /9.81    s2/m
66      real(kind=real8), SAVE  ::  GravF2                   !    (Gravitational  Force) ** 2 !    9.81    m2/s4
67      real(kind=real8), SAVE  ::  vonKrm =    0.4          ! von Karman Constant            !    0.4
68      real(kind=real8), SAVE  ::  A_Stab =    5.8          ! Stability  Coefficient Moment  !    5.8
69      real(kind=real8), SAVE  ::  AhStab =    5.4          ! Stability  Coefficient Heat    !    5.4
70      real(kind=real8), SAVE  ::  AsStab =    4.0          ! Stability  Coefficient Blown * !    4.0
71      real(kind=real8), SAVE  ::  r_Stab =    3.0          ! Turbul.Diffusivit.Ratio K*/Km  !       
72
73      real(kind=real8), SAVE  ::  EarthR = 6371.229e3      ! Earth   Radius                 !            m
74      real(kind=real8), SAVE  ::  DirAxX                   ! x-Axis  Direction              !   90       degrees is the most natural choice
75      real(kind=real8), SAVE  ::  sh_MAX                   ! Highest Domain  Grid  Point    !       
76      real(kind=real8), SAVE  ::  dzaMIN                   ! Thinest Atmosph.Layer Thickness!       
77
78      character(len=3), dimension(0:12)  ::  LabMon  !
79      data LabMon /'---','Jan','Feb','Mar','Apr','May','Jun'           &
80     &                  ,'Jul','Aug','Sep','Oct','Nov','Dec'/
81
82      integer, SAVE         , dimension(0:12)  ::  njYear  ! Nb of Days since Begin of Year !
83      data njYear / 0   ,   0 ,  31 ,  59 ,  90 , 120 , 151            &
84     &                  , 181 , 212 , 243 , 273 , 304 , 334/
85
86      integer, SAVE         , dimension(0:12)  ::  njLeap  ! Nb of added Days for Leap Year !
87      data njLeap / 0   ,   0 ,   0 ,   1 ,   1 ,   1 ,   1            &
88     &                  ,   1 ,   1 ,   1 ,   1 ,   1 ,   1/
89
90
91!     CAUTION: values in the 3rd column are purely indicative
92
93      end module Mod_PHY____dat
Note: See TracBrowser for help on using the repository browser.