Changeset 1027 for palm


Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/advec_ws.f90

    r1020 r1027  
    44! Current revisions:
    55! ------------------
    6 !
     6! Bugfix in calculation indices k_mm, k_pp in accelerator version
    77!
    88! Former revisions:
     
    25252525                ibit6 = IBITS(wall_flags_0(k-1,j,i),6,1)
    25262526
    2527                 k_pp  = k + 2 * ( 1 - ibit6 )
    2528                 k_mm  = k - 2 * ( 1 - ibit6 )
     2527                k_pp  = k + 2 * ibit8
     2528                k_mm  = k - 2 * ( ibit7 + ibit8 )
    25292529                k_mmm = k - 3 * ibit8
    25302530
     
    32243224                ibit15 = IBITS(wall_flags_0(k-1,j,i),15,1)
    32253225
    3226                 k_pp  = k + 2 * ( 1 - ibit15 )
    3227                 k_mm  = k - 2 * ( 1 - ibit15 )
     3226                k_pp  = k + 2 * ibit17
     3227                k_mm  = k - 2 * ( ibit16 + ibit17 )
    32283228                k_mmm = k - 3 * ibit17
    32293229
     
    39413941                ibit24 = IBITS(wall_flags_0(k-1,j,i),24,1)
    39423942
    3943                 k_pp  = k + 2 * ( 1 - ibit24 )
    3944                 k_mm  = k - 2 * ( 1 - ibit24 )
     3943                k_pp  = k + 2 * ibit26
     3944                k_mm  = k - 2 * ( ibit25 + ibit26 )
    39453945                k_mmm = k - 3 * ibit26
    39463946
     
    46374637                ibit33 = IBITS(wall_flags_0(k-1,j,i),33,1)
    46384638
    4639                 k_pp  = k + 2 * ( 1 - ibit33 )
    4640                 k_mm  = k - 2 * ( 1 - ibit33 )
     4639                k_pp  = k + 2 * ibit35
     4640                k_mm  = k - 2 * ( ibit34 + ibit35 )
    46414641                k_mmm = k - 3 * ibit35
    46424642
Note: See TracChangeset for help on using the changeset viewer.