Ignore:
Timestamp:
May 9, 2019 9:48:32 AM (2 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:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/nesting_offl_mod.f90

    r3937 r3964  
    2525! -----------------
    2626! $Id$
     27! Ensure that veloctiy term in calculation of bulk Richardson number does not
     28! become zero
     29!
     30! 3937 2019-04-29 15:09:07Z suehring
    2731! Set boundary conditon on upper-left and upper-south grid point for the u- and
    2832! v-component, respectively.
     
    900904!--          Note, no interpolation of u- and v-component is made, as both
    901905!--          are mainly mean inflow profiles with very small spatial variation.
     906!--          Add a safety factor in case the velocity term becomes zero. This
     907!--          may happen if overhanging 3D structures are directly located at
     908!--          the boundary, where velocity inside the building is zero
     909!--          (k_surface is the index of the lowest upward-facing surface).
    902910             zi_local = 0.0_wp
    903911             DO  k = k_surface+1, nzt
     
    906914                ri_bulk = zu(k) * g / vpt_surface *                            &
    907915                          ( vpt_col(k) - vpt_surface ) /                       &
    908                           ( u_comp * u_comp + v_comp * v_comp )
     916                          ( u_comp * u_comp + v_comp * v_comp + 1E-5_wp )
    909917                       
    910918                IF ( zi_local == 0.0_wp  .AND.  ri_bulk > ri_bulk_crit )       &
     
    943951                ri_bulk = zu(k) * g / vpt_surface *                            &
    944952                       ( vpt_col(k) - vpt_surface ) /                          &
    945                        ( u_comp * u_comp + v_comp * v_comp )
     953                       ( u_comp * u_comp + v_comp * v_comp + 1E-5_wp )
    946954                       
    947955                IF ( zi_local == 0.0_wp  .AND.  ri_bulk > 0.25_wp )            &
Note: See TracChangeset for help on using the changeset viewer.