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

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

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

revision 1.3 by cnh, Thu May 21 18:26:36 1998 UTC revision 1.4 by cnh, Mon May 25 16:17:36 1998 UTC
# Line 31  C     === Global data === Line 31  C     === Global data ===
31  #include "SIZE.h"  #include "SIZE.h"
32  #include "EEPARAMS.h"  #include "EEPARAMS.h"
33  #include "PARAMS.h"  #include "PARAMS.h"
34    #include "GRID.h"
35  #include "CG2D.h"  #include "CG2D.h"
36    
37  C     === Routine arguments ===  C     === Routine arguments ===
# Line 122  C--   Initial residual calculation Line 123  C--   Initial residual calculation
123       &    -aW2d(I  ,J  ,bi,bj)*cg2d_x(I  ,J  ,bi,bj)       &    -aW2d(I  ,J  ,bi,bj)*cg2d_x(I  ,J  ,bi,bj)
124       &    -aW2d(I+1,J  ,bi,bj)*cg2d_x(I  ,J  ,bi,bj)       &    -aW2d(I+1,J  ,bi,bj)*cg2d_x(I  ,J  ,bi,bj)
125       &    -aS2d(I  ,J  ,bi,bj)*cg2d_x(I  ,J  ,bi,bj)       &    -aS2d(I  ,J  ,bi,bj)*cg2d_x(I  ,J  ,bi,bj)
126       &    -aS2d(I  ,J+1,bi,bj)*cg2d_x(I  ,J  ,bi,bj))       &    -aS2d(I  ,J+1,bi,bj)*cg2d_x(I  ,J  ,bi,bj)
127         &    -freeSurfFac*zA(i,j,bi,bj)*
128         &     cg2d_x(I  ,J  ,bi,bj)/deltaTMom/deltaTMom*cg2dNorm
129         &    )
130            err            = err            +            err            = err            +
131       &     cg2d_r(I,J,bi,bj)*cg2d_r(I,J,bi,bj)       &     cg2d_r(I,J,bi,bj)*cg2d_r(I,J,bi,bj)
132            sumRHS            = sumRHS            +            sumRHS            = sumRHS            +
# Line 169  C--    conjugate direction vector "s". Line 173  C--    conjugate direction vector "s".
173       &     +pW(I+1,J  ,bi,bj)*cg2d_r(I+1,J  ,bi,bj)       &     +pW(I+1,J  ,bi,bj)*cg2d_r(I+1,J  ,bi,bj)
174       &     +pS(I  ,J  ,bi,bj)*cg2d_r(I  ,J-1,bi,bj)       &     +pS(I  ,J  ,bi,bj)*cg2d_r(I  ,J-1,bi,bj)
175       &     +pS(I  ,J+1,bi,bj)*cg2d_r(I  ,J+1,bi,bj)       &     +pS(I  ,J+1,bi,bj)*cg2d_r(I  ,J+1,bi,bj)
176  C    &                        cg2d_r(I  ,J  ,bi,bj)  CcnhDebugStarts
177    C          cg2d_q(I,J,bi,bj) = cg2d_r(I  ,J  ,bi,bj)
178    CcnhDebugEnds
179             etaN = etaN             etaN = etaN
180       &     +cg2d_q(I,J,bi,bj)*cg2d_r(I,J,bi,bj)       &     +cg2d_q(I,J,bi,bj)*cg2d_r(I,J,bi,bj)
181            ENDDO            ENDDO
# Line 219  C==    q = A.s Line 225  C==    q = A.s
225       &    -aW2d(I+1,J  ,bi,bj)*cg2d_s(I  ,J  ,bi,bj)       &    -aW2d(I+1,J  ,bi,bj)*cg2d_s(I  ,J  ,bi,bj)
226       &    -aS2d(I  ,J  ,bi,bj)*cg2d_s(I  ,J  ,bi,bj)       &    -aS2d(I  ,J  ,bi,bj)*cg2d_s(I  ,J  ,bi,bj)
227       &    -aS2d(I  ,J+1,bi,bj)*cg2d_s(I  ,J  ,bi,bj)       &    -aS2d(I  ,J+1,bi,bj)*cg2d_s(I  ,J  ,bi,bj)
228         &    -freeSurfFac*zA(i,j,bi,bj)*
229         &     cg2d_s(I  ,J  ,bi,bj)/deltaTMom/deltaTMom*cg2dNorm
230            alpha = alpha+cg2d_s(I,J,bi,bj)*cg2d_q(I,J,bi,bj)            alpha = alpha+cg2d_s(I,J,bi,bj)*cg2d_q(I,J,bi,bj)
231            ENDDO            ENDDO
232           ENDDO           ENDDO
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22