source: LMDZ5/branches/LF-private/libf/phylmd/sisvat/VARphy.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: 7.6 KB
Line 
1MODULE VARphy
2
3IMPLICIT NONE
4
5! MAR 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    ::  zero = 0.0e+0, demi = 0.5e+0, unun = 1.0e+0,     &
13     &                       epsi = 1.0e-6, eps9 = 1.0e-9, eps12= 1.0e-12,    &
14     &                       third=1/3,     thous=1.e3             
15      REAL               ::  argmin,argmax      !hj valeurs??
16!C +...                    zero = 0.0e+0                                       
17!C +                            demi = 0.5e+0                                 
18!C +                                 unun = 1.0e+0                             
19!C +                                      epsi = 1.0e-6                       
20!C +                                           eps9 = 1.0e-9                   
21!C +                                                eps12= 1.0e-12             
22!C +                                                      third=1/3           
23!C +                                                            thous=1.e3     
24!C +                       argmin:        Function exp(x) Minimum Argument     
25!C +                              argmax: Function exp(x) Maximum Argument     
26!C +                                      (Min/Max are Machine Dependant)     
27                                                                             
28      REAL, PARAMETER    ::  pi = 3.141592653589793238462643e0,               &
29     &                       degrad = pi / 180.d0, hourad = pi /  12.d0
30!C +...                pi    = 3.141592653589793238462643e0                   
31!C +                   degrad= pi / 180.d0                                     
32!C +                   hourad= pi /  12.d0                                   
33                                                                               
34      REAL, PARAMETER    ::  earthr = 6371.229e+3, earthv = 729.217e-7,       &
35     &                       gravit = 9.81e0,      gravi2 = gravit**2 ,       &
36     &                       grvinv = 1./gravit                       
37!C +...                earthr: Earth Radius                  =6371.d+3    m   
38!C +...                earthv: Earth Angular Velocity        = 7.29d-5  s-1   
39!C +                   gravit: Earth Gravity Acceleration    = 9.81    m/s2   
40!C +                   gravi2: idem (squared)                                 
41!C +                   grvinv: idem (inverse)                                 
42                                                                             
43      REAL, PARAMETER    ::  akmol = 1.35e-5, vonkar = 0.40e0, A_Turb = 5.8,  &
44     &                       AhTurb = 5.4,    AsTurb = 4.0,    r_Turb = 3.0           
45!C +...                akmol : Air Viscosity                 = 1.35d-5 m2/s   
46!C +                   vonkar: von Karman constant           = 0.4           
47!C +                   A_Turb: Stability  Coefficient Moment = 5.8             
48!C +                   AhTurb: Stability  Coefficient Heat   = 5.4             
49!C +                   AsTurb: Stability  Coefficient BLOW * = 4.0           
50!C +                   r_Turb: Turbulent Diffusivities Ratio K*/Km             
51                                                                             
52!      REAL    ::  RDryAi,cp,CvDrya,racv,cap,pcap,stefan,qv_MIN           
53      REAL, PARAMETER    ::  RDryAi=287.05967, ra=287.05967, cp=1004.708845,  &
54     &                       CvDrya=Cp-RDryAi, racv=RDryAi/CvDryA,            &
55     &                       cap=0.28586e0,    pcap=3.73010e0,                &
56     &                       stefan=5.67e-8,   qv_MIN=3.00e-6         
57      REAL, PARAMETER    ::  RVapor=461.e0,    epsq=3.00e-6, Lv_H2O=2.5008e+6,&
58     &                       Ls_H2O=2.8345e+6, r_LvCp=2490.04,r_LcCp=332.27,  &
59     &                       r_LsCp=2822.31           
60!C +...                RDryAi: perfect gas law constant for dry air (287 J/kg/K)
61!C +                   ra    : perfect gas law constant for dry air (287 J/kg/K)
62!C +                   cp    : dry air specific heat at constant p (1004 J/kg/K)
63!C +                   CvDryA: dry air specific heat at constant V ( 717 J/kg/K)
64!C +                   racv  = ra / cv      = ra /(cp-ra)                       
65!C +                   cap   = ra / cp                                           
66!C +                   pcap  = 100.[kPa]**cap                                   
67!C +...                RVapor: perfect gas law constant Water Vapor (461 J/kg/K)
68!C +                   epsq  : Minimum Water Vapor Content    =    3.00d-6 kg/kg
69!C +                   Lv_H2O: Latent Heat Vaporisation/Water = 2500.00d+3 J/kg 
70!C +                   Ls_H2O: Latent Heat Sublimation /Ice   = 2833.60d+3 J/kg 
71!C +                   r_LvCp [Lv=2500000J/kg]/[Cp=1004J/kg/K]= 2490.04    K/[kg/
72!C +                   r_LvCp [Lc= 333600J/kg]/[Cp=1004J/kg/K]=  332.27    K/[kg/
73!C +                   r_LvCp [Ls=2833600J/kg]/[Cp=1004J/kg/K]= 2822.31    K/[kg/
74!C +                   stefan: Stefan Constant                =    5.67d-8 W/m2/K
75!C +                   qv_MIN: Minimum Specific Humidity      =    3.00d-6 kg/kg
76                                                                               
77      REAL, PARAMETER    ::  tfrwat=271.20e0,  siclf=302.00e6,cdsice=2.04e0,  &
78     &                       hic0=1.00e0, ro_Wat=1000.00e0, C__Wat=4186.00e0, &
79     &                       fracoh=0.75e0             
80!C +...                tfrwat: Sea-Water freezing Temperature =  271.20d+0 K     
81!C +                   siclf : Sea-Ice Heat of Fusion         =  302.00d+6 J/m3 
82!C +                   cdsice: Sea-Ice Thermal Conductivity   =    2.04d+0 W/mK 
83!C +                   hic0  : Sea-Ice Initial Thickness      =    2.00d+0 m     
84!C +                   ro_Wat: Density       of Water         = 1000.00d+0 kg/m3
85!C +                   C__Wat: Heat Capacity of Water         = 4186.00d+0 J/kg/K
86!C +                   fracoh: 25% of cooling = Oceanic Heat Flux               
87!C +***                        (Hibler 1984)    (ANTARCTIC Ocean)               
88                                                                               
89      REAL, PARAMETER    ::  ro_Ice=920.e0, cdice=2.10e0                                 
90!C +...                ro_Ice:  Density      of Pure Ice       =  920.00d+0 kg/m3
91!C +                   cdice :  Conductivity of Pure Ice       =    2.51d+0 W/m/K
92                                                                               
93      REAL, PARAMETER    ::  TfSnow=273.15e+0,csnow=2105.00e+0,r0sno=3.00e+1, &
94     &                       blsno=3.30e+2,   Lf_H2O=3.337e+5                         
95!C +...                TfSnow:        Snow melting Temperature=  273.15d+0 K     
96!C +                    csnow:Heat Capacity of Snow             2105      J/kg/K
97!C +                         (Loth et al. 1993, JGR 98 D6, 2.2.2 2e para p.10453)
98!C +                   r0sno : Fresh  Snow Density                50.00d+0 kg/m3
99!C +                   blsno : Blowed Snow Density               330.00d+0 kg/m3
100!C +                   Lf_H2O: Latent Heat of Fusion of Snow  =    3.34d+5 J/kg 
101
102! MAR Constants
103                     
104                                                                               
105
106END MODULE VARphy
Note: See TracBrowser for help on using the repository browser.