/[MITgcm]/MITgcm/pkg/diagnostics/diag_cg2d.F
ViewVC logotype

Diff of /MITgcm/pkg/diagnostics/diag_cg2d.F

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

revision 1.3 by jmc, Mon Aug 1 20:40:43 2011 UTC revision 1.4 by jmc, Wed Aug 24 02:23:41 2011 UTC
# Line 101  CEOP Line 101  CEOP
101        _RL  s2d(1-1:sNx+1,1-1:sNy+1,nSx,nSy)        _RL  s2d(1-1:sNx+1,1-1:sNy+1,nSx,nSy)
102        _RL  x2dm(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL  x2dm(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
103    
104    #ifdef ALLOW_DEBUG
105          IF (debugMode) CALL DEBUG_ENTER('DIAG_CG2D',myThid)
106    #endif
107    
108  C--   Set matrice main diagnonal:  C--   Set matrice main diagnonal:
109        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
110         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
111    C-    Initialise overlap regions (in case EXCH call do not fill them all)
112            DO j = 1-Oly,sNy+Oly
113             DO i = 1-Olx,sNx+Olx
114               aC2d(i,j,bi,bj) = 0.
115             ENDDO
116            ENDDO
117          DO j=1,sNy          DO j=1,sNy
118           DO i=1,sNx           DO i=1,sNx
119             aC2d(i,j,bi,bj) = -( ( aW2d(i,j,bi,bj)+aW2d(i+1,j,bi,bj) )             aC2d(i,j,bi,bj) = -( ( aW2d(i,j,bi,bj)+aW2d(i+1,j,bi,bj) )
# Line 353  C-    use the lowest residual solution ( Line 363  C-    use the lowest residual solution (
363        ENDIF        ENDIF
364  c     CALL EXCH_XY_RL( x2d, myThid )  c     CALL EXCH_XY_RL( x2d, myThid )
365    
366    #ifdef ALLOW_DEBUG
367          IF (debugMode) CALL DEBUG_LEAVE('DIAG_CG2D',myThid)
368    #endif
369    
370        RETURN        RETURN
371        END        END
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22