!> \file noahmp_tables.f90 !! This file contains Fortran versions of the data tables included with NoahMP in mptable.tbl, soilparm.tbl, and genparm.tbl. !> \ingroup NoahMP_LSM !! \brief Data from MPTABLE.TBL, SOILPARM.TBL, GENPARM.TBL for NoahMP !! !! Note that a subset of the data in the *.TBL files is represented in this file. For example, !! only the data in the noah_mp_modis_parameters section of MPTABLE.TBL and the STAS section of !! SOILPARM.TBL are included in this module. module noahmp_tables implicit none integer, private, parameter :: mvt = 30 ! use 30 instead of 27 integer, private, parameter :: mband = 2 integer, private, parameter :: msc = 8 integer, private, parameter :: max_soiltyp = 30 integer, private, parameter :: ncrop = 5 integer, private, parameter :: nstage = 8 integer :: i integer, private, parameter :: slcats = 30 ! mptable.tbl vegetation parameters integer :: isurban_table = 13 integer :: iswater_table = 17 integer :: isbarren_table = 16 integer :: isice_table = 15 integer :: iscrop_table = 12 integer :: eblforest_table = 2 integer :: natural_table = 14 integer :: low_density_residential_table = 31 integer :: high_density_residential_table = 32 integer :: high_intensity_industrial_table = 33 ! real :: ch2op_table(mvt) !< maximum intercepted h2o per unit lai+sai (mm) data ( ch2op_table(i),i=1,mvt) / 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, & & 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, & & 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, & & 0.1, 0.1, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / real :: dleaf_table(mvt) !< characteristic leaf dimension (m) data ( dleaf_table(i),i=1,mvt) / 0.04, 0.04, 0.04, 0.04, 0.04, 0.04, & & 0.04, 0.04, 0.04, 0.04, 0.04, 0.04, & & 0.04, 0.04, 0.04, 0.04, 0.04, 0.04, & & 0.04, 0.04, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / real :: z0mvt_table(mvt) !< momentum roughness length (m) data ( z0mvt_table(i),i=1,mvt) / 1.09, 1.10, 0.85, 0.80, 0.80, 0.20, & & 0.06, 0.60, 0.50, 0.12, 0.30, 0.15, & & 1.00, 0.14, 0.00, 0.00, 0.00, 0.30, & & 0.20, 0.03, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / ! real :: hvt_table(mvt) !< top of canopy (m) data ( hvt_table(i),i=1,mvt) / 20.0, 20.0, 18.0, 16.0, 16.0, 1.10, & & 1.10, 13.0, 10.0, 1.00, 5.00, 2.00, & & 15.0, 1.50, 0.00, 0.00, 0.00, 4.00, & & 2.00, 0.50, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / real :: hvb_table(mvt) !< bottom of canopy (m) data ( hvb_table(i),i=1,mvt) / 8.50, 8.00, 7.00, 11.5, 10.0, 0.10, & & 0.10, 0.10, 0.10, 0.05, 0.10, 0.10, & & 1.00, 0.10, 0.00, 0.00, 0.00, 0.30, & & 0.20, 0.10, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / real :: den_table(mvt) !< tree density (no. of trunks per m2) data ( den_table (i),i=1,mvt) / 0.28, 0.02, 0.28, 0.10, 0.10, 10.0, & & 10.0, 10.0, 0.02, 100., 5.05, 25.0, & & 0.01, 25.0, 0.00, 0.01, 0.01, 1.00, & & 1.00, 1.00, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / ! real :: rc_table(mvt) !< tree crown radius (m) data ( rc_table (i),i=1,mvt) / 1.20, 3.60, 1.20, 1.40, 1.40, 0.12, & & 0.12, 0.12, 3.00, 0.03, 0.75, 0.08, & & 1.00, 0.08, 0.00, 0.01, 0.01, 0.30, & & 0.30, 0.30, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / real :: mfsno_table(mvt) !< snowmelt curve parameter () data ( mfsno_table(i),i=1,mvt) / 1.00, 1.00, 1.00, 1.00, 1.00, 2.00, & & 2.00, 2.00, 2.00, 2.00, 3.00, 3.00, & & 4.00, 4.00, 2.50, 3.00, 3.00, 3.50, & & 3.50, 3.50, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / real :: scffac_table(mvt) !< snow cover factor (m) data (scffac_table(i),i=1,mvt) / 0.005, 0.005, 0.005, 0.005, 0.005, & & 0.008, 0.008, 0.010, 0.010, 0.010, & & 0.010, 0.007, 0.021, 0.013, 0.015, & & 0.008, 0.015, 0.015, 0.015, 0.015, & & 0.000, 0.000, 0.000, 0.000, 0.000, & & 0.000, 0.000, 0.000, 0.000, 0.000 / ! real :: saim_table(mvt,12) !< monthly stem area index, one-sided data (saim_table (i,1),i=1,mvt) / 0.4, 0.5, 0.3, 0.4, 0.4, 0.3, & & 0.2, 0.4, 0.3, 0.3, 0.3, 0.3, & & 0.0, 0.3, 0.0, 0.0, 0.0, 0.2, & & 0.1, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / ! &_______________________________________________________________________& data (saim_table (i,2),i=1,mvt) / 0.4, 0.5, 0.3, 0.4, 0.4, 0.3, & & 0.2, 0.4, 0.3, 0.3, 0.3, 0.3, & & 0.0, 0.3, 0.0, 0.0, 0.0, 0.2, & & 0.1, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / data (saim_table (i,3),i=1,mvt) / 0.4, 0.5, 0.3, 0.4, 0.4, 0.3, & & 0.2, 0.4, 0.3, 0.3, 0.3, 0.3, & & 0.0, 0.3, 0.0, 0.0, 0.0, 0.2, & & 0.1, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / data (saim_table (i,4),i=1,mvt) / 0.3, 0.5, 0.4, 0.4, 0.4, 0.3, & & 0.2, 0.4, 0.3, 0.3, 0.3, 0.3, & & 0.0, 0.3, 0.0, 0.0, 0.0, 0.2, & & 0.1, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / ! &_______________________________________________________________________& data (saim_table (i,5),i=1,mvt) / 0.4, 0.5, 0.4, 0.4, 0.4, 0.3, & & 0.2, 0.4, 0.3, 0.3, 0.3, 0.3, & & 0.0, 0.3, 0.0, 0.0, 0.0, 0.2, & & 0.1, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / data (saim_table (i,6),i=1,mvt) / 0.5, 0.5, 0.7, 0.4, 0.4, 0.3, & & 0.2, 0.4, 0.4, 0.4, 0.4, 0.3, & & 0.0, 0.4, 0.0, 0.0, 0.0, 0.2, & & 0.2, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / data (saim_table (i,7),i=1,mvt) / 0.5, 0.5, 1.3, 0.9, 0.7, 0.6, & & 0.4, 0.7, 0.8, 0.8, 0.6, 0.4, & & 0.0, 0.6, 0.0, 0.0, 0.0, 0.4, & & 0.4, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / ! &_______________________________________________________________________& data (saim_table (i,8),i=1,mvt) / 0.6, 0.5, 1.2, 1.2, 0.8, 0.9, & & 0.6, 1.2, 1.2, 1.3, 0.9, 0.5, & & 0.0, 0.9, 0.0, 0.0, 0.0, 0.6, & & 0.6, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / data (saim_table (i,9),i=1,mvt) / 0.6, 0.5, 1.0, 1.6, 1.0, 1.2, & & 0.8, 1.4, 1.3, 1.1, 0.9, 0.4, & & 0.0, 0.7, 0.0, 0.0, 0.0, 0.8, & & 0.7, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / data (saim_table (i,10),i=1,mvt) / 0.7, 0.5, 0.8, 1.4, 1.0, 0.9, & & 0.7, 1.1, 0.7, 0.4, 0.6, 0.3, & & 0.0, 0.3, 0.0, 0.0, 0.0, 0.7, & & 0.5, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / ! &_______________________________________________________________________& data (saim_table (i,11),i=1,mvt) / 0.6, 0.5, 0.6, 0.6, 0.5, 0.4, & & 0.3, 0.5, 0.4, 0.4, 0.4, 0.3, & & 0.0, 0.3, 0.0, 0.0, 0.0, 0.3, & & 0.3, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / data (saim_table (i,12),i=1,mvt) / 0.5, 0.5, 0.5, 0.4, 0.4, 0.3, & & 0.2, 0.4, 0.4, 0.4, 0.3, 0.3, & & 0.0, 0.3, 0.0, 0.0, 0.0, 0.2, & & 0.2, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / !! lai real :: laim_table(mvt,12) !< monthly leaf area index, one-sided data (laim_table (i,1),i=1,mvt) / 4.0, 4.5, 0.0, 0.0, 2.0, 0.0, & & 0.0, 0.2, 0.3, 0.4, 0.2, 0.0, & & 0.0, 0.2, 0.0, 0.0, 0.0, 1.0, & & 0.6, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / ! &_______________________________________________________________________& data (laim_table (i,2),i=1,mvt) / 4.0, 4.5, 0.0, 0.0, 2.0, 0.0, & & 0.0, 0.2, 0.3, 0.5, 0.3, 0.0, & & 0.0, 0.3, 0.0, 0.0, 0.0, 1.0, & & 0.6, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / data (laim_table (i,3),i=1,mvt) / 4.0, 4.5, 0.0, 0.3, 2.2, 0.3, & & 0.2, 0.4, 0.5, 0.6, 0.3, 0.0, & & 0.0, 0.3, 0.0, 0.0, 0.0, 1.1, & & 0.7, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / data (laim_table (i,4),i=1,mvt) / 4.0, 4.5, 0.6, 1.2, 2.6, 0.9, & & 0.6, 1.0, 0.8, 0.7, 0.5, 0.0, & & 0.0, 0.4, 0.0, 0.0, 0.0, 1.3, & & 0.8, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / ! &_______________________________________________________________________& data (laim_table (i,5),i=1,mvt) / 4.0, 4.5, 1.2, 3.0, 3.5, 2.2, & & 1.5, 2.4, 1.8, 1.2, 1.5, 1.0, & & 0.0, 1.1, 0.0, 0.0, 0.0, 1.7, & & 1.2, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / data (laim_table (i,6),i=1,mvt) / 4.0, 4.5, 2.0, 4.7, 4.3, 3.5, & & 2.3, 4.1, 3.6, 3.0, 2.9, 2.0, & & 0.0, 2.5, 0.0, 0.0, 0.0, 2.1, & & 1.8, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / data (laim_table (i,7),i=1,mvt) / 4.0, 4.5, 2.6, 4.5, 4.3, 3.5, & & 2.3, 4.1, 3.8, 3.5, 3.5, 3.0, & & 0.0, 3.2, 0.0, 0.0, 0.0, 2.1, & & 1.8, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / ! &_______________________________________________________________________& data (laim_table (i,8),i=1,mvt) / 4.0, 4.5, 1.7, 3.4, 3.7, 2.5, & & 1.7, 2.7, 2.1, 1.5, 2.7, 3.0, & & 0.0, 2.2, 0.0, 0.0, 0.0, 1.8, & & 1.3, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / data (laim_table (i,9),i=1,mvt) / 4.0, 4.5, 1.0, 1.2, 2.6, 0.9, & & 0.6, 1.0, 0.9, 0.7, 1.2, 1.5, & & 0.0, 1.1, 0.0, 0.0, 0.0, 1.3, & & 0.8, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / data (laim_table (i,10),i=1,mvt) / 4.0, 4.5, 0.5, 0.3, 2.2, 0.3, & & 0.2, 0.4, 0.5, 0.6, 0.3, 0.0, & & 0.0, 0.3, 0.0, 0.0, 0.0, 1.1, & & 0.7, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / ! &_______________________________________________________________________& data (laim_table (i,11),i=1,mvt) / 4.0, 4.5, 0.2, 0.0, 2.0, 0.0, & & 0.0, 0.2, 0.3, 0.5, 0.3, 0.0, & & 0.0, 0.3, 0.0, 0.0, 0.0, 1.0, & & 0.6, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / data (laim_table (i,12),i=1,mvt) / 4.0, 4.5, 0.0, 0.0, 2.0, 0.0, & & 0.0, 0.2, 0.3, 0.4, 0.2, 0.0, & & 0.0, 0.2, 0.0, 0.0, 0.0, 1.0, & & 0.6, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / real :: sla_table(mvt) !< single-side leaf area per kg [m2/kg] data ( sla_table (i),i=1,mvt) / 80, 80, 80, 80, 80, 60, & & 60, 60, 50, 60, 80, 80, & & 60, 80, 0, 0, 0, 80, & & 80, 80, 0, 0, 0, 0, & & 0, 0, 0, 0, 0, 0 / real :: dilefc_table(mvt) !< coeficient for leaf stress death [1/s] data (dilefc_table (i),i=1,mvt) / 1.20, 0.50, 1.80, 0.60, 0.80, 0.20, & & 0.20, 0.20, 0.50, 0.20, 0.4, 0.50, & & 0.00, 0.35, 0.00, 0.00, 0.00, 0.30, & & 0.40, 0.30, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / real :: dilefw_table(mvt) !< coeficient for leaf stress death [1/s] data (dilefw_table(i),i=1,mvt) / 0.20, 4.00, 0.20, 0.20, 0.20, 0.20, & & 0.20, 0.20, 0.50, 0.10, 0.2, 0.20, & & 0.00, 0.20, 0.00, 0.00, 0.00, 0.20, & & 0.20, 0.20, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / real :: fragr_table(mvt) !< fraction of growth respiration !original was 0.3 data ( fragr_table(i),i=1,mvt) / 0.10, 0.20, 0.10, 0.20, 0.10, 0.20, & & 0.20, 0.20, 0.20, 0.20, 0.1, 0.20, & & 0.00, 0.20, 0.00, 0.10, 0.00, 0.10, & & 0.10, 0.10, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / real :: ltovrc_table(mvt) !< leaf turnover [1/s] data ( ltovrc_table(i),i=1,mvt) / 0.5, 0.55, 0.2, 0.55, 0.5, 0.65, & & 0.65, 0.65, 0.65, 0.50, 1.4, 1.6, & & 0.0, 1.2, 0.0, 0.0, 0.0, 1.3, & & 1.4, 1.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / ! real :: c3psn_table(mvt) !< photosynthetic pathway: 0. = c4, 1. = c3 data ( c3psn_table (i),i=1,mvt) / 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, & & 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, & & 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, & & 1.0, 1.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / real :: kc25_table(mvt) !< co2 michaelis-menten constant at 25c (pa) data ( kc25_table (i),i=1,mvt) / 30.0, 30.0, 30.0, 30.0, 30.0, 30.0, & & 30.0, 30.0, 30.0, 30.0, 30.0, 30.0, & & 30.0, 30.0, 30.0, 30.0, 30.0, 30.0, & & 30.0, 30.0, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / real :: akc_table(mvt) !< q10 for kc25 data ( akc_table (i),i=1,mvt) / 2.1, 2.1, 2.1, 2.1, 2.1, 2.1, & & 2.1, 2.1, 2.1, 2.1, 2.1, 2.1, & & 2.1, 2.1, 2.1, 2.1, 2.1, 2.1, & & 2.1, 2.1, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / real :: ko25_table(mvt) !< o2 michaelis-menten constant at 25c (pa) data ( ko25_table (i),i=1,mvt) / 3.e4, 3.e4, 3.e4, 3.e4, 3.e4, 3.e4, & & 3.e4, 3.e4, 3.e4, 3.e4, 3.e4, 3.e4, & & 3.e4, 3.e4, 3.e4, 3.e4, 3.e4, 3.e4, & & 3.e4, 3.e4, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / real :: ako_table(mvt) !< q10 for ko25 data ( ako_table (i),i=1,mvt) / 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, & & 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, & & 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, & & 1.2, 1.2, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / real :: vcmx25_table(mvt) !< maximum rate of carboxylation at 25c (umol co2/m**2/s) data ( vcmx25_table(i),i=1,mvt) / 50.0, 60.0, 60.0, 60.0, 55.0, 40.0, & & 40.0, 40.0, 40.0, 40.0, 50.0, 80.0, & & 0.00, 60.0, 0.00, 0.00, 0.00, 50.0, & & 50.0, 50.0, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / real :: avcmx_table(mvt) !< q10 for vcmx25 data ( avcmx_table (i),i=1,mvt) / 2.4, 2.4, 2.4, 2.4, 2.4, 2.4, & & 2.4, 2.4, 2.4, 2.4, 2.4, 2.4, & & 2.4, 2.4, 2.4, 2.4, 2.4, 2.4, & & 2.4, 2.4, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / real :: bp_table(mvt) !< minimum leaf conductance (umol/m**2/s) data ( bp_table (i),i=1,mvt) / 2.e3, 2.e3, 2.e3, 2.e3, 2.e3, 2.e3, & & 2.e3, 2.e3, 2.e3, 2.e3, 2.e3, 2.e3, & & 1.e15, 2.e3,1.e15, 2.e3,1.e15, 2.e3, & & 2.e3, 2.e3, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / real :: mp_table(mvt) !< slope of conductance-to-photosynthesis relationship data ( mp_table (i),i=1,mvt) / 6., 9., 6., 9., 9., 9., & & 9., 9., 9., 9., 9., 9., & & 9., 9., 9., 9., 9., 9., & & 9., 9., 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / real :: qe25_table(mvt) !< quantum efficiency at 25c (umol co2 / umo photon) data ( qe25_table (i),i=1,mvt) / 0.06, 0.06, 0.06, 0.06, 0.06, 0.06, & & 0.06, 0.06, 0.06, 0.06, 0.06, 0.06, & & 0.00, 0.06, 0.00, 0.06, 0.00, 0.06, & & 0.06, 0.06, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / real :: aqe_table(mvt) !< q10 for qe25 data ( aqe_table (i),i=1,mvt) / 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, & & 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, & & 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, & & 1.0, 1.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / real :: rmf25_table(mvt) !< leaf maintenance respiration at 25c (umol co2/m**2/s) data ( rmf25_table (i),i=1,mvt) / 3.00, 0.65, 4.00, 3.00, 3.00, 0.26, & & 0.26, 0.26, 0.80, 1.80, 3.2, 1.00, & & 0.00, 1.45, 0.00, 0.00, 0.00, 3.00, & & 3.00, 3.00, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / real :: rms25_table(mvt) !< stem maintenance respiration at 25c (umol co2/kg bio/s) data ( rms25_table (i),i=1,mvt) / 0.90, 0.30, 0.64, 0.10, 0.80, 0.10, & & 0.10, 0.10, 0.32, 0.10, 0.10, 0.10, & & 0.00, 0.10, 0.00, 0.00, 0.00, 0.10, & & 0.10, 0.00, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / real :: rmr25_table(mvt) !< root maintenance respiration at 25c (umol co2/kg bio/s) data ( rmr25_table (i),i=1,mvt) / 0.36, 0.05, 0.05, 0.01, 0.03, 0.00, & & 0.00, 0.00, 0.01, 1.20, 0.0, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 2.11, & & 2.11, 0.00, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / real :: arm_table(mvt) !< q10 for maintenance respiration data ( arm_table (i),i=1,mvt) / 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, & & 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, & & 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, & & 2.0, 2.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / real :: folnmx_table(mvt) !< foliage nitrogen concentration when f(n)=1 (%) data (folnmx_table (i),i=1,mvt) / 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, & & 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, & & 0.00, 1.5, 0.00, 1.5, 0.00, 1.5, & & 1.5, 1.5, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / real :: tmin_table(mvt) !< minimum temperature for photosynthesis (k) data ( tmin_table (i),i=1,mvt) / 265, 273, 268, 273, 268, 273, & & 273, 273, 273, 273, 268, 273, & & 0, 273, 0, 0, 0, 268, & & 268, 268, 0, 0, 0, 0, & & 0, 0, 0, 0, 0, 0 / ! real :: xl_table(mvt) !< leaf/stem orientation index data ( xl_table (i),i=1,mvt) / 0.010,0.010,0.010,0.250,0.250,0.010, & & 0.010, 0.010, 0.010, -0.30, -0.025, -0.30, & & 0.000, -0.30, 0.000, 0.000, 0.000, 0.250, & & 0.250, 0.250, 0.000, 0.000, 0.000, 0.000, & & 0.001, 0.001, 0.000, 0.000, 0.000, 0.000 / ! real :: rhol_table(mvt,mband) !< leaf reflectance: 1=vis, 2=nir data ( rhol_table (i,1),i=1,mvt) / 0.07, 0.10, 0.07, 0.10, 0.10, 0.07, & & 0.07, 0.07, 0.10, 0.11, 0.105, 0.11, & & 0.00, 0.11, 0.00, 0.00, 0.00, 0.10, & & 0.10, 0.10, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / ! &_______________________________________________________________________& data ( rhol_table (i,2),i=1,mvt) / 0.35, 0.45, 0.35, 0.45, 0.45, 0.35, & & 0.35, 0.35, 0.45, 0.58, 0.515, 0.58, & & 0.00, 0.58, 0.00, 0.00, 0.00, 0.45, & & 0.45, 0.45, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / real :: rhos_table(mvt,mband) !< stem reflectance: 1=vis, 2=nir data ( rhos_table (i,1),i=1,mvt) / 0.16, 0.16, 0.16, 0.16, 0.16, 0.16, & & 0.16, 0.16, 0.16, 0.36, 0.26, 0.36, & & 0.00, 0.36, 0.00, 0.00, 0.00, 0.16, & & 0.16,0.16, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / data ( rhos_table (i,2),i=1,mvt) / 0.39, 0.39, 0.39, 0.39, 0.39, 0.39, & & 0.39, 0.39, 0.39, 0.58, 0.485, 0.58, & & 0.00, 0.58, 0.00, 0.00, 0.00, 0.39, & & 0.39, 0.39, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / ! &_______________________________________________________________________& real :: taul_table(mvt,mband) !< leaf transmittance: 1=vis, 2=nir ! data ( taul_table (i,1),i=1,mvt) / 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, & & 0.05, 0.05, 0.05, 0.07, 0.06, 0.07, & & 0.00, 0.07, 0.00, 0.00, 0.00, 0.05, & & 0.05, 0.05,0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / data ( taul_table (i,2),i=1,mvt) / 0.10, 0.25, 0.10, 0.25, 0.25, 0.10, & & 0.10, 0.10, 0.25, 0.25, 0.25, 0.25, & & 0.00, 0.25, 0.00, 0.00, 0.00, 0.25, & & 0.25, 0.25, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / real :: taus_table(mvt,mband) !< stem transmittance: 1=vis, 2=nir data(taus_table (i,1),i=1,mvt) / 0.001,0.001,0.001,0.001,0.001, 0.001, & & 0.001, 0.001, 0.001, 0.220, 0.1105,0.220, & & 0.000, 0.220, 0.000, 0.000, 0.000, 0.001, & & 0.001, 0.001, 0.000, 0.000, 0.000, 0.000, & & 0.001, 0.001, 0.000, 0.000, 0.000, 0.000 / data(taus_table (i,2),i=1,mvt) / 0.001,0.001,0.001,0.001,0.001, 0.001, & & 0.001, 0.001, 0.001, 0.380, 0.1905, 0.380, & & 0.000, 0.380, 0.000, 0.000, 0.000, 0.001, & & 0.001, 0.001, 0.000, 0.000, 0.000, 0.000, & & 0.001, 0.001, 0.000, 0.000, 0.000, 0.000 / real :: mrp_table(mvt) !< microbial respiration parameter (umol co2 /kg c/ s) data ( mrp_table (i),i=1,mvt) / 0.37, 0.23, 0.37, 0.40, 0.30, 0.19, & & 0.19, 0.19, 0.40, 0.17,0.285, 0.23, & & 0.00, 0.23, 0.00, 0.00, 0.00, 0.23, & & 0.20, 0.00, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / ! real :: cwpvt_table(mvt) !< empirical canopy wind parameter data ( cwpvt_table (i),i=1,mvt) / 0.09, 0.335, 0.09, 0.335, 0.145, 0.50, & & 1.00, 0.65, 0.50, 2.50, 0.585, 0.835, & & 0.835, 0.835, 0.09, 0.09, 0.09, 0.335, & & 0.50, 0.09, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / real :: wrrat_table(mvt) !< wood to non-wood ratio data ( wrrat_table (i),i=1,mvt) / 30.0, 30.0, 30.0, 30.0, 30.0, 3.00, & & 3.00, 3.00, 3.00, 0.00, 15.0, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 3.00, & & 3.00, 0.00, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / real :: wdpool_table(mvt) !< wood pool (switch 1 or 0) depending on woody or not [-] data ( wdpool_table(i),i=1,mvt) / 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, & & 1.00, 1.00, 1.00, 0.00, 0.5, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 1.00, & & 1.00, 0.00, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / real :: tdlef_table(mvt) !< characteristic t for leaf freezing [k] data ( tdlef_table (i),i=1,mvt) / 278, 278, 268, 278, 268, 278, & & 278, 278, 278, 278, 268, 278, & & 278, 278, 0, 0, 0, 268, & & 268, 268, 0, 0, 0, 0, & & 0, 0, 0, 0, 0, 0 / real :: nroot_table(mvt) !< number of soil layers with root present data ( nroot_table (i),i=1,mvt) / 4, 4, 4, 4, 4, 3, & & 3, 3, 3, 3, 2, 3, & & 1, 3, 1, 1, 0, 3, & & 3, 2, 0, 0, 0, 0, & & 0, 0, 0, 0, 0, 0 / real :: rgl_table(mvt) !< parameter used in radiation stress function data ( rgl_table (i),i=1,mvt) / 30.0, 30.0, 30.0, 30.0, 30.0, 100.0,& & 100.0, 100.0, 65.0, 100.0, 65.0, 100.0, & & 999.0, 100.0, 999.0, 999.0, 30.0, 100.0, & & 100.0, 100.0, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / real :: rs_table(mvt) !< minimum stomatal resistance [s m-1] data ( rs_table (i),i=1,mvt) / 125.0, 150.0,150.0,100.0,125.0,300.0,& & 170.0,300.0, 70.0, 40.0, 70.0, 40.0, & & 200.0, 40.0, 999.0,999.0,100.0,150.0, & & 150.0, 200.0,0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / real :: hs_table(mvt) !< parameter used in vapor pressure deficit function data ( hs_table (i),i=1,mvt) / 47.35,41.69,47.35,54.53,51.93,42.00, & & 39.18, 42.00, 54.53, 36.35, 55.97, 36.25, & & 999.0, 36.25, 999.0, 999.0, 51.75, 42.00, & & 42.00, 42.00, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / real :: topt_table(mvt) !< optimum transpiration air temperature [k] data ( topt_table (i),i=1,mvt) / 298.0,298.0,298.0,298.0,298.0,298.0, & & 298.0, 298.0, 298.0, 298.0, 298.0, 298.0, & & 298.0, 298.0, 298.0, 298.0, 298.0, 298.0, & & 298.0, 298.0, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / real :: rsmax_table(mvt) !< maximal stomatal resistance [s m-1] data ( rsmax_table (i),i=1,mvt) / 5000., 5000.,5000.,5000.,5000.,5000.,& & 5000., 5000., 5000., 5000., 5000., 5000., & & 5000., 5000., 5000., 5000., 5000., 5000., & & 5000., 5000., 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / !!!!!!!!!!!!!! Wield not defined but read !!!!!!!!!!!!!!!!1 real :: slarea_table(mvt) data (slarea_table (i),i=1,mvt) / 0.0090,0.0200,0.0200,0.0258,0.0223, & & 0.0227, 0.0188, 0.0227, 0.0236, 0.0060, & & 0.0295, 0.0200, 0.0228, 0.0223, 0.02, & & 0.02, 0.0422, 0.02, 0.02, 0.02, & & 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0 / ! &_______________________________________________________________________& real :: esp1_table(mvt) data (esp1_table (i),i=1,mvt) / 0.46, 0.00, 0.00,46.86,30.98, 21.62, & & 0.11, 21.62, 22.80, 0.02, 0.815, 0.00, & & 41.87, 0.04, 0.0, 0.0, 2.31, 0.0, & & 0.0, 0.0,0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / real :: esp2_table(mvt) data (esp2_table (i),i=1,mvt) / 3.34, 0.00, 0.00, 0.38, 0.96, 0.92, & & 0.22, 0.92, 0.59, 0.05, 0.535, 0.00, & & 0.98, 0.09, 0.0, 0.0, 1.47, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / real :: esp3_table(mvt) data (esp3_table (i),i=1,mvt) / 1.85, 0.00, 0.00, 1.84, 1.84, 1.73, & & 1.26, 1.73, 1.37, 0.03, 0.605, 0.00, & & 1.82, 0.05, 0.0, 0.0, 1.70, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / ! &_______________________________________________________________________& real :: esp4_table(mvt) data (esp4_table (i),i=1,mvt) / 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / real :: esp5_table(mvt) data (esp5_table (i),i=1,mvt) / 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, & & 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 / !!!!!!!!!!!!!!!!!!! what are the tables used for !!!!!!!!!!!!!! ! soilparm.tbl parameters real :: bexp_table(max_soiltyp) data (bexp_table(i), i=1,slcats) /2.79, 4.26, 4.74, 5.33, 3.86, 5.25,& & 6.77, 8.72, 8.17, 10.73, 10.39, 11.55, & & 5.25, 0.0, 2.79, 4.26, 11.55, 2.79, & & 2.79, 0.00, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / real :: smcdry_table(max_soiltyp) data (smcdry_table(i), i=1,slcats) /0.010, 0.028, 0.047, 0.084, 0.061,& & 0.066, 0.069, 0.120, 0.103, 0.100, 0.126, 0.138, & & 0.066, 0.0, 0.006, 0.028, 0.030, 0.006, & & 0.010, 0.000, 0.000, 0.000, 0.000, 0.000, & & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000 / real :: f1_table(max_soiltyp) data (f1_table(i), i=1,slcats) /-0.472, -1.044, -0.569, 0.162, 0.162, & & -0.327, -1.491, -1.118, -1.297, -3.209, -1.916, -2.138, & & -0.327, 0.000, -1.111, -1.044, -10.472, -0.472, & & -0.472, 0.000, 0.000, 0.000, 0.000, 0.000, & & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000 / real :: smcmax_table(max_soiltyp) data (smcmax_table(i), i=1,slcats) /0.339, 0.421, 0.434, 0.476, 0.484,& & 0.439, 0.404, 0.464, 0.465, 0.406, 0.468, 0.468, & & 0.439, 1.000, 0.200, 0.421, 0.468, 0.200, & & 0.339, 0.339, 0.000, 0.000, 0.000, 0.000, & & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000 / real :: smcref_table(max_soiltyp) data (smcref_table(i), i=1,slcats) /0.192, 0.283, 0.312, 0.360, 0.347, & & 0.329, 0.315, 0.387, 0.382, 0.338, 0.404, 0.412, & & 0.329, 0.000, 0.170, 0.283, 0.454, 0.170, & & 0.192, 0.000, 0.000, 0.000, 0.000, 0.000, & & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000 / real :: psisat_table(max_soiltyp) data (psisat_table(i), i=1,slcats) /0.069, 0.036, 0.141, 0.759, 0.955, & & 0.355, 0.135, 0.617, 0.263, 0.098, 0.324, 0.468, & & 0.355, 0.00, 0.069, 0.036, 0.468, 0.069, & & 0.069, 0.00, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / real :: dksat_table(max_soiltyp) data (dksat_table(i), i=1,slcats) /4.66e-5, 1.41e-5, 5.23e-6, 2.81e-6, & & 2.18e-6, 3.38e-6, 4.45e-6, 2.03e-6, 2.45e-6,7.22e-6, & & 1.34e-6, 9.74e-7, 3.38e-6, 0.00, 1.41e-4, & & 1.41e-5, 9.74e-7, 1.41e-4, 4.66e-5,0.0, & & 0.00, 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00 / real :: dwsat_table(max_soiltyp) data (dwsat_table(i), i=1,slcats) / 2.65e-5, 5.14e-6, 8.05e-6, & & 2.39e-5, 1.66e-5, 1.43e-5, 1.01e-5, 2.35e-5, 1.13e-5, 1.87e-5, & & 9.64e-6, 1.12e-5, 1.43e-5, 0.00, 1.36e-4, 5.14e-6, & & 1.12e-5, 1.36e-4, 2.65e-5, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00 / real :: smcwlt_table(max_soiltyp) data (smcwlt_table(i), i=1,slcats) /0.010, 0.028, 0.047, 0.084, 0.061,& & 0.066, 0.069, 0.120, 0.103, 0.100, 0.126, 0.138, & & 0.066, 0.000, 0.006, 0.028, 0.030, 0.006, & & 0.010, 0.000, 0.000, 0.000, 0.000, 0.000, & & 0.000, 0.000, 0.000, 0.000, 0.000, 0.000 / real :: quartz_table(max_soiltyp) data (quartz_table(i), i=1,slcats) /0.92, 0.82, 0.60, 0.25, 0.10, & & 0.40, 0.60, 0.10, 0.35, 0.52, 0.10, & & 0.25, 0.05, 0.60, 0.07, 0.25, 0.60, & & 0.52, 0.92, 0.00, 0.00, 0.00, 0.00,0.00, & & 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 / ! genparm.tbl parameters real :: slope_table(9) !< slope factor for soil drainage data (slope_table(i), i=1,9) /0.1, 0.6, 1.0, 0.35, 0.55, 0.8, & & 0.63, 0.0, 0.0 / real :: csoil_table = 2.00e+6 !< soil heat capacity [j m-3 k-1] real :: refdk_table = 2.0e-6 !< parameter in the surface runoff parameterization real :: refkdt_table = 3.0 !< parameter in the surface runoff parameterization real :: frzk_table =0.15 !< frozen ground parameter real :: zbot_table = -8.0 !< depth [m] of lower boundary soil temperature real :: czil_table = 0.1 !< parameter used in the calculation of the roughness length for heat ! mptable.tbl radiation parameters ! &_______________________________________________________________________& real :: albsat_table(msc,mband) !< saturated soil albedos: 1=vis, 2=nir data(albsat_table(i,1),i=1,8)/0.15,0.11,0.10,0.09,0.08,0.07,0.06,0.05/ data(albsat_table(i,2),i=1,8)/0.30,0.22,0.20,0.18,0.16,0.14,0.12,0.10/ real :: albdry_table(msc,mband) !< dry soil albedos: 1=vis, 2=nir data(albdry_table(i,1),i=1,8)/0.27,0.22,0.20,0.18,0.16,0.14,0.12,0.10/ data(albdry_table(i,2),i=1,8)/0.54,0.44,0.40,0.36,0.32,0.28,0.24,0.20/ real :: albice_table(mband) !< albedo land ice: 1=vis, 2=nir data (albice_table(i),i=1,mband) /0.80, 0.55/ real :: alblak_table(mband) !< albedo frozen lakes: 1=vis, 2=nir data (alblak_table(i),i=1,mband) /0.60, 0.40/ real :: omegas_table(mband) !< two-stream parameter omega for snow data (omegas_table(i),i=1,mband) /0.8, 0.4/ real :: betads_table = 0.5 !< two-stream parameter betad for snow real :: betais_table = 0.5 !< two-stream parameter betad for snow real :: eg_table(2) !< emissivity data eg_table /0.97, 0.98 / real :: betads, betais data betads, betais /0.5, 0.5/ ! mptable.tbl global parameters real :: co2_table = 395.e-06 !< co2 partial pressure real :: o2_table = 0.209 !< o2 partial pressure real :: timean_table = 10.5 !< gridcell mean topgraphic index (global mean) real :: fsatmx_table = 0.38 !< maximum surface saturated fraction (global mean) real :: z0sno_table = 0.002 !< snow surface roughness length (m) (0.002) real :: ssi_table = 0.03 !< liquid water holding capacity for snowpack (m3/m3) (0.03) real :: snow_ret_fac_table = 5.e-5 !< snowpack water release timescale factor (1/s) real :: swemx_table = 1.00 !< new snow mass to fully cover old snow (mm) real :: tau0_table = 1.e6 !< tau0 from yang97 eqn. 10a real :: grain_growth_table = 5000. !< growth from vapor diffusion yang97 eqn. 10b real :: extra_growth_table = 10. !< extra growth near freezing yang97 eqn. 10c real :: dirt_soot_table = 0.3 !< dirt and soot term yang97 eqn. 10d real :: bats_cosz_table = 2.0 !< zenith angle snow albedo adjustment; b in yang97 eqn. 15 real :: bats_vis_new_table = 0.95 !< new snow visible albedo real :: bats_nir_new_table = 0.65 !< new snow nir albedo real :: bats_vis_age_table = 0.2 !< age factor for diffuse visible snow albedo yang97 eqn. 17 real :: bats_nir_age_table = 0.5 !< age factor for diffuse nir snow albedo yang97 eqn. 18 real :: bats_vis_dir_table = 0.4 !< cosz factor for direct visible snow albedo yang97 eqn. 15 real :: bats_nir_dir_table = 0.4 !< cosz factor for direct nir snow albedo yang97 eqn. 16 real :: rsurf_snow_table = 50.0 !< surface resistance for snow(s/m) real :: rsurf_exp_table = 5.0 !< exponent in the shape parameter for soil resistance option 1 real :: snow_emis_table = 0.95 !< surface emissivity ! Noah mp crops ! mptable.tbl crop parameters ! ! NCROP = 5 ! 1: Corn ! 2: Soybean ! 3: Sorghum ! 4: Rice ! 5: Winter wheat ! &_______________________________________________________________________& integer :: pltday_table(ncrop) !< planting date data (pltday_table(i), i=1,5) /130,111,111,111,111/ integer :: hsday_table(ncrop) !< harvest date data (hsday_table(i),i=1,5) /280,300,300,300,300/ real :: plantpop_table(ncrop) !< plant density [per ha] - used? data (plantpop_table(i),i=1,5) /78.0,78.0,78.0,78.0,78.0/ real :: irri_table(ncrop) !< irrigation strategy 0= non-irrigation 1=irrigation (no water-stress) data (irri_table(i),i=1,5) /0.0,0.0,0.0,0.0,0.0/ real :: gddtbase_table(ncrop) !< base temperature for gdd accumulation [c] data (gddtbase_table(i),i=1,5) /10.0,10.0,10.0,10.0,10.0/ real :: gddtcut_table(ncrop) !< upper temperature for gdd accumulation [c] data (gddtcut_table(i),i=1,5) /30.0,30.0,30.0,30.0,30.0/ real :: gdds1_table(ncrop) !< gdd from seeding to emergence data (gdds1_table(i),i=1,5) /60.0,50.0,50.0,50.0,50.0/ real :: gdds2_table(ncrop) !< gdd from seeding to initial vegetative data (gdds2_table(i),i=1,5) /675.0,718.0,718.0,718.0,718.0/ real :: gdds3_table(ncrop) !< gdd from seeding to post vegetative data (gdds3_table(i),i=1,5) /1183.0,933.0,933.0,933.0,933.0/ real :: gdds4_table(ncrop) !< gdd from seeding to intial reproductive data (gdds4_table(i),i=1,5) /1253.0,1103.0,1103.0,1103.0,1103.0/ real :: gdds5_table(ncrop) !< gdd from seeding to pysical maturity data (gdds5_table(i),i=1,5) /1605.0,1555.0,1555.0,1555.0,1555.0/ integer :: c3c4_table(ncrop) !< photosynthetic pathway: 1. = c3 2. = c4 data (c3c4_table(i),i=1,5) /2.0,1.0,2.0,2.0,2.0/ real :: aref_table(ncrop) !< reference maximum co2 assimulation rate data (aref_table(i),i=1,5) /7.0,7.0,7.0,7.0,7.0/ real :: psnrf_table(ncrop) !< co2 assimulation reduction factor(0-1) (caused by non-modeling part,e.g.pest,weeds) data (psnrf_table(i),i=1,5) /0.85,0.85,0.85,0.85,0.85/ real :: i2par_table(ncrop) !< fraction of incoming solar radiation to photosynthetically active radiation data (i2par_table(i),i=1,5) / 0.5,0.5,0.5,0.5,0.5/ real :: tassim0_table(ncrop) !< minimum temperature for co2 assimulation [c] data (tassim0_table(i),i=1,5) /8.0,8.0,8.0,8.0,8.0/ real :: tassim1_table(ncrop) !< co2 assimulation linearly increasing until temperature reaches t1 [c] data (tassim1_table(i),i=1,5) /18.0,18.0,18.0,18.0,18.0/ real :: tassim2_table(ncrop) !< co2 assmilation rate remain at aref until temperature reaches t2 [c] data (tassim2_table(i),i=1,5) /30.0,30.0,30.0,30.0,30.0/ real :: k_table(ncrop) !< light extinction coefficient data ( k_table(i),i=1,5) /0.55,0.55,0.55,0.55,0.55/ real :: epsi_table(ncrop) !< initial light use efficiency data (epsi_table(i),i=1,5) /12.5,12.5,12.5,12.5,12.5/ real :: q10mr_table(ncrop) !< q10 for maintainance respiration data (q10mr_table(i),i=1,5) /2.0,2.0,2.0,2.0,2.0/ real :: foln_mx_table(ncrop) !< foliage nitrogen concentration when f(n)=1 (%) data (foln_mx_table(i),i=1,5) /1.5,1.5,1.5,1.5,1.5/ real :: lefreez_table(ncrop) !< characteristic t for leaf freezing [k] data (lefreez_table(i),i=1,5) /268,268,268,268,268/ real :: dile_fc_table(ncrop,nstage) !< coeficient for temperature leaf stress death [1/s] data (dile_fc_table(i,1),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (dile_fc_table(i,2),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (dile_fc_table(i,3),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (dile_fc_table(i,4),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (dile_fc_table(i,5),i=1,5) /0.5,0.5,0.5,0.5,0.5/ data (dile_fc_table(i,6),i=1,5) /0.5,0.5,0.5,0.5,0.5/ data (dile_fc_table(i,7),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (dile_fc_table(i,8),i=1,5) /0.0,0.0,0.0,0.0,0.0/ real :: dile_fw_table(ncrop,nstage) !< coeficient for water leaf stress death [1/s] data (dile_fw_table(i,1),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (dile_fw_table(i,2),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (dile_fw_table(i,3),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (dile_fw_table(i,4),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (dile_fw_table(i,5),i=1,5) /0.2,0.2,0.2,0.2,0.2/ data (dile_fw_table(i,6),i=1,5) /0.2,0.2,0.2,0.2,0.2/ data (dile_fw_table(i,7),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (dile_fw_table(i,8),i=1,5) /0.0,0.0,0.0,0.0,0.0/ real :: fra_gr_table(ncrop) !< fraction of growth respiration data (fra_gr_table(i),i=1,5) /0.2,0.2,0.2,0.2,0.2/ real :: lf_ovrc_table(ncrop,nstage) !< fraction of leaf turnover [1/s] data (lf_ovrc_table(i,1),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (lf_ovrc_table(i,2),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (lf_ovrc_table(i,3),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (lf_ovrc_table(i,4),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (lf_ovrc_table(i,5),i=1,5) /0.2,0.48,0.48,0.48,0.48/ data (lf_ovrc_table(i,6),i=1,5) /0.3,0.48,0.48,0.48,0.48/ data (lf_ovrc_table(i,7),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (lf_ovrc_table(i,8),i=1,5) /0.0,0.0,0.0,0.0,0.0/ real :: st_ovrc_table(ncrop,nstage) !< fraction of stem turnover [1/s] data (st_ovrc_table(i,1),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (st_ovrc_table(i,2),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (st_ovrc_table(i,3),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (st_ovrc_table(i,4),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (st_ovrc_table(i,5),i=1,5) /0.12,0.12,0.12,0.12,0.12/ data (st_ovrc_table(i,6),i=1,5) /0.06,0.06,0.06,0.06,0.06/ data (st_ovrc_table(i,7),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (st_ovrc_table(i,8),i=1,5) /0.0,0.0,0.0,0.0,0.0/ real :: rt_ovrc_table(ncrop,nstage) !< fraction of root tunrover [1/s] data (rt_ovrc_table(i,1),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (rt_ovrc_table(i,2),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (rt_ovrc_table(i,3),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (rt_ovrc_table(i,4),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (rt_ovrc_table(i,5),i=1,5) /0.12,0.12,0.12,0.12,0.12/ data (rt_ovrc_table(i,6),i=1,5) /0.06,0.06,0.06,0.06,0.06/ data (rt_ovrc_table(i,7),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (rt_ovrc_table(i,8),i=1,5) /0.0,0.0,0.0,0.0,0.0/ real :: lfmr25_table(ncrop) !< leaf maintenance respiration at 25c [umol co2/m**2 /s] data (lfmr25_table(i),i=1,5) /1.0,1.0,1.0,1.0,1.0/ real :: stmr25_table(ncrop) !< stem maintenance respiration at 25c [umol co2/kg bio/s] data (stmr25_table(i),i=1,5) /0.05,0.1,0.1,0.1,0.1/ real :: rtmr25_table(ncrop) !< root maintenance respiration at 25c [umol co2/kg bio/s] data (rtmr25_table(i),i=1,5) /0.05,0.0,0.0,0.0,0.0/ real :: grainmr25_table(ncrop) !< grain maintenance respiration at 25c [umol co2/kg bio/s] data (grainmr25_table(i),i=1,5) /0.0,0.1,0.1,0.1,0.1/ real :: lfpt_table(ncrop,nstage) !< fraction of carbohydrate flux to leaf data (lfpt_table(i,1),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (lfpt_table(i,2),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (lfpt_table(i,3),i=1,5) /0.4,0.4,0.4,0.4,0.4/ data (lfpt_table(i,4),i=1,5) /0.2,0.2,0.2,0.2,0.2/ data (lfpt_table(i,5),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (lfpt_table(i,6),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (lfpt_table(i,7),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (lfpt_table(i,8),i=1,5) /0.0,0.0,0.0,0.0,0.0/ real :: stpt_table(ncrop,nstage) !< fraction of carbohydrate flux to stem data (stpt_table(i,1),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (stpt_table(i,2),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (stpt_table(i,3),i=1,5) /0.2,0.2,0.2,0.2,0.2/ data (stpt_table(i,4),i=1,5) /0.5,0.5,0.5,0.5,0.5/ data (stpt_table(i,5),i=1,5) /0.0,0.15,0.15,0.15,0.15/ data (stpt_table(i,6),i=1,5) /0.0,0.05,0.05,0.05,0.05/ data (stpt_table(i,7),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (stpt_table(i,8),i=1,5) /0.0,0.0,0.0,0.0,0.0/ real :: rtpt_table(ncrop,nstage) !< fraction of carbohydrate flux to root data (rtpt_table(i,1),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (rtpt_table(i,2),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (rtpt_table(i,3),i=1,5) /0.34,0.4,0.4,0.4,0.4/ data (rtpt_table(i,4),i=1,5) /0.3,0.3,0.3,0.3,0.3/ data (rtpt_table(i,5),i=1,5) /0.05,0.05,0.05,0.05,0.05/ data (rtpt_table(i,6),i=1,5) /0.0,0.05,0.05,0.05,0.05/ data (rtpt_table(i,7),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (rtpt_table(i,8),i=1,5) /0.0,0.0,0.0,0.0,0.0/ real :: grainpt_table(ncrop,nstage) !< fraction of carbohydrate flux to grain data (grainpt_table(i,1),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (grainpt_table(i,2),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (grainpt_table(i,3),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (grainpt_table(i,4),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (grainpt_table(i,5),i=1,5) /0.95,0.8,0.8,0.8,0.8/ data (grainpt_table(i,6),i=1,5) /1.0,0.9,0.9,0.9,0.9/ data (grainpt_table(i,7),i=1,5) /0.0,0.0,0.0,0.0,0.0/ data (grainpt_table(i,8),i=1,5) /0.0,0.0,0.0,0.0,0.0/ real :: bio2lai_table(ncrop) !< leaf are per living leaf biomass [m^2/kg] data (bio2lai_table(i),i=1,5) /0.035,0.015,0.015,0.015,0.015/ ! mptable.tbl optional parameters !------------------------------------------------------------------------------ ! Saxton and Rawls 2006 Pedo-transfer function coefficients !------------------------------------------------------------------------------ real :: sr2006_theta_1500t_a = -0.024 !< sand coefficient real :: sr2006_theta_1500t_b = 0.487 !< clay coefficient real :: sr2006_theta_1500t_c = 0.006 !< orgm coefficient real :: sr2006_theta_1500t_d = 0.005 !< sand*orgm coefficient real :: sr2006_theta_1500t_e = -0.013 !< clay*orgm coefficient real :: sr2006_theta_1500t_f = 0.068 !< sand*clay coefficient real :: sr2006_theta_1500t_g = 0.031 !< constant adjustment real :: sr2006_theta_1500_a = 0.14 !< theta_1500t coefficient real :: sr2006_theta_1500_b = -0.02 !< constant adjustment real :: sr2006_theta_33t_a = -0.251 !< sand coefficient real :: sr2006_theta_33t_b = 0.195 !< clay coefficient real :: sr2006_theta_33t_c = 0.011 !< orgm coefficient real :: sr2006_theta_33t_d = 0.006 !< sand*orgm coefficient real :: sr2006_theta_33t_e = -0.027 !< clay*orgm coefficient real :: sr2006_theta_33t_f = 0.452 !< sand*clay coefficient real :: sr2006_theta_33t_g = 0.299 !< constant adjustment real :: sr2006_theta_33_a = 1.283 !< theta_33t*theta_33t coefficient real :: sr2006_theta_33_b = -0.374 !< theta_33t coefficient real :: sr2006_theta_33_c = -0.015 !< constant adjustment real :: sr2006_theta_s33t_a = 0.278 !< sand coefficient real :: sr2006_theta_s33t_b = 0.034 !< clay coefficient real :: sr2006_theta_s33t_c = 0.022 !< orgm coefficient real :: sr2006_theta_s33t_d = -0.018 !< sand*orgm coefficient real :: sr2006_theta_s33t_e = -0.027 !< clay*orgm coefficient real :: sr2006_theta_s33t_f = -0.584 !< sand*clay coefficient real :: sr2006_theta_s33t_g = 0.078 !< constant adjustment real :: sr2006_theta_s33_a = 0.636 !< theta_s33t coefficient real :: sr2006_theta_s33_b = -0.107 !< constant adjustment real :: sr2006_psi_et_a = -21.67 !< sand coefficient real :: sr2006_psi_et_b = -27.93 !< clay coefficient real :: sr2006_psi_et_c = -81.97 !< theta_s33 coefficient real :: sr2006_psi_et_d = 71.12 !< sand*theta_s33 coefficient real :: sr2006_psi_et_e = 8.29 !< clay*theta_s33 coefficient real :: sr2006_psi_et_f = 14.05 !< sand*clay coefficient real :: sr2006_psi_et_g = 27.16 !< constant adjustment real :: sr2006_psi_e_a = 0.02 !< psi_et*psi_et coefficient real :: sr2006_psi_e_b = -0.113 !< psi_et coefficient real :: sr2006_psi_e_c = -0.7 !< constant adjustment real :: sr2006_smcmax_a = -0.097 !< sand adjustment real :: sr2006_smcmax_b = 0.043 !< constant adjustment end module noahmp_tables