set bulk urban roughness to table value

wrf-model · Mar 19, 2019 · a5ca65f · a5ca65f
1 parent 9d45381
commit a5ca65f
Showing 1 changed file with 9 additions and 0 deletions.
diff --git a/phys/module_sf_noahmplsm.F b/phys/module_sf_noahmplsm.F
@@ -1857,6 +1857,15 @@ SUBROUTINE ENERGY (parameters,ICE    ,VEGTYP ,IST    ,NSNOW  ,NSOIL  , & !in
         ZPD  = ZPDG
      END IF
 
+! special case for urban
+
+     IF (parameters%urban_flag) THEN
+       Z0MG = parameters%Z0MVT
+       ZPDG  = 0.65 * parameters%HVT
+       Z0M  = Z0MG
+       ZPD  = ZPDG
+     END IF
+
      ZLVL = MAX(ZPD,parameters%HVT) + ZREF
      IF(ZPDG >= ZLVL) ZLVL = ZPDG + ZREF
 !     UR   = UR*LOG(ZLVL/Z0M)/LOG(10./Z0M)       !input UR is at 10m