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 |
|
|
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 |
|
|