Changeset 2628 for palm/trunk/SOURCE/lpm_pack_arrays.f90

Timestamp:

Nov 20, 2017 12:40:38 PM (6 years ago)

Author:

schwenkel

Message:

enable particle advection with grid stretching and some formatation changes in lpm

File:

• palm/trunk/SOURCE/lpm_pack_arrays.f90 (modified) (4 diffs)

palm/trunk/SOURCE/lpm_pack_arrays.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Enabled particle advection with grid stretching.
+! 
+! 2609 2017-11-14 14:14:44Z schwenkel
 ! Integrated subroutine pack_and_sort into lpm_sort_in_subboxes
 ! 
@@ -112 +115 @@
 
        USE cpulog,                                                              &
-           ONLY:  cpu_log, log_point_s
+          ONLY:  cpu_log, log_point_s
 
        USE indices,                                                             &
-           ONLY:  nxl, nxr, nys, nyn, nzb, nzt
+          ONLY:  nxl, nxr, nys, nyn, nzb, nzt
 
        USE kinds
 
        USE control_parameters,                                                  &
-           ONLY: dz
+          ONLY: dz
        
        USE grid_variables,                                                      &
-           ONLY: dx,dy,ddx, ddy
-          
+          ONLY: dx,dy,ddx, ddy
+           
+       USE arrays_3d,                                                           &
+          ONLY:  zu
        IMPLICIT NONE
 
@@ -167 +172 @@
                       i = (particles(n)%x + 0.5_wp * dx) * ddx
                       j = (particles(n)%y + 0.5_wp * dy) * ddy
-                      k = ( particles(n)%z+ 0.5_wp *dz ) / dz+1 + offset_ocean_nzt
                       
                       IF ( i == ip )  sort_index = sort_index + 4 
-                      IF ( j == jp )  sort_index = sort_index + 2
-                      IF ( k == kp )  sort_index = sort_index + 1
+                      IF ( j == jp )  sort_index = sort_index + 2                      
+                      IF ( zu(kp) > particles(n)%z ) sort_index = sort_index + 1
                       
                       sort_count(sort_index) = sort_count(sort_index) + 1
@@ -343 +347 @@
 
 
- END module lpm_pack_and_sort_mod
+ END MODULE lpm_pack_and_sort_mod