! ! !################################################################## !################################################################## !###### ###### !###### SOILCST.INC ###### !###### ###### !###### Developed by ###### !###### Center for Analysis and Prediction of Storms ###### !###### University of Oklahoma ###### !###### ###### !################################################################## !################################################################## ! ! !----------------------------------------------------------------------- ! ! PURPOSE: ! ! Include file 'soilcst.inc' for ARPS soil model. ! ! This file defines the parameters and common blocks to be used ! in the soil-vegetation energy and moisture model ! !----------------------------------------------------------------------- ! ! AUTHOR: Yuhe Liu ! 09/25/1995 ! ! MODIFICATION HISTORY: ! ! 02/25/02 (J. Brotzge) ! Added several new variables to specify soil moisure/temp ! profiles as a function of soil type. ! Data taken from Chen and Dudhia, MWR, 2001; ! Cosby et al. (1984), Wat Res Res. ! and Rawls et al (1983), J. of Hydraulic Eng. ! !----------------------------------------------------------------------- ! REAL :: bslope(13) ! Slope of the retention curve DATA bslope/ 4.05, 4.38, 4.90, 5.30, 5.39, 7.12, & 7.75, 8.52, 10.40, 10.40, 11.40, 11.55, 0.00/ REAL :: cgsat(13) ! Thermal coefficient for bare ground ! at saturation DATA cgsat/ 3.222E-6, 3.057E-6, 3.560E-6, 4.418E-6, & 4.111E-6, 3.670E-6, 3.593E-6, 3.995E-6, & 3.058E-6, 3.729E-6, 3.600E-6, 1.00E-30, 8.0116E-6/ REAL :: psisat(13) ! Saturated matric potential (m) DATA psisat/ .121, .090, .218, .786, .478, .299, & .356, .630, .153, .490, .405, .036, .000/ REAL :: kns(13) ! Saturated hydraulic conductivity (m/s) DATA kns/ 1.07E-6, 1.41E-5, 5.23E-6, 2.81E-6, 3.38E-6, 4.45E-6, & 2.04E-6, 2.45E-6, 7.22E-6, 1.34E-6, 9.74E-7, 1.34E-6, & 0.0/ INTEGER :: pwgeq(13) ! Exponent in Wgeq formula DATA pwgeq/ 4, 4, 4, 6, 6, 6, & 8, 10, 8, 10, 12, 0, 0/ REAL :: awgeq(13) ! Coefficient in Wgeq formula DATA awgeq/ .387, .404, .219, .105, .148, .135, & .127, .084, .139, .075, .083, .000, .000/ REAL :: c1sat(13) ! Value of c1 at saturation (when d1=10cm) DATA c1sat/ .082, .098, .132, .153, .191, .213, & .385, .227, .421, .375, .342, .000, .000/ REAL :: c2ref(13) ! Value of c2 at wetdp = 0.5*wsat DATA c2ref/ 3.9, 3.7, 1.8, .8, .8, .8, & .4, .6, .3, .3, .3, .0, .0/ REAL :: wsat(13) ! Saturated volumetric moisture content.(ETA) DATA wsat / .421, .421, .434, .434, .439, .404, & .464, .465, .406, .406, .468, 1.0E-20, 1.00/ REAL :: wfc(13) ! Field capacity moisture. (ETA) DATA wfc / .283, .283, .312, .312, .329, .315, & .387, .382, .338, .338, .412, 1.0E-25, 1.00/ REAL :: wwlt(13) ! Wilting volumetric moisture content. (ETA) DATA wwlt / .029, .029, .047, .047, .066, .069, & .119, .103, .020, .020, .139, 1.0E-30, 1.00/ ! REAL :: wsat(13) ! Saturated volumetric moisture content. ! DATA wsat / .395, .410, .435, .485, .451, .420, & ! .477, .476, .426, .492, .482, 1.0E-20, 1.00/ ! ! REAL :: wfc(13) ! Field capacity moisture. ! DATA wfc / .135, .150, .195, .255, .240, .255, & ! .322, .325, .310, .370, .367, 1.0E-25, 1.00/ ! ! REAL :: wwlt(13) ! Wilting volumetric moisture content. ! DATA wwlt / .068, .075, .114, .179, .155, .175, & ! .218, .250, .219, .283, .286, 1.0E-30, 1.00/ REAL :: rgl(14) ! Species-dependent threshold value of ! incoming solar radiation DATA rgl / 100.0, 100.0, 100.0, 30.0, 30.0, 30.0, 30.0, & 30.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0/ REAL :: rsmin(14) ! Minimum of surface resistance DATA rsmin/ 999.0, 500.0, 40.0, 80.0, 80.0, 100.0, 150.0, & 100.0, 999.0, 40.0, 5.0, 80.0, 300.0, 1.0E-30/ REAL :: alpha(7) ! Polynomial fit used for estimating des/dT DATA alpha/ 610.7799961, 44.36518521,1.428945805, & 2.650648471E-2, 3.031240396E-4,2.034080948E-6, & 6.136820929E-9/ REAL :: quartz(13) !Quartz content DATA quartz/ 0.92, 0.82, 0.60, 0.25, 0.40, 0.60, 0.10, & 0.35, 0.52, 0.10, 0.25, 0.00, 0.00/ REAL :: porosity(13) ! Effective porosity (cm3/cm3) DATA porosity/ 0.417, 0.421, 0.434, 0.486, 0.439, 0.404, & 0.464, 0.465, 0.406, 0.423, 0.468, 0.000, & 0.000/ REAL :: rootzone(14) ! Effective root zone depth (m) DATA rootzone/ 1.0, 0.40, 1.0, 1.0, 1.0, 2.0, 2.0, & 2.0, 0.40, 1.0, 1.0, 1.0, 1.0, 2.0/ REAL :: snup(14) ! Maximum snow depth threshold at max albedo DATA snup/ 0.040, 0.025, 0.040, 0.040, 0.040, 0.080, 0.080, & 0.080, 0.025, 0.040, 0.040, 0.040, 0.025, 0.000/ REAL :: hsf2(14) ! Variable used for F2 in Rc estimate DATA hsf2/ 0.00, 42.0, 36.35, 42.0, 54.53, 54.53, 41.69, & 60.0, 999.0, 36.35, 60.0, 42.0, 200.0, 0.0/ REAL :: cgv ! Thermal coefficient for vegetation DATA cgv/1.5E-05/ REAL :: dirfc ! Directional factor of Rayleigh ! scattering ! and absorption transmission function, ! above cloud: dirfc = dirf (not for sfc.) DATA dirfc/1.6666667/ ! below & within cloud: dirfc = 5/3 REAL :: d1 ! Top layer depth of soil column: 10 cm DATA d1/ 0.10/ ! Chaning d1 requires changing C1sat. The ! latter is propotional to the inverse of d1 REAL :: d2 ! Deep layer depth of soil column: 1 m DATA d2/ 1.0/ ! Deep layer depth of soil column: 1 m REAL :: rsmax ! Maximum surface resistance in s/m DATA rsmax/5000./ ! Maximum surface resistence REAL :: wetsml ! NP, Eq. 18, small value to limit DATA wetsml/1.0E-30/ ! Small value to limit c2wg at saturation ! w2 = wsat ! DATA roufns/ 0.002, 0.020, 0.010, 0.100, 0.150, 0.750, 1.000, ! : 2.000, 0.005, 0.010, 0.020, 0.060, 0.040, 0.001/ ! DATA veg / 0.10, 0.10, 0.60, 0.40, 0.40, 0.90, 0.99, ! : 0.99, 0.01, 0.30, 0.99, 0.40, 0.20, 0.00/ REAL :: snow_albedo ! Albedo of snow covered ground (was .65, now .6) REAL :: snowdepth_crit ! Depth (m) at which snow cover is considered to ! to be 100%. REAL :: snowflxfac ! Factor by which to modify snow fluxes REAL :: cg_snow ! heat capacity of snow PARAMETER (snow_albedo=0.60,snowdepth_crit=0.025, & cg_snow=1.0/60610.0,snowflxfac=1.0)