Ignore:
Timestamp:
Jun 19, 2018 3:30:55 PM (3 years ago)
Author:
gronemeier
Message:

bugfix: limit increase of dt_3d only in case of RANS mode

File:
1 edited

Legend:

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

    r3083 r3084  
    2525! -----------------
    2626! $Id$
     27! limit increase of dt_3d only in case of RANS mode
     28!
     29! 3083 2018-06-19 14:03:12Z gronemeier
    2730! limit dt_3d to be at maximum 2*old_dt; define old_dt at beginning of routine
    2831! Add km/kh_max
     
    128131    USE control_parameters,                                                    &
    129132        ONLY:  cfl_factor, coupling_mode, dt_3d, dt_fixed, dt_max,             &
    130                galilei_transformation, old_dt, message_string,                 &
     133               galilei_transformation, old_dt, message_string, rans_mode,      &
    131134               stop_dt, terminate_coupled, terminate_coupled_remote,           &
    132135               timestep_reason, u_gtrans, use_ug_for_galilei_tr, v_gtrans
     
    328331!--    The time step is the minimum of the 3-4 components and the diffusion time
    329332!--    step minus a reduction (cfl_factor) to be on the safe side.
    330 !--    The time step must not exceed the maximum allowed value and must not
    331 !--    increase by more than a factor of 10.
     333!--    The time step must not exceed the maximum allowed value.
    332334       dt_3d = cfl_factor * MIN( dt_diff, dt_u, dt_v, dt_w, dt_precipitation )
    333        dt_3d = MIN( dt_3d, dt_max, 2.0_wp * old_dt )
     335       dt_3d = MIN( dt_3d, dt_max )
     336!
     337!--    In RANS mode, the time step must not increase by more than a factor of 2
     338       IF ( rans_mode )  dt_3d = MIN( dt_3d, dt_max, 2.0_wp * old_dt )
    334339
    335340!
Note: See TracChangeset for help on using the changeset viewer.