From a5ca65fe4a3312316091c3eec0d7036641eb5f97 Mon Sep 17 00:00:00 2001
From: Michael Barlage <barlage@ucar.edu>
Date: Tue, 19 Mar 2019 10:08:19 -0600
Subject: [PATCH] set bulk urban roughness to table value

---
 phys/module_sf_noahmplsm.F | 9 +++++++++
 1 file changed, 9 insertions(+)

diff --git a/phys/module_sf_noahmplsm.F b/phys/module_sf_noahmplsm.F
index 5dadbe8728..ba994d601c 100644
--- 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