/[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.8 by cnh, Fri Aug 21 12:29:55 1998 UTC revision 1.9 by cnh, Sun Sep 6 14:45:11 1998 UTC
# Line 124  C--   Initial residual calculation Line 124  C--   Initial residual calculation
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*_rA(i,j,bi,bj)*       &    -freeSurfFac*_rA(i,j,bi,bj)* Gravity*rhoConst*
128       &     cg2d_x(I  ,J  ,bi,bj)/deltaTMom/deltaTMom*cg2dNorm       &     cg2d_x(I  ,J  ,bi,bj)/deltaTMom/deltaTMom*cg2dNorm
129       &    )       &    )
130            err            = err            +            err            = err            +
# Line 144  C--   Initial residual calculation Line 144  C--   Initial residual calculation
144  C     WRITE(6,*) ' mythid, err = ', mythid, SQRT(err)  C     WRITE(6,*) ' mythid, err = ', mythid, SQRT(err)
145        _GLOBAL_SUM_R8( errBuf    , err   , myThid )        _GLOBAL_SUM_R8( errBuf    , err   , myThid )
146        err    = errBuf(1,1)        err    = errBuf(1,1)
147  C     write(0,*) 'cg2d: Sum(rhs) = ',sumRHS        write(0,*) 'cg2d: Sum(rhs) = ',sumRHS
148    
149        actualIts      = 0        actualIts      = 0
150        actualResidual = SQRT(err)        actualResidual = SQRT(err)
# Line 157  C     >>>>>>>>>>>>>>> BEGIN SOLVER <<<<< Line 157  C     >>>>>>>>>>>>>>> BEGIN SOLVER <<<<<
157        DO 10 it2d=1, cg2dMaxIters        DO 10 it2d=1, cg2dMaxIters
158    
159  CcnhDebugStarts  CcnhDebugStarts
160  C      WRITE(0,*) ' CG2D: Iteration ',it2d-1,' residual = ',actualResidual         WRITE(0,*) ' CG2D: Iteration ',it2d-1,' residual = ',actualResidual
161  CcnhDebugEnds  CcnhDebugEnds
162         IF ( err .LT. cg2dTargetResidual ) GOTO 11         IF ( err .LT. cg2dTargetResidual ) GOTO 11
163  C--    Solve preconditioning equation and update  C--    Solve preconditioning equation and update
# Line 225  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*_rA(i,j,bi,bj)*       &    -freeSurfFac*_rA(i,j,bi,bj)* Gravity*rhoConst*
229       &     cg2d_s(I  ,J  ,bi,bj)/deltaTMom/deltaTMom*cg2dNorm       &     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
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22