Changeset 3065 for palm/trunk/SOURCE/lpm_set_attributes.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_set_attributes.f90 (modified) (5 diffs)

palm/trunk/SOURCE/lpm_set_attributes.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! dz was replaced by dzw to allow for right vertical stretching
+! 
+! 2718 2018-01-02 08:49:38Z maronga
 ! Corrected "Former revisions" section
 ! 
@@ -80 +83 @@
 
     USE arrays_3d,                                                             &
-        ONLY:  pt, u, v, w, zu, zw
+        ONLY:  dzw, pt, u, v, w, zu, zw
 
     USE control_parameters,                                                    &
-        ONLY:  u_gtrans, v_gtrans, dz
+        ONLY:  u_gtrans, v_gtrans
 
     USE cpulog,                                                                &
@@ -207 +210 @@
                                      ( 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
@@ -240 +243 @@
                                       ( 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
@@ -344 +347 @@
                                  ) / ( 3.0_wp * gg ) - sums(k,4)
 
-                      pt_int = pt_int_l + ( zv(n) - zu(k) ) / dz *    &
+                      pt_int = pt_int_l + ( zv(n) - zu(k) ) / dzw(k) *          &
                                           ( pt_int_u - pt_int_l )