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

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

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

revision 1.1 by dgoldberg, Thu Jul 19 18:52:31 2012 UTC revision 1.4 by dgoldberg, Thu Mar 7 15:23:19 2013 UTC
# Line 77  C     LOCAL VARIABLES Line 77  C     LOCAL VARIABLES
77        _RL dot_p1_tile (nSx,nSy)        _RL dot_p1_tile (nSx,nSy)
78        _RL dot_p2_tile (nSx,nSy)        _RL dot_p2_tile (nSx,nSy)
79        _RL ad_tolerance        _RL ad_tolerance
80          CHARACTER*(MAX_LEN_MBUF) msgBuf
81    
82  !       iters = streamice_max_cg_iter  !       iters = streamice_max_cg_iter
83    
84  #ifdef ALLOW_STREAMICE  #ifdef ALLOW_STREAMICE
85    
86          WRITE(msgBuf,'(A)') 'CALLING MANUAL CG ADJOINT'
87        PRINT *, "CALLING MANUAL CG ADJOINT"         CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
88         &                     SQUEEZE_RIGHT , 1)
89          
90    
91        conv_flag = 0        conv_flag = 0
92        ad_tolerance = 1.e-14        ad_tolerance = 1.e-14
# Line 105  C     LOCAL VARIABLES Line 108  C     LOCAL VARIABLES
108         ENDDO         ENDDO
109        ENDDO        ENDDO
110    
111    
112    
113        CALL STREAMICE_CG_SOLVE(        CALL STREAMICE_CG_SOLVE(
114       &  U_state,       &  U_state,
115       &  V_state,       &  V_state,
# Line 120  C     LOCAL VARIABLES Line 125  C     LOCAL VARIABLES
125    
126        tmpiter = 0        tmpiter = 0
127    
128          _EXCH_XY_RL( cg_Bu, myThid )
129          _EXCH_XY_RL( cg_Bv, myThid )
130    
131        CALL STREAMICE_CG_SOLVE(        CALL STREAMICE_CG_SOLVE(
132       &  cg_Uin,       &  cg_Uin,
133       &  cg_Vin,       &  cg_Vin,
# Line 133  C     LOCAL VARIABLES Line 141  C     LOCAL VARIABLES
141       &  tmpiter,       &  tmpiter,
142       &  myThid )       &  myThid )
143    
144          _EXCH_XY_RL( cg_Uin, myThid )
145          _EXCH_XY_RL( cg_Vin, myThid )
146          _EXCH_XY_RL( U_state, myThid )
147          _EXCH_XY_RL( V_state, myThid )  
148    
149  !       DO bj = myByLo(myThid), myByHi(myThid)  !       DO bj = myByLo(myThid), myByHi(myThid)
150  !        DO bi = myBxLo(myThid), myBxHi(myThid)  !        DO bi = myBxLo(myThid), myBxHi(myThid)
151  !         DO j=1,sNy  !         DO j=1,sNy
# Line 180  C     LOCAL VARIABLES Line 193  C     LOCAL VARIABLES
193  !         ENDDO  !         ENDDO
194  !  !
195  !  !
 ! !         print *, "GOT HERE 2"  
196  !  !
197  ! ! #else  ! ! #else
198  ! !  ! !
# Line 634  c     if iters has reached max_iters the Line 646  c     if iters has reached max_iters the
646         ENDDO         ENDDO
647        ENDDO        ENDDO
648    
         
649    
650        PRINT *, "DONE WITH MANUAL CG ADJOINT:",tmpiter,"ITERS"        WRITE(msgBuf,'(A,I5,A)') 'DONE WITH MANUAL CG ADJOINT:',tmpiter,
651         & 'ITERS'
652           CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
653         &                     SQUEEZE_RIGHT , 1)
654          
655    
656  #endif  #endif
657        RETURN        RETURN
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22