/[MITgcm]/MITgcm_contrib/heimbach/OpenAD/code_shallow_openad/cg2d.F
ViewVC logotype

Diff of /MITgcm_contrib/heimbach/OpenAD/code_shallow_openad/cg2d.F

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

revision 1.4 by utke, Tue May 8 14:04:56 2007 UTC revision 1.5 by utke, Tue May 8 19:17:23 2007 UTC
# Line 92  c     INTEGER ks Line 92  c     INTEGER ks
92        _RL    sumRHS,sumRHStile        _RL    sumRHS,sumRHStile
93        _RL    rhsMax        _RL    rhsMax
94        _RL    rhsNorm        _RL    rhsNorm
95          CHARACTER*(MAX_LEN_MBUF) msgBuf
96  CEOP  CEOP
97    
98    
# Line 190  c    &    ) Line 191  c    &    )
191          err    = err    + errTile          err    = err    + errTile
192         ENDDO         ENDDO
193        ENDDO        ENDDO
 C     _EXCH_XY_R8( cg2d_r, myThid )  
194  #ifdef LETS_MAKE_JAM  #ifdef LETS_MAKE_JAM
195        CALL EXCH_XY_O1_R8_JAM( cg2d_r )        CALL EXCH_XY_O1_R8_JAM( cg2d_r )
196  #else  #else
197        CALL EXCH_XY_RL( cg2d_r, myThid )        CALL EXCH_XY_RL( cg2d_r, myThid )
198  #endif  #endif
 C     _EXCH_XY_R8( cg2d_s, myThid )  
199  #ifdef LETS_MAKE_JAM  #ifdef LETS_MAKE_JAM
200        CALL EXCH_XY_O1_R8_JAM( cg2d_s )        CALL EXCH_XY_O1_R8_JAM( cg2d_s )
201  #else  #else
# Line 346  C      Now compute "interior" points. Line 345  C      Now compute "interior" points.
345         err = SQRT(err)         err = SQRT(err)
346         actualIts      = it2d         actualIts      = it2d
347         actualResidual = err         actualResidual = err
348  C      _EXCH_XY_R8(cg2d_r, myThid )         IF ( debugLevel.GT.debLevB ) THEN
349    c       IF ( DIFFERENT_MULTIPLE(monitorFreq,myTime,deltaTClock)
350    c    &     ) THEN
351            _BEGIN_MASTER( myThid )
352            WRITE(msgBuf,'(A,I6,A,1PE21.14)')
353         &    ' cg2d: iter=', actualIts, ' ; resid.= ', actualResidual
354            CALL PRINT_MESSAGE(msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
355            _END_MASTER( myThid )
356    c       ENDIF
357           ENDIF
358  #ifdef LETS_MAKE_JAM  #ifdef LETS_MAKE_JAM
359        CALL EXCH_XY_O1_R8_JAM( cg2d_r )        CALL EXCH_XY_O1_R8_JAM( cg2d_r )
360  #else  #else
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22