/[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.15 by jmc, Wed Feb 25 00:56:47 2004 UTC revision 1.19 by adcroft, Wed May 26 14:50:10 2004 UTC
# Line 233  C note (jmc) : Dissipation and Vort3 adv Line 233  C note (jmc) : Dissipation and Vort3 adv
233  C              use the same maskZ (and hFacZ)  => needs 2 call(s)  C              use the same maskZ (and hFacZ)  => needs 2 call(s)
234  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)
235    
236        CALL MOM_VI_CALC_KE(bi,bj,k,uFld,vFld,KE,myThid)        CALL MOM_CALC_KE(bi,bj,k,2,uFld,vFld,KE,myThid)
237    
238        CALL MOM_VI_CALC_HDIV(bi,bj,k,uFld,vFld,hDiv,myThid)        CALL MOM_CALC_HDIV(bi,bj,k,2,uFld,vFld,hDiv,myThid)
239    
240        CALL MOM_VI_CALC_RELVORT3(bi,bj,k,uFld,vFld,hFacZ,vort3,myThid)        CALL MOM_CALC_RELVORT3(bi,bj,k,uFld,vFld,hFacZ,vort3,myThid)
241    
242  c     CALL MOM_VI_CALC_ABSVORT3(bi,bj,k,vort3,omega3,myThid)  c     CALL MOM_CALC_ABSVORT3(bi,bj,k,vort3,omega3,myThid)
243    
244        IF (momViscosity) THEN        IF (momViscosity) THEN
245  C      Calculate del^2 u and del^2 v for bi-harmonic term  C      Calculate del^2 u and del^2 v for bi-harmonic term
246         IF (viscA4.NE.0. .OR. viscA4Grid.NE.0.) THEN         IF (viscA4.NE.0.
247         &     .OR. viscA4Grid.NE.0.
248         &     .OR. viscC4leith.NE.0.
249         &    ) THEN
250           CALL MOM_VI_DEL2UV(bi,bj,k,hDiv,vort3,hFacZ,           CALL MOM_VI_DEL2UV(bi,bj,k,hDiv,vort3,hFacZ,
251       O                      del2u,del2v,       O                      del2u,del2v,
252       &                      myThid)       &                      myThid)
253           CALL MOM_VI_CALC_HDIV(bi,bj,k,del2u,del2v,dStar,myThid)           CALL MOM_CALC_HDIV(bi,bj,k,2,del2u,del2v,dStar,myThid)
254           CALL MOM_VI_CALC_RELVORT3(           CALL MOM_CALC_RELVORT3(
255       &                         bi,bj,k,del2u,del2v,hFacZ,zStar,myThid)       &                         bi,bj,k,del2u,del2v,hFacZ,zStar,myThid)
256         ENDIF         ENDIF
257  C      Calculate dissipation terms for U and V equations  C      Calculate dissipation terms for U and V equations
258  C      in terms of vorticity and divergence  C      in terms of vorticity and divergence
259         IF (viscAh.NE.0. .OR. viscA4.NE.0. .OR.         IF (viscAh.NE.0. .OR. viscA4.NE.0.
260       &      viscAhGrid.NE.0. .OR. viscA4Grid.NE.0. ) THEN       &    .OR.  viscAhGrid.NE.0. .OR. viscA4Grid.NE.0.
261         &    .OR.  viscC2leith.NE.0. .OR. viscC4leith.NE.0.
262         &    ) THEN
263           CALL MOM_VI_HDISSIP(bi,bj,k,hDiv,vort3,hFacZ,dStar,zStar,           CALL MOM_VI_HDISSIP(bi,bj,k,hDiv,vort3,hFacZ,dStar,zStar,
264       O                       uDiss,vDiss,       O                       uDiss,vDiss,
265       &                       myThid)       &                       myThid)
Legend:
Removed from v.1.15  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.22