Changeset 1027 for palm

Timestamp:

Oct 15, 2012 5:18:39 PM (12 years ago)

Author:

suehring

Message:

Bugfix in calculation indices k_mm, k_pp in accelerator version

File:

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

palm/trunk/SOURCE/advec_ws.f90

@@ -4 +4 @@
 ! Current revisions:
 ! ------------------
-!
+! Bugfix in calculation indices k_mm, k_pp in accelerator version
 !
 ! Former revisions:
@@ -2525 +2525 @@
                 ibit6 = IBITS(wall_flags_0(k-1,j,i),6,1)
 
-                k_pp  = k + 2 * ( 1 - ibit6 )
-                k_mm  = k - 2 * ( 1 - ibit6 )
+                k_pp  = k + 2 * ibit8
+                k_mm  = k - 2 * ( ibit7 + ibit8 )
                 k_mmm = k - 3 * ibit8
 
@@ -3224 +3224 @@
                 ibit15 = IBITS(wall_flags_0(k-1,j,i),15,1)
 
-                k_pp  = k + 2 * ( 1 - ibit15 )
-                k_mm  = k - 2 * ( 1 - ibit15 )
+                k_pp  = k + 2 * ibit17
+                k_mm  = k - 2 * ( ibit16 + ibit17 )
                 k_mmm = k - 3 * ibit17
 
@@ -3941 +3941 @@
                 ibit24 = IBITS(wall_flags_0(k-1,j,i),24,1)
 
-                k_pp  = k + 2 * ( 1 - ibit24 )
-                k_mm  = k - 2 * ( 1 - ibit24 )
+                k_pp  = k + 2 * ibit26
+                k_mm  = k - 2 * ( ibit25 + ibit26 )
                 k_mmm = k - 3 * ibit26
 
@@ -4637 +4637 @@
                 ibit33 = IBITS(wall_flags_0(k-1,j,i),33,1)
 
-                k_pp  = k + 2 * ( 1 - ibit33 )
-                k_mm  = k - 2 * ( 1 - ibit33 )
+                k_pp  = k + 2 * ibit35
+                k_mm  = k - 2 * ( ibit34 + ibit35 )
                 k_mmm = k - 3 * ibit35