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