Changeset 5087 for LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis
- Timestamp:
- Jul 20, 2024, 12:00:23 PM (4 months ago)
- Location:
- LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis
- Files:
-
- 3 edited
-
VARdSV.F90 (modified) (1 diff)
-
VARphy.F90 (modified) (4 diffs)
-
surf_inlandsis_mod.F90 (modified) (5 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/VARdSV.F90
r3792 r5087 57 57 58 58 REAL,PARAMETER :: etadSV(0:nsot) = (/ 1.000,0.395,0.410, & 59 & 0.435,0.485,0.451,0.420,0.477,0.476,0.426,0.492,0.482,0.001 /)59 0.435,0.485,0.451,0.420,0.477,0.476,0.426,0.492,0.482,0.001 /) 60 60 ! Water Content at Saturation [m3/m3] 61 61 62 62 REAL,PARAMETER :: psidSV(0:nsot) = (/ 1.000,0.121,0.090, & 63 &0.218,0.786,0.478,0.299,0.356,0.630,0.153,0.490,0.405,0.001 /)63 0.218,0.786,0.478,0.299,0.356,0.630,0.153,0.490,0.405,0.001 /) 64 64 ! Water Succion at Saturation [m] 65 65 66 66 REAL,PARAMETER :: Ks_dSV(0:nsot) = (/ 0.e00, 176.0e-6, & 67 &156.3e-6, 34.1e-6, 7.2e-6, 7.0e-6, 6.3e-6, &68 &1.7e-6, 2.5e-6, 2.2e-6, 1.0e-6, 1.3e-6,0.0e0 /)67 156.3e-6, 34.1e-6, 7.2e-6, 7.0e-6, 6.3e-6, & 68 1.7e-6, 2.5e-6, 2.2e-6, 1.0e-6, 1.3e-6,0.0e0 /) 69 69 ! Hydraulic Conductivity 70 70 ! at Saturation [m/s] 71 71 REAL,PARAMETER :: bCHdSV(0:nsot) = (/ 1.00, 4.05, 4.38, & 72 &4.90, 5.30, 5.39, 7.12, 7.75, 8.52,10.40,10.40,11.40, 0.02 /)72 4.90, 5.30, 5.39, 7.12, 7.75, 8.52,10.40,10.40,11.40, 0.02 /) 73 73 ! Clapp-Hornberger Coefficient b [-] 74 74 -
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/VARphy.F90
r3900 r5087 26 26 INTEGER, PARAMETER :: iun=1 27 27 REAL, PARAMETER :: zer0 = 0.0e+0, half = 0.5e+0, un_1 = 1.0e+0, & 28 & eps6 = 1.0e-6, R_1000=1.e328 eps6 = 1.0e-6, R_1000=1.e3 29 29 REAL, PARAMETER :: zero = 0.0e+0, demi = 0.5e+0, unun = 1.0e+0, & 30 & epsi = 1.0e-6, eps9 = 1.0e-930 epsi = 1.0e-6, eps9 = 1.0e-9 31 31 REAL :: ea_MAX,ea_MIN 32 32 REAL, PARAMETER :: pi = 3.141592653589793238462643e0 … … 40 40 41 41 REAL, PARAMETER :: gravit = 9.81e0, gravi2 = gravit**2 , & 42 & grvinv = 1./gravit42 grvinv = 1./gravit 43 43 !C + gravit: Earth Gravity Acceleration = 9.81 m/s2 44 44 !C + gravi2: idem (squared) … … 55 55 56 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.3157 Ls_H2O=2.8345e+6, r_LvCp=2490.04,r_LcCp=332.27, & 58 r_LsCp=2822.31 59 59 !C + cp : dry air specific heat at constant p (1004 J/kg/K) 60 60 61 61 REAL, PARAMETER :: LhfH2O=3.34e+5, LhvH2O=2.5008e+6,& 62 & LhsH2O=2.8345e+662 LhsH2O=2.8345e+6 63 63 REAL, PARAMETER :: rhoWat=1000.00e0 64 64 !C + rhoWat: Water Specific Mass = 1000.00d+0 kg/m3 … … 81 81 82 82 REAL, PARAMETER :: TfSnow=273.15e+0,csnow=2105.00e+0,r0sno=3.00e+1, & 83 & blsno=3.30e+2, Lf_H2O=3.337e+583 blsno=3.30e+2, Lf_H2O=3.337e+5 84 84 !C +... TfSnow: Snow melting Temperature= 273.15d+0 K 85 85 !C + csnow:Heat Capacity of Snow 2105 J/kg/K -
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/surf_inlandsis_mod.F90
r5086 r5087 501 501 write(un_outfor, *) fn_outfor, ikl, dt__SV, rlon(ikl2i(ikl)), rlat(ikl2i(ikl)) 502 502 write(un_outfor, *) 'nsnow - albedo - z0m - z0h , dz [m,30], temp [K,41], rho [kg/m3,41], eta [kg/kg,41] & 503 &G1 [-,30], G2 [-,30], agesnow [d,30], history [-,30], DOP [m,30]'503 G1 [-,30], G2 [-,30], agesnow [d,30], history [-,30], DOP [m,30]' 504 504 END IF 505 505 … … 963 963 dz_8SV(isl) = 0.125 * dz_dSV(isl) 964 964 dzAvSV(isl) = 0.125 * dz_dSV(islmSV(isl)) & 965 &+ 0.750 * dz_dSV(isl) &966 &+ 0.125 * dz_dSV(islpSV(isl))965 + 0.750 * dz_dSV(isl) & 966 + 0.125 * dz_dSV(islpSV(isl)) 967 967 zz_dSV = zz_dSV + dz_dSV(isl) 968 968 END DO … … 986 986 rocsSV(ist) = (1.0 - etadSV(ist)) * 1.2E+6 ! Soil Contrib. to (ro c)_s 987 987 s1__SV(ist) = bCHdSV(ist) & ! Factor of (eta)**(b+2) 988 &* psidSV(ist) * Ks_dSV(ist) & ! in DR97, Eqn.(3.36)989 &/ (etadSV(ist)**(bCHdSV(ist) + 3.)) !988 * psidSV(ist) * Ks_dSV(ist) & ! in DR97, Eqn.(3.36) 989 / (etadSV(ist)**(bCHdSV(ist) + 3.)) ! 990 990 s2__SV(ist) = Ks_dSV(ist) & ! Factor of (eta)**(2b+3) 991 &/ (etadSV(ist)**(2. * bCHdSV(ist) + 3.)) ! in DR97, Eqn.(3.35)991 / (etadSV(ist)**(2. * bCHdSV(ist) + 3.)) ! in DR97, Eqn.(3.35) 992 992 993 993 !C +--Soil Minimum Humidity (from a prescribed minimum relative Humidity) … … 995 995 Psimax = -(log(RHsMin)) / 7.2E-5 ! DR97, Eqn 3.15 Inversion 996 996 etamSV(ist) = etadSV(ist) & 997 &* (PsiMax / psidSV(ist))**(-min(10., 1. / bCHdSV(ist)))997 * (PsiMax / psidSV(ist))**(-min(10., 1. / bCHdSV(ist))) 998 998 END DO 999 999 etamSV(12) = 0. … … 1008 1008 DO ikh = 0, nkhy 1009 1009 Khyd_1 = s2__SV(ist) & ! DR97, Eqn.(3.35) 1010 &* (eta__1 **(2. * bCHdSV(ist) + 3.)) !1010 * (eta__1 **(2. * bCHdSV(ist) + 3.)) ! 1011 1011 Khyd_2 = s2__SV(ist) &! 1012 &* (eta__2 **(2. * bCHdSV(ist) + 3.)) !1012 * (eta__2 **(2. * bCHdSV(ist) + 3.)) ! 1013 1013 1014 1014 a_Khyd = (Khyd_2 - Khyd_1) / d__eta !
Note: See TracChangeset
for help on using the changeset viewer.
