source: LMDZ5/branches/LF-private/libf/phylmd/sisvat/PHY_SV.F90

Last change on this file was 1865, checked in by Laurent Fairhead, 11 years ago

Inclusion de la bibliothèque SISVAT/MAR à LMDZ pour le traitement des surfaces
"land ice"

  1. Menegoz

Integration of the SISVAT/MAR library to LMDZ to model the land ice surfaces

  1. Menegoz
File size: 4.7 KB
Line 
1MODULE PHY_SV
2
3IMPLICIT NONE
4
5! Constants
6                                                                             
7      CHARACTER(LEN=1), PARAMETER :: labnum(0:9) =                            &
8              (/'0','1','2','3','4','5','6','7','8','9'/)             
9!C +...                labnum: Alphanumeric Character                         
10                                                                               
11      INTEGER, PARAMETER ::  izr=0,iun=1                                             
12      REAL, PARAMETER    ::  zer0 = 0.0e+0, half = 0.5e+0, un_1 = 1.0e+0,     &
13     &                       eps6 = 1.0e-6, epsn = 1.0e-9, R_1000=1.e3             
14      REAL               ::  ea_MAX,ea_MIN
15!C +...                    zer0 = 0.0e+0                                       
16!C +                       half = 0.5e+0                                 
17!C +                       un_1 = 1.0e+0                             
18!C +                       eps6 = 1.0e-6  ! Arbirary Small Value           !    1.e-6     
19!C +                       epsn = 1.0e-9  ! Arbirary Small Value           !    1.e-6
20!C +                       R_1000=1.e3     
21!C +                       ea_MIN   ! MAX allowed exponential Argum. ! computed by HOST     
22!C +                       ea_MIN   ! MIN allowed exponential Argum. ! computed by HOST     
23                                                                             
24      REAL, PARAMETER    ::  piNmbr = 3.141592653589793238462643e0
25      REAL, SAVE         ::  Dg2Rad = piNmbr / 180.
26!C +...                    piNmbr   ! pi    = 3.141592653589793238462643e0             
27!C +...                    Dg2rad   ! pi / 180                                       
28      REAL, PARAMETER    ::  Grav_F = 9.81e0, Grav_I = 1./Grav_F                       
29!C +                   Grav_F:  Gravitational  Force         =  9.81    m/s2           
30!C +                   Grav_I: 1 / Grav_F                    = 1 /9.81    s2/m   
31                                                                             
32      REAL, PARAMETER    ::  A_MolV = 1.35e-5, vonKrm = 0.40e0, A_Stab = 5.8,  &
33     &                       AhStab = 5.4,     AsStab = 4.0,    r_Stab = 3.0           
34!C +...                A_MolV: Air Viscosity                 = 1.35d-5 m2/s   
35!C +                   vonKrm: von Karman constant           = 0.4           
36!C +                   A_Stab: Stability  Coefficient Moment = 5.8             
37!C +                   AhStab: Stability  Coefficient Heat   = 5.4             
38!C +                   AsStab: Stability  Coefficient Blown * = 4.0           
39!C +                   r_Stab: Turbulent Diffusivities Ratio K*/Km             
40                                                                             
41      REAL, PARAMETER    ::  R_DAir=287.05967, CpdAir=1004.708845,             &
42     &                       RCp=R_DAir/CpdAir,p0_kap=3.73010e0,               &
43     &                       StefBo=5.67e-8         
44      REAL, PARAMETER    ::  LhfH2O=3.34e+5, LhvH2O=2.5008e+6,&
45     &                       LhsH2O=2.8345e+6     
46!C +...                R_DAir: perfect gas law constant for dry air (287 J/kg/K)
47!C +                   CpdAir: dry air specific heat at constant p (1004 J/kg/K)
48!C +                   RCp   : R / Cp                     
49!C +                   p0_kap: 100 kPa ** (R/Cp)                                   
50!C +                   LhfH2O: Latent Heat of Fusion of Snow  =    3.34d+5 J/kg
51!C +                   LhvH2O: Latent Heat Vaporisation Snow  = 2500.00d+3 J/kg 
52!C +                   LhsH2O: Latent Heat Sublimation  Snow  = 2833.60d+3 J/kg 
53!C +                   StefBo: Stefan-Boltzman Constant       =    5.67d-8 W/m2/K4
54                                                                               
55      REAL, PARAMETER    ::  rhoWat=1000.00e0, hC_Wat=4186.00e0       
56!C +                   rhoWat: Water Specific Mass            = 1000.00d+0 kg/m3
57!C +                   hC_Wat: Water Heat Capacity            = 4186.00d+0 J/kg/K
58                                                                               
59      REAL, PARAMETER    ::  rhoIce=920.e0, BSnoRo=3.30e+2                             
60!C +...                rhoIce:  Density      of Pure Ice      =  920.00d+0 kg/m3
61!C +                   BSnoRo : Blowed Snow Density           =  255.00d+0 kg/m3
62                                                                               
63      REAL, PARAMETER    ::  Tf_Sno=273.16e+0,Tf_Sea=271.2e+0             
64!C +...                Tf_Sno: Snow  Melting  Point           =  273.16    K     
65!C +...                Tf_Sea: Sea   Melting  Point           =  271.2     K
66
67                     
68                                                                               
69
70END MODULE PHY_SV
Note: See TracBrowser for help on using the repository browser.