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

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

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

revision 1.23 by adcroft, Wed Dec 9 16:11:52 1998 UTC revision 1.24 by adcroft, Mon Mar 22 15:54:04 1999 UTC
# Line 126  CcnhDebugEnds Line 126  CcnhDebugEnds
126        _EXCH_XY_R4(aW2d, myThid)        _EXCH_XY_R4(aW2d, myThid)
127        _EXCH_XY_R4(aS2d, myThid)        _EXCH_XY_R4(aS2d, myThid)
128  CcnhDebugStarts  CcnhDebugStarts
129        CALL PLOT_FIELD_XYRS( aW2d, 'AW2D INI_CG2D.2' , 1, myThid )  C     CALL PLOT_FIELD_XYRS( aW2d, 'AW2D INI_CG2D.2' , 1, myThid )
130        CALL PLOT_FIELD_XYRS( aS2d, 'AS2D INI_CG2D.2' , 1, myThid )  C     CALL PLOT_FIELD_XYRS( aS2d, 'AS2D INI_CG2D.2' , 1, myThid )
131  CcnhDebugEnds  CcnhDebugEnds
132    
133  C--   Initialise preconditioner  C--   Initialise preconditioner
# Line 151  C           defaults to 0.51 but can be Line 151  C           defaults to 0.51 but can be
151            pC(I,J,bi,bj) = 1. _d 0            pC(I,J,bi,bj) = 1. _d 0
152            aC = -(            aC = -(
153       &     aW2d(I,J,bi,bj) + aW2d(I+1,J  ,bi,bj)       &     aW2d(I,J,bi,bj) + aW2d(I+1,J  ,bi,bj)
154       &    +aS2d(I,J,bi,bj) + aS2D(I  ,J+1,bi,bj)       &    +aS2d(I,J,bi,bj) + aS2d(I  ,J+1,bi,bj)
155       &    +freeSurfFac*myNorm* horiVertRatio*       &    +freeSurfFac*myNorm* horiVertRatio*
156       &     rA(I,J,bi,bj)/deltaTMom/deltaTMom       &     rA(I,J,bi,bj)/deltaTMom/deltaTMom
157       &    )       &    )
# Line 196  C--   Update overlap regions Line 196  C--   Update overlap regions
196        _EXCH_XY_R4(pW, myThid)        _EXCH_XY_R4(pW, myThid)
197        _EXCH_XY_R4(pS, myThid)        _EXCH_XY_R4(pS, myThid)
198  CcnhDebugStarts  CcnhDebugStarts
199        CALL PLOT_FIELD_XYRS( pC, 'pC   INI_CG2D.2' , 1, myThid )  C     CALL PLOT_FIELD_XYRS( pC, 'pC   INI_CG2D.2' , 1, myThid )
200        CALL PLOT_FIELD_XYRS( pW, 'pW   INI_CG2D.2' , 1, myThid )  C     CALL PLOT_FIELD_XYRS( pW, 'pW   INI_CG2D.2' , 1, myThid )
201        CALL PLOT_FIELD_XYRS( pS, 'pS   INI_CG2D.2' , 1, myThid )  C     CALL PLOT_FIELD_XYRS( pS, 'pS   INI_CG2D.2' , 1, myThid )
202  CcnhDebugEnds  CcnhDebugEnds
203    
204  C--   Set default values for initial guess and RHS  C--   Set default values for initial guess and RHS
Legend:
Removed from v.1.23  
changed lines
  Added in v.1.24
  ViewVC Help
Powered by ViewVC 1.1.22