/[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.46 by jmc, Sun Sep 4 19:29:03 2005 UTC revision 1.47 by baylor, Fri Sep 16 19:32:20 2005 UTC
# Line 188  c     CALL MOM_VI_HFACZ_DISS(bi,bj,k,hFa Line 188  c     CALL MOM_VI_HFACZ_DISS(bi,bj,k,hFa
188    
189        CALL MOM_CALC_RELVORT3(bi,bj,k,uFld,vFld,hFacZ,vort3,myThid)        CALL MOM_CALC_RELVORT3(bi,bj,k,uFld,vFld,hFacZ,vort3,myThid)
190    
191          CALL MOM_CALC_TENSION(bi,bj,k,uFld,vFld,tension,myThid)
192    
193          CALL MOM_CALC_STRAIN(bi,bj,k,uFld,vFld,hFacZ,strain,myThid)
194    
195        IF (useAbsVorticity)        IF (useAbsVorticity)
196       & CALL MOM_CALC_ABSVORT3(bi,bj,k,vort3,omega3,myThid)       & CALL MOM_CALC_ABSVORT3(bi,bj,k,vort3,omega3,myThid)
197    
# Line 198  C      Calculate del^2 u and del^2 v for Line 202  C      Calculate del^2 u and del^2 v for
202       &     .OR. viscA4Grid.NE.0.       &     .OR. viscA4Grid.NE.0.
203       &     .OR. viscC4leith.NE.0.       &     .OR. viscC4leith.NE.0.
204       &     .OR. viscC4leithD.NE.0.       &     .OR. viscC4leithD.NE.0.
205         &     .OR. viscC4smag.NE.0.
206       &    ) THEN       &    ) THEN
207           CALL MOM_VI_DEL2UV(bi,bj,k,hDiv,vort3,hFacZ,           CALL MOM_VI_DEL2UV(bi,bj,k,hDiv,vort3,hFacZ,
208       O                      del2u,del2v,       O                      del2u,del2v,
209       &                      myThid)       &                      myThid)
          CALL MOM_CALC_HDIV(bi,bj,k,2,del2u,del2v,dStar,myThid)  
          CALL MOM_CALC_RELVORT3(  
      &                         bi,bj,k,del2u,del2v,hFacZ,zStar,myThid)  
210         ENDIF         ENDIF
211    
212  C      Calculate dissipation terms for U and V equations  C      Calculate dissipation terms for U and V equations
213  C      in terms of vorticity and divergence  
214         IF (    viscAhD.NE.0. .OR. viscAhZ.NE.0.  C      in terms of tension and strain
215       &    .OR. viscA4D.NE.0. .OR. viscA4Z.NE.0.         IF (useStrainTensionVisc) THEN
216       &    .OR. viscAhGrid.NE.0. .OR. viscA4Grid.NE.0.           CALL MOM_HDISSIP(bi,bj,k,hDiv,vort3,tension,strain,KE,
217       &    .OR. viscC2leith.NE.0. .OR. viscC4leith.NE.0.       I                    hFacZ,
      &    .OR. viscC2leithD.NE.0. .OR. viscC4leithD.NE.0.  
      &    ) THEN  
          CALL MOM_VI_HDISSIP(bi,bj,k,hDiv,vort3,hFacZ,dStar,zStar,  
      O                       guDiss,gvDiss,  
      &                       myThid)  
        ENDIF  
 C      or in terms of tension and strain  
        IF (viscAstrain.NE.0. .OR. viscAtension.NE.0.  
      O      .OR. viscC2smag.ne.0) THEN  
          CALL MOM_CALC_TENSION(bi,bj,k,uFld,vFld,  
      O                         tension,  
      I                         myThid)  
          CALL MOM_CALC_STRAIN(bi,bj,k,uFld,vFld,hFacZ,  
      O                        strain,  
      I                        myThid)  
          CALL MOM_HDISSIP(bi,bj,k,  
      I                    tension,strain,hFacZ,viscAtension,viscAstrain,  
218       O                    guDiss,gvDiss,       O                    guDiss,gvDiss,
219       I                    myThid)       I                    myThid)
220           ELSE
221    C      in terms of vorticity and divergence
222             CALL MOM_VI_HDISSIP(bi,bj,k,hDiv,vort3,tension,strain,KE,
223         I                       hFacZ,dStar,zStar,
224         O                       guDiss,gvDiss,
225         &                       myThid)        
226         ENDIF         ENDIF
227        ENDIF        ENDIF
228    
Legend:
Removed from v.1.46  
changed lines
  Added in v.1.47
  ViewVC Help
Powered by ViewVC 1.1.22