Ignore:
Timestamp:
Oct 26, 2007 1:11:14 PM (16 years ago)
Author:
raasch
Message:

bugfix in poisfft

File:
1 edited

Legend:

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

    r90 r128  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Bugfix: wavenumber calculation for even nx in routines maketri
    77!
    88! Former revisions:
     
    343343#if defined( __parallel )
    344344          DO  i = nxl_z, nxr_z
    345              IF ( i >= 0 .AND. i < nnxh )  THEN
     345             IF ( i >= 0 .AND. i <= nnxh )  THEN
    346346                ll(i) = 2.0 * ( 1.0 - COS( ( 2.0 * pi * i ) / &
    347                                           FLOAT( nx+1 ) ) ) / ( dx * dx ) + &
    348                         2.0 * ( 1.0 - COS( ( 2.0 * pi * j ) / &
    349                                           FLOAT( ny+1 ) ) ) / ( dy * dy )
    350              ELSEIF ( i == nnxh )  THEN
    351                 ll(i) = 2.0 * ( 1.0 - COS( ( 2.0 * pi * ( nx+1-i ) ) / &
    352347                                          FLOAT( nx+1 ) ) ) / ( dx * dx ) + &
    353348                        2.0 * ( 1.0 - COS( ( 2.0 * pi * j ) / &
     
    14131408!--       before it is Fourier-transformed
    14141409          DO  i = 0, nx
    1415              IF ( i >= 0 .AND. i < nnxh ) THEN
     1410             IF ( i >= 0 .AND. i <= nnxh ) THEN
    14161411                l(i) = 2.0 * ( 1.0 - COS( ( 2.0 * pi * i ) / &
    1417                                          FLOAT( nx+1 ) ) ) * ddx2 + &
    1418                        2.0 * ( 1.0 - COS( ( 2.0 * pi * j ) / &
    1419                                          FLOAT( ny+1 ) ) ) * ddy2
    1420              ELSEIF ( i == nnxh ) THEN
    1421                 l(i) = 2.0 * ( 1.0 - COS( ( 2.0 * pi * ( nx+1-i ) ) / &
    14221412                                         FLOAT( nx+1 ) ) ) * ddx2 + &
    14231413                       2.0 * ( 1.0 - COS( ( 2.0 * pi * j ) / &
Note: See TracChangeset for help on using the changeset viewer.