/[MITgcm]/MITgcm/pkg/streamice/streamice_get_vel_resid_err_oad.F
ViewVC logotype

Diff of /MITgcm/pkg/streamice/streamice_get_vel_resid_err_oad.F

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

revision 1.1 by dgoldberg, Mon Mar 2 18:11:58 2015 UTC revision 1.2 by dgoldberg, Mon Mar 23 14:07:16 2015 UTC
# Line 29  C     === Global variables === Line 29  C     === Global variables ===
29  C     !INPUT/OUTPUT ARGUMENTS  C     !INPUT/OUTPUT ARGUMENTS
30        _RL err_max        _RL err_max
31        INTEGER myThid        INTEGER myThid
32          Real*8 u_dummy (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
33          Real*8 v_dummy (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
34    
35  #ifdef ALLOW_STREAMICE  #ifdef ALLOW_STREAMICE
36  #if (defined (ALLOW_OPENAD) && defined (ALLOW_STREAMICE_OAD_FP))  #if (defined (ALLOW_STREAMICE_OAD_FP))
37    
38        INTEGER conv_flag, i, j, bi, bj, myIter        INTEGER conv_flag, i, j, bi, bj, myIter
39    
40    
41        INTEGER ikey_nl        INTEGER ikey_nl
42        _RL err_tempu, err_tempv        _RL err_tempu, err_tempv
43        _RL max_vel, tempu, tempv, err_lastchange, cgtol        _RL max_vel, tempu, tempv, err_lastchange, cgtol
44        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
45        LOGICAL isTape        LOGICAL isTape
46    
47  !#ifdef ALLOW_AUTODIFF_TAMC        err_max = 0. _d 0
 !!$TAF STORE U_streamice = comlev1_stream_nl, key=ikey_nl  
 !!$TAF STORE V_streamice = comlev1_stream_nl, key=ikey_nl  
 !#endif  
   
         err_max = 0. _d 0  
48    
49          DO bj = myByLo(myThid), myByHi(myThid)          DO bj = myByLo(myThid), myByHi(myThid)
50           DO bi = myBxLo(myThid), myBxHi(myThid)           DO bi = myBxLo(myThid), myBxHi(myThid)
51            DO j=1,sNy            DO j=1-oly,sNy+oly
52             DO i=1,sNx             DO i=1-olx,sNx+olx
53              Au_SI (i,j,bi,bj) = 0. _d 0              Au_SI (i,j,bi,bj) = 0. _d 0
54              Av_SI (i,j,bi,bj) = 0. _d 0              Av_SI (i,j,bi,bj) = 0. _d 0
55                u_dummy (i,j,bi,bj) = u_streamice(i,j,bi,bj)
56                v_dummy (i,j,bi,bj) = v_streamice(i,j,bi,bj)
57              ubd_SI (i,j,bi,bj) = 0. _d 0              ubd_SI (i,j,bi,bj) = 0. _d 0
58              vbd_SI (i,j,bi,bj) = 0. _d 0              vbd_SI (i,j,bi,bj) = 0. _d 0
59             ENDDO             ENDDO
# Line 65  C     !INPUT/OUTPUT ARGUMENTS Line 65  C     !INPUT/OUTPUT ARGUMENTS
65       O    ubd_SI,       O    ubd_SI,
66       O    vbd_SI)       O    vbd_SI)
67    
 !#ifdef ALLOW_AUTODIFF_TAMC  
 !!$TAF STORE U_streamice = comlev1_stream_nl, key=ikey_nl  
 !!$TAF STORE V_streamice = comlev1_stream_nl, key=ikey_nl  
 !#endif  
   
68          CALL STREAMICE_CG_ACTION( myThid,          CALL STREAMICE_CG_ACTION( myThid,
69       O    Au_SI,       O    Au_SI,
70       O    Av_SI,       O    Av_SI,
71       I    U_streamice,       I    U_dummy,
72       I    V_streamice,       I    V_dummy,
73       I    0, sNx+1, 0, sNy+1 )       I    0, sNx+1, 0, sNy+1 )
74    
 !#ifdef ALLOW_AUTODIFF_TAMC  
 !!$TAF STORE U_streamice = comlev1_stream_nl, key=ikey_nl  
 !!$TAF STORE V_streamice = comlev1_stream_nl, key=ikey_nl  
 !#endif  
75    
76          DO bj = myByLo(myThid), myByHi(myThid)          DO bj = myByLo(myThid), myByHi(myThid)
77           DO bi = myBxLo(myThid), myBxHi(myThid)           DO bi = myBxLo(myThid), myBxHi(myThid)
# Line 106  C     !INPUT/OUTPUT ARGUMENTS Line 97  C     !INPUT/OUTPUT ARGUMENTS
97           ENDDO           ENDDO
98          ENDDO          ENDDO
99    
100    
101    
102          CALL GLOBAL_MAX_R8 (err_max, myThid)          CALL GLOBAL_MAX_R8 (err_max, myThid)
 !#ifdef ALLOW_AUTODIFF_TAMC  
 !!$TAF STORE err_max = comlev1_stream_nl, key=ikey_dynamics  
 !#endif  
103    
104  #endif  #endif
105  #endif  #endif
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22