Changeset 150 for palm/trunk/SOURCE/particle_boundary_conds.f90

Timestamp:

Feb 29, 2008 8:19:58 AM (16 years ago)

Author:

raasch

Message:

particle advection allowed for ocean runs

File:

• palm/trunk/SOURCE/particle_boundary_conds.f90 (modified) (8 diffs)

palm/trunk/SOURCE/particle_boundary_conds.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! Vertical index calculations adjusted for ocean runs.
 !
 ! Former revisions:
@@ -18 +18 @@
 !        reflection conditions is basically the same for all four cases, so it
 !        should be possible to further simplify/shorten it.
+! THIS ROUTINE HAS NOT BEEN TESTED FOR OCEAN RUNS SO FAR! (see offset_ocean_*)
 !------------------------------------------------------------------------------!
 
@@ -54 +55 @@
        i2 = ( particles(n)%x + 0.5 * dx ) * ddx
        j2 = ( particles(n)%y + 0.5 * dy ) * ddy
-       k2 = particles(n)%z / dz + 1
+       k2 = particles(n)%z / dz + 1 + offset_ocean_nzt_m1
 
        prt_x = particles(n)%x
@@ -69 +70 @@
           i1 = ( pos_x_old + 0.5 * dx ) * ddx
           j1 = ( pos_y_old + 0.5 * dy ) * ddy
-          k1 = pos_z_old / dz
+          k1 = pos_z_old / dz + offset_ocean_nzt_m1
 
 !
@@ -132 +133 @@
                 i3 = ( pos_x + 0.5 * dx ) * ddx   
                 j3 = ( pos_y + 0.5 * dy ) * ddy
-                k3 = pos_z / dz
+                k3 = pos_z / dz + offset_ocean_nzt_m1
 
                 i5 = pos_x * ddx
                 j5 = pos_y * ddy
-                k5 = pos_z / dz
+                k5 = pos_z / dz + offset_ocean_nzt_m1
 
                 IF ( k5 <= nzb_s_inner(j5,i3)  .AND. &
@@ -246 +247 @@
                 i3 = ( pos_x + 0.5 * dx ) * ddx
                 j3 = ( pos_y + 0.5 * dy ) * ddy
-                k3 = pos_z / dz
+                k3 = pos_z / dz + offset_ocean_nzt_m1
 
                 i5 = pos_x * ddx
                 j5 = pos_y * ddy
-                k5 = pos_z / dz
+                k5 = pos_z / dz + offset_ocean_nzt_m1
 
                 IF ( k5 <= nzb_s_inner(j3,i5)  .AND. &
@@ -361 +362 @@
                 i3 = ( pos_x + 0.5 * dx ) * ddx
                 j3 = ( pos_y + 0.5 * dy ) * ddy
-                k3 =  pos_z / dz
+                k3 = pos_z / dz + offset_ocean_nzt_m1
 
                 i5 = pos_x * ddx
                 j5 = pos_y * ddy
-                k5 = pos_z / dz
+                k5 = pos_z / dz + offset_ocean_nzt_m1
 
                 IF ( k5 <= nzb_s_inner(j5,i3)  .AND. &
@@ -476 +477 @@
                 i3 = ( pos_x + 0.5 * dx ) * ddx   
                 j3 = ( pos_y + 0.5 * dy ) * ddy
-                k3 = pos_z / dz
+                k3 = pos_z / dz + offset_ocean_nzt_m1
 
                 i5 = pos_x * ddx
                 j5 = pos_y * ddy
-                k5 = pos_z / dz
+                k5 = pos_z / dz + offset_ocean_nzt_m1
 
                 IF ( k3 <= nzb_s_inner(j3,i3)  .AND. &