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

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

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

revision 1.49 by jmc, Mon Dec 12 15:50:51 2005 UTC revision 1.50 by jmc, Mon Dec 19 21:10:34 2005 UTC
# Line 253  C     see CG2D.h for the interface to th Line 253  C     see CG2D.h for the interface to th
253        firstResidual=0.        firstResidual=0.
254        lastResidual=0.        lastResidual=0.
255        numIters=cg2dMaxIters        numIters=cg2dMaxIters
256    c     CALL TIMER_START('CG2D   [SOLVE_FOR_PRESSURE]',myThid)
257        CALL CG2D(        CALL CG2D(
258       U           cg2d_b,       U           cg2d_b,
259       U           cg2d_x,       U           cg2d_x,
# Line 261  C     see CG2D.h for the interface to th Line 262  C     see CG2D.h for the interface to th
262       U           numIters,       U           numIters,
263       I           myThid )       I           myThid )
264        _EXCH_XY_R8(cg2d_x, myThid )        _EXCH_XY_R8(cg2d_x, myThid )
265    c     CALL TIMER_STOP ('CG2D   [SOLVE_FOR_PRESSURE]',myThid)
266    
267  #ifdef ALLOW_DEBUG  #ifdef ALLOW_DEBUG
268        IF ( debugLevel .GE. debLevB ) THEN        IF ( debugLevel .GE. debLevB ) THEN
# Line 411  C-    end bi,bj loops Line 413  C-    end bi,bj loops
413        firstResidual=0.        firstResidual=0.
414        lastResidual=0.        lastResidual=0.
415        numIters=cg3dMaxIters        numIters=cg3dMaxIters
416          CALL TIMER_START('CG3D   [SOLVE_FOR_PRESSURE]',myThid)
417        CALL CG3D(        CALL CG3D(
418       U           cg3d_b,       U           cg3d_b,
419       U           phi_nh,       U           phi_nh,
# Line 419  C-    end bi,bj loops Line 422  C-    end bi,bj loops
422       U           numIters,       U           numIters,
423       I           myThid )       I           myThid )
424        _EXCH_XYZ_R8(phi_nh, myThid )        _EXCH_XYZ_R8(phi_nh, myThid )
425          CALL TIMER_STOP ('CG3D   [SOLVE_FOR_PRESSURE]',myThid)
426    
427        IF ( DIFFERENT_MULTIPLE(monitorFreq,myTime,deltaTClock)        IF ( DIFFERENT_MULTIPLE(monitorFreq,myTime,deltaTClock)
428       &   ) THEN       &   ) THEN
Legend:
Removed from v.1.49  
changed lines
  Added in v.1.50
  ViewVC Help
Powered by ViewVC 1.1.22