Changeset 1685 for palm/trunk/SOURCE/lpm_init.f90

Timestamp:

Oct 8, 2015 7:32:13 AM (9 years ago)

Author:

raasch

Message:

bugfix concerning vertical index calculation for particles in case of ocean runs

File:

• palm/trunk/SOURCE/lpm_init.f90 (modified) (3 diffs)

palm/trunk/SOURCE/lpm_init.f90

@@ -19 +19 @@
 ! Current revisions:
 ! -----------------
-! 
+! bugfix concerning vertical index offset in case of ocean
 ! 
 ! Former revisions:
@@ -128 +128 @@
                 number_of_initial_tails, number_of_particles,                  &
                 number_of_particle_groups, number_of_sublayers,                &
-                number_of_tails, offset_ocean_nzt, offset_ocean_nzt_m1,        & 
+                number_of_tails, offset_ocean_nzt, offset_ocean_nzt_m1,        &
                 particles, particle_advection_start, particle_groups,          &
                 particle_groups_type, particles_per_point,                     &
@@ -760 +760 @@
                             ip = ( tmp_particle%x + 0.5_wp * dx ) * ddx
                             jp = ( tmp_particle%y + 0.5_wp * dy ) * ddy
-                            kp = tmp_particle%z / dz + 1 + offset_ocean_nzt_m1
+                            kp = tmp_particle%z / dz + 1 + offset_ocean_nzt
 
                             IF ( seed_follows_topography )  THEN