/[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.22 by adcroft, Tue May 29 14:01:37 2001 UTC revision 1.23 by adcroft, Wed Jun 6 14:55:45 2001 UTC
# Line 135  C Western boundary Line 135  C Western boundary
135         ENDDO         ENDDO
136        ENDDO        ENDDO
137    
138    #ifndef EXCLUDE_DEBUGMODE
139           CALL DEBUG_STATS_RL(1,cg2d_b,'cg2d_b (SOLVE_FOR_PRESSURE)',
140         &                        myThid)
141    #endif
142    
143  C--   Find the surface pressure using a two-dimensional conjugate  C--   Find the surface pressure using a two-dimensional conjugate
144  C--   gradient solver.  C--   gradient solver.
# Line 151  C     see CG2D.h for the interface to th Line 155  C     see CG2D.h for the interface to th
155       I           myThid )       I           myThid )
156        _EXCH_XY_R8(cg2d_x, myThid )        _EXCH_XY_R8(cg2d_x, myThid )
157    
158    #ifndef EXCLUDE_DEBUGMODE
159           CALL DEBUG_STATS_RL(1,cg2d_x,'cg2d_x (SOLVE_FOR_PRESSURE)',
160         &                        myThid)
161    #endif
162    
163        _BEGIN_MASTER( myThid )        _BEGIN_MASTER( myThid )
164         WRITE(*,'(A,I6,1PE30.14)') ' CG2D iters, err = ',         WRITE(*,'(A,I6,1PE30.14)') ' CG2D iters, err = ',
165       &                               0, firstResidual       &                               0, firstResidual
Legend:
Removed from v.1.22  
changed lines
  Added in v.1.23
  ViewVC Help
Powered by ViewVC 1.1.22