/[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.4 by jmc, Wed Aug 24 02:23:41 2011 UTC revision 1.5 by jmc, Mon Mar 19 18:22:09 2012 UTC
# Line 91  CEOP Line 91  CEOP
91        _RS  pW  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS  pW  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
92        _RS  pS  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS  pS  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
93        _RS  pC  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS  pC  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
94        _RL  q2d(1-1:sNx+1,1-1:sNy+1,nSx,nSy)        _RL  q2d(0:sNx+1,0:sNy+1,nSx,nSy)
95  #ifdef DEBUG_DIAG_CG2D  #ifdef DEBUG_DIAG_CG2D
96        CHARACTER*(10) sufx        CHARACTER*(10) sufx
97        _RL  r2d(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL  r2d(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
98  #else  #else
99        _RL  r2d(1-1:sNx+1,1-1:sNy+1,nSx,nSy)        _RL  r2d(0:sNx+1,0:sNy+1,nSx,nSy)
100  #endif  #endif
101        _RL  s2d(1-1:sNx+1,1-1:sNy+1,nSx,nSy)        _RL  s2d(0:sNx+1,0: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  #ifdef ALLOW_DEBUG
# Line 109  C--   Set matrice main diagnonal: Line 109  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)  C-    Initialise overlap regions (in case EXCH call do not fill them all)
112          DO j = 1-Oly,sNy+Oly          DO j = 1-OLy,sNy+OLy
113           DO i = 1-Olx,sNx+Olx           DO i = 1-OLx,sNx+OLx
114             aC2d(i,j,bi,bj) = 0.             aC2d(i,j,bi,bj) = 0.
115           ENDDO           ENDDO
116          ENDDO          ENDDO
# Line 166  C--   Initialise inverter Line 166  C--   Initialise inverter
166  C--   Initial residual calculation  C--   Initial residual calculation
167        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
168         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
169          DO j=1-1,sNy+1          DO j=0,sNy+1
170           DO i=1-1,sNx+1           DO i=0,sNx+1
171              r2d(i,j,bi,bj) = 0.
172            s2d(i,j,bi,bj) = 0.            s2d(i,j,bi,bj) = 0.
173            x2dm(i,j,bi,bj) = x2d(i,j,bi,bj)            x2dm(i,j,bi,bj) = x2d(i,j,bi,bj)
174           ENDDO           ENDDO
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22