Changeset 3964 for palm/trunk/SOURCE/nesting_offl_mod.f90

Timestamp:

May 9, 2019 9:48:32 AM (5 years ago)

Author:

suehring

Message:

In a nested child domain, distinguish between soil moisture and temperature initialization in case the parent domain is initialized via the dynamic input file; in the offline nesting, add a safety factor for the calculation of bulk Richardson number in order to avoid division by zero which can potentially happen if 3D buildings are located directly at the lateral model boundaries

File:

• palm/trunk/SOURCE/nesting_offl_mod.f90 (modified) (4 diffs)

palm/trunk/SOURCE/nesting_offl_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Ensure that veloctiy term in calculation of bulk Richardson number does not
+! become zero
+! 
+! 3937 2019-04-29 15:09:07Z suehring
 ! Set boundary conditon on upper-left and upper-south grid point for the u- and
 ! v-component, respectively.
@@ -900 +904 @@
 !--          Note, no interpolation of u- and v-component is made, as both 
 !--          are mainly mean inflow profiles with very small spatial variation. 
+!--          Add a safety factor in case the velocity term becomes zero. This
+!--          may happen if overhanging 3D structures are directly located at
+!--          the boundary, where velocity inside the building is zero 
+!--          (k_surface is the index of the lowest upward-facing surface).
              zi_local = 0.0_wp
              DO  k = k_surface+1, nzt
@@ -906 +914 @@
                 ri_bulk = zu(k) * g / vpt_surface *                            &
                           ( vpt_col(k) - vpt_surface ) /                       &
-                          ( u_comp * u_comp + v_comp * v_comp ) 
+                          ( u_comp * u_comp + v_comp * v_comp + 1E-5_wp ) 
                        
                 IF ( zi_local == 0.0_wp  .AND.  ri_bulk > ri_bulk_crit )       &
@@ -943 +951 @@
                 ri_bulk = zu(k) * g / vpt_surface *                            &
                        ( vpt_col(k) - vpt_surface ) /                          &
-                       ( u_comp * u_comp + v_comp * v_comp ) 
+                       ( u_comp * u_comp + v_comp * v_comp + 1E-5_wp ) 
                        
                 IF ( zi_local == 0.0_wp  .AND.  ri_bulk > 0.25_wp )            &