/[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.37.4.1 by heimbach, Mon Apr 8 20:10:37 2002 UTC revision 1.41 by jmc, Mon Jan 13 01:31:22 2003 UTC
# Line 194  c       WRITE(*,*) ' mythid, sumArea = ' Line 194  c       WRITE(*,*) ' mythid, sumArea = '
194          _GLOBAL_SUM_R8( sumArea, myThid )          _GLOBAL_SUM_R8( sumArea, myThid )
195  C-  convert Target-Residual (in W unit) to cg2d-solver residual unit [m^2/s^2]  C-  convert Target-Residual (in W unit) to cg2d-solver residual unit [m^2/s^2]
196          cg2dTolerance = cg2dNorm * cg2dTargetResWunit          cg2dTolerance = cg2dNorm * cg2dTargetResWunit
197       &                           * sumArea / deltaTMom       &                           * sumArea / deltaTmom
198          WRITE(*,'(2A,1P2E22.14)') ' ini_cg2d: ',          _BEGIN_MASTER( myThid )
199            WRITE(standardMessageUnit,'(2A,1P2E22.14)') ' ini_cg2d: ',
200       &          'sumArea,cg2dTolerance =', sumArea,cg2dTolerance       &          'sumArea,cg2dTolerance =', sumArea,cg2dTolerance
201            _END_MASTER( myThid )
202        ENDIF        ENDIF
203            
204  C--   Initialise preconditioner  C--   Initialise preconditioner
# Line 222  C           defaults to 0.51 but can be Line 224  C           defaults to 0.51 but can be
224       &     aW2d(I,J,bi,bj) + aW2d(I+1,J  ,bi,bj)       &     aW2d(I,J,bi,bj) + aW2d(I+1,J  ,bi,bj)
225       &    +aS2d(I,J,bi,bj) + aS2d(I  ,J+1,bi,bj)       &    +aS2d(I,J,bi,bj) + aS2d(I  ,J+1,bi,bj)
226       &    +freeSurfFac*myNorm*recip_Bo(I,J,bi,bj)*       &    +freeSurfFac*myNorm*recip_Bo(I,J,bi,bj)*
227       &     rA(I,J,bi,bj)/deltaTMom/deltaTMom       &     rA(I,J,bi,bj)/deltaTMom/deltaTfreesurf
228       &    )       &    )
229            aCs = -(            aCs = -(
230       &     aW2d(I,J-1,bi,bj) + aW2d(I+1,J-1,bi,bj)       &     aW2d(I,J-1,bi,bj) + aW2d(I+1,J-1,bi,bj)
231       &    +aS2d(I,J-1,bi,bj) + aS2d(I  ,J  ,bi,bj)       &    +aS2d(I,J-1,bi,bj) + aS2d(I  ,J  ,bi,bj)
232       &    +freeSurfFac*myNorm*recip_Bo(I,J-1,bi,bj)*       &    +freeSurfFac*myNorm*recip_Bo(I,J-1,bi,bj)*
233       &     rA(I,J-1,bi,bj)/deltaTMom/deltaTMom       &     rA(I,J-1,bi,bj)/deltaTMom/deltaTfreesurf
234       &    )       &    )
235            aCw = -(            aCw = -(
236       &     aW2d(I-1,J,bi,bj) + aW2d(I  ,J  ,bi,bj)       &     aW2d(I-1,J,bi,bj) + aW2d(I  ,J  ,bi,bj)
237       &    +aS2d(I-1,J,bi,bj) + aS2d(I-1,J+1,bi,bj)       &    +aS2d(I-1,J,bi,bj) + aS2d(I-1,J+1,bi,bj)
238       &    +freeSurfFac*myNorm*recip_Bo(I-1,J,bi,bj)*       &    +freeSurfFac*myNorm*recip_Bo(I-1,J,bi,bj)*
239       &     rA(I-1,J,bi,bj)/deltaTMom/deltaTMom       &     rA(I-1,J,bi,bj)/deltaTMom/deltaTfreesurf
240       &    )       &    )
241            IF ( aC .EQ. 0. ) THEN            IF ( aC .EQ. 0. ) THEN
242              pC(I,J,bi,bj) = 1. _d 0              pC(I,J,bi,bj) = 1. _d 0
Legend:
Removed from v.1.37.4.1  
changed lines
  Added in v.1.41
  ViewVC Help
Powered by ViewVC 1.1.22