PageRenderTime 73ms CodeModel.GetById 27ms RepoModel.GetById 0ms app.codeStats 3ms

/wrfv2_fire/phys/module_sf_noahmplsm.F

http://github.com/jbeezley/wrf-fire
FORTRAN Legacy | 8875 lines | 5630 code | 1321 blank | 1924 comment | 120 complexity | 9bdab95b1e3a0cc9ff0a5e1b9fd0b8d3 MD5 | raw file
Possible License(s): AGPL-1.0

Large files files are truncated, but you can click here to view the full file

    

  1. module noahmp_globals
    2. 

  2. 

  3.   ! Maybe most of these can be moved to a REDPRM use statement?
    4. 

  4.   use module_sf_noahlsm, only: &
    5. 

  5.        &                       SLCATS,     &
    6. 

  6.        &                       LUCATS,     &
    7. 

  7.        &                       CSOIL_DATA, & 
    8. 

  8.        &                       BB,         &
    9. 

  9.        &                       SATDK,      &
    10. 

  10.        &                       SATDW,      &
    11. 

  11.        &                       F11,        &
    12. 

  12.        &                       SATPSI,     &
    13. 

  13.        &                       QTZ,        &
    14. 

  14.        &                       DRYSMC,     &
    15. 

  15.        &                       MAXSMC,     &
    16. 

  16.        &                       REFSMC,     &
    17. 

  17.        &                       WLTSMC,     &
    18. 

  18.        &                       RSTBL,      &
    19. 

  19.        &                       RGLTBL,     &
    20. 

  20.        &                       HSTBL,      &
    21. 

  21.        &                       NROTBL,     &
    22. 

  22.        &                       TOPT_DATA,  &
    23. 

  23.        &                       RSMAX_DATA, &
    24. 

  24.        &                       ZBOT_DATA,  &
    25. 

  25.        &                       CZIL_DATA,  &
    26. 

  26.        &                       FRZK_DATA,  &
    27. 

  27.        &                       SLOPE_DATA, &
    28. 

  28.        &                       REFDK_DATA, &
    29. 

  29.        &                       REFKDT_DATA
    30. 

  30.        
    31. 

  31.   implicit none
    32. 

  32. 

  33. ! ==================================================================================================
    34. 

  34. !------------------------------------------------------------------------------------------!
    35. 

  35. ! Physical Constants:                                                                      !
    36. 

  36. !------------------------------------------------------------------------------------------!
    37. 

  37. 

  38.   REAL, PARAMETER :: GRAV   = 9.80616   !acceleration due to gravity (m/s2)
    39. 

  39.   REAL, PARAMETER :: SB     = 5.67E-08  !Stefan-Boltzmann constant (w/m2/k4)
    40. 

  40.   REAL, PARAMETER :: VKC    = 0.40      !von Karman constant
    41. 

  41.   REAL, PARAMETER :: TFRZ   = 273.16    !freezing/melting point (k)
    42. 

  42.   REAL, PARAMETER :: HSUB   = 2.8440E06 !latent heat of sublimation (j/kg)
    43. 

  43.   REAL, PARAMETER :: HVAP   = 2.5104E06 !latent heat of vaporization (j/kg)
    44. 

  44.   REAL, PARAMETER :: HFUS   = 0.3336E06 !latent heat of fusion (j/kg)
    45. 

  45.   REAL, PARAMETER :: CWAT   = 4.188E06  !specific heat capacity of water (j/m3/k)
    46. 

  46.   REAL, PARAMETER :: CICE   = 2.094E06  !specific heat capacity of ice (j/m3/k)
    47. 

  47.   REAL, PARAMETER :: CPAIR  = 1004.64   !heat capacity dry air at const pres (j/kg/k)
    48. 

  48.   REAL, PARAMETER :: TKWAT  = 0.6       !thermal conductivity of water (w/m/k)
    49. 

  49.   REAL, PARAMETER :: TKICE  = 2.2       !thermal conductivity of ice (w/m/k)
    50. 

  50.   REAL, PARAMETER :: TKAIR  = 0.023     !thermal conductivity of air (w/m/k)
    51. 

  51.   REAL, PARAMETER :: RAIR   = 287.04    !gas constant for dry air (j/kg/k)
    52. 

  52.   REAL, PARAMETER :: RW     = 461.269   !gas constant for  water vapor (j/kg/k)
    53. 

  53.   REAL, PARAMETER :: DENH2O = 1000.     !density of water (kg/m3)
    54. 

  54.   REAL, PARAMETER :: DENICE = 917.      !density of ice (kg/m3)
    55. 

  55. 

  56. !------------------------------------------------------------------------------------------!
    57. 

  57. ! From the VEGPARM.TBL tables, as functions of vegetation category.
    58. 

  58. !------------------------------------------------------------------------------------------!
    59. 

  59.   INTEGER :: NROOT        !rooting depth [as the number of layers] ( Assigned in REDPRM )
    60. 

  60.   REAL    :: RGL          !parameter used in radiation stress function ( Assigned in REDPRM )
    61. 

  61.   REAL    :: RSMIN        !minimum Canopy Resistance [s/m] ( Assigned in REDPRM )
    62. 

  62.   REAL    :: HS           !parameter used in vapor pressure deficit function ( Assigned in REDPRM )
    63. 

  63.   REAL    :: RSMAX        !maximum stomatal resistance ( Assigned in REDPRM )
    64. 

  64.   REAL    :: TOPT         !optimum transpiration air temperature.
    65. 

  65. 

  66. !KWM  CHARACTER(LEN=256) ::  LUTYPE
    67. 

  67. !KWM  INTEGER                    :: LUCATS, BARE
    68. 

  68. !KWM  INTEGER, PARAMETER         :: NLUS=50
    69. 

  69. !KWM  INTEGER, DIMENSION(1:NLUS) :: NROTBL
    70. 

  70. !KWM  REAL,    DIMENSION(1:NLUS) :: RSTBL, RGLTBL, HSTBL
    71. 

  71. !KWM  REAL                       :: TOPT_DATA,RSMAX_DATA
    72. 

  72. 

  73. ! not further used in this version (niu):
    74. 

  74. 

  75. !KWM  REAL,    DIMENSION(1:NLUS) ::  SNUPTBL, LAITBL,        &
    76. 

  76. !KWM                                 ALBTBL, SHDTBL, MAXALB
    77. 

  77. !KWM  REAL                       ::  CMCMAX_DATA,CFACTR_DATA,SBETA_DATA,&
    78. 

  78. !KWM                                 SALP_DATA  ,SMLOW_DATA ,SMHIGH_DATA
    79. 

  79. 

  80. !KWM  REAL,    DIMENSION(NLUS)   ::  LAIMINTBL    !KWM
    81. 

  81. !KWM  REAL,    DIMENSION(NLUS)   ::  LAIMAXTBL    !KWM
    82. 

  82. !KWM  REAL,    DIMENSION(NLUS)   ::  EMISSMINTBL  !KWM
    83. 

  83. !KWM  REAL,    DIMENSION(NLUS)   ::  EMISSMAXTBL  !KWM
    84. 

  84. !KWM  REAL,    DIMENSION(NLUS)   ::  ALBEDOMINTBL !KWM
    85. 

  85. !KWM  REAL,    DIMENSION(NLUS)   ::  ALBEDOMAXTBL !KWM
    86. 

  86. !KWM  REAL,    DIMENSION(NLUS)   ::  Z0MINTBL     !KWM
    87. 

  87. !KWM  REAL,    DIMENSION(NLUS)   ::  Z0MAXTBL     !KWM
    88. 

  88. 

  89. 

  90. !------------------------------------------------------------------------------------------!
    91. 

  91. ! From the SOILPARM.TBL tables, as functions of soil category.
    92. 

  92. !------------------------------------------------------------------------------------------!
    93. 

  93.   REAL    :: BEXP         !B parameter ( Assigned in REDPRM )
    94. 

  94.   REAL    :: SMCDRY       !dry soil moisture threshold where direct evap from top
    95. 

  95.                           !layer ends (volumetric) ( Assigned in REDPRM )
    96. 

  96.   REAL    :: F1           !soil thermal diffusivity/conductivity coef ( Assigned in REDPRM )
    97. 

  97.   REAL    :: SMCMAX       !porosity, saturated value of soil moisture (volumetric)
    98. 

  98.   REAL    :: SMCREF       !reference soil moisture (field capacity) (volumetric) ( Assigned in REDPRM )
    99. 

  99.   REAL    :: PSISAT       !saturated soil matric potential ( Assigned in REDPRM )
    100. 

  100.   REAL    :: DKSAT        !saturated soil hydraulic conductivity ( Assigned in REDPRM )
    101. 

  101.   REAL    :: DWSAT        !saturated soil hydraulic diffusivity ( Assigned in REDPRM )
    102. 

  102.   REAL    :: SMCWLT       !wilting point soil moisture (volumetric) ( Assigned in REDPRM )
    103. 

  103.   REAL    :: QUARTZ       !soil quartz content ( Assigned in REDPRM )
    104. 

  104. 

  105. !KWM  CHARACTER*4 SLTYPE
    106. 

  106. !KWM  INTEGER                     :: SLCATS
    107. 

  107. !KWM  INTEGER, PARAMETER          :: NSLTYPE=30
    108. 

  108. !KWM  REAL, DIMENSION (1:NSLTYPE) :: BB,DRYSMC,F11,                           &
    109. 

  109. !KWM        MAXSMC, REFSMC,SATPSI,SATDK,SATDW, WLTSMC,QTZ
    110. 

  110. !------------------------------------------------------------------------------------------!
    111. 

  111. ! From the GENPARM.TBL file
    112. 

  112. !------------------------------------------------------------------------------------------!
    113. 

  113.   REAL    :: SLOPE       !slope index (0 - 1) ( Assigned in REDPRM )
    114. 

  114.   REAL    :: CSOIL       !vol. soil heat capacity [j/m3/K] ( Assigned in REDPRM )
    115. 

  115.   REAL    :: ZBOT        !Depth (m) of lower boundary soil temperature ( Assigned in REDPRM )
    116. 

  116.   REAL    :: CZIL        !Calculate roughness length of heat ( Assigned in REDPRM )
    117. 

  117. 

  118.   REAL    :: KDT         !used in compute maximum infiltration rate (in INFIL) ( Assigned in REDPRM )
    119. 

  119.   REAL    :: FRZX        !used in compute maximum infiltration rate (in INFIL) ( Assigned in REDPRM )
    120. 

  120. 

  121. ! LSM GENERAL PARAMETERS
    122. 

  122. 

  123. !KWM  INTEGER :: SLPCATS
    124. 

  124. !KWM  INTEGER, PARAMETER :: NSLOPE=30
    125. 

  125. !KWM  REAL, DIMENSION (1:NSLOPE) :: SLOPE_DATA
    126. 

  126. !KWM  REAL ::  FXEXP_DATA,CSOIL_DATA,REFDK_DATA ,         &
    127. 

  127. !KWM           REFKDT_DATA,FRZK_DATA ,ZBOT_DATA ,CZIL_DATA
    128. 

  128. 

  129. ! =====================================options for different schemes================================
    130. 

  130. ! options for dynamic vegetation: 
    131. 

  131. ! 1 -> off (use table LAI; use FVEG = SHDFAC from input)
    132. 

  132. ! 2 -> on (together with OPT_CRS = 1)
    133. 

  133. ! 3 -> off (use table LAI; calculate FVEG)
    134. 

  134. ! 4 -> off (use table LAI; use maximum vegetation fraction)
    135. 

  135. 

  136.   INTEGER :: DVEG    != 2   !
    137. 

  137. 

  138. ! options for canopy stomatal resistance
    139. 

  139. ! 1-> Ball-Berry; 2->Jarvis
    140. 

  140. 

  141.   INTEGER :: OPT_CRS != 1    !(must 1 when DVEG = 2)
    142. 

  142. 

  143. ! options for soil moisture factor for stomatal resistance
    144. 

  144. ! 1-> Noah (soil moisture) 
    145. 

  145. ! 2-> CLM  (matric potential)
    146. 

  146. ! 3-> SSiB (matric potential)
    147. 

  147. 

  148.   INTEGER :: OPT_BTR != 1    !(suggested 1)
    149. 

  149. 

  150. ! options for runoff and groundwater
    151. 

  151. ! 1 -> TOPMODEL with groundwater (Niu et al. 2007 JGR) ;
    152. 

  152. ! 2 -> TOPMODEL with an equilibrium water table (Niu et al. 2005 JGR) ;
    153. 

  153. ! 3 -> original surface and subsurface runoff (free drainage)
    154. 

  154. ! 4 -> BATS surface and subsurface runoff (free drainage)
    155. 

  155. 

  156.   INTEGER :: OPT_RUN != 1    !(suggested 1)
    157. 

  157. 

  158. ! options for surface layer drag coeff (CH & CM)
    159. 

  159. ! 1->M-O ; 2->original Noah (Chen97); 3->MYJ consistent; 4->YSU consistent. 
    160. 

  160. 

  161.   INTEGER :: OPT_SFC != 1    !(1 or 2 or 3 or 4)
    162. 

  162. 

  163. ! options for supercooled liquid water (or ice fraction)
    164. 

  164. ! 1-> no iteration (Niu and Yang, 2006 JHM); 2: Koren's iteration 
    165. 

    166. 

  166.   INTEGER :: OPT_FRZ != 1    !(1 or 2)
    167. 

  167. 

  168. ! options for frozen soil permeability
    169. 

  169. ! 1 -> linear effects, more permeable (Niu and Yang, 2006, JHM)
    170. 

  170. ! 2 -> nonlinear effects, less permeable (old)
    171. 

  171. 

  172.   INTEGER :: OPT_INF != 1    !(suggested 1)
    173. 

  173. 

  174. ! options for radiation transfer
    175. 

  175. ! 1 -> modified two-stream (gap = F(solar angle, 3D structure ...)<1-FVEG)
    176. 

  176. ! 2 -> two-stream applied to grid-cell (gap = 0)
    177. 

  177. ! 3 -> two-stream applied to vegetated fraction (gap=1-FVEG)
    178. 

  178. 

  179.   INTEGER :: OPT_RAD != 1    !(suggested 1)
    180. 

  180. 

  181. ! options for ground snow surface albedo
    182. 

  182. ! 1-> BATS; 2 -> CLASS
    183. 

  183. 

  184.   INTEGER :: OPT_ALB != 2    !(suggested 2)
    185. 

  185. 

  186. ! options for partitioning  precipitation into rainfall & snowfall
    187. 

  187. ! 1 -> Jordan (1991); 2 -> BATS: when SFCTMP<TFRZ+2.2 ; 3-> SFCTMP<TFRZ
    188. 

  188. 

  189.   INTEGER :: OPT_SNF != 1    !(suggested 1)
    190. 

  190. 

  191. ! options for lower boundary condition of soil temperature
    192. 

  192. ! 1 -> zero heat flux from bottom (ZBOT and TBOT not used)
    193. 

  193. ! 2 -> TBOT at ZBOT (8m) read from a file (original Noah)
    194. 

  194. 

  195.   INTEGER :: OPT_TBOT != 2   !(suggested 2)
    196. 

  196. 

  197. ! options for snow/soil temperature time scheme (only layer 1)
    198. 

  198. ! 1 -> semi-implicit; 2 -> full implicit (original Noah)
    199. 

  199. 

  200.   INTEGER :: OPT_STC != 1    !(suggested 1)
    201. 

  201. ! ==================================================================================================
    202. 

  202. ! runoff parameters used for SIMTOP and SIMGM:
    203. 

  203.   REAL, PARAMETER :: TIMEAN = 10.5   !gridcell mean topgraphic index (global mean)
    204. 

  204.   REAL, PARAMETER :: FSATMX = 0.38   !maximum surface saturated fraction (global mean)
    205. 

  205. 

  206. ! adjustable parameters for snow processes
    207. 

  207. 

  208.   REAL, PARAMETER :: M      = 1.0 ! 2.50   !melting factor (-) 
    209. 

  209.   REAL, PARAMETER :: Z0SNO  = 0.002  !snow surface roughness length (m) (0.002)
    210. 

  210.   REAL, PARAMETER :: SSI    = 0.03   !liquid water holding capacity for snowpack (m3/m3) (0.03)
    211. 

  211.   REAL, PARAMETER :: SWEMX  = 1.00   !new snow mass to fully cover old snow (mm)
    212. 

  212.                                      !equivalent to 10mm depth (density = 100 kg/m3)
    213. 

  213. 

  214. ! NOTES: things to add or improve
    215. 

  215. ! 1. lake model: explicit representation of lake water storage, sunlight through lake
    216. 

  216. !    with different purity, turbulent mixing of surface laker water, snow on frozen lake, etc.
    217. 

  217. ! 2. shallow snow wihtout a layer: melting energy
    218. 

  218. ! 3. urban model to be added.
    219. 

  219. ! 4. irrigation
    220. 

  220. !------------------------------------------------------------------------------------------!
    221. 

  221. END MODULE NOAHMP_GLOBALS
    222. 

  222. !------------------------------------------------------------------------------------------!
    223. 

  223. !------------------------------------------------------------------------------------------!
    224. 

  224. MODULE NOAHMP_VEG_PARAMETERS
    225. 

  225. 

  226.     IMPLICIT NONE
    227. 

  227. 

  228.     INTEGER, PARAMETER :: MAX_VEG_PARAMS = 33
    229. 

  229.     INTEGER, PARAMETER :: MVT   = 27
    230. 

  230.     INTEGER, PARAMETER :: MBAND = 2
    231. 

  231. 

  232.     INTEGER, PRIVATE :: ISURBAN
    233. 

  233.     INTEGER :: ISWATER
    234. 

  234.     INTEGER :: ISBARREN
    235. 

  235.     INTEGER :: ISSNOW
    236. 

  236.     INTEGER :: EBLFOREST
    237. 

  237. 

  238.     REAL :: CH2OP(MVT)       !maximum intercepted h2o per unit lai+sai (mm)
    239. 

  239.     REAL :: DLEAF(MVT)       !characteristic leaf dimension (m)
    240. 

  240.     REAL :: Z0MVT(MVT)       !momentum roughness length (m)
    241. 

  241.     REAL :: HVT(MVT)         !top of canopy (m)
    242. 

  242.     REAL :: HVB(MVT)         !bottom of canopy (m)
    243. 

  243.     REAL :: DEN(MVT)         !tree density (no. of trunks per m2)
    244. 

  244.     REAL :: RC(MVT)          !tree crown radius (m)
    245. 

  245.     REAL :: SAIM(MVT,12)     !monthly stem area index, one-sided
    246. 

  246.     REAL :: LAIM(MVT,12)     !monthly leaf area index, one-sided
    247. 

  247.     REAL :: SLA(MVT)         !single-side leaf area per Kg [m2/kg]
    248. 

  248.     REAL :: DILEFC(MVT)      !coeficient for leaf stress death [1/s]
    249. 

  249.     REAL :: DILEFW(MVT)      !coeficient for leaf stress death [1/s]
    250. 

  250.     REAL :: FRAGR(MVT)       !fraction of growth respiration  !original was 0.3 
    251. 

  251.     REAL :: LTOVRC(MVT)      !leaf turnover [1/s]
    252. 

  252. 

  253.     REAL :: C3PSN(MVT)       !photosynthetic pathway: 0. = c4, 1. = c3
    254. 

  254.     REAL :: KC25(MVT)        !co2 michaelis-menten constant at 25c (pa)
    255. 

  255.     REAL :: AKC(MVT)         !q10 for kc25
    256. 

  256.     REAL :: KO25(MVT)        !o2 michaelis-menten constant at 25c (pa)
    257. 

  257.     REAL :: AKO(MVT)         !q10 for ko25
    258. 

  258.     REAL :: VCMX25(MVT)      !maximum rate of carboxylation at 25c (umol co2/m**2/s)
    259. 

  259.     REAL :: AVCMX(MVT)       !q10 for vcmx25
    260. 

  260.     REAL :: BP(MVT)          !minimum leaf conductance (umol/m**2/s)
    261. 

  261.     REAL :: MP(MVT)          !slope of conductance-to-photosynthesis relationship
    262. 

  262.     REAL :: QE25(MVT)        !quantum efficiency at 25c (umol co2 / umol photon)
    263. 

  263.     REAL :: AQE(MVT)         !q10 for qe25
    264. 

  264.     REAL :: RMF25(MVT)       !leaf maintenance respiration at 25c (umol co2/m**2/s)
    265. 

  265.     REAL :: RMS25(MVT)       !stem maintenance respiration at 25c (umol co2/kg bio/s)
    266. 

  266.     REAL :: RMR25(MVT)       !root maintenance respiration at 25c (umol co2/kg bio/s)
    267. 

  267.     REAL :: ARM(MVT)         !q10 for maintenance respiration
    268. 

  268.     REAL :: FOLNMX(MVT)      !foliage nitrogen concentration when f(n)=1 (%)
    269. 

  269.     REAL :: TMIN(MVT)        !minimum temperature for photosynthesis (k)
    270. 

  270. 

  271.     REAL :: XL(MVT)          !leaf/stem orientation index
    272. 

  272.     REAL :: RHOL(MVT,MBAND)  !leaf reflectance: 1=vis, 2=nir
    273. 

  273.     REAL :: RHOS(MVT,MBAND)  !stem reflectance: 1=vis, 2=nir
    274. 

  274.     REAL :: TAUL(MVT,MBAND)  !leaf transmittance: 1=vis, 2=nir
    275. 

  275.     REAL :: TAUS(MVT,MBAND)  !stem transmittance: 1=vis, 2=nir
    276. 

  276. 

  277.     REAL :: MRP(MVT)         !microbial respiration parameter (umol co2 /kg c/ s)
    278. 

  278.     REAL :: CWPVT(MVT)       !empirical canopy wind parameter
    279. 

  279. 

  280.     REAL :: WRRAT(MVT)       !wood to non-wood ratio
    281. 

  281.     REAL :: WDPOOL(MVT)      !wood pool (switch 1 or 0) depending on woody or not [-]
    282. 

  282.     REAL :: TDLEF(MVT)       !characteristic T for leaf freezing [K]
    283. 

  283. 

  284.     INTEGER :: IK,IM
    285. 

  285.     REAL :: TMP10(MVT*MBAND)
    286. 

  286.     REAL :: TMP11(MVT*MBAND)
    287. 

  287.     REAL :: TMP12(MVT*MBAND)
    288. 

  288.     REAL :: TMP13(MVT*MBAND)
    289. 

  289.     REAL :: TMP14(MVT*12)
    290. 

  290.     REAL :: TMP15(MVT*12)
    291. 

  291.     REAL :: TMP16(MVT*5)
    292. 

  292. 

  293.     real slarea(MVT)
    294. 

  294.     real eps(MVT,5)
    295. 

  295. 

  296. CONTAINS
    297. 

  297.   subroutine read_mp_veg_parameters(DATASET_IDENTIFIER)
    298. 

  298.     implicit none
    299. 

  299.     character(len=*), intent(in) :: DATASET_IDENTIFIER
    300. 

  300.     integer :: ierr
    301. 

  301. 

  302.     ! Temporary arrays used in reshaping namelist arrays
    303. 

  303.     REAL :: TMP10(MVT*MBAND)
    304. 

  304.     REAL :: TMP11(MVT*MBAND)
    305. 

  305.     REAL :: TMP12(MVT*MBAND)
    306. 

  306.     REAL :: TMP13(MVT*MBAND)
    307. 

  307.     REAL :: TMP14(MVT*12)
    308. 

  308.     REAL :: TMP15(MVT*12)
    309. 

  309.     REAL :: TMP16(MVT*5)
    310. 

  310. 

  311.     integer :: NVEG
    312. 

  312.     character(len=256) :: VEG_DATASET_DESCRIPTION
    313. 

  313. 

  314.     NAMELIST / noah_mp_usgs_veg_categories / VEG_DATASET_DESCRIPTION, NVEG
    315. 

  315.     NAMELIST / noah_mp_usgs_parameters / ISURBAN, ISWATER, ISBARREN, ISSNOW, EBLFOREST, &
    316. 

  316.          CH2OP, DLEAF, Z0MVT, HVT, HVB, DEN, RC, RHOL,  RHOS, TAUL, TAUS, XL, CWPVT, C3PSN, KC25, AKC, KO25, AKO, AVCMX, AQE, &
    317. 

  317.          LTOVRC,  DILEFC,  DILEFW,  RMF25 ,  SLA   ,  FRAGR ,  TMIN  ,  VCMX25,  TDLEF ,  BP, MP, QE25, RMS25, RMR25, ARM, FOLNMX, WDPOOL, WRRAT, MRP,   &
    318. 

  318.          SAIM,  LAIM,  SLAREA, EPS
    319. 

  319. 

  320.     NAMELIST / noah_mp_modis_veg_categories / VEG_DATASET_DESCRIPTION, NVEG
    321. 

  321.     NAMELIST / noah_mp_modis_parameters / ISURBAN, ISWATER, ISBARREN, ISSNOW, EBLFOREST, &
    322. 

  322.          CH2OP, DLEAF, Z0MVT, HVT, HVB, DEN, RC, RHOL,  RHOS, TAUL, TAUS, XL, CWPVT, C3PSN, KC25, AKC, KO25, AKO, AVCMX, AQE, &
    323. 

  323.          LTOVRC,  DILEFC,  DILEFW,  RMF25 ,  SLA   ,  FRAGR ,  TMIN  ,  VCMX25,  TDLEF ,  BP, MP, QE25, RMS25, RMR25, ARM, FOLNMX, WDPOOL, WRRAT, MRP,   &
    324. 

  324.          SAIM,  LAIM,  SLAREA, EPS
    325. 

  325. 

  326.     ! Initialize our variables to bad values, so that if the namelist read fails, we come to a screeching halt as soon as we try to use anything.
    327. 

  327.     CH2OP  = -1.E36
    328. 

  328.     DLEAF  = -1.E36
    329. 

  329.     Z0MVT  = -1.E36
    330. 

  330.     HVT    = -1.E36
    331. 

  331.     HVB    = -1.E36
    332. 

  332.     DEN    = -1.E36
    333. 

  333.     RC     = -1.E36
    334. 

  334.     RHOL   = -1.E36
    335. 

  335.     RHOS   = -1.E36
    336. 

  336.     TAUL   = -1.E36
    337. 

  337.     TAUS   = -1.E36
    338. 

  338.     XL     = -1.E36
    339. 

  339.     CWPVT  = -1.E36
    340. 

  340.     C3PSN  = -1.E36
    341. 

  341.     KC25   = -1.E36
    342. 

  342.     AKC    = -1.E36
    343. 

  343.     KO25   = -1.E36
    344. 

  344.     AKO    = -1.E36
    345. 

  345.     AVCMX  = -1.E36
    346. 

  346.     AQE    = -1.E36
    347. 

  347.     LTOVRC = -1.E36
    348. 

  348.     DILEFC = -1.E36
    349. 

  349.     DILEFW = -1.E36
    350. 

  350.     RMF25  = -1.E36
    351. 

  351.     SLA    = -1.E36
    352. 

  352.     FRAGR  = -1.E36
    353. 

  353.     TMIN   = -1.E36
    354. 

  354.     VCMX25 = -1.E36
    355. 

  355.     TDLEF  = -1.E36
    356. 

  356.     BP     = -1.E36
    357. 

  357.     MP     = -1.E36
    358. 

  358.     QE25   = -1.E36
    359. 

  359.     RMS25  = -1.E36
    360. 

  360.     RMR25  = -1.E36
    361. 

  361.     ARM    = -1.E36
    362. 

  362.     FOLNMX = -1.E36
    363. 

  363.     WDPOOL = -1.E36
    364. 

  364.     WRRAT  = -1.E36
    365. 

  365.     MRP    = -1.E36
    366. 

  366.     SAIM   = -1.E36
    367. 

  367.     LAIM   = -1.E36
    368. 

  368.     SLAREA = -1.E36
    369. 

  369.     EPS    = -1.E36
    370. 

  370. 

  371.     open(15, file="MPTABLE.TBL", status='old', form='formatted', action='read', iostat=ierr)
    372. 

  372.     if (ierr /= 0) then
    373. 

  373.        write(*,'("****** Error ******************************************************")')
    374. 

  374.        write(*,'("Cannot find file MPTABLE.TBL")')
    375. 

  375.        write(*,'("STOP")')
    376. 

  376.        write(*,'("*******************************************************************")')
    377. 

  377.        call wrf_error_fatal("STOP in Noah-MP read_mp_veg_parameters")
    378. 

  378.     endif
    379. 

  379. 

  380.     if ( trim(DATASET_IDENTIFIER) == "USGS" ) then
    381. 

  381.        read(15,noah_mp_usgs_veg_categories)
    382. 

  382.        read(15,noah_mp_usgs_parameters)
    383. 

  383.     else if ( trim(DATASET_IDENTIFIER) == "MODIFIED_IGBP_MODIS_NOAH" ) then
    384. 

  384.        read(15,noah_mp_modis_veg_categories)
    385. 

  385.        read(15,noah_mp_modis_parameters)
    386. 

  386.     else
    387. 

  387.        write(*,'("Unrecognized DATASET_IDENTIFIER in subroutine READ_MP_VEG_PARAMETERS")')
    388. 

  388.        write(*,'("DATASET_IDENTIFIER = ''", A, "''")') trim(DATASET_IDENTIFIER)
    389. 

  389.        call wrf_error_fatal("STOP in Noah-MP read_mp_veg_parameters")
    390. 

  390.     endif
    391. 

  391.     close(15)
    392. 

  392. 

  393.     ! Problem.  Namelist reading of 2-d arrays doesn't work well when the arrays are declared with larger dimension than the
    394. 

  394.     ! variables in the provided namelist.  So we need to reshape the 2-d arrays after we've read them.
    395. 

    396. 

  396.     if ( MVT > NVEG ) then
    397. 

  397. 

  398.        ! 
    399. 

  399.        ! Reshape the 2-d arrays:
    400. 

  400.        ! 
    401. 

  401. 

  402.        TMP10 = reshape( RHOL, (/ MVT*size(RHOL,2) /))
    403. 

  403.        TMP11 = reshape( RHOS, (/ MVT*size(RHOS,2) /))
    404. 

  404.        TMP12 = reshape( TAUL, (/ MVT*size(TAUL,2) /))
    405. 

  405.        TMP13 = reshape( TAUS, (/ MVT*size(TAUS,2) /))
    406. 

  406.        TMP14 = reshape( SAIM, (/ MVT*size(SAIM,2) /))
    407. 

  407.        TMP15 = reshape( LAIM, (/ MVT*size(LAIM,2) /))
    408. 

  408.        TMP16 = reshape( EPS,  (/ MVT*size(EPS ,2) /))
    409. 

  409. 

  410.        RHOL(1:NVEG,:) = reshape( TMP10, (/ NVEG, size(RHOL,2) /))
    411. 

  411.        RHOS(1:NVEG,:) = reshape( TMP11, (/ NVEG, size(RHOS,2) /))
    412. 

  412.        TAUL(1:NVEG,:) = reshape( TMP12, (/ NVEG, size(TAUL,2) /))
    413. 

  413.        TAUS(1:NVEG,:) = reshape( TMP13, (/ NVEG, size(TAUS,2) /))
    414. 

  414.        SAIM(1:NVEG,:) = reshape( TMP14, (/ NVEG, size(SAIM,2) /))
    415. 

  415.        LAIM(1:NVEG,:) = reshape( TMP15, (/ NVEG, size(LAIM,2) /))
    416. 

  416.        EPS(1:NVEG,:)  = reshape( TMP16, (/ NVEG, size(EPS,2)  /))
    417. 

  417. 

  418.        RHOL(NVEG+1:MVT,:) = -1.E36
    419. 

  419.        RHOS(NVEG+1:MVT,:) = -1.E36
    420. 

  420.        TAUL(NVEG+1:MVT,:) = -1.E36
    421. 

  421.        TAUS(NVEG+1:MVT,:) = -1.E36
    422. 

  422.        SAIM(NVEG+1:MVT,:) = -1.E36
    423. 

  423.        LAIM(NVEG+1:MVT,:) = -1.E36
    424. 

  424.        EPS( NVEG+1:MVT,:) = -1.E36
    425. 

  425.     endif
    426. 

  426. 

  427.   end subroutine read_mp_veg_parameters
    428. 

  428. 

  429. END MODULE NOAHMP_VEG_PARAMETERS
    430. 

  430. ! ==================================================================================================
    431. 

  431. ! ==================================================================================================
    432. 

  432. MODULE NOAHMP_RAD_PARAMETERS
    433. 

  433. 

  434.     IMPLICIT NONE
    435. 

  435.  
    436. 

  436.     INTEGER I                ! loop index
    437. 

  437.     INTEGER, PARAMETER :: MSC   = 9
    438. 

  438.     INTEGER, PARAMETER :: MBAND = 2
    439. 

  439. 

  440.     REAL :: ALBSAT(MSC,MBAND)   !saturated soil albedos: 1=vis, 2=nir
    441. 

  441.     REAL :: ALBDRY(MSC,MBAND)   !dry soil albedos: 1=vis, 2=nir
    442. 

  442.     REAL :: ALBICE(MBAND)       !albedo land ice: 1=vis, 2=nir
    443. 

  443.     REAL :: ALBLAK(MBAND)       !albedo frozen lakes: 1=vis, 2=nir
    444. 

  444.     REAL :: OMEGAS(MBAND)       !two-stream parameter omega for snow
    445. 

  445.     REAL :: BETADS              !two-stream parameter betad for snow
    446. 

  446.     REAL :: BETAIS              !two-stream parameter betad for snow
    447. 

  447.     REAL :: EG(2)               !emissivity
    448. 

  448. 

  449. ! saturated soil albedos: 1=vis, 2=nir
    450. 

  450.     DATA(ALBSAT(I,1),I=1,8)/0.15,0.11,0.10,0.09,0.08,0.07,0.06,0.05/
    451. 

  451.     DATA(ALBSAT(I,2),I=1,8)/0.30,0.22,0.20,0.18,0.16,0.14,0.12,0.10/
    452. 

  452. 

  453. ! dry soil albedos: 1=vis, 2=nir
    454. 

  454.     DATA(ALBDRY(I,1),I=1,8)/0.27,0.22,0.20,0.18,0.16,0.14,0.12,0.10/
    455. 

  455.     DATA(ALBDRY(I,2),I=1,8)/0.54,0.44,0.40,0.36,0.32,0.28,0.24,0.20/
    456. 

  456. 

  457. ! albedo land ice: 1=vis, 2=nir
    458. 

  458.     DATA (ALBICE(I),I=1,MBAND) /0.80, 0.55/
    459. 

  459. 

  460. ! albedo frozen lakes: 1=vis, 2=nir
    461. 

  461.     DATA (ALBLAK(I),I=1,MBAND) /0.60, 0.40/
    462. 

  462. 

  463. ! omega,betad,betai for snow
    464. 

  464.     DATA (OMEGAS(I),I=1,MBAND) /0.8, 0.4/
    465. 

  465.     DATA BETADS, BETAIS /0.5, 0.5/
    466. 

  466. 

  467. ! emissivity ground surface    
    468. 

  468.       DATA EG /0.97, 0.98/ ! 1-soil;2-lake
    469. 

  469. 

  470. END MODULE NOAHMP_RAD_PARAMETERS
    471. 

  471. ! ==================================================================================================
    472. 

  472. 

  473. MODULE NOAHMP_ROUTINES
    474. 

  474.   USE NOAHMP_GLOBALS
    475. 

  475.   IMPLICIT NONE
    476. 

  476. 

  477.   public  :: noahmp_options
    478. 

  478.   public  :: NOAHMP_SFLX
    479. 

  479.   public  :: REDPRM
    480. 

  480.   public  :: FRH2O           
    481. 

  481. 

  482.   private :: ATM
    483. 

  483.   private :: PHENOLOGY
    484. 

  484.   private :: ENERGY
    485. 

  485.   private ::       THERMOPROP
    486. 

  486.   private ::               CSNOW
    487. 

  487.   private ::               TDFCND
    488. 

  488.   private ::       RADIATION
    489. 

  489.   private ::               ALBEDO
    490. 

  490.   private ::                         SNOW_AGE
    491. 

  491.   private ::                         SNOWALB_BATS  
    492. 

  492.   private ::                         SNOWALB_CLASS
    493. 

  493.   private ::                         GROUNDALB
    494. 

  494.   private ::                         TWOSTREAM
    495. 

  495.   private ::               SURRAD
    496. 

  496.   private ::       VEGE_FLUX
    497. 

  497.   private ::               SFCDIF1                  
    498. 

  498.   private ::               SFCDIF2                
    499. 

  499.   private ::               STOMATA                  
    500. 

  500.   private ::               CANRES                  
    501. 

  501.   private ::               ESAT
    502. 

  502.   private ::               RAGRB
    503. 

  503.   private ::       BARE_FLUX
    504. 

  504.   private ::       TSNOSOI
    505. 

  505.   private ::               HRT
    506. 

  506.   private ::               HSTEP   
    507. 

  507.   private ::                         ROSR12
    508. 

  508.   private ::       PHASECHANGE
    509. 

  509. 

  510.   private :: WATER
    511. 

  511.   private ::       CANWATER
    512. 

  512.   private ::       SNOWWATER
    513. 

  513.   private ::               SNOWFALL
    514. 

  514.   private ::               COMBINE
    515. 

  515.   private ::               DIVIDE
    516. 

  516.   private ::                         COMBO
    517. 

  517.   private ::               COMPACT
    518. 

  518.   private ::               SNOWH2O
    519. 

  519.   private ::       SOILWATER
    520. 

  520.   private ::               ZWTEQ
    521. 

  521.   private ::               INFIL
    522. 

  522.   private ::               SRT
    523. 

  523.   private ::                         WDFCND1        
    524. 

  524.   private ::                         WDFCND2       
    525. 

  525. !  private ::                         INFIL       
    526. 

  526.   private ::               SSTEP
    527. 

  527.   private ::       GROUNDWATER
    528. 

  528. 

  529.   private :: CARBON
    530. 

  530.   private ::       CO2FLUX
    531. 

  531. !  private ::       BVOCFLUX
    532. 

  532. !  private ::       CH4FLUX
    533. 

  533. 

  534.   private :: ERROR
    535. 

  535. 

  536. contains
    537. 

  537. !
    538. 

  538. ! ==================================================================================================
    539. 

  539. 

  540.   SUBROUTINE NOAHMP_SFLX (&
    541. 

  541.        &           ILOC    , JLOC    , LAT     , YEARLEN , JULIAN  , COSZ    , & ! IN : Time/Space-related
    542. 

  542.        &           DT      , DX      , DZ8W    , NSOIL   , ZSOIL   , NSNOW   , & ! IN : Model configuration 
    543. 

  543.        &           SHDFAC  , SHDMAX  , VEGTYP  , ISURBAN , ICE     , IST     , & ! IN : Vegetation/Soil characteristics
    544. 

  544.        &           ISC     ,                                                   & ! IN : Vegetation/Soil characteristics
    545. 

  545.        &           IZ0TLND ,                                                   & ! IN : User options
    546. 

  546.        &           SFCTMP  , SFCPRS  , PSFC    , UU      , VV      , Q2      , & ! IN : Forcing
    547. 

  547.        &           QC      , SOLDN   , LWDN    , PRCP    , TBOT    , CO2AIR  , & ! IN : Forcing
    548. 

  548.        &           O2AIR   , FOLN    , FICEOLD , PBLH    ,                     & ! IN : Forcing
    549. 

  549.        &           ZLVL    , ALBOLD  , SNEQVO  ,                               & ! IN/OUT : 
    550. 

  550.        &           STC     , SH2O    , SMC     , TAH     , EAH     , FWET    , & ! IN/OUT : 
    551. 

  551.        &           CANLIQ  , CANICE  , TV      , TG      , QSFC    , QSNOW   , & ! IN/OUT : 
    552. 

  552.        &           ISNOW   , ZSNSO   , SNOWH   , SNEQV   , SNICE   , SNLIQ   , & ! IN/OUT : 
    553. 

  553.        &           ZWT     , WA      , WT      , WSLAKE  , LFMASS  , RTMASS  , & ! IN/OUT : 
    554. 

  554.        &           STMASS  , WOOD    , STBLCP  , FASTCP  , LAI     , SAI     , & ! IN/OUT : 
    555. 

  555.        &           CM      , CH      , CHSTAR  ,                               & ! IN/OUT : 
    556. 

  556.        &           FSA     , FSR     , FIRA    , FSH     , SSOIL   , FCEV    , & ! OUT : 
    557. 

  557.        &           FGEV    , FCTR    , ECAN    , ETRAN   , EDIR    , TRAD    , & ! OUT :
    558. 

  558.        &           TS      , TGB     , TGV     , T2MV    , T2MB    , TSTAR   , & ! OUT :
    559. 

  559.        &           Q1      , Q2V     , Q2B     , RUNSRF  , RUNSUB  , APAR    , & ! OUT :
    560. 

  560.        &           PSN     , SAV     , SAG     , FSNO    , NEE     , GPP     , & ! OUT :
    561. 

  561.        &           NPP     , FVEG    , ALBEDO  , QMELT   , PONDING , PONDING1, & ! OUT :
    562. 

  562.        &           PONDING2, RSSUN   , RSSHA   , BGAP    , WGAP    , GAP     , & ! OUT :
    563. 

  563.        &           ERRWAT  , CHV     , CHB     , EMISSI)                         ! OUT :
    564. 

  564. 

  565. ! --------------------------------------------------------------------------------------------------
    566. 

  566. ! Initial code: Guo-Yue Niu, Oct. 2007
    567. 

  567. ! --------------------------------------------------------------------------------------------------
    568. 

  568.   USE NOAHMP_VEG_PARAMETERS
    569. 

  569.   USE NOAHMP_RAD_PARAMETERS
    570. 

  570. ! --------------------------------------------------------------------------------------------------
    571. 

  571.   implicit none
    572. 

  572. ! --------------------------------------------------------------------------------------------------
    573. 

  573. ! input
    574. 

  574.   INTEGER                        , INTENT(IN)    :: ICE    !ice (ice = 1)
    575. 

  575.   INTEGER                        , INTENT(IN)    :: IST    !surface type 1->soil; 2->lake
    576. 

  576.   INTEGER                        , INTENT(IN)    :: VEGTYP !vegetation type 
    577. 

  577.   INTEGER                        , INTENT(IN)    :: ISC    !soil color type (1-lighest; 8-darkest)
    578. 

  578.   INTEGER                        , INTENT(IN)    :: NSNOW  !maximum no. of snow layers        
    579. 

  579.   INTEGER                        , INTENT(IN)    :: NSOIL  !no. of soil layers        
    580. 

  580.   INTEGER                        , INTENT(IN)    :: ILOC   !grid index
    581. 

  581.   INTEGER                        , INTENT(IN)    :: JLOC   !grid index
    582. 

  582.   REAL                           , INTENT(IN)    :: DT     !time step [sec]
    583. 

  583.   REAL, DIMENSION(       1:NSOIL), INTENT(IN)    :: ZSOIL  !layer-bottom depth from soil surf (m)
    584. 

  584.   REAL                           , INTENT(IN)    :: Q2     !mixing ratio (kg/kg) lowest model layer
    585. 

  585.   REAL                           , INTENT(IN)    :: SFCTMP !surface air temperature [K]
    586. 

  586.   REAL                           , INTENT(IN)    :: UU     !wind speed in eastward dir (m/s)
    587. 

  587.   REAL                           , INTENT(IN)    :: VV     !wind speed in northward dir (m/s)
    588. 

  588.   REAL                           , INTENT(IN)    :: SOLDN  !downward shortwave radiation (w/m2)
    589. 

  589.   REAL                           , INTENT(IN)    :: PRCP   !precipitation rate (kg m-2 s-1)
    590. 

  590.   REAL                           , INTENT(IN)    :: LWDN   !downward longwave radiation (w/m2)
    591. 

  591.   REAL                           , INTENT(IN)    :: SFCPRS !pressure (pa)
    592. 

  592.   REAL                           , INTENT(INOUT) :: ZLVL   !reference height (m)
    593. 

  593.   REAL                           , INTENT(IN)    :: COSZ   !cosine solar zenith angle [0-1]
    594. 

  594.   REAL                           , INTENT(IN)    :: TBOT   !bottom condition for soil temp. [K]
    595. 

  595.   REAL                           , INTENT(IN)    :: FOLN   !foliage nitrogen (%) [1-saturated]
    596. 

  596.   REAL                           , INTENT(IN)    :: SHDFAC !green vegetation fraction [0.0-1.0]
    597. 

  597.   INTEGER                        , INTENT(IN)    :: YEARLEN!Number of days in the particular year.
    598. 

  598.   REAL                           , INTENT(IN)    :: JULIAN !Julian day of year (floating point)
    599. 

  599.   REAL                           , INTENT(IN)    :: LAT    !latitude (radians)
    600. 

  600.   REAL, DIMENSION(-NSNOW+1:    0), INTENT(IN)    :: FICEOLD!ice fraction at last timestep
    601. 

  601. 

  602. !jref:start; in 
    603. 

  603.   INTEGER                        , INTENT(IN)    :: ISURBAN
    604. 

  604.   INTEGER                        , INTENT(IN)    :: IZ0TLND
    605. 

  605.   REAL                           , INTENT(IN)    :: QC     !cloud water mixing ratio
    606. 

  606.   REAL                           , INTENT(IN)    :: PBLH   !planetary boundary layer height
    607. 

  607.   REAL                           , INTENT(INOUT)    :: QSFC   !mixing ratio at lowest model layer
    608. 

  608.   REAL                           , INTENT(IN)    :: PSFC   !pressure at lowest model layer
    609. 

  609.   REAL                           , INTENT(IN)    :: DZ8W   !thickness of lowest layer
    610. 

  610.   REAL                           , INTENT(IN)    :: DX
    611. 

  611.   REAL                           , INTENT(IN)    :: SHDMAX  !yearly max vegetation fraction
    612. 

  612. !jref:end
    613. 

  613. 

  614. ! input/output : need arbitary intial values
    615. 

  615.   REAL                           , INTENT(INOUT) :: QSNOW  !snowfall [mm/s]
    616. 

  616.   REAL                           , INTENT(INOUT) :: FWET   !wetted or snowed fraction of canopy (-)
    617. 

  617.   REAL                           , INTENT(INOUT) :: SNEQVO !snow mass at last time step (mm)
    618. 

  618.   REAL                           , INTENT(INOUT) :: EAH    !canopy air vapor pressure (pa)
    619. 

  619.   REAL                           , INTENT(INOUT) :: TAH    !canopy air tmeperature (k)
    620. 

  620.   REAL                           , INTENT(INOUT) :: ALBOLD !snow albedo at last time step (CLASS type)
    621. 

  621.   REAL                           , INTENT(INOUT) :: CM     !momentum drag coefficient
    622. 

  622.   REAL                           , INTENT(INOUT) :: CH     !sensible heat exchange coefficient
    623. 

  623. 

  624. ! prognostic variables
    625. 

  625.   INTEGER                        , INTENT(INOUT) :: ISNOW  !actual no. of snow layers [-]
    626. 

  626.   REAL                           , INTENT(INOUT) :: CANLIQ !intercepted liquid water (mm)
    627. 

  627.   REAL                           , INTENT(INOUT) :: CANICE !intercepted ice mass (mm)
    628. 

  628.   REAL                           , INTENT(INOUT) :: SNEQV  !snow water eqv. [mm]
    629. 

  629.   REAL, DIMENSION(       1:NSOIL), INTENT(INOUT) :: SMC    !soil moisture (ice + liq.) [m3/m3]
    630. 

  630.   REAL, DIMENSION(-NSNOW+1:NSOIL), INTENT(INOUT) :: ZSNSO  !layer-bottom depth from snow surf [m]
    631. 

  631.   REAL                           , INTENT(INOUT) :: SNOWH  !snow height [m]
    632. 

  632.   REAL, DIMENSION(-NSNOW+1:    0), INTENT(INOUT) :: SNICE  !snow layer ice [mm]
    633. 

  633.   REAL, DIMENSION(-NSNOW+1:    0), INTENT(INOUT) :: SNLIQ  !snow layer liquid water [mm]
    634. 

  634.   REAL                           , INTENT(INOUT) :: TV     !vegetation temperature (k)
    635. 

  635.   REAL                           , INTENT(INOUT) :: TG     !ground temperature (k)
    636. 

  636.   REAL, DIMENSION(-NSNOW+1:NSOIL), INTENT(INOUT) :: STC    !snow/soil temperature [k]
    637. 

  637.   REAL, DIMENSION(       1:NSOIL), INTENT(INOUT) :: SH2O   !liquid soil moisture [m3/m3]
    638. 

  638.   REAL                           , INTENT(INOUT) :: ZWT    !depth to water table [m]
    639. 

  639.   REAL                           , INTENT(INOUT) :: WA     !water storage in aquifer [mm]
    640. 

  640.   REAL                           , INTENT(INOUT) :: WT     !water in aquifer&saturated soil [mm]
    641. 

  641.   REAL                           , INTENT(INOUT) :: WSLAKE !lake water storage (can be neg.) (mm)
    642. 

  642. 

  643. ! output
    644. 

  644.   REAL                           , INTENT(OUT)   :: FSA    !total absorbed solar radiation (w/m2)
    645. 

  645.   REAL                           , INTENT(OUT)   :: FSR    !total reflected solar radiation (w/m2)
    646. 

  646.   REAL                           , INTENT(OUT)   :: FIRA   !total net LW rad (w/m2)  [+ to atm]
    647. 

  647.   REAL                           , INTENT(OUT)   :: FSH    !total sensible heat (w/m2) [+ to atm]
    648. 

  648.   REAL                           , INTENT(OUT)   :: FCEV   !canopy evap heat (w/m2) [+ to atm]
    649. 

  649.   REAL                           , INTENT(OUT)   :: FGEV   !ground evap heat (w/m2) [+ to atm]
    650. 

  650.   REAL                           , INTENT(OUT)   :: FCTR   !transpiration heat (w/m2) [+ to atm]
    651. 

  651.   REAL                           , INTENT(OUT)   :: SSOIL  !ground heat flux (w/m2)   [+ to soil]
    652. 

  652.   REAL                           , INTENT(OUT)   :: TRAD   !surface radiative temperature (k)
    653. 

  653.   REAL                           , INTENT(OUT)   :: TS     !surface temperature (k)
    654. 

  654.   REAL                           , INTENT(OUT)   :: ECAN   !evaporation of intercepted water (mm/s)
    655. 

  655.   REAL                           , INTENT(OUT)   :: ETRAN  !transpiration rate (mm/s)
    656. 

  656.   REAL                           , INTENT(OUT)   :: EDIR   !soil surface evaporation rate (mm/s]
    657. 

  657.   REAL                           , INTENT(OUT)   :: RUNSRF !surface runoff [mm/s] 
    658. 

  658.   REAL                           , INTENT(OUT)   :: RUNSUB !baseflow (saturation excess) [mm/s]
    659. 

  659.   REAL                           , INTENT(OUT)   :: PSN    !total photosynthesis (umol co2/m2/s) [+]
    660. 

  660.   REAL                           , INTENT(OUT)   :: APAR   !photosyn active energy by canopy (w/m2)
    661. 

  661.   REAL                           , INTENT(OUT)   :: SAV    !solar rad absorbed by veg. (w/m2)
    662. 

  662.   REAL                           , INTENT(OUT)   :: SAG    !solar rad absorbed by ground (w/m2)
    663. 

  663.   REAL                           , INTENT(OUT)   :: FSNO   !snow cover fraction on the ground (-)
    664. 

  664.   REAL                           , INTENT(OUT)   :: FVEG   !green vegetation fraction [0.0-1.0]
    665. 

  665.   REAL                           , INTENT(OUT)   :: ALBEDO !surface albedo [-]
    666. 

  666.   REAL                           , INTENT(OUT)   :: ERRWAT !water error [kg m{-2}]
    667. 

  667.   REAL                           , INTENT(OUT)   :: QMELT  !snowmelt [mm/s]
    668. 

  668.   REAL                           , INTENT(OUT)   :: PONDING!surface ponding [mm]
    669. 

  669.   REAL                           , INTENT(OUT)   :: PONDING1!surface ponding [mm]
    670. 

  670.   REAL                           , INTENT(OUT)   :: PONDING2!surface ponding [mm]
    671. 

  671. 

  672. !jref:start; output
    673. 

  673.   REAL                           , INTENT(OUT)     :: CHSTAR !effective exchange coefficient
    674. 

  674.   REAL                           , INTENT(OUT)     :: TSTAR  !effective skin temperature
    675. 

  675.   REAL                           , INTENT(OUT)     :: T2MV   !2-m air temperature over vegetated part [k]
    676. 

  676.   REAL                           , INTENT(OUT)     :: T2MB   !2-m air temperature over bare ground part [k]
    677. 

  677.   REAL, INTENT(OUT) :: RSSUN        !sunlit leaf stomatal resistance (s/m)
    678. 

  678.   REAL, INTENT(OUT) :: RSSHA        !shaded leaf stomatal resistance (s/m)
    679. 

  679.   REAL, INTENT(OUT) :: BGAP
    680. 

  680.   REAL, INTENT(OUT) :: WGAP
    681. 

  681.   REAL, INTENT(OUT) :: GAP
    682. 

  682.   REAL, INTENT(OUT) :: TGV
    683. 

  683.   REAL, INTENT(OUT) :: TGB
    684. 

  684.   REAL, INTENT(OUT) :: Q1
    685. 

  685.   REAL, INTENT(OUT) :: EMISSI
    686. 

  686. !jref:end
    687. 

  687. 

  688. ! local
    689. 

  689.   INTEGER                                        :: IZ     !do-loop index
    690. 

  690.   INTEGER, DIMENSION(-NSNOW+1:NSOIL)             :: IMELT  !phase change index [1-melt; 2-freeze]
    691. 

  691.   REAL                                           :: CMC    !intercepted water (CANICE+CANLIQ) (mm)
    692. 

  692.   REAL                                           :: TAUX   !wind stress: e-w (n/m2)
    693. 

  693.   REAL                                           :: TAUY   !wind stress: n-s (n/m2)
    694. 

  694.   REAL                                           :: RHOAIR !density air (kg/m3)
    695. 

  695. !  REAL, DIMENSION(       1:    5)                :: VOCFLX !voc fluxes [ug C m-2 h-1]
    696. 

  696.   REAL, DIMENSION(-NSNOW+1:NSOIL)                :: DZSNSO !snow/soil layer thickness [m]
    697. 

  697.   REAL                                           :: THAIR  !potential temperature (k)
    698. 

  698.   REAL                                           :: QAIR   !specific humidity (kg/kg) (q2/(1+q2))
    699. 

  699.   REAL                                           :: EAIR   !vapor pressure air (pa)
    700. 

  700.   REAL, DIMENSION(       1:    2)                :: SOLAD  !incoming direct solar rad (w/m2)
    701. 

  701.   REAL, DIMENSION(       1:    2)                :: SOLAI  !incoming diffuse solar rad (w/m2)
    702. 

  702.   REAL                                           :: QPRECC !convective precipitation (mm/s)
    703. 

  703.   REAL                                           :: QPRECL !large-scale precipitation (mm/s)
    704. 

  704.   REAL                                           :: IGS    !growing season index (0=off, 1=on)
    705. 

  705.   REAL                                           :: ELAI   !leaf area index, after burying by snow
    706. 

  706.   REAL                                           :: ESAI   !stem area index, after burying by snow
    707. 

  707.   REAL                                           :: BEVAP  !soil water evaporation factor (0 - 1)
    708. 

  708.   REAL, DIMENSION(       1:NSOIL)                :: BTRANI !Soil water transpiration factor (0 - 1)
    709. 

  709.   REAL                                           :: BTRAN  !soil water transpiration factor (0 - 1)
    710. 

  710.   REAL                                           :: HTOP   !top of canopy layer (m)
    711. 

  711.   REAL                                           :: QIN    !groundwater recharge [mm/s]
    712. 

  712.   REAL                                           :: QDIS   !groundwater discharge [mm/s]
    713. 

  713.   REAL, DIMENSION(       1:NSOIL)                :: SICE   !soil ice content (m3/m3)
    714. 

  714.   REAL, DIMENSION(-NSNOW+1:    0)                :: SNICEV !partial volume ice of snow [m3/m3]
    715. 

  715.   REAL, DIMENSION(-NSNOW+1:    0)                :: SNLIQV !partial volume liq of snow [m3/m3]
    716. 

  716.   REAL, DIMENSION(-NSNOW+1:    0)                :: EPORE  !effective porosity [m3/m3]
    717. 

  717.   REAL                                           :: TOTSC  !total soil carbon (g/m2)
    718. 

  718.   REAL                                           :: TOTLB  !total living carbon (g/m2)
    719. 

  719.   REAL                                           :: T2M    !2-meter air temperature (k)
    720. 

  720.   REAL                                           :: QDEW   !ground surface dew rate [mm/s]
    721. 

  721.   REAL                                           :: QVAP   !ground surface evap. rate [mm/s]
    722. 

  722.   REAL                                           :: LATHEA !latent heat [j/kg]
    723. 

  723.   REAL                                           :: SWDOWN !downward solar [w/m2]
    724. 

  724.   REAL                                           :: BEG_WB !water storage at begin of a step [mm]
    725. 

  725. 

  726. !jref:start 
    727. 

  727.   REAL                                           :: FSRV
    728. 

  728.   REAL                                           :: FSRG
    729. 

  729.   REAL,INTENT(OUT)                               :: Q2V
    730. 

  730.   REAL,INTENT(OUT)                               :: Q2B
    731. 

  731.   REAL :: Q2E
    732. 

  732.   REAL :: QFX
    733. 

  733.   REAL,INTENT(OUT)                               :: CHV    !sensible heat exchange coefficient over vegetated fraction
    734. 

  734.   REAL,INTENT(OUT)                               :: CHB    !sensible heat exchange coefficient over bare-ground
    735. 

  735. !jref:end  
    736. 

  736. 

  737. ! carbon
    738. 

  738. ! inputs
    739. 

  739.   REAL                           , INTENT(IN)    :: CO2AIR !atmospheric co2 concentration (pa)
    740. 

  740.   REAL                           , INTENT(IN)    :: O2AIR  !atmospheric o2 concentration (pa)
    741. 

  741. 

  742. ! inputs and outputs : prognostic variables
    743. 

  743.   REAL                        , INTENT(INOUT)    :: LFMASS !leaf mass [g/m2]
    744. 

  744.   REAL                        , INTENT(INOUT)    :: RTMASS !mass of fine roots [g/m2]
    745. 

  745.   REAL                        , INTENT(INOUT)    :: STMASS !stem mass [g/m2]
    746. 

  746.   REAL                        , INTENT(INOUT)    :: WOOD   !mass of wood (incl. woody roots) [g/m2]
    747. 

  747.   REAL                        , INTENT(INOUT)    :: STBLCP !stable carbon in deep soil [g/m2]
    748. 

  748.   REAL                        , INTENT(INOUT)    :: FASTCP !short-lived carbon, shallow soil [g/m2]
    749. 

  749.   REAL                        , INTENT(INOUT)    :: LAI    !leaf area index [-]
    750. 

  750.   REAL                        , INTENT(INOUT)    :: SAI    !stem area index [-]
    751. 

  751. 

  752. ! outputs
    753. 

  753.   REAL                          , INTENT(OUT)    :: NEE    !net ecosys exchange (g/m2/s CO2)
    754. 

  754.   REAL                          , INTENT(OUT)    :: GPP    !net instantaneous assimilation [g/m2/s C]
    755. 

  755.   REAL                          , INTENT(OUT)    :: NPP    !net primary productivity [g/m2/s C]
    756. 

  756.   REAL                                           :: AUTORS !net ecosystem respiration (g/m2/s C)
    757. 

  757.   REAL                                           :: HETERS !organic respiration (g/m2/s C)
    758. 

  758.   REAL                                           :: TROOT  !root-zone averaged temperature (k)
    759. 

  759. 

  760.   ! INTENT (OUT) variables need to be assigned a value.  These normally get assigned values
    761. 

  761.   ! only if DVEG == 2.
    762. 

  762.   nee = 0.0
    763. 

  763.   npp = 0.0
    764. 

  764.   gpp = 0.0
    765. 

  765. 

  766. ! --------------------------------------------------------------------------------------------------
    767. 

  767. ! re-process atmospheric forcing
    768. 

  768. 

  769.    CALL ATM (SFCPRS ,SFCTMP ,Q2     ,PRCP   ,SOLDN  ,COSZ   ,THAIR  , & 
    770. 

  770.              QAIR   ,EAIR   ,RHOAIR ,QPRECC ,QPRECL ,SOLAD  ,SOLAI  , &
    771. 

  771.              SWDOWN )
    772. 

  772. 

  773. ! snow/soil layer thickness (m)
    774. 

  774. 

  775.      DO IZ = ISNOW+1, NSOIL
    776. 

  776.          IF(IZ == ISNOW+1) THEN
    777. 

  777.            DZSNSO(IZ) = - ZSNSO(IZ)
    778. 

  778.          ELSE
    779. 

  779.            DZSNSO(IZ) = ZSNSO(IZ-1) - ZSNSO(IZ)
    780. 

  780.          END IF
    781. 

  781.      END DO
    782. 

  782. 

  783. ! root-zone temperature
    784. 

  784. 

  785.      TROOT  = 0.
    786. 

  786.      DO IZ=1,NROOT
    787. 

  787.         TROOT = TROOT + STC(IZ)*DZSNSO(IZ)/(-ZSOIL(NROOT))
    788. 

  788.      ENDDO
    789. 

  789. 

  790. ! total water storage for water balance check
    791. 

  791.     
    792. 

  792.      IF(IST == 1) THEN
    793. 

  793.      BEG_WB = CANLIQ + CANICE + SNEQV + WA
    794. 

  794.      DO IZ = 1,NSOIL
    795. 

  795.         BEG_WB = BEG_WB + SMC(IZ) * DZSNSO(IZ) * 1000.
    796. 

  796.      END DO
    797. 

  797.      END IF
    798. 

  798. 

  799. ! vegetation phenology
    800. 

  800. 

  801.      CALL PHENOLOGY (VEGTYP , SNOWH  , TV     , LAT   , YEARLEN , JULIAN , & !in
    802. 

  802.                      LAI    , SAI    , TROOT  , HTOP  , ELAI    , ESAI   ,IGS)
    803. 

  803. 

  804. !input GVF should be consistent with LAI
    805. 

  805.      IF(DVEG == 1) THEN
    806. 

  806.         FVEG = SHDFAC
    807. 

  807.         IF(FVEG <= 0.05) FVEG = 0.05
    808. 

  808.      ELSE IF (DVEG == 2 .or. DVEG == 3) THEN
    809. 

  809.         FVEG = 1.-EXP(-0.52*(LAI+SAI))
    810. 

  810.         IF(FVEG <= 0.05) FVEG = 0.05
    811. 

  811.      ELSE IF (DVEG == 4) THEN
    812. 

  812.         FVEG = SHDMAX
    813. 

  813.         IF(FVEG <= 0.05) FVEG = 0.05
    814. 

  814.      ELSE
    815. 

  815.         WRITE(*,*) "-------- FATAL CALLED IN SFLX -----------"
    816. 

  816.         CALL wrf_error_fatal("Namelist parameter DVEG unknown") 
    817. 

  817.      ENDIF
    818. 

  818. 

  819. !     CALL PHENOLOGY (VEGTYP,IMONTH ,IDAY   ,SNOWH  ,TV     ,LAT   , & !in
    820. 

  820. !                     LAI   ,SAI    ,TROOT  ,                & !in
    821. 

  821. !                     HTOP  ,ELAI   ,ESAI   ,IGS    )   !out
    822. 

  822. 

  823. ! compute energy budget (momentum & energy fluxes and phase changes) 
    824. 

  824. 

  825.     CALL ENERGY (ICE    ,VEGTYP ,IST    ,ISC    ,NSNOW  ,NSOIL  , & !in
    826. 

  826.                  ISNOW  ,NROOT  ,DT     ,RHOAIR ,SFCPRS ,QAIR   , & !in
    827. 

  827.                  SFCTMP ,THAIR  ,LWDN   ,UU     ,VV     ,ZLVL   , & !in
    828. 

  828.                  CO2AIR ,O2AIR  ,SOLAD  ,SOLAI  ,COSZ   ,IGS    , & !in
    829. 

  829.                  EAIR   ,HTOP   ,TBOT   ,ZBOT   ,ZSNSO  ,ZSOIL  , & !in
    830. 

  830.                  ELAI   ,ESAI   ,CSOIL  ,FWET   ,FOLN   ,         & !in
    831. 

  831.                  FVEG   ,                                         & !in
    832. 

  832.                  QSNOW  ,DZSNSO ,LAT    ,CANLIQ ,CANICE ,iloc, jloc , & !in
    833. 

  833.                  IMELT  ,SNICEV ,SNLIQV ,EPORE  ,T2M    ,FSNO   , & !out
    834. 

  834.                  SAV    ,SAG    ,QMELT  ,FSA    ,FSR    ,TAUX   , & !out
    835. 

  835.                  TAUY   ,FIRA   ,FSH    ,FCEV   ,FGEV   ,FCTR   , & !out
    836. 

  836.                  TRAD   ,PSN    ,APAR   ,SSOIL  ,BTRANI ,BTRAN  , & !out
    837. 

  837.                  PONDING,TS     ,LATHEA ,                         & !out
    838. 

  838.                  TV     ,TG     ,STC    ,SNOWH  ,EAH    ,TAH    , & !inout
    839. 

  839.                  SNEQVO ,SNEQV  ,SH2O   ,SMC    ,SNICE  ,SNLIQ  , & !inout
    840. 

  840.                  ALBOLD ,CM     ,CH     ,DX     ,DZ8W   ,Q2     , & !inout
    841. 

  841. !jref:start
    842. 

  842.                  QC     ,PBLH   ,QSFC   ,PSFC   ,ISURBAN,IZ0TLND, & !in 
    843. 

  843.                  CHSTAR ,TSTAR  ,T2MV   ,T2MB  ,FSRV   , &
    844. 

  844.                  FSRG   ,RSSUN   ,RSSHA ,BGAP   ,WGAP, GAP,TGV,TGB,&
    845. 

  845.                  Q1     ,Q2V    ,Q2B    ,Q2E    ,CHV   ,CHB     , & !out
    846. 

  846.                  EMISSI)                                            !out
    847. 

  847. !jref:end                            
    848. 

  848. 

  849.     SICE(:) = MAX(0.0, SMC(:) - SH2O(:))   
    850. 

  850.     SNEQVO  = SNEQV
    851. 

  851. 

  852.     QVAP = MAX( FGEV/LATHEA, 0.)       ! positive part of fgev
    853. 

  853.     QDEW = ABS( MIN(FGEV/LATHEA, 0.))  ! negative part of fgev
    854. 

  854.     EDIR = QVAP - QDEW
    855. 

  855. 

  856. ! compute water budgets (water storages, ET components, and runoff)
    857. 

  857. 

  858.      CALL WATER (VEGTYP ,NSNOW  ,NSOIL  ,IMELT  ,DT     ,UU     , & !in
    859. 

  859.                  VV     ,FCEV   ,FCTR   ,QPRECC ,QPRECL ,ELAI   , & !in
    860. 

  860.                  ESAI   ,SFCTMP ,QVAP   ,QDEW   ,ZSOIL  ,BTRANI , & !in
    861. 

  861.                  FICEOLD,PONDING,TG     ,IST    ,FVEG   ,iloc,jloc , & !in
    862. 

  862.                  ISNOW  ,CANLIQ ,CANICE ,TV     ,SNOWH  ,SNEQV  , & !inout
    863. 

  863.                  SNICE  ,SNLIQ  ,STC    ,ZSNSO  ,SH2O   ,SMC    , & !inout
    864. 

  864.                  SICE   ,ZWT    ,WA     ,WT     ,DZSNSO ,WSLAKE , & !inout
    865. 

  865.                  CMC    ,ECAN   ,ETRAN  ,FWET   ,RUNSRF ,RUNSUB , & !out
    866. 

  866.                  QIN    ,QDIS   ,QSNOW  ,PONDING1       ,PONDING2,&
    867. 

  867.                  ISURBAN)  !out
    868. 

  868. 

  869. !     write(*,'(a20,10F15.5)') 'SFLX:RUNOFF=',RUNSRF*DT,RUNSUB*DT,EDIR*DT
    870. 

  870. 

  871. ! compute carbon budgets (carbon storages and co2 & bvoc fluxes)
    872. 

  872. 

  873.    IF (DVEG == 2) THEN
    874. 

  874.     CALL CARBON (NSNOW  ,NSOIL  ,VEGTYP ,NROOT  ,DT     ,ZSOIL  , & !in
    875. 

  875.                  DZSNSO ,STC    ,SMC    ,TV     ,TG     ,PSN    , & !in
    876. 

  876.                  FOLN   ,SMCMAX ,BTRAN  ,APAR   ,FVEG   ,IGS    , & !in
    877. 

  877.                  TROOT  ,IST    ,LAT    ,iloc   ,jloc   ,         & !in
    878. 

  878.                  LFMASS ,RTMASS ,STMASS ,WOOD   ,STBLCP ,FASTCP , & !inout
    879. 

  879.                  GPP    ,NPP    ,NEE    ,AUTORS ,HETERS ,TOTSC  , & !out
    880. 

  880.                  TOTLB  ,LAI    ,SAI    )                   !out
    881. 

  881.    END IF
    882. 

  882. 

  883. ! water and energy balance check
    884. 

  884. 

  885.      CALL ERROR (SWDOWN ,FSA    ,FSR    ,FIRA   ,FSH    ,FCEV   , & !in
    886. 

  886.                  FGEV   ,FCTR   ,SSOIL  ,BEG_WB ,CANLIQ ,CANICE , & !in
    887. 

  887.                  SNEQV  ,WA     ,SMC    ,DZSNSO ,PRCP   ,ECAN   , & !in
    888. 

  888.                  ETRAN  ,EDIR   ,RUNSRF ,RUNSUB ,DT     ,NSOIL  , & !in
    889. 

  889.                  NSNOW  ,IST    ,ERRWAT ,ILOC   , JLOC  ,FVEG   , &
    890. 

  890.                  SAV    ,SAG    ,FSRV   ,FSRG)   !in ( Except ERRWAT, which is out )
    891. 

  891. 

  892. ! urban - jref
    893. 

  893.     QFX = ETRAN + ECAN + EDIR
    894. 

  894.     IF ( VEGTYP == ISURBAN ) THEN
    895. 

  895.        QSFC = (QFX/RHOAIR*CHSTAR) + QAIR
    896. 

  896.        Q2B = QSFC
    897. 

  897.     END IF
    898. 

  898. 

  899.     IF(SWDOWN.NE.0.) THEN
    900. 

  900.       ALBEDO = FSR / SWDOWN
    901. 

  901.     ELSE
    902. 

  902.       ALBEDO = -999.9
    903. 

  903.     END IF
    904. 

  904.     
    905. 

  905. 

  906.   END SUBROUTINE NOAHMP_SFLX
    907. 

  907. ! ==================================================================================================
    908. 

  908.   SUBROUTINE ATM (SFCPRS ,SFCTMP ,Q2     ,PRCP   ,SOLDN  ,COSZ   ,THAIR  , &
    909. 

  909.                   QAIR   ,EAIR   ,RHOAIR ,QPRECC ,QPRECL ,SOLAD  ,SOLAI  , &
    910. 

  910.                   SWDOWN )     
    911. 

  911. ! --------------------------------------------------------------------------------------------------
    912. 

  912. ! re-process atmospheric forcing
    913. 

  913. ! --------------------------------------------------------------------------------------------------
    914. 

  914.   IMPLICIT NONE
    915. 

  915. ! --------------------------------------------------------------------------------------------------
    916. 

  916. ! inputs
    917. 

  917. 

  918.   REAL                          , INTENT(IN)  :: SFCPRS !pressure (pa)
    919. 

  919.   REAL                          , INTENT(IN)  :: SFCTMP !surface air temperature [k]
    920. 

  920.   REAL                          , INTENT(IN)  :: Q2     !mixing ratio (kg/kg)
    921. 

  921.   REAL                          , INTENT(IN)  :: SOLDN  !downward shortwave radiation (w/m2)
    922. 

  922.   REAL                          , INTENT(IN)  :: PRCP   !precipitation rate (kg m-2 s-1)
    923. 

  923.   REAL                          , INTENT(IN)  :: COSZ   !cosine solar zenith angle [0-1]
    924. 

  924. 

  925. ! outputs
    926. 

  926. 

  927.   REAL                          , INTENT(OUT) :: THAIR  !potential temperature (k)
    928. 

  928.   REAL                          , INTENT(OUT) :: QAIR   !specific humidity (kg/kg) (q2/(1+q2))
    929. 

  929.   REAL                          , INTENT(OUT) :: EAIR   !vapor pressure air (pa)
    930. 

  930.   REAL, DIMENSION(       1:   2), INTENT(OUT) :: SOLAD  !incoming direct solar radiation (w/m2)
    931. 

  931.   REAL, DIMENSION(       1:   2), INTENT(OUT) :: SOLAI  !incoming diffuse solar radiation (w/m2)
    932. 

  932.   REAL                          , INTENT(OUT) :: QPRECC !convective precipitation (mm/s)
    933. 

  933.   REAL                          , INTENT(OUT) :: QPRECL !large-scale precipitation (mm/s)
    934. 

  934.   REAL                          , INTENT(OUT) :: RHOAIR !density air (kg/m3)
    935. 

  935.   REAL                          , INTENT(OUT) :: SWDOWN !downward solar filtered by sun angle [w/m2]
    936. 

  936. 

  937. !locals
    938. 

  938. 

  939.   REAL                                        :: PAIR   !atm bottom level pressure (pa)
    940. 

  940. ! --------------------------------------------------------------------------------------------------
    941. 

  941. 

  942. !jref: seems like PAIR should be P1000mb??
    943. 

  943.        PAIR   = SFCPRS                   ! atm bottom level pressure (pa)
    944. 

  944.        THAIR  = SFCTMP * (SFCPRS/PAIR)**(RAIR/CPAIR) 
    945. 

  945. !jref: mixing ratio to specific       
    946. 

  946.        QAIR   = Q2 / (1.0+Q2)           ! mixing ratio to specific humidity [kg/kg]
    947. 

  947. !       QAIR   = Q2                       ! GLDAS forcing: Q2 = specific humidity [kg/kg]
    948. 

  948. 

  949.        EAIR   = QAIR*SFCPRS / (0.622+0.378*QAIR)
    950. 

  950.        RHOAIR = (SFCPRS-0.378*EAIR) / (RAIR*SFCTMP)
    951. 

  951. 

  952.        QPRECC = 0.10 * PRCP          ! should be from the atmospheric model
    953. 

  953.        QPRECL = 0.90 * PRCP          ! should be from the atmospheric model
    954. 

  954. 

  955.        IF(COSZ <= 0.) THEN 
    956. 

  956.           SWDOWN = 0.
    957. 

  957.        ELSE
    958. 

  958.           SWDOWN = SOLDN
    959. 

  959.        END IF 
    960. 

  960. 

  961.        SOLAD(1) = SWDOWN*0.7*0.5     ! direct  vis
    962. 

  962.        SOLAD(2) = SWDOWN*0.7*0.5     ! direct  nir
    963. 

  963.        SOLAI(1) = SWDOWN*0.3*0.5     ! diffuse vis
    964. 

  964.        SOLAI(2) = SWDOWN*0.3*0.5     ! diffuse nir
    965. 

  965. 

  966.   END SUBROUTINE ATM
    967. 

  967. ! ==================================================================================================
    968. 

  968. ! --------------------------------------------------------------------------------------------------
    969. 

  969.   SUBROUTINE PHENOLOGY (VEGTYP , SNOWH  , TV     , LAT   , YEARLEN , JULIAN , & !in
    970. 

  970.                         LAI    , SAI    , TROOT  , HTOP  , ELAI    , ESAI   , IGS)
    971. 

  971. 

  972. ! --------------------------------------------------------------------------------------------------
    973. 

  973. ! vegetation phenology considering vegeation canopy being buries by snow and evolution in time
    974. 

  974. ! --------------------------------------------------------------------------------------------------
    975. 

  975.   USE NOAHMP_VEG_PARAMETERS
    976. 

  976. ! --------------------------------------------------------------------------------------------------
    977. 

  977.   IMPLICIT NONE
    978. 

  978. ! --------------------------------------------------------------------------------------------------
    979. 

  979. ! inputs
    980. 

  980.   INTEGER                , INTENT(IN   ) :: VEGTYP !vegetation type 
    981. 

  981.   REAL                   , INTENT(IN   ) :: SNOWH  !snow height [m]
    982. 

  982.   REAL                   , INTENT(IN   ) :: TV  …
    983.

Large files files are truncated, but you can click here to view the full file