Changeset 3256 for palm/trunk/SOURCE

Timestamp:

Sep 17, 2018 12:20:07 PM (6 years ago)

Author:

suehring

Message:

Enable initialization of z0q for vegetation, pavement and water surfaces via namelist input

File:

• palm/trunk/SOURCE/land_surface_model_mod.f90 (modified) (10 diffs)

palm/trunk/SOURCE/land_surface_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Enable initialization of z0q for vegetation, pavement and water surfaces via
+! namelist input. 
+! 
+! 3248 2018-09-14 09:42:06Z sward
 ! Minor formating changes
 ! 
@@ -3140 +3144 @@
           IF ( z0h_vegetation == 9999999.9_wp )  THEN
              z0h_vegetation = vegetation_pars(ind_v_z0qh,vegetation_type)
-          ENDIF   
+          ENDIF 
+          
+          IF ( z0q_vegetation == 9999999.9_wp )  THEN
+             z0q_vegetation = vegetation_pars(ind_v_z0qh,vegetation_type)
+          ENDIF
           
           IF ( lambda_surface_stable == 9999999.9_wp )  THEN
@@ -3177 +3185 @@
              surf_lsm_h%z0(m)               = z0_vegetation
              surf_lsm_h%z0h(m)              = z0h_vegetation
-             surf_lsm_h%z0q(m)              = z0h_vegetation
+             surf_lsm_h%z0q(m)              = z0q_vegetation
              surf_lsm_h%lambda_surface_s(m) = lambda_surface_stable
              surf_lsm_h%lambda_surface_u(m) = lambda_surface_unstable
@@ -3206 +3214 @@
                 surf_lsm_v(l)%z0(m)               = z0_vegetation
                 surf_lsm_v(l)%z0h(m)              = z0h_vegetation
-                surf_lsm_v(l)%z0q(m)              = z0h_vegetation
+                surf_lsm_v(l)%z0q(m)              = z0q_vegetation
                 surf_lsm_v(l)%lambda_surface_s(m) = lambda_surface_stable
                 surf_lsm_v(l)%lambda_surface_u(m) = lambda_surface_unstable
@@ -3461 +3469 @@
              z0h_water = water_pars(ind_w_z0h,water_type)       
           ENDIF  
+          
+          IF ( z0q_water == 9999999.9_wp )  THEN
+             z0q_water = water_pars(ind_w_z0h,water_type)       
+          ENDIF 
 
           IF ( albedo_type == 9999999  .AND.  albedo == 9999999.9_wp )  THEN
@@ -3479 +3491 @@
              surf_lsm_h%z0(m)               = z0_water
              surf_lsm_h%z0h(m)              = z0h_water
-             surf_lsm_h%z0q(m)              = z0h_water
+             surf_lsm_h%z0q(m)              = z0q_water
              surf_lsm_h%lambda_surface_s(m) = 1.0E10_wp
              surf_lsm_h%lambda_surface_u(m) = 1.0E10_wp               
@@ -3500 +3512 @@
                 surf_lsm_v(l)%z0(m)               = z0_water
                 surf_lsm_v(l)%z0h(m)              = z0h_water
-                surf_lsm_v(l)%z0q(m)              = z0h_water
+                surf_lsm_v(l)%z0q(m)              = z0q_water
                 surf_lsm_v(l)%lambda_surface_s(m) = 1.0E10_wp
                 surf_lsm_v(l)%lambda_surface_u(m) = 1.0E10_wp               
@@ -3712 +3724 @@
              z0h_pavement = pavement_pars(ind_p_z0h,pavement_type)
           ENDIF
+          
+          IF ( z0q_pavement == 9999999.9_wp )  THEN
+             z0q_pavement = pavement_pars(ind_p_z0h,pavement_type)
+          ENDIF
 
           IF ( pavement_heat_conduct == 9999999.9_wp )  THEN
@@ -3755 +3771 @@
              surf_lsm_h%z0(m)                  = z0_pavement
              surf_lsm_h%z0h(m)                 = z0h_pavement
-             surf_lsm_h%z0q(m)                 = z0h_pavement
+             surf_lsm_h%z0q(m)                 = z0q_pavement
              surf_lsm_h%lambda_surface_s(m)    = pavement_heat_conduct         &
                                                   * ddz_soil(nzb_soil)         &
@@ -3793 +3809 @@
                 surf_lsm_v(l)%z0(m)                  = z0_pavement
                 surf_lsm_v(l)%z0h(m)                 = z0h_pavement
-                surf_lsm_v(l)%z0q(m)                 = z0h_pavement
+                surf_lsm_v(l)%z0q(m)                 = z0q_pavement
                 surf_lsm_v(l)%lambda_surface_s(m)    = pavement_heat_conduct   &
                                                   * ddz_soil(nzb_soil)         &