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/parin.f90

    r4022 r4079  
    2525! -----------------
    2626! $Id$
     27! +monotonic_limiter_z
     29! 4022 2019-06-12 11:52:39Z suehring
    2730! Change default top boundary condition for pressure to Neumann in offline
    2831! nesting case
    590593             loop_optimization, lsf_exception, masking_method, mg_cycles,      &
    591594             mg_switch_to_pe0_level, mixing_length_1d, momentum_advec,         &
     595             monotonic_limiter_z,                                              &
    592596             netcdf_precision, neutral, ngsrb,                                 &
    593597             nsor, nsor_ini, nudging, nx, ny, nz, ocean_mode, omega,           &
    660664             loop_optimization, lsf_exception, masking_method, mg_cycles,      &
    661665             mg_switch_to_pe0_level, mixing_length_1d, momentum_advec,         &
     666             monotonic_limiter_z,                                              &
    662667             netcdf_precision, neutral, ngsrb,                                 &
    663668             nsor, nsor_ini, nudging, nx, ny, nz, ocean_mode, omega,           &
Note: See TracChangeset for help on using the changeset viewer.