Changeset 4640 for palm


Ignore:
Timestamp:
Aug 11, 2020 4:28:32 PM (4 years ago)
Author:
suehring
Message:

synthetic turbulence generator: to avoid that the correction term in r11/r22 computation becomes unrealistically high, limit Obukhov length (term is not valid for near neutral conditions); to avoid unrealistically large perturbations, change computation of r21 so that this resembles the vertical transport of horizontal momentum

File:
1 edited

Legend:

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

    r4629 r4640  
    2525! -----------------
    2626! $Id$
     27! - to avoid that the correction term in r11/r22 computation becomes unrealistically high, limit
     28!   Obukhov length (term is not valid for near neutral conditions)
     29! - to avoid unrealistically large perturbations, change computation of r21 so that this resembles
     30!   the vertical transport of horizontal momentum
     31!
     32! 4629 2020-07-29 09:37:56Z raasch
    2733! support for MPI Fortran77 interface (mpif.h) removed
    2834!
     
    21302136                                       )**( 2.0_wp / 3.0_wp ),                                     &
    21312137                          0.0_wp,                                                                  &
    2132                           scale_l < 0.0_wp )
     2138                          scale_l < -5.0_wp )
    21332139!
    21342140!--    Determine normalized height coordinate
     
    21442150       r22(k) = ABS( - sin_phi(k) * rlon1 + cos_phi(k) * rlat1 )
    21452151!
    2146 !--    For u'v' no parametrization exist so far - ?. For simplicity assume a similar profile as
    2147 !--    for u'w'.
    2148        r21(k) = 0.5_wp * ( r11(k) + r22(k) )
    2149 !
    21502152!--    w'w'
    21512153       r33(k) = scale_wm**2 * (                                                                    &
     
    21642166       r32(k) = SIGN( ABS( - sin_phi(k) * rlon2 + cos_phi(k) * rlat2 ),                            &
    21652167                   - ( v_init(k+1) - v_init(k-1) ) )
     2168!
     2169!--    For u'v' no parametrization exist so far. For simplicity assume a similar profile as
     2170!--    for the vertical transport.
     2171       r21(k) = 0.5_wp * ( r31(k) + r32(k) )
    21662172!
    21672173!--    Compute turbulent time scales according to Brost et al. (1982). Note, time scales are
     
    22172223       r31(k) = r31_zi * blend
    22182224       r32(k) = r32_zi * blend
    2219        r21(k) = 0.5_wp * ( r11(k) + r22(k) )
     2225       r21(k) = 0.5_wp * ( r31(k) + r32(k) )
    22202226
    22212227!
Note: See TracChangeset for help on using the changeset viewer.