[Mitgcm-support] rigidlid

[mitgcm-support at dev.mitgcm.org]

J-M

Doesn't the following code (calc_gt.F)

      ELSEIF ( rigidLid ) THEN
       DO j=1-Oly,sNy+Oly
        DO i=1-Olx,sNx+Olx
         wT(i,j) = rTrans(i,j)*
     &       (Tracer(i,j,k,bi,bj)+Tracer(i,j,kM1,bi,bj))*0.5 _d 0
        ENDDO
       ENDDO

do the same as

      ELSE
C-  include "free-surface correction" :
       DO j=1-Oly,sNy+Oly
        DO i=1-Olx,sNx+Olx
         wT(i,j) = rTrans(i,j)*(
     &      maskC(i,j,kM1,bi,bj)*
     &       (Tracer(i,j,k,bi,bj)+Tracer(i,j,kM1,bi,bj))*0.5 _d 0
     &    +(maskC(i,j,k,bi,bj)-maskC(i,j,kM1,bi,bj))*
     &        Tracer(i,j,k,bi,bj) )
        ENDDO
       ENDDO

because rTrans(i,j)=0 if rigidLid?

A.