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