Changeset 2969 for palm/trunk/SOURCE

Timestamp:

Apr 13, 2018 11:55:09 AM (7 years ago)

Author:

thiele

Message:

Bugfix in Interpolation indices

File:

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

palm/trunk/SOURCE/lpm_advec.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix in Interpolation indices.
+! 
+! 2886 2018-03-14 11:51:53Z thiele
 ! Bugfix in passive particle SGS Model:
 ! Sometimes the added SGS velocities would lead to a violation of the CFL 
@@ -388 +391 @@
           ELSE
 
-             x  = xv(n) + ( 0.5_wp - i ) * dx
-             y  = yv(n) - j * dy
+             x  = xv(n) - i * dx
+             y  = yv(n) + ( 0.5_wp - j ) * dy
              aa = x**2          + y**2
              bb = ( dx - x )**2 + y**2
@@ -482 +485 @@
 
           ELSE
-             x  = xv(n) - i * dx
-             y  = yv(n) + ( 0.5_wp - j ) * dy
+             x  = xv(n) + ( 0.5_wp - i ) * dx
+             y  = yv(n) - j * dy
              aa = x**2          + y**2
              bb = ( dx - x )**2 + y**2
@@ -517 +520 @@
           IF ( vertical_particle_advection(particles(n)%group) )  THEN
 
-             x  = xv(n) - i * dx
-             y  = yv(n) - j * dy
+             x  = xv(n) + ( 0.5_wp - i ) * dx
+             y  = yv(n) + ( 0.5_wp - j ) * dy
              aa = x**2          + y**2
              bb = ( dx - x )**2 + y**2
@@ -588 +591 @@
 !
 !--             Interpolate TKE
-                x  = xv(n) - i * dx
-                y  = yv(n) - j * dy
+                x  = xv(n) + ( 0.5_wp - i ) * dx
+                y  = yv(n) + ( 0.5_wp - j ) * dy
                 aa = x**2          + y**2
                 bb = ( dx - x )**2 + y**2