[MITgcm-support] zonal filter

[Gianmaria Sannino]

Hi there,

 

compiling the code with the check bound (on IBM pwr5) it appears an error on
the zonal filter; in particular

 

while ampFactorV is defined:

_RL ampFactorV( 1-Olx:sNx+Olx, 1-Oly:sNy+Oly, nSx, nSy )

 

 

It is used as:

 

C     o Set amplitude scale factor as function of latitude and mode number

      DO bj=1,nSy

       DO bi=1,nSx

        DO j=1-oLy,sNy+Oly

         ampFactor(1,J,bi,bj) = one

         ampFactorV(1,J,bi,bj) = one

         DO i=1,Nx/2-1

          ampFactor(2*I,J,bi,bj) = ampfact( yC(1,J,bi,bj) , I )

C         IF (ampFactor(2*I,J,bi,bj).LE..9) ampFactor(2*I,J,bi,bj)=0.

c-->          ampFactor(2*I+1,J,bi,bj) = ampFactor(2*I,J,bi,bj)

c-->         ampFactorV(2*I,J,bi,bj) = ampfact( yG(1,J,bi,bj) , I )

C         IF (ampFactorV(2*I,J,bi,bj).LE..9) ampFactorV(2*I,J,bi,bj)=0.

          ampFactorV(2*I+1,J,bi,bj) = ampFactorV(2*I,J,bi,bj)

         ENDDO

 

         I=Nx/2

         IF ( zonal_filt_mode2dx.EQ.0 ) THEN

           ampFactor(Nx,J,bi,bj) = ampfact( yC(1,J,bi,bj) , I )

c-->           ampFactorV(Nx,J,bi,bj) = ampfact( yG(1,J,bi,bj) , I )

         ELSE

           ampFactor(Nx,J,bi,bj) = 0.

c-->           ampFactorV(Nx,J,bi,bj) = 0.

         ENDIF

 

        ENDDO

 

 

Is it correct??

 

A presto

Gm

 

 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20061006/0324a48e/attachment.htm>