source: LMDZ6/branches/Ocean_skin/libf/phylmd/inlandsis/VARphy.F90 @ 4096

Last change on this file since 4096 was 4013, checked in by lguez, 3 years ago

Sync latest trunk changes to Ocean_skin

File size: 5.3 KB
Line 
1MODULE VARphy
2
3
4IMPLICIT NONE
5
6
7
8! Etienne, 08112020: merging YOMCST_SISVAT,  PHY_SV and VARphy
9! but redundancies remain
10!
11!  ATTENTION!!!!: ce fichier include est compatible format fixe/format libre
12!                 veillez à n'utiliser que des ! pour les commentaires
13!                 et à bien positionner les & des lignes de continuation
14!                 (les placer en colonne 6 et en colonne 73)
15!
16!
17! A1.0 Fundamental constants
18!----------------------------
19! Constants
20                                                                             
21      CHARACTER(LEN=1), PARAMETER :: labnum(0:9) =                            &
22              (/'0','1','2','3','4','5','6','7','8','9'/)             
23!C +...                labnum: Alphanumeric Character     
24     
25                                                                               
26      INTEGER, PARAMETER ::  iun=1                                             
27      REAL, PARAMETER    ::  zer0 = 0.0e+0, half = 0.5e+0, un_1 = 1.0e+0,     &
28     &                       eps6 = 1.0e-6, R_1000=1.e3   
29      REAL, PARAMETER    ::  zero = 0.0e+0, demi = 0.5e+0, unun = 1.0e+0,     &
30     &                       epsi = 1.0e-6, eps9 = 1.0e-9         
31      REAL               ::  ea_MAX,ea_MIN                                                                               
32      REAL, PARAMETER    ::  pi = 3.141592653589793238462643e0
33     
34! A1.2 Geoide
35!-------------
36      REAL, PARAMETER    ::  Grav_F = 9.81e0, Grav_I = 1./Grav_F                       
37!C +                   Grav_F:  Gravitational  Force         =  9.81    m/s2           
38!C +                   Grav_I: 1 / Grav_F                    = 1 /9.81    s2/m   
39
40                                                                               
41      REAL, PARAMETER    ::  gravit = 9.81e0,      gravi2 = gravit**2 ,       &
42     &                       grvinv = 1./gravit                       
43!C +                   gravit: Earth Gravity Acceleration    = 9.81    m/s2   
44!C +                   gravi2: idem (squared)                                 
45!C +                   grvinv: idem (inverse)                 
46! A1.3 Radiation
47!----------------
48      REAL,PARAMETER :: StefBo=5.67e-8 
49!C +                   StefBo: Stefan-Boltzman Constant       =    5.67d-8 W/m2/K4
50
51! A1.4 Thermodynamic
52!-----------------------------
53     
54      REAL, PARAMETER    ::  cp=1004.708845
55                                                       
56      REAL, PARAMETER    ::  RVapor=461.e0, Lv_H2O=2.5008e+6,&
57     &                       Ls_H2O=2.8345e+6, r_LvCp=2490.04,r_LcCp=332.27,  &
58     &                       r_LsCp=2822.31           
59!C +                   cp    : dry air specific heat at constant p (1004 J/kg/K)
60
61      REAL, PARAMETER    ::  LhfH2O=3.34e+5, LhvH2O=2.5008e+6,&
62     &                       LhsH2O=2.8345e+6     
63      REAL, PARAMETER    ::  rhoWat=1000.00e0   
64!C +                   rhoWat: Water Specific Mass            = 1000.00d+0 kg/m3
65
66      REAL, PARAMETER    ::  rhoIce=920.e0, BSnoRo=3.30e+2, frsno=300.0                             
67!C +...                rhoIce:  Density      of Pure Ice      =  920.00d+0 kg/m3
68!C +                   BSnoRo : Blowed Snow Density           =  255.00d+0 kg/m3
69                                                                               
70      REAL, PARAMETER    ::  Tf_Sno=273.16e+0
71!C +...                Tf_Sno: Snow  Melting  Point           =  273.16    K     
72
73
74      REAL, PARAMETER    ::  ro_Wat=1000.00e0, C__Wat=4186.00e0
75!C +                   ro_Wat: Density       of Water         = 1000.00d+0 kg/m3
76!C +                   C__Wat: Heat Capacity of Water         = 4186.00d+0 J/kg/K
77                                                                               
78      REAL, PARAMETER    ::  ro_Ice=920.e0, cdice=2.10e0                                 
79!C +...                ro_Ice:  Density      of Pure Ice       =  920.00d+0 kg/m3
80!C +                   cdice :  Conductivity of Pure Ice       =    2.51d+0 W/m/K
81                                                                               
82      REAL, PARAMETER    ::  TfSnow=273.15e+0,csnow=2105.00e+0,r0sno=3.00e+1, &
83     &                       blsno=3.30e+2,   Lf_H2O=3.337e+5                         
84!C +...                TfSnow:        Snow melting Temperature=  273.15d+0 K     
85!C +                    csnow:Heat Capacity of Snow             2105      J/kg/K
86!C +                         (Loth et al. 1993, JGR 98 D6, 2.2.2 2e para p.10453)
87!C +                   r0sno : Fresh  Snow Density                50.00d+0 kg/m3
88!C +                   blsno : Blowed Snow Density               330.00d+0 kg/m3
89!C +                   Lf_H2O: Latent Heat of Fusion of Snow  =    3.34d+5 J/kg 
90
91! A1.6 Turbulent and molecular diffusion
92!----------------------------------------
93      REAL, PARAMETER    ::  A_MolV = 1.35e-5, vonKrm = 0.40e0, r_turb=3.0
94      REAL, PARAMETER    ::  A_turb=5.8, akmol=1.35e-5
95!C +...                A_MolV: Air Viscosity                 = 1.35d-5 m2/s   
96!C +                   vonKrm: von Karman constant           = 0.4           
97!C +                   r_turb:   Turbulent Diffusivities Ratio K*/Km       
98!C +                   A_turb:    Stability  Coefficient Moment                                 
99!C +                   Air Viscosity                 = 1.35d-5 m2/s                 
100
101
102
103END MODULE VARphy
Note: See TracBrowser for help on using the repository browser.