Ignore:
Timestamp:
Apr 13, 2018 11:55:09 AM (4 years ago)
Author:
thiele
Message:

Bugfix in Interpolation indices

File:
1 edited

Legend:

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

    r2886 r2969  
    2525! -----------------
    2626! $Id$
     27! Bugfix in Interpolation indices.
     28!
     29! 2886 2018-03-14 11:51:53Z thiele
    2730! Bugfix in passive particle SGS Model:
    2831! Sometimes the added SGS velocities would lead to a violation of the CFL
     
    388391          ELSE
    389392
    390              x  = xv(n) + ( 0.5_wp - i ) * dx
    391              y  = yv(n) - j * dy
     393             x  = xv(n) - i * dx
     394             y  = yv(n) + ( 0.5_wp - j ) * dy
    392395             aa = x**2          + y**2
    393396             bb = ( dx - x )**2 + y**2
     
    482485
    483486          ELSE
    484              x  = xv(n) - i * dx
    485              y  = yv(n) + ( 0.5_wp - j ) * dy
     487             x  = xv(n) + ( 0.5_wp - i ) * dx
     488             y  = yv(n) - j * dy
    486489             aa = x**2          + y**2
    487490             bb = ( dx - x )**2 + y**2
     
    517520          IF ( vertical_particle_advection(particles(n)%group) )  THEN
    518521
    519              x  = xv(n) - i * dx
    520              y  = yv(n) - j * dy
     522             x  = xv(n) + ( 0.5_wp - i ) * dx
     523             y  = yv(n) + ( 0.5_wp - j ) * dy
    521524             aa = x**2          + y**2
    522525             bb = ( dx - x )**2 + y**2
     
    588591!
    589592!--             Interpolate TKE
    590                 x  = xv(n) - i * dx
    591                 y  = yv(n) - j * dy
     593                x  = xv(n) + ( 0.5_wp - i ) * dx
     594                y  = yv(n) + ( 0.5_wp - j ) * dy
    592595                aa = x**2          + y**2
    593596                bb = ( dx - x )**2 + y**2
Note: See TracChangeset for help on using the changeset viewer.