diff --git a/drivers/wrf/module_sf_noahmpdrv.F b/drivers/wrf/module_sf_noahmpdrv.F index a750be24..9fc16db9 100644 --- a/drivers/wrf/module_sf_noahmpdrv.F +++ b/drivers/wrf/module_sf_noahmpdrv.F @@ -3024,7 +3024,9 @@ SUBROUTINE noahmp_urban(sf_urban_physics, NSOIL, IVGTYP, ITIMESTEP, USE module_sf_bep, only: bep USE module_sf_bep_bem, only: bep_bem USE module_ra_gfdleta, only: cal_mon_day - USE NOAHMP_TABLES, ONLY: ISURBAN_TABLE + USE NOAHMP_TABLES, ONLY: ISURBAN_TABLE,LCZ_1_TABLE,LCZ_2_TABLE,LCZ_3_TABLE,LCZ_4_TABLE, & + LCZ_5_TABLE,LCZ_6_TABLE,LCZ_7_TABLE,LCZ_8_TABLE, & + LCZ_9_TABLE,LCZ_10_TABLE,LCZ_11_TABLE USE module_model_constants, only: KARMAN, CP, XLV !---------------------------------------------------------------- IMPLICIT NONE @@ -3345,8 +3347,10 @@ SUBROUTINE noahmp_urban(sf_urban_physics, NSOIL, IVGTYP, ITIMESTEP, ILOOP : DO I = its, ite - IF( IVGTYP(I,J) == ISURBAN_TABLE .or. IVGTYP(I,J) == 31 .or. & - IVGTYP(I,J) == 32 .or. IVGTYP(I,J) == 33 ) THEN + IF( IVGTYP(I,J) == ISURBAN_TABLE .or. IVGTYP(I,J) == LCZ_1_TABLE .or. IVGTYP(I,J) == LCZ_2_TABLE .or. & + IVGTYP(I,J) == LCZ_3_TABLE .or. IVGTYP(I,J) == LCZ_4_TABLE .or. IVGTYP(I,J) == LCZ_5_TABLE .or. & + IVGTYP(I,J) == LCZ_6_TABLE .or. IVGTYP(I,J) == LCZ_7_TABLE .or. IVGTYP(I,J) == LCZ_8_TABLE .or. & + IVGTYP(I,J) == LCZ_9_TABLE .or. IVGTYP(I,J) == LCZ_10_TABLE .or. IVGTYP(I,J) == LCZ_11_TABLE ) THEN UTYPE_URB = UTYPE_URB2D(I,J) !urban type (low, high or industrial) diff --git a/src/module_sf_noahmplsm.F b/src/module_sf_noahmplsm.F index 3cd25c49..2df8e0d1 100644 --- a/src/module_sf_noahmplsm.F +++ b/src/module_sf_noahmplsm.F @@ -7696,7 +7696,7 @@ SUBROUTINE INFIL (parameters,NSOIL ,DT ,ZSOIL ,SH2O ,SICE , & !in CALL WDFCND2 (parameters,WDF,WCND,SH2O(1),SICEMAX,1) INFMAX = MAX (INFMAX,WCND) - INFMAX = MIN (INFMAX,PX) + INFMAX = MIN (INFMAX,PX/DT) RUNSRF= MAX(0.0, QINSUR - INFMAX) PDDUM = QINSUR - RUNSRF