/[MITgcm]/MITgcm/pkg/gmredi/gmredi_xtransport.F
ViewVC logotype

Diff of /MITgcm/pkg/gmredi/gmredi_xtransport.F

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

revision 1.9 by heimbach, Thu Nov 14 22:43:49 2002 UTC revision 1.10 by heimbach, Fri Jan 10 00:48:39 2003 UTC
# Line 90  C--   Area integrated zonal flux Line 90  C--   Area integrated zonal flux
90  #else  #else
91       &    *(GM_isopycK       &    *(GM_isopycK
92  #ifdef GM_VISBECK_VARIABLE_K  #ifdef GM_VISBECK_VARIABLE_K
93  cph     &    +0.5*(VisbeckK(i,j,bi,bj)+VisbeckK(i-1,j,bi,bj))  cph     &    +op5*(VisbeckK(i,j,bi,bj)+VisbeckK(i-1,j,bi,bj))
94  #endif  #endif
95       &    )       &    )
96  #endif /* GM_NON_UNITY_DIAGONAL */  #endif /* GM_NON_UNITY_DIAGONAL */
# Line 108  cph     &    +0.5*(VisbeckK(i,j,bi,bj)+V Line 108  cph     &    +0.5*(VisbeckK(i,j,bi,bj)+V
108          DO i=iMin,iMax          DO i=iMin,iMax
109    
110  C-      Vertical gradients interpolated to U points  C-      Vertical gradients interpolated to U points
111          dTdz =  0.5*(          dTdz =  op5*(
112       &   +0.5*recip_drC(k)*       &   +op5*recip_drC(k)*
113       &       ( maskC(i-1,j,k,bi,bj)*       &       ( maskC(i-1,j,k,bi,bj)*
114       &           (Tracer(i-1,j,km1,bi,bj)-Tracer(i-1,j,k,bi,bj))       &           (Tracer(i-1,j,km1,bi,bj)-Tracer(i-1,j,k,bi,bj))
115       &        +maskC( i ,j,k,bi,bj)*       &        +maskC( i ,j,k,bi,bj)*
116       &           (Tracer( i ,j,km1,bi,bj)-Tracer( i ,j,k,bi,bj))       &           (Tracer( i ,j,km1,bi,bj)-Tracer( i ,j,k,bi,bj))
117       &       )       &       )
118       &   +0.5*recip_drC(kp1)*       &   +op5*recip_drC(kp1)*
119       &       ( maskC(i-1,j,kp1,bi,bj)*       &       ( maskC(i-1,j,kp1,bi,bj)*
120       &           (Tracer(i-1,j,k,bi,bj)-Tracer(i-1,j,kp1,bi,bj))       &           (Tracer(i-1,j,k,bi,bj)-Tracer(i-1,j,kp1,bi,bj))
121       &        +maskC( i ,j,kp1,bi,bj)*       &        +maskC( i ,j,kp1,bi,bj)*
# Line 141  C-      Off-diagonal components of horiz Line 141  C-      Off-diagonal components of horiz
141       &                            -GM_PsiX(i,j,k,bi,bj) )       &                            -GM_PsiX(i,j,k,bi,bj) )
142       &                          *maskW(i,j,k,bi,bj)       &                          *maskW(i,j,k,bi,bj)
143           df(i,j) = df(i,j)           df(i,j) = df(i,j)
144       &    +uTrans*0.5*(Tracer(i,j,k,bi,bj)+Tracer(i-1,j,k,bi,bj))       &    +uTrans*op5*(Tracer(i,j,k,bi,bj)+Tracer(i-1,j,k,bi,bj))
145          ENDDO          ENDDO
146         ENDDO         ENDDO
147        ENDIF        ENDIF
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22