/[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.10 by cnh, Sun Sep 6 17:35:19 1998 UTC
# Line 72  C--   Initialise inverter Line 72  C--   Initialise inverter
72        sumRHSBuf(1,myThid) = 0. D0        sumRHSBuf(1,myThid) = 0. D0
73        etaNM1              = 1. D0        etaNM1              = 1. D0
74    
75    CcnhDebugStarts
76          _EXCH_XY_R8( cg2d_b, myThid )
77          CALL PLOT_FIELD_XYRL( cg2d_b, 'CG2D.0 CG2D_B' , 1, myThid )
78    CcnhDebugEnds
79    
80  C--   Normalise RHS  C--   Normalise RHS
81        rhsMax = 0. _d 0        rhsMax = 0. _d 0
82        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
# Line 104  C--   Update overlaps Line 109  C--   Update overlaps
109        _EXCH_XY_R8( cg2d_b, myThid )        _EXCH_XY_R8( cg2d_b, myThid )
110        _EXCH_XY_R8( cg2d_x, myThid )        _EXCH_XY_R8( cg2d_x, myThid )
111  CcnhDebugStarts  CcnhDebugStarts
112  C     CALL PLOT_FIELD_XYRL( cg2d_b, 'CG2D.1 CG2D_B' , 1, myThid )        CALL PLOT_FIELD_XYRL( cg2d_b, 'CG2D.1 CG2D_B' , 1, myThid )
113  CcnhDebugEnds  CcnhDebugEnds
114    
115  C--   Initial residual calculation  C--   Initial residual calculation
# Line 124  C--   Initial residual calculation Line 129  C--   Initial residual calculation
129       &    -aW2d(I+1,J  ,bi,bj)*cg2d_x(I  ,J  ,bi,bj)       &    -aW2d(I+1,J  ,bi,bj)*cg2d_x(I  ,J  ,bi,bj)
130       &    -aS2d(I  ,J  ,bi,bj)*cg2d_x(I  ,J  ,bi,bj)       &    -aS2d(I  ,J  ,bi,bj)*cg2d_x(I  ,J  ,bi,bj)
131       &    -aS2d(I  ,J+1,bi,bj)*cg2d_x(I  ,J  ,bi,bj)       &    -aS2d(I  ,J+1,bi,bj)*cg2d_x(I  ,J  ,bi,bj)
132       &    -freeSurfFac*_rA(i,j,bi,bj)*       &    -freeSurfFac*_rA(i,j,bi,bj)* horiVertRatio*
133       &     cg2d_x(I  ,J  ,bi,bj)/deltaTMom/deltaTMom*cg2dNorm       &     cg2d_x(I  ,J  ,bi,bj)/deltaTMom/deltaTMom*cg2dNorm
134       &    )       &    )
135            err            = err            +            err            = err            +
# Line 144  C--   Initial residual calculation Line 149  C--   Initial residual calculation
149  C     WRITE(6,*) ' mythid, err = ', mythid, SQRT(err)  C     WRITE(6,*) ' mythid, err = ', mythid, SQRT(err)
150        _GLOBAL_SUM_R8( errBuf    , err   , myThid )        _GLOBAL_SUM_R8( errBuf    , err   , myThid )
151        err    = errBuf(1,1)        err    = errBuf(1,1)
152  C     write(0,*) 'cg2d: Sum(rhs) = ',sumRHS        write(0,*) 'cg2d: Sum(rhs) = ',sumRHS
153    
154        actualIts      = 0        actualIts      = 0
155        actualResidual = SQRT(err)        actualResidual = SQRT(err)
# Line 225  C==    q = A.s Line 230  C==    q = A.s
230       &    -aW2d(I+1,J  ,bi,bj)*cg2d_s(I  ,J  ,bi,bj)       &    -aW2d(I+1,J  ,bi,bj)*cg2d_s(I  ,J  ,bi,bj)
231       &    -aS2d(I  ,J  ,bi,bj)*cg2d_s(I  ,J  ,bi,bj)       &    -aS2d(I  ,J  ,bi,bj)*cg2d_s(I  ,J  ,bi,bj)
232       &    -aS2d(I  ,J+1,bi,bj)*cg2d_s(I  ,J  ,bi,bj)       &    -aS2d(I  ,J+1,bi,bj)*cg2d_s(I  ,J  ,bi,bj)
233       &    -freeSurfFac*_rA(i,j,bi,bj)*       &    -freeSurfFac*_rA(i,j,bi,bj)* horiVertRatio*
234       &     cg2d_s(I  ,J  ,bi,bj)/deltaTMom/deltaTMom*cg2dNorm       &     cg2d_s(I  ,J  ,bi,bj)/deltaTMom/deltaTMom*cg2dNorm
235            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)
236            ENDDO            ENDDO
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22