Changeset 3065 for palm/trunk/SOURCE/lpm_init.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_init.f90 (modified) (4 diffs)

palm/trunk/SOURCE/lpm_init.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! dz was replaced by dzw or dz(1) to allow for right vertical stretching
+! 
+! 3049 2018-05-29 13:52:36Z Giersch
 ! Error messages revised
 ! 
@@ -794 +797 @@
                                ip = tmp_particle%x * ddx
                                jp = tmp_particle%y * ddy
-                               kp = tmp_particle%z / dz + 1 + offset_ocean_nzt                               
+                               kp = tmp_particle%z / dz(1) + 1 + offset_ocean_nzt                               
                                DO WHILE( zw(kp) < tmp_particle%z ) 
                                   kp = kp + 1
@@ -969 +972 @@
                                      pdz(particles(n)%group)
                       particles(n)%z = particles(n)%z +                        &
-                              MERGE( rand_contr, SIGN( dz, rand_contr ),       &
-                                     ABS( rand_contr ) < dz                    &
+                              MERGE( rand_contr, SIGN( dzw(kp), rand_contr ),  &
+                                     ABS( rand_contr ) < dzw(kp)               &
                                    )
                    ENDIF
@@ -987 +990 @@
                    i = particles(n)%x * ddx
                    j = particles(n)%y * ddy
-                   k = particles(n)%z / dz + 1 + offset_ocean_nzt
+                   k = particles(n)%z / dz(1) + 1 + offset_ocean_nzt
                    DO WHILE( zw(k) < particles(n)%z )
                       k = k + 1