Skip to content

Noah mp snow #1403

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
davegill merged 11 commits into from
Feb 23, 2021
Merged

Noah mp snow #1403

Show file tree
Hide file tree
Changes from 10 commits
Commits
Show all changes
11 commits
Select commit Hold shift + click to select a range
45d1380
Merge pull request #1 from wrf-model/develop
cenlinhe Dec 17, 2020
4329b9f
Add snow retention at the snowpack bottom to be consistent with NWM u…
cenlinhe Dec 17, 2020
575a450
Modify CWPVT parameter values to be vegetation dependent in MPTABLE b…
cenlinhe Dec 17, 2020
fb13ed8
Bring hard-coded snow emissivity of 1.0 to a global tunable MPTABLE p…
cenlinhe Dec 18, 2020
d883131
Bring hard-coded 2.5*z0 (z0=0.002m) constant in snow cover formulatio…
cenlinhe Dec 18, 2020
de1e035
Update MFSNO with optimized vegetation-dependent values based on SNOT…
cenlinhe Dec 18, 2020
7f799a6
Add wet-bulb temperature snow-rain partitioning scheme for OPT_SNF=5b…
cenlinhe Dec 18, 2020
004eeda
Fix RB issue of value range in RAGRB for leaf resistence, now using m…
cenlinhe Dec 23, 2020
8ff1bd3
Some formatting changes
cenlinhe Dec 23, 2020
9b09e33
Bug Fix for the maximum SWE threshold in noahmpdrv file to be consist…
cenlinhe Feb 15, 2021
143e43c
Bug fix for SCFFAC values in MPTABLE & some format changes
cenlinhe Feb 22, 2021
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
13 changes: 8 additions & 5 deletions phys/module_sf_noahmpdrv.F
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1086,6 +1086,7 @@ SUBROUTINE TRANSFER_MP_PARAMETERS(VEGTYPE,SOILTYPE,SLOPETYPE,SOILCOLOR,CROPTYPE,
parameters%DEN = DEN_TABLE(VEGTYPE) !tree density (no. of trunks per m2)
parameters%RC = RC_TABLE(VEGTYPE) !tree crown radius (m)
parameters%MFSNO = MFSNO_TABLE(VEGTYPE) !snowmelt m parameter ()
parameters%SCFFAC = SCFFAC_TABLE(VEGTYPE) !snow cover factor (m) (originally hard-coded 2.5*z0 in SCF formulation)
parameters%SAIM = SAIM_TABLE(VEGTYPE,:) !monthly stem area index, one-sided
parameters%LAIM = LAIM_TABLE(VEGTYPE,:) !monthly leaf area index, one-sided
parameters%SLA = SLA_TABLE(VEGTYPE) !single-side leaf area per Kg [m2/kg]
Expand Down Expand Up @@ -1200,7 +1201,9 @@ SUBROUTINE TRANSFER_MP_PARAMETERS(VEGTYPE,SOILTYPE,SLOPETYPE,SOILCOLOR,CROPTYPE,
parameters%FSATMX = FSATMX_TABLE
parameters%Z0SNO = Z0SNO_TABLE
parameters%SSI = SSI_TABLE
parameters%SWEMX = SWEMX_TABLE
parameters%SNOW_RET_FAC = SNOW_RET_FAC_TABLE
parameters%SNOW_EMIS = SNOW_EMIS_TABLE
parameters%SWEMX = SWEMX_TABLE
parameters%TAU0 = TAU0_TABLE
parameters%GRAIN_GROWTH = GRAIN_GROWTH_TABLE
parameters%EXTRA_GROWTH = EXTRA_GROWTH_TABLE
Expand Down Expand Up @@ -1559,7 +1562,7 @@ SUBROUTINE NOAHMP_INIT ( MMINLU, SNOW , SNOWH , CANWAT , ISLTYP , IVGTYP, XLAT
ENDIF


! Check if snow/snowh are consistent and cap SWE at 2000mm;
! Check if snow/snowh are consistent and cap SWE at 5000mm;
! the Noah-MP code does it internally but if we don't do it here, problems ensue
DO J = jts,jtf
DO I = its,itf
Expand All @@ -1568,9 +1571,9 @@ SUBROUTINE NOAHMP_INIT ( MMINLU, SNOW , SNOWH , CANWAT , ISLTYP , IVGTYP, XLAT
,i,j,snow(i,j),snowh(i,j)
CALL wrf_message(err_message)
ENDIF
IF ( SNOW( i,j ) > 2000. ) THEN
SNOWH(I,J) = SNOWH(I,J) * 2000. / SNOW(I,J) ! SNOW in mm and SNOWH in m
SNOW (I,J) = 2000. ! cap SNOW at 2000, maintain density
IF ( SNOW( i,j ) > 5000. ) THEN
SNOWH(I,J) = SNOWH(I,J) * 5000. / SNOW(I,J) ! SNOW in mm and SNOWH in m
SNOW (I,J) = 5000. ! cap SNOW at 5000, maintain density
ENDIF
ENDDO
ENDDO
Expand Down
82 changes: 70 additions & 12 deletions phys/module_sf_noahmplsm.F
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -142,7 +142,8 @@ MODULE MODULE_SF_NOAHMPLSM
! 2 -> BATS: when SFCTMP<TFRZ+2.2
! 3 -> SFCTMP < TFRZ
! 4 -> Use WRF microphysics output

! 5 -> Use wetbulb temperature (Wang et al., 2019 GRL) C.He, 12/18/2020
Comment thread
davegill marked this conversation as resolved.
Outdated

INTEGER :: OPT_TBOT ! options for lower boundary condition of soil temperature
! 1 -> zero heat flux from bottom (ZBOT and TBOT not used)
! **2 -> TBOT at ZBOT (8m) read from a file (original Noah)
Expand Down Expand Up @@ -219,6 +220,7 @@ MODULE MODULE_SF_NOAHMPLSM
REAL :: DEN !tree density (no. of trunks per m2)
REAL :: RC !tree crown radius (m)
REAL :: MFSNO !snowmelt m parameter ()
REAL :: SCFFAC !snow cover factor (m) (originally hard-coded 2.5*z0 in SCF formulation)
REAL :: SAIM(12) !monthly stem area index, one-sided
REAL :: LAIM(12) !monthly leaf area index, one-sided
REAL :: SLA !single-side leaf area per Kg [m2/kg]
Expand Down Expand Up @@ -291,6 +293,8 @@ MODULE MODULE_SF_NOAHMPLSM
REAL :: FSATMX !maximum surface saturated fraction (global mean)
REAL :: Z0SNO !snow surface roughness length (m) (0.002)
REAL :: SSI !liquid water holding capacity for snowpack (m3/m3)
REAL :: SNOW_RET_FAC !snowpack water release timescale factor (1/s)
REAL :: SNOW_EMIS !snow emissivity
REAL :: SWEMX !new snow mass to fully cover old snow (mm)
REAL :: TAU0 !tau0 from Yang97 eqn. 10a
REAL :: GRAIN_GROWTH !growth from vapor diffusion Yang97 eqn. 10b
Expand Down Expand Up @@ -925,6 +929,15 @@ SUBROUTINE ATM (parameters,SFCPRS ,SFCTMP ,Q2 ,
REAL :: PRCP_FROZEN !total frozen precipitation [mm/s] ! MB/AN : v3.7
REAL, PARAMETER :: RHO_GRPL = 500.0 ! graupel bulk density [kg/m3] ! MB/AN : v3.7
REAL, PARAMETER :: RHO_HAIL = 917.0 ! hail bulk density [kg/m3] ! MB/AN : v3.7
! wet-bulb scheme Wang et al., 2019 GRL, C.He, 12/18/2020
REAL :: ESATAIR ! saturated vapor pressure of air
REAL :: LATHEA ! latent heat of vapor/sublimation
REAL :: GAMMA_b ! (cp*p)/(eps*L)
REAL :: TDC ! air temperature [C]
REAL :: TWET ! wetbulb temperature
INTEGER :: ITER
INTEGER, PARAMETER :: NITER = 10 ! iterations for Twet calculation

! --------------------------------------------------------------------------------------------------

!jref: seems like PAIR should be P1000mb??
Expand Down Expand Up @@ -1014,6 +1027,23 @@ SUBROUTINE ATM (parameters,SFCPRS ,SFCTMP ,Q2 ,

ENDIF

! wet-bulb scheme (Wang et al., 2019 GRL), C.He, 12/18/2020
IF(OPT_SNF == 5) THEN
TDC = MIN( 50., MAX(-50.,(SFCTMP-TFRZ)) ) !Kelvin to degree Celsius with limit -50 to +50
IF (SFCTMP > TFRZ) THEN
LATHEA = HVAP
ELSE
LATHEA = HSUB
END IF
GAMMA_b = CPAIR*SFCPRS/(0.622*LATHEA)
TWET = TDC - 5. ! first guess wetbulb temperature
DO ITER = 1, NITER
ESATAIR = 610.8 * EXP((17.27*TWET)/(237.3+TWET))
TWET = TWET - (ESATAIR-EAIR)/ GAMMA_b ! Wang et al., 2019 GRL Eq.2
END DO
FPICE = 1.0/(1.0+6.99E-5*exp(2.0*(TWET+3.97))) ! Wang et al., 2019 GRL Eq. 1
ENDIF

RAIN = PRCP * (1.-FPICE)
SNOW = PRCP * FPICE

Expand Down Expand Up @@ -1831,7 +1861,8 @@ SUBROUTINE ENERGY (parameters,ICE ,VEGTYP ,IST ,NSNOW ,NSOIL , & !in
IF(SNOWH.GT.0.) THEN
BDSNO = SNEQV / SNOWH
FMELT = (BDSNO/100.)**parameters%MFSNO
FSNO = TANH( SNOWH /(2.5* Z0 * FMELT))
!FSNO = TANH( SNOWH /(2.5* Z0 * FMELT))
Comment thread
davegill marked this conversation as resolved.
FSNO = TANH( SNOWH /(parameters%SCFFAC * FMELT)) ! C.He: bring hard-coded 2.5*z0 to MPTABLE tunable parameter SCFFAC
ENDIF

! ground roughness length
Expand Down Expand Up @@ -1900,9 +1931,9 @@ SUBROUTINE ENERGY (parameters,ICE ,VEGTYP ,IST ,NSNOW ,NSOIL , & !in

EMV = 1. - EXP(-(ELAI+ESAI)/1.0)
IF (ICE == 1) THEN
EMG = 0.98*(1.-FSNO) + 1.0*FSNO
EMG = 0.98*(1.-FSNO) + parameters%SNOW_EMIS*FSNO ! move hard-coded snow emissivity as a global parameter to MPTABLE
ELSE
EMG = parameters%EG(IST)*(1.-FSNO) + 1.0*FSNO
EMG = parameters%EG(IST)*(1.-FSNO) + parameters%SNOW_EMIS*FSNO
END IF

! soil moisture factor controlling stomatal resistance
Expand Down Expand Up @@ -4333,7 +4364,8 @@ SUBROUTINE RAGRB(parameters,ITER ,VAI ,RHOAIR ,HG ,TAH , & !in

TMPRB = CWPC*50. / (1. - EXP(-CWPC/2.))
RB = TMPRB * SQRT(parameters%DLEAF/UC)
RB = MAX(RB,100.0)
RB = MIN(MAX(RB, 5.0),50.0) ! limit RB to 5~50
! RB = MAX(RB,100.0) ! typically RB<50,so this leads to constant RB of 100
! RB = 200
Comment thread
davegill marked this conversation as resolved.
Outdated

END SUBROUTINE RAGRB
Expand Down Expand Up @@ -5494,6 +5526,7 @@ SUBROUTINE PHASECHANGE (parameters,NSNOW ,NSOIL ,ISNOW ,DT ,FACT ,
SNEQV = MAX(0.,TEMP1-XM(1))
PROPOR = SNEQV/TEMP1
SNOWH = MAX(0.,PROPOR * SNOWH)
SNOWH = MIN(MAX(SNOWH,SNEQV/500.0),SNEQV/50.0) ! limit adjustment to a reasonable density
HEATR = HM(1) - HFUS*(TEMP1-SNEQV)/DT
IF (HEATR > 0.) THEN
XM(1) = HEATR*DT/HFUS
Expand Down Expand Up @@ -6330,10 +6363,12 @@ SUBROUTINE COMBINE (parameters,NSNOW ,NSOIL ,ILOC ,JLOC , & !in
IF(J /= 0) THEN
SNLIQ(J+1) = SNLIQ(J+1) + SNLIQ(J)
SNICE(J+1) = SNICE(J+1) + SNICE(J)
DZSNSO(J+1) = DZSNSO(J+1) + DZSNSO(J)
ELSE
IF (ISNOW_OLD < -1) THEN ! MB/KM: change to ISNOW
SNLIQ(J-1) = SNLIQ(J-1) + SNLIQ(J)
SNICE(J-1) = SNICE(J-1) + SNICE(J)
DZSNSO(J-1) = DZSNSO(J-1) + DZSNSO(J)
ELSE
IF(SNICE(J) >= 0.) THEN
PONDING1 = SNLIQ(J) ! ISNOW WILL GET SET TO ZERO BELOW; PONDING1 WILL GET
Expand Down Expand Up @@ -6741,6 +6776,7 @@ SUBROUTINE COMPACT (parameters,NSNOW ,NSOIL ,DT ,STC ,SNICE , & !in
! The change in DZ due to compaction

DZSNSO(J) = DZSNSO(J)*(1.+PDZDTC)
DZSNSO(J) = max(DZSNSO(J),SNICE(J)/DENICE + SNLIQ(J)/DENH2O)
END IF

! Pressure of overlying snow
Expand Down Expand Up @@ -6803,6 +6839,7 @@ SUBROUTINE SNOWH2O (parameters,NSNOW ,NSOIL ,DT ,QSNFRO ,QSNSUB , & !in
REAL, DIMENSION(-NSNOW+1:0) :: EPORE !effective porosity = porosity - VOL_ICE
REAL :: PROPOR, TEMP
REAL :: PONDING1, PONDING2
REAL, PARAMETER :: max_liq_mass_fraction = 0.4
! ----------------------------------------------------------------------

!for the case when SNEQV becomes '0' after 'COMBINE'
Expand All @@ -6825,6 +6862,7 @@ SUBROUTINE SNOWH2O (parameters,NSNOW ,NSOIL ,DT ,QSNFRO ,QSNSUB , & !in
SNEQV = SNEQV - QSNSUB*DT + QSNFRO*DT
PROPOR = SNEQV/TEMP
SNOWH = MAX(0.,PROPOR * SNOWH)
SNOWH = MIN(MAX(SNOWH,SNEQV/500.0),SNEQV/50.0) ! limit adjustment to a reasonable density

IF(SNEQV < 0.) THEN
SICE(1) = SICE(1) + SNEQV/(DZSNSO(1)*1000.)
Expand Down Expand Up @@ -6876,11 +6914,22 @@ SUBROUTINE SNOWH2O (parameters,NSNOW ,NSOIL ,DT ,QSNFRO ,QSNSUB , & !in
SNLIQ(J) = SNLIQ(J) + QIN
VOL_LIQ(J) = SNLIQ(J)/(DZSNSO(J)*DENH2O)
QOUT = MAX(0.,(VOL_LIQ(J)-parameters%SSI*EPORE(J))*DZSNSO(J))
IF(J == 0) THEN
QOUT = MAX((VOL_LIQ(J)- EPORE(J))*DZSNSO(J) , parameters%SNOW_RET_FAC*DT*QOUT)
END IF
QOUT = QOUT*DENH2O
SNLIQ(J) = SNLIQ(J) - QOUT
IF((SNLIQ(J)/(SNICE(J)+SNLIQ(J))) > max_liq_mass_fraction) THEN
QOUT = QOUT + (SNLIQ(J) - max_liq_mass_fraction/(1.0 - max_liq_mass_fraction)*SNICE(J))
SNLIQ(J) = max_liq_mass_fraction/(1.0 - max_liq_mass_fraction)*SNICE(J)
ENDIF
QIN = QOUT
END DO

DO J = ISNOW+1, 0
DZSNSO(J) = MAX(DZSNSO(J),SNLIQ(J)/DENH2O + SNICE(J)/DENICE)
END DO

! Liquid water from snow bottom to soil

QSNBOT = QOUT / DT ! mm/s
Expand Down Expand Up @@ -9119,7 +9168,7 @@ subroutine noahmp_options(idveg ,iopt_crs ,iopt_btr ,iopt_run ,iopt_sfc
! 1 -> semi-implicit; 2 -> full implicit (original Noah)
INTEGER, INTENT(IN) :: iopt_rsf !surface resistance (1->Sakaguchi/Zeng; 2->Seller; 3->mod Sellers; 4->1+snow)
INTEGER, INTENT(IN) :: iopt_soil !soil parameters set-up option
INTEGER, INTENT(IN) :: iopt_pedo !pedo-transfer function
INTEGER, INTENT(IN) :: iopt_pedo !pedo-transfer function (1->Saxton and Rawls)
INTEGER, INTENT(IN) :: iopt_crop !crop model option (0->none; 1->Liu et al.; 2->Gecros)

! -------------------------------------------------------------------------------------------------
Expand Down Expand Up @@ -9178,6 +9227,7 @@ MODULE NOAHMP_TABLES
REAL :: DEN_TABLE(MVT) !tree density (no. of trunks per m2)
REAL :: RC_TABLE(MVT) !tree crown radius (m)
REAL :: MFSNO_TABLE(MVT) !snowmelt curve parameter ()
REAL :: SCFFAC_TABLE(MVT) !snow cover factor (m) (replace original hard-coded 2.5*z0 in SCF formulation)
REAL :: SAIM_TABLE(MVT,12) !monthly stem area index, one-sided
REAL :: LAIM_TABLE(MVT,12) !monthly leaf area index, one-sided
REAL :: SLA_TABLE(MVT) !single-side leaf area per Kg [m2/kg]
Expand Down Expand Up @@ -9269,6 +9319,8 @@ MODULE NOAHMP_TABLES
REAL :: FSATMX_TABLE !maximum surface saturated fraction (global mean)
REAL :: Z0SNO_TABLE !snow surface roughness length (m) (0.002)
REAL :: SSI_TABLE !liquid water holding capacity for snowpack (m3/m3) (0.03)
REAL :: SNOW_RET_FAC_TABLE !snowpack water release timescale factor (1/s)
REAL :: SNOW_EMIS_TABLE!snow emissivity
REAL :: SWEMX_TABLE !new snow mass to fully cover old snow (mm)
REAL :: TAU0_TABLE !tau0 from Yang97 eqn. 10a
REAL :: GRAIN_GROWTH_TABLE !growth from vapor diffusion Yang97 eqn. 10b
Expand Down Expand Up @@ -9412,15 +9464,15 @@ subroutine read_mp_veg_parameters(DATASET_IDENTIFIER)
LAI_JUL,LAI_AUG,LAI_SEP,LAI_OCT,LAI_NOV,LAI_DEC
REAL, DIMENSION(MVT) :: RHOL_VIS, RHOL_NIR, RHOS_VIS, RHOS_NIR, &
TAUL_VIS, TAUL_NIR, TAUS_VIS, TAUS_NIR
REAL, DIMENSION(MVT) :: CH2OP, DLEAF, Z0MVT, HVT, HVB, DEN, RC, MFSNO, XL, CWPVT, C3PSN, KC25, AKC, KO25, AKO, &
REAL, DIMENSION(MVT) :: CH2OP, DLEAF, Z0MVT, HVT, HVB, DEN, RC, MFSNO, SCFFAC, XL, CWPVT, C3PSN, KC25, AKC, KO25, AKO, &
AVCMX, AQE, LTOVRC, DILEFC, DILEFW, RMF25 , SLA , FRAGR , TMIN , VCMX25, TDLEF , &
BP, MP, QE25, RMS25, RMR25, ARM, FOLNMX, WDPOOL, WRRAT, MRP, NROOT, RGL, RS, HS, TOPT, RSMAX, &
SLAREA, EPS1, EPS2, EPS3, EPS4, EPS5

NAMELIST / noahmp_usgs_veg_categories / VEG_DATASET_DESCRIPTION, NVEG
NAMELIST / noahmp_usgs_parameters / ISURBAN, ISWATER, ISBARREN, ISICE, ISCROP, EBLFOREST, NATURAL, &
LOW_DENSITY_RESIDENTIAL, HIGH_DENSITY_RESIDENTIAL, HIGH_INTENSITY_INDUSTRIAL, &
CH2OP, DLEAF, Z0MVT, HVT, HVB, DEN, RC, MFSNO, XL, CWPVT, C3PSN, KC25, AKC, KO25, AKO, AVCMX, AQE, &
CH2OP, DLEAF, Z0MVT, HVT, HVB, DEN, RC, MFSNO, SCFFAC, XL, CWPVT, C3PSN, KC25, AKC, KO25, AKO, AVCMX, AQE, &
LTOVRC, DILEFC, DILEFW, RMF25 , SLA , FRAGR , TMIN , VCMX25, TDLEF , BP, MP, QE25, RMS25, RMR25, ARM, &
FOLNMX, WDPOOL, WRRAT, MRP, NROOT, RGL, RS, HS, TOPT, RSMAX, &
SAI_JAN, SAI_FEB, SAI_MAR, SAI_APR, SAI_MAY, SAI_JUN,SAI_JUL,SAI_AUG,SAI_SEP,SAI_OCT,SAI_NOV,SAI_DEC, &
Expand All @@ -9430,7 +9482,7 @@ subroutine read_mp_veg_parameters(DATASET_IDENTIFIER)
NAMELIST / noahmp_modis_veg_categories / VEG_DATASET_DESCRIPTION, NVEG
NAMELIST / noahmp_modis_parameters / ISURBAN, ISWATER, ISBARREN, ISICE, ISCROP, EBLFOREST, NATURAL, &
LOW_DENSITY_RESIDENTIAL, HIGH_DENSITY_RESIDENTIAL, HIGH_INTENSITY_INDUSTRIAL, &
CH2OP, DLEAF, Z0MVT, HVT, HVB, DEN, RC, MFSNO, XL, CWPVT, C3PSN, KC25, AKC, KO25, AKO, AVCMX, AQE, &
CH2OP, DLEAF, Z0MVT, HVT, HVB, DEN, RC, MFSNO, SCFFAC, XL, CWPVT, C3PSN, KC25, AKC, KO25, AKO, AVCMX, AQE, &
LTOVRC, DILEFC, DILEFW, RMF25 , SLA , FRAGR , TMIN , VCMX25, TDLEF , BP, MP, QE25, RMS25, RMR25, ARM, &
FOLNMX, WDPOOL, WRRAT, MRP, NROOT, RGL, RS, HS, TOPT, RSMAX, &
SAI_JAN, SAI_FEB, SAI_MAR, SAI_APR, SAI_MAY, SAI_JUN,SAI_JUL,SAI_AUG,SAI_SEP,SAI_OCT,SAI_NOV,SAI_DEC, &
Expand All @@ -9446,6 +9498,7 @@ subroutine read_mp_veg_parameters(DATASET_IDENTIFIER)
DEN_TABLE = -1.E36
RC_TABLE = -1.E36
MFSNO_TABLE = -1.E36
SCFFAC_TABLE = -1.E36
RHOL_TABLE = -1.E36
RHOS_TABLE = -1.E36
TAUL_TABLE = -1.E36
Expand Down Expand Up @@ -9541,6 +9594,7 @@ subroutine read_mp_veg_parameters(DATASET_IDENTIFIER)
DEN_TABLE(1:NVEG) = DEN(1:NVEG)
RC_TABLE(1:NVEG) = RC(1:NVEG)
MFSNO_TABLE(1:NVEG) = MFSNO(1:NVEG)
SCFFAC_TABLE(1:NVEG) = SCFFAC(1:NVEG)
XL_TABLE(1:NVEG) = XL(1:NVEG)
CWPVT_TABLE(1:NVEG) = CWPVT(1:NVEG)
C3PSN_TABLE(1:NVEG) = C3PSN(1:NVEG)
Expand Down Expand Up @@ -9780,12 +9834,12 @@ subroutine read_mp_global_parameters()
integer :: ierr
logical :: file_named

REAL :: CO2,O2,TIMEAN,FSATMX,Z0SNO,SSI, &
REAL :: CO2,O2,TIMEAN,FSATMX,Z0SNO,SSI,SNOW_RET_FAC,SNOW_EMIS,&
SWEMX,TAU0,GRAIN_GROWTH,EXTRA_GROWTH,DIRT_SOOT,&
BATS_COSZ,BATS_VIS_NEW,BATS_NIR_NEW,BATS_VIS_AGE,BATS_NIR_AGE,BATS_VIS_DIR,BATS_NIR_DIR,&
RSURF_SNOW,RSURF_EXP

NAMELIST / noahmp_global_parameters / CO2,O2,TIMEAN,FSATMX,Z0SNO,SSI, &
NAMELIST / noahmp_global_parameters / CO2,O2,TIMEAN,FSATMX,Z0SNO,SSI,SNOW_RET_FAC,SNOW_EMIS,&
SWEMX,TAU0,GRAIN_GROWTH,EXTRA_GROWTH,DIRT_SOOT,&
BATS_COSZ,BATS_VIS_NEW,BATS_NIR_NEW,BATS_VIS_AGE,BATS_NIR_AGE,BATS_VIS_DIR,BATS_NIR_DIR,&
RSURF_SNOW,RSURF_EXP
Expand All @@ -9798,7 +9852,9 @@ subroutine read_mp_global_parameters()
FSATMX_TABLE = -1.E36
Z0SNO_TABLE = -1.E36
SSI_TABLE = -1.E36
SWEMX_TABLE = -1.E36
SNOW_RET_FAC_TABLE = -1.E36
SNOW_EMIS_TABLE = -1.E36
SWEMX_TABLE = -1.E36
TAU0_TABLE = -1.E36
GRAIN_GROWTH_TABLE = -1.E36
EXTRA_GROWTH_TABLE = -1.E36
Expand Down Expand Up @@ -9834,6 +9890,8 @@ subroutine read_mp_global_parameters()
FSATMX_TABLE = FSATMX
Z0SNO_TABLE = Z0SNO
SSI_TABLE = SSI
SNOW_RET_FAC_TABLE = SNOW_RET_FAC
SNOW_EMIS_TABLE = SNOW_EMIS
SWEMX_TABLE = SWEMX
TAU0_TABLE = TAU0
GRAIN_GROWTH_TABLE = GRAIN_GROWTH
Expand Down
Loading