/[MITgcm]/MITgcm/model/src/calc_gw.F
ViewVC logotype

Diff of /MITgcm/model/src/calc_gw.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.43 by jmc, Sat Jan 23 00:04:03 2010 UTC revision 1.45 by jmc, Tue May 3 19:26:03 2011 UTC
# Line 289  C     to compute d/dx(W), fill corners w Line 289  C     to compute d/dx(W), fill corners w
289  #ifdef COSINEMETH_III  #ifdef COSINEMETH_III
290       &              *sqCosFacU(j,bi,bj)       &              *sqCosFacU(j,bi,bj)
291  #endif  #endif
292    #ifdef ALLOW_OBCS
293         &              *maskInW(i,j,bi,bj)
294    #endif
295             ENDDO             ENDDO
296            ENDDO            ENDDO
297    
# Line 311  C     to compute d/dy(W), fill corners w Line 314  C     to compute d/dy(W), fill corners w
314       &              *sqCosFacV(j,bi,bj)       &              *sqCosFacV(j,bi,bj)
315  #endif  #endif
316  #endif  #endif
317    #ifdef ALLOW_OBCS
318         &              *maskInS(i,j,bi,bj)
319    #endif
320             ENDDO             ENDDO
321            ENDDO            ENDDO
322    
# Line 509  C     Tendency is minus divergence of ad Line 515  C     Tendency is minus divergence of ad
515  C     anelastic: all transports & advect. fluxes are scaled by rhoFac  C     anelastic: all transports & advect. fluxes are scaled by rhoFac
516            DO j=jMin,jMax            DO j=jMin,jMax
517             DO i=iMin,iMax             DO i=iMin,iMax
518    C to recover old (before 2009/11/30) results (since flxAdvUp(k=2) was zero)
519    c            IF (k.EQ.2) flxAdvUp(i,j) = 0.
520               gW(i,j,k,bi,bj) =               gW(i,j,k,bi,bj) =
521       &        -(   ( flx_EW(i+1,j)-flx_EW(i,j) )       &        -(   ( flx_EW(i+1,j)-flx_EW(i,j) )
522       &           + ( flx_NS(i,j+1)-flx_NS(i,j) )       &           + ( flx_NS(i,j+1)-flx_NS(i,j) )
Legend:
Removed from v.1.43  
changed lines
  Added in v.1.45
  ViewVC Help
Powered by ViewVC 1.1.22