| 1 | MODULE VARphy |
|---|
| 2 | |
|---|
| 3 | |
|---|
| 4 | IMPLICIT 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 | |
|---|
| 103 | END MODULE VARphy |
|---|
