Changeset 3065 for palm/trunk/SOURCE/lpm_advec.f90

Timestamp:

Jun 12, 2018 7:03:02 AM (6 years ago)

Author:

Giersch

Message:

New vertical stretching procedure has been introduced

File:

• palm/trunk/SOURCE/lpm_advec.f90 (modified) (11 diffs)

palm/trunk/SOURCE/lpm_advec.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! dz values were replaced by dzw or dz(1) to allow for right vertical stretching
+! 
+! 2969 2018-04-13 11:55:09Z thiele
 ! Bugfix in Interpolation indices.
 ! 
@@ -158 +161 @@
 
     USE arrays_3d,                                                             &
-        ONLY:  de_dx, de_dy, de_dz, diss, e, km, u, v, w, zu, zw
+        ONLY:  de_dx, de_dy, de_dz, diss, dzw, e, km, u, v, w, zu, zw
 
     USE cpulog
@@ -409 +412 @@
                             + ( gg-cc ) * u(k+1,j+1,i) + ( gg-dd ) *           &
                             u(k+1,j+1,i+1) ) / ( 3.0_wp * gg ) - u_gtrans
-                u_int(n) = u_int_l + ( zv(n) - zu(k) ) / dz *                  &
+                u_int(n) = u_int_l + ( zv(n) - zu(k) ) / dzw(k) *               &
                            ( u_int_u - u_int_l )
              ENDIF
@@ -503 +506 @@
                           + ( gg-cc ) * v(k+1,j+1,i) + ( gg-dd ) * v(k+1,j+1,i+1) &
                           ) / ( 3.0_wp * gg ) - v_gtrans
-                v_int(n) = v_int_l + ( zv(n) - zu(k) ) / dz * &
+                v_int(n) = v_int_l + ( zv(n) - zu(k) ) / dzw(k) *               &
                                   ( v_int_u - v_int_l )
              ENDIF
@@ -540 +543 @@
                             ( gg-dd ) * w(k+1,j+1,i+1) &
                           ) / ( 3.0_wp * gg )
-                w_int(n) = w_int_l + ( zv(n) - zw(k) ) / dz * &
+                w_int(n) = w_int_l + ( zv(n) - zw(k) ) / dzw(k) *               &
                            ( w_int_u - w_int_l )
              ENDIF
@@ -611 +614 @@
                                ( gg - dd ) * e(k+1,j+1,i+1) &
                             ) / ( 3.0_wp * gg )
-                   e_int(n) = e_int_l + ( zv(n) - zu(k) ) / dz * &
+                   e_int(n) = e_int_l + ( zv(n) - zu(k) ) / dzw(k) *            &
                                      ( e_int_u - e_int_l )
                 ENDIF
@@ -637 +640 @@
                                    ( gg - dd ) * de_dx(k+1,j+1,i+1) &
                                   ) / ( 3.0_wp * gg )
-                   de_dx_int(n) = de_dx_int_l + ( zv(n) - zu(k) ) / dz * &
+                   de_dx_int(n) = de_dx_int_l + ( zv(n) - zu(k) ) / dzw(k) *    &
                                               ( de_dx_int_u - de_dx_int_l )
                 ENDIF
@@ -655 +658 @@
                                    ( gg - dd ) * de_dy(k+1,j+1,i+1) &
                                   ) / ( 3.0_wp * gg )
-                      de_dy_int(n) = de_dy_int_l + ( zv(n) - zu(k) ) / dz * &
+                      de_dy_int(n) = de_dy_int_l + ( zv(n) - zu(k) ) / dzw(k) * &
                                                  ( de_dy_int_u - de_dy_int_l )
                 ENDIF
@@ -661 +664 @@
 !
 !--             Interpolate the TKE gradient along z
-                IF ( zv(n) < 0.5_wp * dz )  THEN
+                IF ( zv(n) < 0.5_wp * dz(1) )  THEN
                    de_dz_int(n) = 0.0_wp
                 ELSE
@@ -678 +681 @@
                                       ( gg - dd ) * de_dz(k+1,j+1,i+1) &
                                      ) / ( 3.0_wp * gg )
-                      de_dz_int(n) = de_dz_int_l + ( zv(n) - zu(k) ) / dz * &
+                      de_dz_int(n) = de_dz_int_l + ( zv(n) - zu(k) ) / dzw(k) * &
                                                  ( de_dz_int_u - de_dz_int_l )
                    ENDIF
@@ -699 +702 @@
                                   ( gg - dd ) * diss(k+1,j+1,i+1) &
                                  ) / ( 3.0_wp * gg )
-                   diss_int(n) = diss_int_l + ( zv(n) - zu(k) ) / dz * &
+                   diss_int(n) = diss_int_l + ( zv(n) - zu(k) ) / dzw(k) *      &
                                             ( diss_int_u - diss_int_l )
                 ENDIF