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

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

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

revision 1.2 by heimbach, Wed May 2 02:36:01 2012 UTC revision 1.3 by dgoldberg, Thu Jul 19 18:46:56 2012 UTC
# Line 94  C     LOCAL VARIABLES Line 94  C     LOCAL VARIABLES
94         ENDDO         ENDDO
95        ENDDO        ENDDO
96    
97    
98        CALL STREAMICE_CG_BOUND_VALS( myThid,            CALL STREAMICE_CG_BOUND_VALS( myThid,    
99       O    ubd_SI,       O    ubd_SI,
100       O    vbd_SI)       O    vbd_SI)
# Line 201  C       but after convergence loop does Line 202  C       but after convergence loop does
202  !$TAF STORE zv_old_si = comlev1_stream_nl, key=ikey_nl  !$TAF STORE zv_old_si = comlev1_stream_nl, key=ikey_nl
203  #endif  #endif
204    
205         CALL STREAMICE_CG_SOLVE(         DO bj = myByLo(myThid), myByHi(myThid)
206       &  U_streamice,          DO bi = myBxLo(myThid), myBxHi(myThid)
207       &  V_streamice,           DO j=1,sNy
208       &  taudx_SI,            DO i=1,sNx
209       &  taudy_SI,             U_streamice (i,j,bi,bj) = 0. _d 0
210       &  cgtol,             V_streamice (i,j,bi,bj) = 0. _d 0
211       &  cg_iters,            ENDDO
212       &  myThid )           ENDDO
213            ENDDO
214           ENDDO
215    
216           CALL STREAMICE_CG_WRAPPER(
217         &       U_streamice,
218         &       V_streamice,
219         &       taudx_SI,
220         &       taudy_SI,
221         &       cgtol,
222         &       cg_iters,
223         &       myThid )
224    
225                
226    
# Line 329  C--------------------------------------- Line 341  C---------------------------------------
341       & iter, " iterations"       & iter, " iterations"
342        endif        endif
343    
344    !        DO bj = myByLo(myThid), myByHi(myThid)
345    !         DO bi = myBxLo(myThid), myBxHi(myThid)
346    !          DO j=1,sNy
347    !           DO i=1,sNx
348    !            U_streamice (i,j,bi,bj) = 0. _d 0
349    !            V_streamice (i,j,bi,bj) = 0. _d 0
350    !           ENDDO
351    !          ENDDO
352    !         ENDDO
353    !        ENDDO
354    !
355    !        CALL STREAMICE_CG_WRAPPER(
356    !      &       U_streamice,
357    !      &       V_streamice,
358    !      &       taudx_SI,
359    !      &       taudy_SI,
360    !      &       cgtol,
361    !      &       cg_iters,
362    !      &       myThid )
363    
364    
365        _EXCH_XY_RL (U_streamice, myThid)        _EXCH_XY_RL (U_streamice, myThid)
366        _EXCH_XY_RL (V_streamice, myThid)        _EXCH_XY_RL (V_streamice, myThid)
367    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22