/[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.7 by dgoldberg, Tue Sep 4 21:11:44 2012 UTC revision 1.8 by dgoldberg, Tue Sep 18 17:06:48 2012 UTC
# Line 74  C     LOCAL VARIABLES Line 74  C     LOCAL VARIABLES
74  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
75  !$TAF STORE U_streamice = comlev1, key=ikey_dynamics  !$TAF STORE U_streamice = comlev1, key=ikey_dynamics
76  !$TAF STORE V_streamice = comlev1, key=ikey_dynamics  !$TAF STORE V_streamice = comlev1, key=ikey_dynamics
77    #ifdef STREAMICE_HYBRID_STRESS
78    !$TAF STORE streamice_taubx = comlev1_stream_nl, key=ikey_nl
79    !$TAF STORE streamice_tauby = comlev1_stream_nl, key=ikey_nl
80    !$TAF STORE visc_streamice_full = comlev1_stream_nl, key=ikey_nl
81    #endif
82    #endif
83    
84    #ifdef STREAMICE_HYBRID_STRESS
85           CALL STREAMICE_VISC_BETA_HYBRID ( myThid )
86    #else
87           CALL STREAMICE_VISC_BETA ( myThid )
88  #endif  #endif
89    
90        CALL STREAMICE_VISC_BETA ( myThid )  #ifdef STREAMICE_HYBRID_STRESS
91    !$TAF STORE visc_streamice_full = comlev1_stream_nl, key=ikey_nl
92    #endif
93    
94        _EXCH_XY_RL ( tau_beta_eff_streamice , myThid )        _EXCH_XY_RL ( tau_beta_eff_streamice , myThid )
95        _EXCH_XY_RL ( visc_streamice , myThid )        _EXCH_XY_RL ( visc_streamice , myThid )
# Line 126  C     LOCAL VARIABLES Line 139  C     LOCAL VARIABLES
139             err_tempu =             err_tempu =
140       &      ABS (Au_SI(i,j,bi,bj)+ubd_SI(i,j,bi,bj) -       &      ABS (Au_SI(i,j,bi,bj)+ubd_SI(i,j,bi,bj) -
141       &           taudx_SI(i,j,bi,bj))       &           taudx_SI(i,j,bi,bj))
142    !            print *, "err_temp_u", err_tempu
143            ENDIF            ENDIF
144            IF (STREAMICE_vmask(i,j,bi,bj).eq.1) THEN            IF (STREAMICE_vmask(i,j,bi,bj).eq.1) THEN
145             err_tempv = MAX( err_tempu,             err_tempv = MAX( err_tempu,
# Line 235  C       but after convergence loop does Line 249  C       but after convergence loop does
249       &       cg_iters,       &       cg_iters,
250       &       myThid )       &       myThid )
251    
252          #ifdef STREAMICE_HYBRID_STRESS
253    #ifdef ALLOW_AUTODIFF_TAMC
254    !$TAF STORE U_streamice = comlev1_stream_nl, key=ikey_nl
255    !$TAF STORE V_streamice = comlev1_stream_nl, key=ikey_nl
256    #endif
257    #endif
258    
259    #ifdef STREAMICE_HYBRID_STRESS
260            CALL STREAMICE_TAUB (myThid)
261    #endif
262    
263         WRITE(msgBuf,'(A,I5,A,I4,A)') 'streamice linear solve number',         WRITE(msgBuf,'(A,I5,A,I4,A)') 'streamice linear solve number',
264       &                       iter, " ",       &                       iter, " ",
# Line 247  C       but after convergence loop does Line 270  C       but after convergence loop does
270  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
271  !$TAF STORE U_streamice = comlev1_stream_nl, key=ikey_nl  !$TAF STORE U_streamice = comlev1_stream_nl, key=ikey_nl
272  !$TAF STORE V_streamice = comlev1_stream_nl, key=ikey_nl  !$TAF STORE V_streamice = comlev1_stream_nl, key=ikey_nl
273    #ifdef STREAMICE_HYBRID_STRESS
274    !$TAF STORE streamice_taubx = comlev1_stream_nl, key=ikey_nl
275    !$TAF STORE streamice_tauby = comlev1_stream_nl, key=ikey_nl
276    !$TAF STORE visc_streamice_full = comlev1_stream_nl, key=ikey_nl
277    #endif
278  #endif  #endif
279    
280    #ifdef STREAMICE_HYBRID_STRESS
281           CALL STREAMICE_VISC_BETA_HYBRID ( myThid )
282    #else
283         CALL STREAMICE_VISC_BETA ( myThid )         CALL STREAMICE_VISC_BETA ( myThid )
284    #endif
285    
286    
287  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
# Line 348  C       but after convergence loop does Line 380  C       but after convergence loop does
380         ENDDO         ENDDO
381    
382         CALL GLOBAL_MAX_R8 (err_max_fp, myThid)         CALL GLOBAL_MAX_R8 (err_max_fp, myThid)
383         print *, "err_max_fp", err_max_fp         WRITE(msgBuf,'(A,E15.7)') '||x_i-x_{i-1}||_inf',
384         &                       err_max_fp
385           CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
386         &                    SQUEEZE_RIGHT , 1)
387    
388         DO bj = myByLo(myThid), myByHi(myThid)         DO bj = myByLo(myThid), myByHi(myThid)
389          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
# Line 361  C       but after convergence loop does Line 396  C       but after convergence loop does
396          ENDDO          ENDDO
397         ENDDO         ENDDO
398    
399         WRITE(msgBuf,'(A,F11.7)') 'err/err_init',         WRITE(msgBuf,'(A,E15.7)') 'err/err_init',
400       &                       err_max/err_init       &                       err_max/err_init
401         CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,         CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
402       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22