Changeset 2703 for palm

Timestamp:

Dec 15, 2017 8:12:38 PM (6 years ago)

Author:

maronga

Message:

workaround to keep aerodyn. resistances positive in urban surface model

File:

• palm/trunk/SOURCE/urban_surface_mod.f90 (modified) (2 diffs)

palm/trunk/SOURCE/urban_surface_mod.f90

@@ -21 +21 @@
 ! Current revisions:
 ! ------------------
-! 
+! Workaround for calculation of r_a
 ! 
 ! Former revisions:
@@ -6890 +6890 @@
 !--        pt, us, ts are not available for the prognostic time step,
 !--        data from the last time step is used here.
-               
-           r_a = ( surf_usm_h%pt1(m) - t_surf_h(m) / exn(k) ) /                              &
-                 ( surf_usm_h%ts(m) * surf_usm_h%us(m) + 1.0E-20_wp )
-           r_a_window = ( surf_usm_h%pt1(m) - t_surf_window_h(m) / exn(k) ) /                &
-                 ( surf_usm_h%ts(m) * surf_usm_h%us(m) + 1.0E-20_wp )
-           r_a_green = ( surf_usm_h%pt1(m) - t_surf_green_h(m) / exn(k) ) /                  &
-                 ( surf_usm_h%ts(m) * surf_usm_h%us(m) + 1.0E-20_wp )
+
+!--        Workarounf to keep resistances positive
+           r_a = ( surf_usm_h%pt1(m) - surf_usm_h%pt_surface(m) ) /                            &
+                 ( surf_usm_h%ts(m) * surf_usm_h%us(m) + 1.0E-20_wp )   
+           r_a_window = r_a
+           r_a_green  = r_a
+!            r_a = ( surf_usm_h%pt1(m) - t_surf_h(m) / exn(k) ) /                              &
+!                  ( surf_usm_h%ts(m) * surf_usm_h%us(m) + 1.0E-20_wp )
+!            r_a_window = ( surf_usm_h%pt1(m) - t_surf_window_h(m) / exn(k) ) /                &
+!                  ( surf_usm_h%ts(m) * surf_usm_h%us(m) + 1.0E-20_wp )
+!            r_a_green = ( surf_usm_h%pt1(m) - t_surf_green_h(m) / exn(k) ) /                  &
+!                  ( surf_usm_h%ts(m) * surf_usm_h%us(m) + 1.0E-20_wp )
                 
 !--        make sure that the resistance does not drop to zero