Jul 9, 2019 6:04:41 PM (2 years ago)

Implementation of a monotonic flux limiter for vertical advection term in Wicker-Skamarock scheme. The flux limiter is currently only applied for passive scalars (passive scalar, chemical species, aerosols) within the region up to the highest topography, in order to avoid the built-up of large concentrations within poorly resolved cavities in urban environments. To enable the limiter monotonic_limiter_z = .T. must be set. Note, the limiter is currently only implemented for the cache-optimized version of advec_ws. Further changes in offline nesting: Set boundary condition for w at nzt+1 at all lateral boundaries (even though these won't enter the numerical solution), in order to avoid high vertical velocities in the run-control file which might built-up due to the mass-conservation; bugfix in offline nesting for chemical species

1 edited


  • palm/trunk/SOURCE/salsa_mod.f90

    r4069 r4079  
    2626! -----------------
    2727! $Id$
     28! Application of monotonic flux limiter for the vertical scalar advection
     29! up to the topography top (only for the cache-optimized version at the
     30! moment).
     32! 4069 2019-07-01 14:05:51Z Giersch
    2833! Masked output running index mid has been introduced as a local variable to
    2934! avoid runtime error (Loop variable has been modified) in time_integration
    77737778    IF ( timestep_scheme(1:5) == 'runge' )  THEN
    77747779       IF ( ws_scheme_sca )  THEN
    7775           CALL advec_s_ws( i, j, rs, id, flux_s, diss_s, flux_l, diss_l, i_omp_start, tn )
     7780          CALL advec_s_ws( i, j, rs, id, flux_s, diss_s, flux_l, diss_l, i_omp_start, tn,          &
     7781                           monotonic_limiter_z )
    77767782       ELSE
    77777783          CALL advec_s_pw( i, j, rs )
Note: See TracChangeset for help on using the changeset viewer.