/[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.10 by dgoldberg, Thu Sep 27 20:29:01 2012 UTC revision 1.11 by dgoldberg, Wed Jan 9 21:56:18 2013 UTC
# Line 116  C     LOCAL VARIABLES Line 116  C     LOCAL VARIABLES
116       O    ubd_SI,       O    ubd_SI,
117       O    vbd_SI)       O    vbd_SI)
118    
119    !      CALL WRITE_FLD_XY_RL("u_bound_cont","",ubd_SI,0,myThid)
120    !      CALL WRITE_FLD_XY_RL("v_bound_cont","",vbd_SI,0,myThid)
121    !      CALL WRITE_FLD_XY_RL("taudx_u","",taudx_SI,0,myThid)
122    !      CALL WRITE_FLD_XY_RL("taudx_v","",taudy_SI,0,myThid)
123    
124  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
125  !$TAF STORE U_streamice = comlev1, key=ikey_dynamics  !$TAF STORE U_streamice = comlev1, key=ikey_dynamics
126  !$TAF STORE V_streamice = comlev1, key=ikey_dynamics  !$TAF STORE V_streamice = comlev1, key=ikey_dynamics
# Line 128  C     LOCAL VARIABLES Line 133  C     LOCAL VARIABLES
133       I    V_streamice,       I    V_streamice,
134       I    0, sNx+1, 0, sNy+1 )       I    0, sNx+1, 0, sNy+1 )
135    
136    
137    
138        err_init = 0. _d 0        err_init = 0. _d 0
139    
140        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
# Line 141  C     LOCAL VARIABLES Line 148  C     LOCAL VARIABLES
148            err_tempv = 0. _d 0            err_tempv = 0. _d 0
149            IF (STREAMICE_umask(i,j,bi,bj).eq.1) THEN            IF (STREAMICE_umask(i,j,bi,bj).eq.1) THEN
150             err_tempu =             err_tempu =
151       &      ABS (Au_SI(i,j,bi,bj)+ubd_SI(i,j,bi,bj) -       &      ABS (Au_SI(i,j,bi,bj)+0*ubd_SI(i,j,bi,bj) -
152       &           taudx_SI(i,j,bi,bj))       &           taudx_SI(i,j,bi,bj))
153  !            print *, "err_temp_u", err_tempu  !            print *, "err_temp_u", err_tempu
154            ENDIF            ENDIF
155            IF (STREAMICE_vmask(i,j,bi,bj).eq.1) THEN            IF (STREAMICE_vmask(i,j,bi,bj).eq.1) THEN
156             err_tempv = MAX( err_tempu,             err_tempv = MAX( err_tempu,
157       &      ABS (Av_SI(i,j,bi,bj)+vbd_SI(i,j,bi,bj) -       &      ABS (Av_SI(i,j,bi,bj)+0*vbd_SI(i,j,bi,bj) -
158       &           taudy_SI(i,j,bi,bj)))       &           taudy_SI(i,j,bi,bj)))
159            ENDIF            ENDIF
160            IF (err_tempv .ge. err_init) err_init = err_tempv            IF (err_tempv .ge. err_init) THEN
161                err_init = err_tempv
162              ENDIF
163           ENDDO           ENDDO
164          ENDDO          ENDDO
165         ENDDO         ENDDO
# Line 325  C       but after convergence loop does Line 334  C       but after convergence loop does
334       I    V_streamice,       I    V_streamice,
335       I    0, sNx+1, 0, sNy+1 )       I    0, sNx+1, 0, sNy+1 )
336    
337          if (iter .eq. streamice_max_nl_iter) then
338          CALL WRITE_FLD_XY_RL("u_bound_cont_A","",ubd_SI,0,myThid)
339          CALL WRITE_FLD_XY_RL("v_bound_cont_A","",vbd_SI,0,myThid)
340          CALL WRITE_FLD_XY_RL("u_bound_cont_B","",Au_SI,0,myThid)
341          CALL WRITE_FLD_XY_RL("v_bound_cont_B","",Av_SI,0,myThid)
342          endif
343    
344  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
345  !$TAF STORE U_streamice = comlev1_stream_nl, key=ikey_nl  !$TAF STORE U_streamice = comlev1_stream_nl, key=ikey_nl
346  !$TAF STORE V_streamice = comlev1_stream_nl, key=ikey_nl  !$TAF STORE V_streamice = comlev1_stream_nl, key=ikey_nl
# Line 347  C       but after convergence loop does Line 363  C       but after convergence loop does
363             err_tempv = 0. _d 0             err_tempv = 0. _d 0
364             IF (STREAMICE_umask(i,j,bi,bj).eq.1) THEN             IF (STREAMICE_umask(i,j,bi,bj).eq.1) THEN
365              err_tempu =              err_tempu =
366       &       ABS (Au_SI(i,j,bi,bj)+ubd_SI(i,j,bi,bj) -       &       ABS (Au_SI(i,j,bi,bj)+0*ubd_SI(i,j,bi,bj) -
367       &            taudx_SI(i,j,bi,bj))       &            taudx_SI(i,j,bi,bj))
368             ENDIF             ENDIF
369             IF (STREAMICE_vmask(i,j,bi,bj).eq.1) THEN             IF (STREAMICE_vmask(i,j,bi,bj).eq.1) THEN
370              err_tempv = MAX( err_tempu,              err_tempv = MAX( err_tempu,
371       &       ABS (Av_SI(i,j,bi,bj)+vbd_SI(i,j,bi,bj) -       &       ABS (Av_SI(i,j,bi,bj)+0*vbd_SI(i,j,bi,bj) -
372       &            taudy_SI(i,j,bi,bj)))       &            taudy_SI(i,j,bi,bj)))
373             ENDIF             ENDIF
374             IF (err_tempv .ge. err_max) err_max = err_tempv  !           if (err_tempu.ge.1.e2.or.err_tempv.ge.1.e2) THEN
375    !            print *, "FOUND MAX ", i,j,err_tempu,err_tempv,
376    !     &      ubd_SI(i,j,bi,bj),vbd_SI(i,j,bi,bj)
377    !           endif
378               IF (err_tempv .ge. err_max) THEN
379                err_max = err_tempv
380               ENDIF
381            ENDDO            ENDDO
382           ENDDO           ENDDO
383          ENDDO          ENDDO
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22