/[MITgcm]/MITgcm_contrib/dgoldberg/streamice/streamice_cg_solve.F
ViewVC logotype

Diff of /MITgcm_contrib/dgoldberg/streamice/streamice_cg_solve.F

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

revision 1.4 by dgoldberg, Thu Jul 19 18:46:56 2012 UTC revision 1.6 by dgoldberg, Thu Jul 26 16:22:58 2012 UTC
# Line 56  C     LOCAL VARIABLES Line 56  C     LOCAL VARIABLES
56        _RL dot_p1, dot_p2, alpha_k, beta_k, resid, resid_0        _RL dot_p1, dot_p2, alpha_k, beta_k, resid, resid_0
57        _RL dot_p1_tile (nSx,nSy)        _RL dot_p1_tile (nSx,nSy)
58        _RL dot_p2_tile (nSx,nSy)        _RL dot_p2_tile (nSx,nSy)
59          CHARACTER*(MAX_LEN_MBUF) msgBuf
60    
61        iters = streamice_max_cg_iter        iters = streamice_max_cg_iter
62    
# Line 151  C     FIND INITIAL RESIDUAL, and initial Line 152  C     FIND INITIAL RESIDUAL, and initial
152        CALL GLOBAL_SUM_TILE_RL( dot_p1_tile, dot_p1, myThid )        CALL GLOBAL_SUM_TILE_RL( dot_p1_tile, dot_p1, myThid )
153        resid_0 = sqrt(dot_p1)        resid_0 = sqrt(dot_p1)
154    
155          WRITE(msgBuf,'(A,E14.7)') 'CONJ GRAD INIT RESID, ',
156         &         resid_0
157           CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
158         &                     SQUEEZE_RIGHT , 1)
159    
160  C    CCCCCCCCCCCCCCCCCCCC  C    CCCCCCCCCCCCCCCCCCCC
161    
162        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
# Line 194  c  !!!!!!!!!!!!!!!!!! Line 200  c  !!!!!!!!!!!!!!!!!!
200        
201  c  ! initially, b-grid data is valid up to 3 halo nodes out -- right? (check for MITgcm!!)  c  ! initially, b-grid data is valid up to 3 halo nodes out -- right? (check for MITgcm!!)
202    
203        print *, "BEGINNING MAIN CG LOOP"        WRITE(msgBuf,'(A)') 'BEGINNING MAIN CG LOOP'
204           CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
205         &                     SQUEEZE_RIGHT , 1)
206    
207  !       IF(STREAMICE_construct_matrix) CALL STREAMICE_CG_MAKE_A(myThid)  !       IF(STREAMICE_construct_matrix) CALL STREAMICE_CG_MAKE_A(myThid)
208    
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22