Ignore:
Timestamp:
Jun 12, 2018 7:03:02 AM (3 years ago)
Author:
Giersch
Message:

New vertical stretching procedure has been introduced

File:
1 edited

Legend:

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

    r2718 r3065  
    2525! -----------------
    2626! $Id$
     27! dz was replaced by dzw to allow for right vertical stretching
     28!
     29! 2718 2018-01-02 08:49:38Z maronga
    2730! Corrected "Former revisions" section
    2831!
     
    8083
    8184    USE arrays_3d,                                                             &
    82         ONLY:  pt, u, v, w, zu, zw
     85        ONLY:  dzw, pt, u, v, w, zu, zw
    8386
    8487    USE control_parameters,                                                    &
    85         ONLY:  u_gtrans, v_gtrans, dz
     88        ONLY:  u_gtrans, v_gtrans
    8689
    8790    USE cpulog,                                                                &
     
    207210                                     ( gg - dd ) * u(k+1,j+1,i+1)                &
    208211                                   ) / ( 3.0_wp * gg ) - u_gtrans
    209                          u_int(n) = u_int_l + ( zv(n) - zu(k) ) / dz *         &
     212                         u_int(n) = u_int_l + ( zv(n) - zu(k) ) / dzw(k) *      &
    210213                                           ( u_int_u - u_int_l )
    211214                      ENDIF
     
    240243                                      ( gg - dd ) * v(k+1,j+1,i+1)                &
    241244                                    ) / ( 3.0_wp * gg ) - v_gtrans
    242                          v_int(n) = v_int_l + ( zv(n) - zu(k) ) / dz *         &
     245                         v_int(n) = v_int_l + ( zv(n) - zu(k) ) / dzw(k) *      &
    243246                                           ( v_int_u - v_int_l )
    244247                      ENDIF
     
    344347                                 ) / ( 3.0_wp * gg ) - sums(k,4)
    345348
    346                       pt_int = pt_int_l + ( zv(n) - zu(k) ) / dz *    &
     349                      pt_int = pt_int_l + ( zv(n) - zu(k) ) / dzw(k) *          &
    347350                                          ( pt_int_u - pt_int_l )
    348351
Note: See TracChangeset for help on using the changeset viewer.