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

Timestamp:

Nov 10, 2017 10:36:31 AM (6 years ago)

Author:

schwenkel

Message:

Modified particle box location and further changes in particle model

File:

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

palm/trunk/SOURCE/lpm_advec.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Changed particle box locations: center of particle box now coincides 
+! with scalar grid point of same index.
+! Renamed module and subroutines: lpm_pack_arrays_mod -> lpm_pack_and_sort_mod
+! lpm_pack_all_arrays -> lpm_sort_in_subboxes, lpm_pack_arrays -> lpm_pack
+! lpm_sort -> lpm_sort_timeloop_done 
+! 
+! 2417 2017-09-06 15:22:27Z suehring
 ! Particle loops adapted for sub-box structure, i.e. for each sub-box the 
 ! particle loop runs from start_index up to end_index instead from 1 to 
@@ -273 +280 @@
 
     DO  nb = 0, 7
-       
+!
+!--    Interpolate u velocity-component       
        i = ip
        j = jp + block_offset(nb)%j_off
        k = kp + block_offset(nb)%k_off
-!
-!--    Interpolate u velocity-component
+
        DO  n = start_index(nb), end_index(nb)
 !
@@ -289 +296 @@
 !--       First, check if particle is located below first vertical grid level 
 !--       above topography (Prandtl-layer height)
-          ilog = ( particles(n)%x + 0.5_wp * dx ) * ddx
-          jlog = ( particles(n)%y + 0.5_wp * dy ) * ddy
+          ilog = particles(n)%x * ddx
+          jlog = particles(n)%y * ddy
 !
 !--       Determine vertical index of topography top
@@ -378 +385 @@
 
        ENDDO
-
+!
+!--    Same procedure for interpolation of the v velocity-component 
        i = ip + block_offset(nb)%i_off
        j = jp
        k = kp + block_offset(nb)%k_off
-!
-!--    Same procedure for interpolation of the v velocity-component 
+
        DO  n = start_index(nb), end_index(nb)
 
-          ilog = ( particles(n)%x + 0.5_wp * dx ) * ddx
-          jlog = ( particles(n)%y + 0.5_wp * dy ) * ddy
+          ilog = particles(n)%x * ddx
+          jlog = particles(n)%y * ddy
 !
 !--       Determine vertical index of topography top
@@ -474 +481 @@
 
        ENDDO
-
+!
+!--    Same procedure for interpolation of the w velocity-component
        i = ip + block_offset(nb)%i_off
        j = jp + block_offset(nb)%j_off
        k = kp - 1
-!
-!--    Same procedure for interpolation of the w velocity-component
+
        DO  n = start_index(nb), end_index(nb)
 
@@ -660 +667 @@
                 i = particles(n)%x * ddx
                 j = particles(n)%y * ddy
-                k = ( zv(n) + 0.5_wp * dz * atmos_ocean_sign ) / dz  &
+                k = ( zv(n) + dz * atmos_ocean_sign ) / dz  &
                     + offset_ocean_nzt                      ! only exact if eq.dist
 !