From 8a61c8fbf3d12505f558611dcb1cc5e4e9337f74 Mon Sep 17 00:00:00 2001 From: Prasanth Valayamkunnath Date: Wed, 9 Jun 2021 19:41:25 +0000 Subject: [PATCH] Add Xinanjiang Surface Runoff Scheme (#559) * Xinanjiang runoff scheme added * spatial soil XAJ params added * change NWM tests to use new runoff scheme Co-authored-by: Katelyn FitzGerald --- .../IO_code/module_NoahMP_hrldas_driver.F | 12 ++ .../NoahMP/IO_code/module_hrldas_netcdf_io.F | 63 +++++++++- .../NoahMP/phys/module_sf_noahmpdrv.F | 13 +- .../NoahMP/phys/module_sf_noahmplsm.F | 111 ++++++++++++++++-- .../NDHMS/Land_models/NoahMP/run/SOILPARM.TBL | 80 ++++++------- trunk/NDHMS/hrldas_namelists.json | 10 +- 6 files changed, 234 insertions(+), 55 deletions(-) diff --git a/trunk/NDHMS/Land_models/NoahMP/IO_code/module_NoahMP_hrldas_driver.F b/trunk/NDHMS/Land_models/NoahMP/IO_code/module_NoahMP_hrldas_driver.F index b8f29a0f0..c3cbcd247 100644 --- a/trunk/NDHMS/Land_models/NoahMP/IO_code/module_NoahMP_hrldas_driver.F +++ b/trunk/NDHMS/Land_models/NoahMP/IO_code/module_NoahMP_hrldas_driver.F @@ -156,6 +156,9 @@ module module_NoahMP_hrldas_driver REAL, ALLOCATABLE, DIMENSION(:,:) :: hvt_2D ! Canopy Height REAL, ALLOCATABLE, DIMENSION(:,:) :: mfsno_2D ! Snow cover m parameter REAL, ALLOCATABLE, DIMENSION(:,:) :: rsurfexp_2D! exponent in the shape parameter for soil resistance option 1 + REAL, ALLOCATABLE, DIMENSION(:,:) :: axaj_2D ! Tension water distribution inflection parameter [-] + REAL, ALLOCATABLE, DIMENSION(:,:) :: bxaj_2D ! Tension water distribution shape parameter [-] + REAL, ALLOCATABLE, DIMENSION(:,:) :: xxaj_2D ! Free water distribution shape parameter [-] #endif ! INOUT (with generic LSM equivalent) (as defined in WRF) @@ -741,6 +744,9 @@ subroutine land_driver_ini(NTIME_out, state, wrfits,wrfite,wrfjts,wrfjte) ALLOCATE ( hvt_2D (XSTART:XEND,YSTART:YEND) ) ! Canopy Height ALLOCATE ( mfsno_2D (XSTART:XEND,YSTART:YEND) ) ! Snow cover m parameter ALLOCATE ( rsurfexp_2D(XSTART:XEND,YSTART:YEND) ) ! exponent in the shape parameter for soil resistance option 1 + ALLOCATE ( axaj_2D (XSTART:XEND,YSTART:YEND) ) ! Tension water distribution inflection parameter [-] + ALLOCATE ( bxaj_2D (XSTART:XEND,YSTART:YEND) ) ! Tension water distribution shape parameter [-] + ALLOCATE ( xxaj_2D (XSTART:XEND,YSTART:YEND) ) ! Free water distribution shape parameter [-] #endif ! INOUT (with generic LSM equivalent) (as defined in WRF) @@ -1087,6 +1093,11 @@ subroutine land_driver_ini(NTIME_out, state, wrfits,wrfite,wrfjts,wrfjte) NSOIL,BEXP_3D,SMCDRY_3D,SMCWLT_3D,SMCREF_3D,SMCMAX_3D, & DKSAT_3D,DWSAT_3D,PSISAT_3D,QUARTZ_3D,REFDK_2D,REFKDT_2D,& SLOPE_2D,CWPVT_2D,VCMX25_2D,MP_2D,HVT_2D,MFSNO_2D,RSURFEXP_2D) + + if (noah_lsm%runoff_option == 7) then + CALL READ_XAJ_RUNOFF(noah_lsm%SPATIAL_FILENAME,XSTART, XEND, YSTART, YEND, & + AXAJ_2D,BXAJ_2D,XXAJ_2D) + end if #endif !------------------------------------------------------------------------ @@ -1658,6 +1669,7 @@ subroutine land_driver_exe(itime, state) DKSAT_3D,DWSAT_3D,PSISAT_3D,QUARTZ_3D, & REFDK_2D,REFKDT_2D,SLOPE_2D, & CWPVT_2D,VCMX25_2D,MP_2D,HVT_2D,MFSNO_2D,RSURFEXP_2D, & + AXAJ_2D,BXAJ_2D,XXAJ_2D, & #endif #ifdef WRF_HYDRO sfcheadrt,INFXSRT,soldrain, & !O diff --git a/trunk/NDHMS/Land_models/NoahMP/IO_code/module_hrldas_netcdf_io.F b/trunk/NDHMS/Land_models/NoahMP/IO_code/module_hrldas_netcdf_io.F index 0e4d3b888..5ebdc7d1a 100644 --- a/trunk/NDHMS/Land_models/NoahMP/IO_code/module_hrldas_netcdf_io.F +++ b/trunk/NDHMS/Land_models/NoahMP/IO_code/module_hrldas_netcdf_io.F @@ -771,6 +771,67 @@ subroutine read_mmf_runoff(wrfinput_flnm, & end subroutine read_mmf_runoff +!------------------------------------------------------------------------------------------- + + subroutine read_xaj_runoff(spatial_filename,xstart, xend,ystart, yend, & + axaj_2d,bxaj_2d,xxaj_2d) + + implicit none + character(len=*), intent(in) :: spatial_filename + integer, intent(in) :: xstart, xend, ystart, yend + real, dimension(xstart:xend,ystart:yend), intent(out) :: axaj_2d + real, dimension(xstart:xend,ystart:yend), intent(out) :: bxaj_2d + real, dimension(xstart:xend,ystart:yend), intent(out) :: xxaj_2d + + character(len=24) :: name + character(len=256) :: units + integer :: ierr,iret, varid + integer :: ncid + real, dimension(xstart:xend,ystart:yend) :: xdum +!------------------------------------------------------------------------------------------------------ + + ierr = nf90_open(spatial_filename, NF90_NOWRITE, ncid) + if (ierr /= 0) then + write(*,'("FATAL ERROR: In read_xaj_runoff(): Problem opening soil file: ''", A, "''")') trim(spatial_filename) + stop + endif + + name = "AXAJ" + iret = nf90_inq_varid(ncid, trim(name), varid) + if (iret /= 0) then + print*, 'ncid = ', ncid + write(*,*) "FATAL ERROR: In read_xaj_runoff(): Problem finding variable '"//trim(name)//"' in NetCDF file: " // trim(spatial_filename) + stop + endif + + iret = nf90_get_var(ncid, varid, axaj_2d, start=(/xstart,ystart/), count=(/xend-xstart+1,yend-ystart+1/)) +! + name = "BXAJ" + iret = nf90_inq_varid(ncid, trim(name), varid) + if (iret /= 0) then + print*, 'ncid = ', ncid + write(*,*) "FATAL ERROR: In read_xaj_runoff(): Problem finding variable '"//trim(name)//"' in NetCDF file: " // trim(spatial_filename) + stop + endif + + iret = nf90_get_var(ncid, varid, bxaj_2d, start=(/xstart,ystart/), count=(/xend-xstart+1,yend-ystart+1/)) +! + name = "XXAJ" + iret = nf90_inq_varid(ncid, trim(name), varid) + if (iret /= 0) then + print*, 'ncid = ', ncid + write(*,*) "FATAL ERROR: In read_xaj_runoff(): Problem finding variable '"//trim(name)//"' in NetCDF file: " // trim(spatial_filename) + stop + endif + + iret = nf90_get_var(ncid, varid, xxaj_2d, start=(/xstart,ystart/), count=(/xend-xstart+1,yend-ystart+1/)) +! + ! Close the NetCDF file + ierr = nf90_close(ncid) + if (ierr /= 0) stop "FATAL ERROR: In module_hrldas_netcdf_io.F read_xaj_runoff() - NF90_CLOSE" + + end subroutine read_xaj_runoff + !--------------------------------------------------------------------------------------------------------- subroutine read_3d_soil(spatial_filename,xstart, xend,ystart, yend, & @@ -1027,7 +1088,7 @@ subroutine read_3d_soil(spatial_filename,xstart, xend,ystart, yend, & endif iret = nf90_get_var(ncid, varid, rsurfexp_2d, start=(/xstart,ystart/), count=(/xend-xstart+1,yend-ystart+1/)) - +! ! Close the NetCDF file ierr = nf90_close(ncid) if (ierr /= 0) stop "FATAL ERROR: In module_hrldas_netcdf_io.F read_3d_soil() - NF90_CLOSE" diff --git a/trunk/NDHMS/Land_models/NoahMP/phys/module_sf_noahmpdrv.F b/trunk/NDHMS/Land_models/NoahMP/phys/module_sf_noahmpdrv.F index 66f0e6c01..b16126c51 100644 --- a/trunk/NDHMS/Land_models/NoahMP/phys/module_sf_noahmpdrv.F +++ b/trunk/NDHMS/Land_models/NoahMP/phys/module_sf_noahmpdrv.F @@ -63,6 +63,7 @@ SUBROUTINE noahmplsm(ITIMESTEP, YR, JULIAN, COSZIN, XLATIN, & ! IN DKSAT_3D,DWSAT_3D,PSISAT_3D,QUARTZ_3D, & REFDK_2D,REFKDT_2D,SLOPE_2D, & CWPVT_2D,VCMX25_2D,MP_2D,HVT_2D,MFSNO_2D,RSURFEXP_2D, & + AXAJ_2D,BXAJ_2D,XXAJ_2D, & #endif #ifdef WRF_HYDRO sfcheadrt,INFXSRT,soldrain, & @@ -161,6 +162,9 @@ SUBROUTINE noahmplsm(ITIMESTEP, YR, JULIAN, COSZIN, XLATIN, & ! IN REAL, DIMENSION( ims:ime, jms:jme ), INTENT(IN) :: HVT_2D ! Canopy Height REAL, DIMENSION( ims:ime, jms:jme ), INTENT(IN) :: MFSNO_2D ! Snow cover m parameter REAL, DIMENSION( ims:ime, jms:jme ), INTENT(IN) :: RSURFEXP_2D ! exponent in the shape parameter for soil resistance option 1 + REAL, DIMENSION( ims:ime, jms:jme ), INTENT(IN) :: AXAJ_2D ! Xinanjiang: Tension water distribution inflection parameter [-] + REAL, DIMENSION( ims:ime, jms:jme ), INTENT(IN) :: BXAJ_2D ! Xinanjiang: Tension water distribution shape parameter [-] + REAL, DIMENSION( ims:ime, jms:jme ), INTENT(IN) :: XXAJ_2D ! Xinanjiang: Free water distribution shape parameter [-] #endif ! INOUT (with generic LSM equivalent) @@ -660,6 +664,9 @@ SUBROUTINE noahmplsm(ITIMESTEP, YR, JULIAN, COSZIN, XLATIN, & ! IN parameters%hvt = HVT_2D(I,J) ! Canopy Height parameters%mfsno = MFSNO_2D(I,J) ! Snow cover m parameter parameters%rsurf_exp = RSURFEXP_2D(I,J) ! exponent in the shape parameter for soil resistance option 1 + parameters%axaj = AXAJ_2D(I,J) ! Xinanjiang: Tension water distribution inflection parameter [-] + parameters%bxaj = BXAJ_2D(I,J) ! Xinanjiang: Tension water distribution shape parameter [-] + parameters%xxaj = XXAJ_2D(I,J) ! Xinanjiang: Free water distribution shape parameter [-] #endif CALL TRANSFER_MP_PARAMETERS(VEGTYP,SOILTYP,SLOPETYP,SOILCOLOR,CROPTYPE,parameters) @@ -1193,9 +1200,11 @@ SUBROUTINE TRANSFER_MP_PARAMETERS(VEGTYPE,SOILTYPE,SLOPETYPE,SOILCOLOR,CROPTYPE, parameters%MFSNO = MFSNO_TABLE(VEGTYPE) !snowmelt m parameter () parameters%HVT = HVT_TABLE(VEGTYPE) !top of canopy (m) parameters%MP = MP_TABLE(VEGTYPE) !slope of conductance-to-photosynthesis relationship - parameters%RSURF_EXP = RSURF_EXP_TABLE + parameters%RSURF_EXP = RSURF_EXP_TABLE + parameters%AXAJ = AXAJ_TABLE(SOILTYPE) + parameters%BXAJ = BXAJ_TABLE(SOILTYPE) + parameters%XXAJ = XXAJ_TABLE(SOILTYPE) #endif - ! ---------------------------------------------------------------------- ! Transfer GENPARM parameters ! ---------------------------------------------------------------------- diff --git a/trunk/NDHMS/Land_models/NoahMP/phys/module_sf_noahmplsm.F b/trunk/NDHMS/Land_models/NoahMP/phys/module_sf_noahmplsm.F index 7af8a4395..fdbef4879 100644 --- a/trunk/NDHMS/Land_models/NoahMP/phys/module_sf_noahmplsm.F +++ b/trunk/NDHMS/Land_models/NoahMP/phys/module_sf_noahmplsm.F @@ -361,6 +361,9 @@ MODULE MODULE_SF_NOAHMPLSM REAL :: DWSAT(NSOIL) !saturated soil hydraulic diffusivity REAL :: QUARTZ(NSOIL) !soil quartz content REAL :: F1 !soil thermal diffusivity/conductivity coef (not used MB: 20140718) + REAL :: AXAJ !Xinanjiang: Tension water distribution inflection parameter [-] + REAL :: BXAJ !Xinanjiang: Tension water distribution shape parameter [-] + REAL :: XXAJ !Xinanjiang: Free water distribution shape parameter [-] !------------------------------------------------------------------------------------------! ! From the GENPARM.TBL file !------------------------------------------------------------------------------------------! @@ -5841,7 +5844,7 @@ SUBROUTINE WATER (parameters,VEGTYP ,NSNOW ,NSOIL ,IMELT ,DT ,UU , & RUNSUB = QDIS !mm/s END IF - IF(OPT_RUN == 3 .or. OPT_RUN == 4) THEN + IF(OPT_RUN == 3 .or. OPT_RUN == 4 .or. OPT_RUN == 7) THEN RUNSUB = RUNSUB + QDRAIN !mm/s END IF @@ -6995,7 +6998,9 @@ SUBROUTINE SOILWATER (parameters,NSOIL ,NSNOW ,DT ,ZSOIL ,DZSNSO , & !in PDDUM = QINSUR - RUNSRF ! m/s END IF END IF - + IF (OPT_RUN == 7) THEN + CALL COMPUTE_XAJ_SURFRUNOFF(parameters,DT,FCR,NSOIL,SMC,ZSOIL,QINSUR,RUNSRF,PDDUM) + END IF ! determine iteration times and finer time step NITER = 1 @@ -7020,6 +7025,11 @@ SUBROUTINE SOILWATER (parameters,NSOIL ,NSNOW ,DT ,ZSOIL ,DZSNSO , & !in PDDUM ,RUNSRF ) !out END IF + IF (QINSUR > 0. .AND. OPT_RUN == 7) THEN + CALL COMPUTE_XAJ_SURFRUNOFF(parameters,DTFINE,FCR,NSOIL,SMC,ZSOIL,QINSUR,& ! in + RUNSRF,PDDUM) ! out + END IF + CALL SRT (parameters,NSOIL ,ZSOIL ,DTFINE ,PDDUM ,ETRANI , & !in QSEVA ,SH2O ,SMC ,ZWT ,FCR , & !in SICEMAX,FCRMAX ,ILOC ,JLOC ,SMCWTD , & !in @@ -7343,7 +7353,7 @@ SUBROUTINE SRT (parameters,NSOIL ,ZSOIL ,DT ,PDDUM ,ETRANI , & !in IF(OPT_RUN == 1 .or. OPT_RUN == 2) THEN QDRAIN = 0. END IF - IF(OPT_RUN == 3) THEN + IF(OPT_RUN == 3 .OR. OPT_RUN == 7) THEN QDRAIN = parameters%SLOPE*WCND(K) END IF IF(OPT_RUN == 4) THEN @@ -7511,6 +7521,82 @@ SUBROUTINE SSTEP (parameters,NSOIL ,NSNOW ,DT ,ZSOIL ,DZSNSO , & !in END SUBROUTINE SSTEP +!== begin xinanjiag================================================================================= + SUBROUTINE COMPUTE_XAJ_SURFRUNOFF(parameters,DT,FCR,NSOIL,SMC,ZSOIL,QINSUR,RUNSRF,PDDUM) +! ---------------------------------------------------------------------- +! Calculate the saturated area and runoff based on Xinanjiag runoff scheme. +! Reference: Knoben, W. J., Freer, J. E., Fowler, K. J., Peel, M. C., & Woods, R. A. (2019). +! Modular Assessment of Rainfall-Runoff Models Toolbox (MARRMoT) v1. 2: +! an open-source, extendable framework providing implementations of 46 conceptual +! hydrologic models as continuous state-space formulations. +! ---------------------------------------------------------------------- +! Author: Prasanth Valayamkunnath +! Date: August 03, 2020 +! ---------------------------------------------------------------------- + IMPLICIT NONE +! ---------------------------------------------------------------------- +! Inputs + TYPE (noahmp_parameters), INTENT(IN) :: parameters + INTEGER, INTENT(IN) :: NSOIL + REAL, DIMENSION(1:NSOIL), INTENT(IN) :: SMC + REAL, DIMENSION(1:NSOIL), INTENT(IN) :: ZSOIL + REAL, DIMENSION(1:NSOIL), INTENT(IN) :: FCR !fraction of imperviousness (-) = IMP + REAL , INTENT(IN) :: QINSUR + REAL , INTENT(IN) :: DT +! Output + REAL , INTENT(INOUT):: RUNSRF + REAL , INTENT(INOUT):: PDDUM +! local + REAL :: WM,WM_MAX,SM,SM_MAX,IRUNOFF,PRUNOFF + INTEGER :: IZ +!------------------------------------------------------------------------ +!initialize + WM = 0.0 + WM_MAX = 0.0 + SM = 0.0 + SM_MAX = 0.0 + IRUNOFF = 0.0 + PRUNOFF = 0.0 + RUNSRF = 0.0 + + DO IZ=1,NSOIL-2 + IF ((SMC(IZ)-parameters%SMCREF(IZ)) .GT. 0.) THEN ! soil moisture greater than field capacity + SM = SM + (SMC(IZ) - parameters%SMCREF(IZ) )*-1*ZSOIL(IZ) !m + WM = WM + (parameters%SMCREF(IZ)*-1*ZSOIL(IZ)) !m + ELSE + WM = WM + (SMC(IZ)*-1*ZSOIL(IZ)) + END IF + WM_MAX = WM_MAX + (parameters%SMCREF(IZ)*-1*ZSOIL(IZ)) + SM_MAX = SM_MAX + (parameters%SMCMAX(IZ) - parameters%SMCREF(IZ))*-1*ZSOIL(IZ) + END DO + + WM = MIN(WM,WM_MAX) ! tension water (m) + SM = MIN(SM,SM_MAX) ! free water (m) + +! impervious surface runoff R_IMP + IRUNOFF = FCR(1)*QINSUR*DT + +! solve pervious surface runoff (m) based on Eq. (310) + IF ((WM/WM_MAX) .LE. (0.5-parameters%AXAJ))THEN + PRUNOFF = (1-FCR(1))*QINSUR*DT*((0.5-parameters%AXAJ)**(1-parameters%BXAJ))*((WM/WM_MAX)**parameters%BXAJ) + ELSE + PRUNOFF = (1-FCR(1))*QINSUR*DT*(1-(((0.5+parameters%AXAJ)**(1-parameters%BXAJ))*((1-(WM/WM_MAX))**parameters%BXAJ))) + END IF + +! estimate surface runoff based on Eq. (313) + IF(QINSUR .EQ. 0.0) THEN + RUNSRF = 0.0 + ELSE + RUNSRF = PRUNOFF*(1-((1-(SM/SM_MAX))**parameters%XXAJ))+IRUNOFF + END IF + RUNSRF = RUNSRF/DT !m/s + RUNSRF = MAX(0.0, RUNSRF) + RUNSRF = MIN(QINSUR, RUNSRF) + PDDUM = QINSUR - RUNSRF + + END SUBROUTINE COMPUTE_XAJ_SURFRUNOFF +!== end xinanjiag ================================================================================== + !== begin wdfcnd1 ================================================================================== SUBROUTINE WDFCND1 (parameters,WDF,WCND,SMC,FCR,ISOIL) @@ -9117,17 +9203,19 @@ MODULE NOAHMP_TABLES INTEGER :: SLCATS - REAL :: BEXP_TABLE(MAX_SOILTYP) !maximum intercepted h2o per unit lai+sai (mm) + REAL :: BEXP_TABLE(MAX_SOILTYP) !maximum intercepted h2o per unit lai+sai (mm) REAL :: SMCDRY_TABLE(MAX_SOILTYP) !characteristic leaf dimension (m) - REAL :: F1_TABLE(MAX_SOILTYP) !momentum roughness length (m) + REAL :: F1_TABLE(MAX_SOILTYP) !momentum roughness length (m) REAL :: SMCMAX_TABLE(MAX_SOILTYP) !top of canopy (m) REAL :: SMCREF_TABLE(MAX_SOILTYP) !bottom of canopy (m) REAL :: PSISAT_TABLE(MAX_SOILTYP) !tree density (no. of trunks per m2) REAL :: DKSAT_TABLE(MAX_SOILTYP) !tree crown radius (m) REAL :: DWSAT_TABLE(MAX_SOILTYP) !monthly stem area index, one-sided REAL :: SMCWLT_TABLE(MAX_SOILTYP) !monthly leaf area index, one-sided - REAL :: QUARTZ_TABLE(MAX_SOILTYP) !single-side leaf area per Kg [m2/kg] - + REAL :: QUARTZ_TABLE(MAX_SOILTYP) !single-side leaf area per Kg [m2/kg] + REAL :: AXAJ_TABLE(MAX_SOILTYP) !Xinanjiang: Tension water distribution inflection parameter [-] + REAL :: BXAJ_TABLE(MAX_SOILTYP) !Xinanjiang: Tension water distribution shape parameter [-] + REAL :: XXAJ_TABLE(MAX_SOILTYP) !Xinanjiang: Free water distribution shape parameter [-] ! GENPARM.TBL parameters REAL :: SLOPE_TABLE(9) !slope factor for soil drainage @@ -9532,8 +9620,10 @@ subroutine read_mp_soil_parameters() FRZK_TABLE = -1.E36 ZBOT_TABLE = -1.E36 CZIL_TABLE = -1.E36 - -! + AXAJ_TABLE = -1.E36 + BXAJ_TABLE = -1.E36 + XXAJ_TABLE = -1.E36 +! !-----READ IN SOIL PROPERTIES FROM SOILPARM.TBL ! inquire( file='SOILPARM.TBL', exist=file_named ) @@ -9560,7 +9650,8 @@ subroutine read_mp_soil_parameters() DO LC=1,SLCATS READ (21,*) ITMP,BEXP_TABLE(LC),SMCDRY_TABLE(LC),F1_TABLE(LC),SMCMAX_TABLE(LC), & SMCREF_TABLE(LC),PSISAT_TABLE(LC),DKSAT_TABLE(LC), DWSAT_TABLE(LC), & - SMCWLT_TABLE(LC), QUARTZ_TABLE(LC) + SMCWLT_TABLE(LC), QUARTZ_TABLE(LC), AXAJ_TABLE (LC), BXAJ_TABLE (LC), & + XXAJ_TABLE(LC) ENDDO CLOSE (21) diff --git a/trunk/NDHMS/Land_models/NoahMP/run/SOILPARM.TBL b/trunk/NDHMS/Land_models/NoahMP/run/SOILPARM.TBL index 518ba6f8b..86ce9e2cf 100644 --- a/trunk/NDHMS/Land_models/NoahMP/run/SOILPARM.TBL +++ b/trunk/NDHMS/Land_models/NoahMP/run/SOILPARM.TBL @@ -1,45 +1,45 @@ Soil Parameters STAS -19,1 'BB DRYSMC F11 MAXSMC REFSMC SATPSI SATDK SATDW WLTSMC QTZ ' -1, 2.79, 0.010, -0.472, 0.339, 0.192, 0.069, 4.66E-5, 2.65E-5, 0.010, 0.92, 'SAND' -2, 4.26, 0.028, -1.044, 0.421, 0.283, 0.036, 1.41E-5, 5.14E-6, 0.028, 0.82, 'LOAMY SAND' -3, 4.74, 0.047, -0.569, 0.434, 0.312, 0.141, 5.23E-6, 8.05E-6, 0.047, 0.60, 'SANDY LOAM' -4, 5.33, 0.084, 0.162, 0.476, 0.360, 0.759, 2.81E-6, 2.39E-5, 0.084, 0.25, 'SILT LOAM' -5, 3.86, 0.061, 0.162, 0.484, 0.347, 0.955, 2.18E-6, 1.66E-5, 0.061, 0.10, 'SILT' -6, 5.25, 0.066, -0.327, 0.439, 0.329, 0.355, 3.38E-6, 1.43E-5, 0.066, 0.40, 'LOAM' -7, 6.77, 0.069, -1.491, 0.404, 0.315, 0.135, 4.45E-6, 1.01E-5, 0.069, 0.60, 'SANDY CLAY LOAM' -8, 8.72, 0.120, -1.118, 0.464, 0.387, 0.617, 2.03E-6, 2.35E-5, 0.120, 0.10, 'SILTY CLAY LOAM' -9, 8.17, 0.103, -1.297, 0.465, 0.382, 0.263, 2.45E-6, 1.13E-5, 0.103, 0.35, 'CLAY LOAM' -10, 10.73, 0.100, -3.209, 0.406, 0.338, 0.098, 7.22E-6, 1.87E-5, 0.100, 0.52, 'SANDY CLAY' -11, 10.39, 0.126, -1.916, 0.468, 0.404, 0.324, 1.34E-6, 9.64E-6, 0.126, 0.10, 'SILTY CLAY' -12, 11.55, 0.138, -2.138, 0.468, 0.412, 0.468, 9.74E-7, 1.12E-5, 0.138, 0.25, 'CLAY' -13, 5.25, 0.066, -0.327, 0.439, 0.329, 0.355, 3.38E-6, 1.43E-5, 0.066, 0.05, 'ORGANIC MATERIAL' -14, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.60, 'WATER' -15, 2.79, 0.006, -1.111, 0.20, 0.17, 0.069, 1.41E-4, 1.36E-4, 0.006, 0.07, 'BEDROCK' -16, 4.26, 0.028, -1.044, 0.421, 0.283, 0.036, 1.41E-5, 5.14E-6, 0.028, 0.25, 'OTHER(land-ice)' -17, 11.55, 0.030, -10.472, 0.468, 0.454, 0.468, 9.74E-7, 1.12E-5, 0.030, 0.60, 'PLAYA' -18, 2.79, 0.006, -0.472, 0.200, 0.17, 0.069, 1.41E-4, 1.36E-4, 0.006, 0.52, 'LAVA' -19, 2.79, 0.01, -0.472, 0.339, 0.192, 0.069, 4.66E-5, 2.65E-5, 0.01, 0.92, 'WHITE SAND' +19,1 'BB DRYSMC F11 MAXSMC REFSMC SATPSI SATDK SATDW WLTSMC QTZ AXAJ BXAJ XXAJ ' +1, 2.79, 0.010, -0.472, 0.339, 0.192, 0.069, 4.66E-5, 2.65E-5, 0.010, 0.92, 0.009, 0.05, 0.05, 'SAND' +2, 4.26, 0.028, -1.044, 0.421, 0.283, 0.036, 1.41E-5, 5.14E-6, 0.028, 0.82, 0.010, 0.08, 0.08, 'LOAMY SAND' +3, 4.74, 0.047, -0.569, 0.434, 0.312, 0.141, 5.23E-6, 8.05E-6, 0.047, 0.60, 0.009, 0.09, 0.09, 'SANDY LOAM' +4, 5.33, 0.084, 0.162, 0.476, 0.360, 0.759, 2.81E-6, 2.39E-5, 0.084, 0.25, 0.010, 0.25, 0.25, 'SILT LOAM' +5, 3.86, 0.061, 0.162, 0.484, 0.347, 0.955, 2.18E-6, 1.66E-5, 0.061, 0.10, 0.012, 0.15, 0.15, 'SILT' +6, 5.25, 0.066, -0.327, 0.439, 0.329, 0.355, 3.38E-6, 1.43E-5, 0.066, 0.40, 0.013, 0.18, 0.18, 'LOAM' +7, 6.77, 0.069, -1.491, 0.404, 0.315, 0.135, 4.45E-6, 1.01E-5, 0.069, 0.60, 0.014, 0.20, 0.20, 'SANDY CLAY LOAM' +8, 8.72, 0.120, -1.118, 0.464, 0.387, 0.617, 2.03E-6, 2.35E-5, 0.120, 0.10, 0.015, 0.22, 0.22, 'SILTY CLAY LOAM' +9, 8.17, 0.103, -1.297, 0.465, 0.382, 0.263, 2.45E-6, 1.13E-5, 0.103, 0.35, 0.016, 0.23, 0.23, 'CLAY LOAM' +10, 10.73, 0.100, -3.209, 0.406, 0.338, 0.098, 7.22E-6, 1.87E-5, 0.100, 0.52, 0.015, 0.25, 0.25, 'SANDY CLAY' +11, 10.39, 0.126, -1.916, 0.468, 0.404, 0.324, 1.34E-6, 9.64E-6, 0.126, 0.10, 0.016, 0.28, 0.28, 'SILTY CLAY' +12, 11.55, 0.138, -2.138, 0.468, 0.412, 0.468, 9.74E-7, 1.12E-5, 0.138, 0.25, 0.017, 0.30, 0.30, 'CLAY' +13, 5.25, 0.066, -0.327, 0.439, 0.329, 0.355, 3.38E-6, 1.43E-5, 0.066, 0.05, 0.012, 0.26, 0.26, 'ORGANIC MATERIAL' +14, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.60, 0.001, 0.00, 0.00, 'WATER' +15, 2.79, 0.006, -1.111, 0.20, 0.17, 0.069, 1.41E-4, 1.36E-4, 0.006, 0.07, 0.017, 1.00, 1.00, 'BEDROCK' +16, 4.26, 0.028, -1.044, 0.421, 0.283, 0.036, 1.41E-5, 5.14E-6, 0.028, 0.25, 0.017, 1.00, 1.00, 'OTHER(land-ice)' +17, 11.55, 0.030, -10.472, 0.468, 0.454, 0.468, 9.74E-7, 1.12E-5, 0.030, 0.60, 0.017, 1.00, 1.00, 'PLAYA' +18, 2.79, 0.006, -0.472, 0.200, 0.17, 0.069, 1.41E-4, 1.36E-4, 0.006, 0.52, 0.015, 0.35, 0.35, 'LAVA' +19, 2.79, 0.01, -0.472, 0.339, 0.192, 0.069, 4.66E-5, 2.65E-5, 0.01, 0.92, 0.009, 0.15, 0.15, 'WHITE SAND' Soil Parameters STAS-RUC -19,1 'BB DRYSMC HC MAXSMC REFSMC SATPSI SATDK SATDW WLTSMC QTZ ' -1, 4.05, 0.002, 1.47, 0.395, 0.174, 0.121, 1.76E-4, 0.608E-6, 0.033, 0.92, 'SAND' -2, 4.38, 0.035, 1.41, 0.410, 0.179, 0.090, 1.56E-4, 0.514E-5, 0.055, 0.82, 'LOAMY SAND' -3, 4.90, 0.041, 1.34, 0.435, 0.249, 0.218, 3.47E-5, 0.805E-5, 0.095, 0.60, 'SANDY LOAM' -4, 5.30, 0.034, 1.27, 0.485, 0.369, 0.786, 7.20E-6, 0.239E-4, 0.143, 0.25, 'SILT LOAM' -5, 5.30, 0.034, 1.27, 0.485, 0.369, 0.786, 7.20E-6, 0.239E-4, 0.143, 0.10, 'SILT' -6, 5.39, 0.050, 1.21, 0.451, 0.314, 0.478, 6.95E-6, 0.143E-4, 0.137, 0.40, 'LOAM' -7, 7.12, 0.068, 1.18, 0.420, 0.299, 0.299, 6.30E-6, 0.990E-5, 0.148, 0.60, 'SANDY CLAY LOAM' -8, 7.75, 0.060, 1.32, 0.477, 0.357, 0.356, 1.70E-6, 0.237E-4, 0.208, 0.10, 'SILTY CLAY LOAM' -9, 8.52, 0.085, 1.23, 0.476, 0.391, 0.630, 2.45E-6, 0.113E-4, 0.230, 0.35, 'CLAY LOAM' -10, 10.40, 0.100, 1.18, 0.426, 0.316, 0.153, 2.17E-6, 0.187E-4, 0.210, 0.52, 'SANDY CLAY' -11, 10.40, 0.070, 1.15, 0.492, 0.409, 0.490, 1.03E-6, 0.964E-5, 0.250, 0.10, 'SILTY CLAY' -12, 11.40, 0.068, 1.09, 0.482, 0.400, 0.405, 1.28E-6, 0.112E-4, 0.268, 0.25, 'CLAY' -13, 5.39, 0.027, 1.21, 0.451, 0.314, 0.478, 6.95E-6, 0.143E-4, 0.117, 0.05, 'ORGANIC MATERIAL' -14, 0.0, 0.0, 4.18, 1.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.00, 'WATER' -15, 4.05, 0.004, 2.03, 0.200, 0.10 , 0.121, 1.41E-4, 0.136E-3, 0.006, 0.60, 'BEDROCK' -16, 4.90, 0.065, 2.10, 0.435, 0.249, 0.218, 3.47E-5, 0.514E-5, 0.114, 0.05, 'OTHER(land-ice)' -17, 11.40, 0.030, 1.41, 0.468, 0.454, 0.468, 9.74E-7, 0.112E-4, 0.030, 0.60, 'PLAYA' -18, 4.05, 0.006, 1.41, 0.200, 0.17, 0.069, 1.41E-4, 0.136E-3, 0.006, 0.52, 'LAVA' -19, 4.05, 0.01, 1.47, 0.339, 0.236, 0.069, 1.76E-4, 0.608E-6, 0.060, 0.92, 'WHITE SAND' +19,1 'BB DRYSMC HC MAXSMC REFSMC SATPSI SATDK SATDW WLTSMC QTZ AXAJ BXAJ XXAJ ' +1, 4.05, 0.002, 1.47, 0.395, 0.174, 0.121, 1.76E-4, 0.608E-6, 0.033, 0.92, 0.009, 0.05, 0.05, 'SAND' +2, 4.38, 0.035, 1.41, 0.410, 0.179, 0.090, 1.56E-4, 0.514E-5, 0.055, 0.82, 0.010, 0.08, 0.08, 'LOAMY SAND' +3, 4.90, 0.041, 1.34, 0.435, 0.249, 0.218, 3.47E-5, 0.805E-5, 0.095, 0.60, 0.009, 0.09, 0.09, 'SANDY LOAM' +4, 5.30, 0.034, 1.27, 0.485, 0.369, 0.786, 7.20E-6, 0.239E-4, 0.143, 0.25, 0.010, 0.25, 0.25, 'SILT LOAM' +5, 5.30, 0.034, 1.27, 0.485, 0.369, 0.786, 7.20E-6, 0.239E-4, 0.143, 0.10, 0.012, 0.15, 0.15, 'SILT' +6, 5.39, 0.050, 1.21, 0.451, 0.314, 0.478, 6.95E-6, 0.143E-4, 0.137, 0.40, 0.013, 0.18, 0.18, 'LOAM' +7, 7.12, 0.068, 1.18, 0.420, 0.299, 0.299, 6.30E-6, 0.990E-5, 0.148, 0.60, 0.014, 0.20, 0.20, 'SANDY CLAY LOAM' +8, 7.75, 0.060, 1.32, 0.477, 0.357, 0.356, 1.70E-6, 0.237E-4, 0.208, 0.10, 0.015, 0.22, 0.22, 'SILTY CLAY LOAM' +9, 8.52, 0.085, 1.23, 0.476, 0.391, 0.630, 2.45E-6, 0.113E-4, 0.230, 0.35, 0.016, 0.23, 0.23, 'CLAY LOAM' +10, 10.40, 0.100, 1.18, 0.426, 0.316, 0.153, 2.17E-6, 0.187E-4, 0.210, 0.52, 0.015, 0.25, 0.25, 'SANDY CLAY' +11, 10.40, 0.070, 1.15, 0.492, 0.409, 0.490, 1.03E-6, 0.964E-5, 0.250, 0.10, 0.016, 0.28, 0.28, 'SILTY CLAY' +12, 11.40, 0.068, 1.09, 0.482, 0.400, 0.405, 1.28E-6, 0.112E-4, 0.268, 0.25, 0.017, 0.30, 0.30, 'CLAY' +13, 5.39, 0.027, 1.21, 0.451, 0.314, 0.478, 6.95E-6, 0.143E-4, 0.117, 0.05, 0.012, 0.26, 0.26, 'ORGANIC MATERIAL' +14, 0.0, 0.0, 4.18, 1.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.00, 0.001, 0.00, 0.00, 'WATER' +15, 4.05, 0.004, 2.03, 0.200, 0.10 , 0.121, 1.41E-4, 0.136E-3, 0.006, 0.60, 0.017, 1.00, 1.00, 'BEDROCK' +16, 4.90, 0.065, 2.10, 0.435, 0.249, 0.218, 3.47E-5, 0.514E-5, 0.114, 0.05, 0.017, 1.00, 1.00, 'OTHER(land-ice)' +17, 11.40, 0.030, 1.41, 0.468, 0.454, 0.468, 9.74E-7, 0.112E-4, 0.030, 0.60, 0.017, 1.00, 1.00, 'PLAYA' +18, 4.05, 0.006, 1.41, 0.200, 0.17, 0.069, 1.41E-4, 0.136E-3, 0.006, 0.52, 0.015, 0.35, 0.35, 'LAVA' +19, 4.05, 0.01, 1.47, 0.339, 0.236, 0.069, 1.76E-4, 0.608E-6, 0.060, 0.92, 0.009, 0.15, 0.15, 'WHITE SAND' diff --git a/trunk/NDHMS/hrldas_namelists.json b/trunk/NDHMS/hrldas_namelists.json index 3628b8443..1f496da1d 100644 --- a/trunk/NDHMS/hrldas_namelists.json +++ b/trunk/NDHMS/hrldas_namelists.json @@ -43,12 +43,16 @@ }, "nwm_ana": { - "noahlsm_offline": {}, + "noahlsm_offline": { + "runoff_option": 7 + }, "wrf_hydro_offline": {} }, "nwm_long_range" : { - "noahlsm_offline": {}, + "noahlsm_offline": { + "runoff_option": 7 + }, "wrf_hydro_offline": {} }, @@ -66,6 +70,7 @@ "nwm_output_ana" : { "noahlsm_offline": { + "runoff_option": 7, "khour": 1, "restart_frequency_hours": 1 }, @@ -74,6 +79,7 @@ "nwm_output_long_range" : { "noahlsm_offline": { + "runoff_option": 7, "khour": 1, "restart_frequency_hours": 1 },