/[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.10 by dgoldberg, Thu Sep 27 20:29:01 2012 UTC
# Line 49  C     LOCAL VARIABLES Line 49  C     LOCAL VARIABLES
49  !       _RL taudx (1-OLx:sNx+OLx,1-OLy:sNy+Oly,nSx,nSy)  !       _RL taudx (1-OLx:sNx+OLx,1-OLy:sNy+Oly,nSx,nSy)
50  !       _RL taudy (1-OLx:sNx+OLx,1-OLy:sNy+Oly,nSx,nSy)  !       _RL taudy (1-OLx:sNx+OLx,1-OLy:sNy+Oly,nSx,nSy)
51    
52        CALL STREAMICE_DRIVING_STRESS (myThid)        IF (STREAMICE_ppm_driving_stress) THEN
53           CALL STREAMICE_DRIVING_STRESS_PPM (myThid)
54          ELSE
55           CALL STREAMICE_DRIVING_STRESS (myThid)
56          ENDIF
57    
58        cgtol = streamice_cg_tol        cgtol = streamice_cg_tol
59    
# Line 74  C     LOCAL VARIABLES Line 78  C     LOCAL VARIABLES
78  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
79  !$TAF STORE U_streamice = comlev1, key=ikey_dynamics  !$TAF STORE U_streamice = comlev1, key=ikey_dynamics
80  !$TAF STORE V_streamice = comlev1, key=ikey_dynamics  !$TAF STORE V_streamice = comlev1, key=ikey_dynamics
81    #ifdef STREAMICE_HYBRID_STRESS
82    !$TAF STORE streamice_taubx = comlev1, key=ikey_dynamics
83    !$TAF STORE streamice_tauby = comlev1, key=ikey_dynamics
84    !$TAF STORE visc_streamice_full = comlev1, key=ikey_dynamics
85    #endif
86    #endif
87    
88    #ifdef STREAMICE_HYBRID_STRESS
89           CALL STREAMICE_VISC_BETA_HYBRID ( myThid )
90    #else
91           CALL STREAMICE_VISC_BETA ( myThid )
92  #endif  #endif
93    
94        CALL STREAMICE_VISC_BETA ( myThid )  #ifdef STREAMICE_HYBRID_STRESS
95    !$TAF STORE visc_streamice_full = comlev1, key=ikey_dynamics
96    #endif
97    
98        _EXCH_XY_RL ( tau_beta_eff_streamice , myThid )        _EXCH_XY_RL ( tau_beta_eff_streamice , myThid )
99        _EXCH_XY_RL ( visc_streamice , myThid )        _EXCH_XY_RL ( visc_streamice , myThid )
# Line 126  C     LOCAL VARIABLES Line 143  C     LOCAL VARIABLES
143             err_tempu =             err_tempu =
144       &      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) -
145       &           taudx_SI(i,j,bi,bj))       &           taudx_SI(i,j,bi,bj))
146    !            print *, "err_temp_u", err_tempu
147            ENDIF            ENDIF
148            IF (STREAMICE_vmask(i,j,bi,bj).eq.1) THEN            IF (STREAMICE_vmask(i,j,bi,bj).eq.1) THEN
149             err_tempv = MAX( err_tempu,             err_tempv = MAX( err_tempu,
# Line 235  C       but after convergence loop does Line 253  C       but after convergence loop does
253       &       cg_iters,       &       cg_iters,
254       &       myThid )       &       myThid )
255    
256          #ifdef STREAMICE_HYBRID_STRESS
257    #ifdef ALLOW_AUTODIFF_TAMC
258    !$TAF STORE U_streamice = comlev1_stream_nl, key=ikey_nl
259    !$TAF STORE V_streamice = comlev1_stream_nl, key=ikey_nl
260    #endif
261    #endif
262    
263    #ifdef STREAMICE_HYBRID_STRESS
264            CALL STREAMICE_TAUB (myThid)
265    #endif
266    
267         WRITE(msgBuf,'(A,I5,A,I4,A)') 'streamice linear solve number',         WRITE(msgBuf,'(A,I5,A,I4,A)') 'streamice linear solve number',
268       &                       iter, " ",       &                       iter, " ",
# Line 247  C       but after convergence loop does Line 274  C       but after convergence loop does
274  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
275  !$TAF STORE U_streamice = comlev1_stream_nl, key=ikey_nl  !$TAF STORE U_streamice = comlev1_stream_nl, key=ikey_nl
276  !$TAF STORE V_streamice = comlev1_stream_nl, key=ikey_nl  !$TAF STORE V_streamice = comlev1_stream_nl, key=ikey_nl
277    #ifdef STREAMICE_HYBRID_STRESS
278    !$TAF STORE streamice_taubx = comlev1_stream_nl, key=ikey_nl
279    !$TAF STORE streamice_tauby = comlev1_stream_nl, key=ikey_nl
280    !$TAF STORE visc_streamice_full = comlev1_stream_nl, key=ikey_nl
281    #endif
282  #endif  #endif
283    
284    #ifdef STREAMICE_HYBRID_STRESS
285           CALL STREAMICE_VISC_BETA_HYBRID ( myThid )
286    #else
287         CALL STREAMICE_VISC_BETA ( myThid )         CALL STREAMICE_VISC_BETA ( myThid )
288    #endif
289    
290    
291  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
# Line 348  C       but after convergence loop does Line 384  C       but after convergence loop does
384         ENDDO         ENDDO
385    
386         CALL GLOBAL_MAX_R8 (err_max_fp, myThid)         CALL GLOBAL_MAX_R8 (err_max_fp, myThid)
387         print *, "err_max_fp", err_max_fp         WRITE(msgBuf,'(A,E15.7)') '||x_i-x_{i-1}||_inf',
388         &                       err_max_fp
389           CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
390         &                    SQUEEZE_RIGHT , 1)
391    
392         DO bj = myByLo(myThid), myByHi(myThid)         DO bj = myByLo(myThid), myByHi(myThid)
393          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
# Line 361  C       but after convergence loop does Line 400  C       but after convergence loop does
400          ENDDO          ENDDO
401         ENDDO         ENDDO
402    
403         WRITE(msgBuf,'(A,F11.7)') 'err/err_init',         WRITE(msgBuf,'(A,E15.7)') 'err/err_init',
404       &                       err_max/err_init       &                       err_max/err_init
405         CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,         CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
406       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22