/[MITgcm]/MITgcm/pkg/mom_vecinv/mom_vecinv.F
ViewVC logotype

Diff of /MITgcm/pkg/mom_vecinv/mom_vecinv.F

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

revision 1.7 by jmc, Sun Aug 3 02:49:35 2003 UTC revision 1.8 by heimbach, Fri Oct 10 23:00:01 2003 UTC
# Line 78  C     == Local variables == Line 78  C     == Local variables ==
78        _RS r_hFacZ(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RS r_hFacZ(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
79        _RS xA(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RS xA(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
80        _RS yA(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RS yA(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
       _RL uTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)  
       _RL vTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)  
81        _RL uFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL uFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
82        _RL vFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL vFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
83        _RL dStar(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL dStar(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
# Line 151  C     Initialise intermediate terms Line 149  C     Initialise intermediate terms
149          vort3(i,j) = 0.          vort3(i,j) = 0.
150          omega3(i,j) = 0.          omega3(i,j) = 0.
151          ke(i,j) = 0.          ke(i,j) = 0.
152    #ifdef ALLOW_AUTODIFF_TAMC
153            strain(i,j)  = 0. _d 0
154            tension(i,j) = 0. _d 0
155            fVerU(i,j,1) = 0. _d 0
156            fVerU(i,j,2) = 0. _d 0
157            fVerV(i,j,1) = 0. _d 0
158            fVerV(i,j,2) = 0. _d 0
159    #endif
160         ENDDO         ENDDO
161        ENDDO        ENDDO
162    
# Line 217  C     Make local copies of horizontal fl Line 223  C     Make local copies of horizontal fl
223         ENDDO         ENDDO
224        ENDDO        ENDDO
225    
 C     Calculate velocity field "volume transports" through tracer cell faces.  
       DO j=1-OLy,sNy+OLy  
        DO i=1-OLx,sNx+OLx  
         uTrans(i,j) = uFld(i,j)*xA(i,j)  
         vTrans(i,j) = vFld(i,j)*yA(i,j)  
        ENDDO  
       ENDDO  
   
226  C note (jmc) : Dissipation and Vort3 advection do not necesary  C note (jmc) : Dissipation and Vort3 advection do not necesary
227  C              use the same maskZ (and hFacZ)  => needs 2 call(s)  C              use the same maskZ (and hFacZ)  => needs 2 call(s)
228  c     CALL MOM_VI_HFACZ_DISS(bi,bj,k,hFacZ,r_hFacZ,myThid)  c     CALL MOM_VI_HFACZ_DISS(bi,bj,k,hFacZ,r_hFacZ,myThid)
# Line 310  C-     No-slip BCs impose a drag at wall Line 308  C-     No-slip BCs impose a drag at wall
308          ENDDO          ENDDO
309         ENDDO         ENDDO
310        ENDIF        ENDIF
311    
312  C-    No-slip BCs impose a drag at bottom  C-    No-slip BCs impose a drag at bottom
313        IF (momViscosity.AND.bottomDragTerms) THEN        IF (momViscosity.AND.bottomDragTerms) THEN
314         CALL MOM_U_BOTTOMDRAG(bi,bj,k,uFld,KE,KappaRU,vF,myThid)         CALL MOM_U_BOTTOMDRAG(bi,bj,k,uFld,KE,KappaRU,vF,myThid)
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22