Ignore:
Timestamp:
Oct 5, 2018 12:34:56 PM (3 years ago)
Author:
raasch
Message:

Stokes drift is regarded in timestep calculation, check if ocean mode is used for invalid combinations

File:
1 edited

Legend:

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

    r3274 r3311  
    2525! -----------------
    2626! $Id$
     27! Stokes drift is regarded in timestep calculation
     28!
     29! 3274 2018-09-24 15:42:55Z knoop
    2730! Modularization of all bulk cloud physics code components
    2831!
     
    136139
    137140    USE arrays_3d,                                                             &
    138         ONLY:  dzu, dzw, kh, km, u, v, w
     141        ONLY:  dzu, dzw, kh, km, u, u_stokes_zu, v, v_stokes_zu, w
    139142
    140143    USE control_parameters,                                                    &
     
    278281             DO  k = nzb+1, nzt
    279282                dt_u_l = MIN( dt_u_l, ( dx     /                               &
    280                                  ( ABS( u(k,j,i) - u_gtrans ) + 1.0E-10_wp ) ) )
     283                                 ( ABS( u(k,j,i) - u_gtrans + u_stokes_zu(k) ) &
     284                                   + 1.0E-10_wp ) ) )
    281285                dt_v_l = MIN( dt_v_l, ( dy     /                               &
    282                                  ( ABS( v(k,j,i) - v_gtrans ) + 1.0E-10_wp ) ) )
     286                                 ( ABS( v(k,j,i) - v_gtrans + v_stokes_zu(k) ) &
     287                                   + 1.0E-10_wp ) ) )
    283288                dt_w_l = MIN( dt_w_l, ( dzu(k) /                               &
    284289                                 ( ABS( w(k,j,i) )            + 1.0E-10_wp ) ) )
Note: See TracChangeset for help on using the changeset viewer.