123 |
_RL rTransU(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL rTransU(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
124 |
_RL rTransV(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL rTransV(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
125 |
_RL KE(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL KE(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
126 |
_RL viscAhD(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
c _RL viscAh_D(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
127 |
_RL viscAhZ(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
c _RL viscAh_Z(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
128 |
_RL viscA4D(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
c _RL viscA4_D(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
129 |
_RL viscA4Z(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
c _RL viscA4_Z(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
130 |
_RL vort3(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
c _RL vort3(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
131 |
_RL hDiv(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
c _RL hDiv(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
132 |
_RL strain(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL strain(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
133 |
_RL tension(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL tension(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
134 |
C I,J,K - Loop counters |
C I,J,K - Loop counters |
265 |
|
|
266 |
CALL MOM_CALC_KE(bi,bj,k,3,uFld,vFld,KE,myThid) |
CALL MOM_CALC_KE(bi,bj,k,3,uFld,vFld,KE,myThid) |
267 |
|
|
268 |
c IF (viscAstrain.NE.0. .OR. viscAtension.NE.0.) THEN |
IF (viscAstrain.NE.0. .OR. viscAtension.NE.0.) THEN |
269 |
CALL MOM_CALC_TENSION(bi,bj,k,uFld,vFld, |
CALL MOM_CALC_TENSION(bi,bj,k,uFld,vFld, |
270 |
O tension, |
O tension, |
271 |
I myThid) |
I myThid) |
272 |
CALL MOM_CALC_STRAIN(bi,bj,k,uFld,vFld,hFacZ, |
CALL MOM_CALC_STRAIN(bi,bj,k,uFld,vFld,hFacZ, |
273 |
O strain, |
O strain, |
274 |
I myThid) |
I myThid) |
275 |
c ENDIF |
ENDIF |
276 |
|
|
277 |
C--- First call (k=1): compute vertical adv. flux fVerU(kUp) & fVerV(kUp) |
C--- First call (k=1): compute vertical adv. flux fVerU(kUp) & fVerV(kUp) |
278 |
IF (momAdvection.AND.k.EQ.1) THEN |
IF (momAdvection.AND.k.EQ.1) THEN |
311 |
c IF (momViscosity) THEN |
c IF (momViscosity) THEN |
312 |
c & CALL MOM_CALC_VISCOSITY(bi,bj,k, |
c & CALL MOM_CALC_VISCOSITY(bi,bj,k, |
313 |
c I uFld,vFld, |
c I uFld,vFld, |
314 |
c O viscAhD,viscAhZ,myThid) |
c O viscAh_D,viscAh_Z,myThid) |
315 |
|
|
316 |
C---- Zonal momentum equation starts here |
C---- Zonal momentum equation starts here |
317 |
|
|