66 |
_RL pF (1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL pF (1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
67 |
_RL del2u(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL del2u(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
68 |
_RL del2v(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL del2v(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
69 |
|
_RL tension(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
70 |
|
_RL strain(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
71 |
_RS hFacZ(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RS hFacZ(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
72 |
_RS r_hFacZ(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RS r_hFacZ(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
73 |
_RS xA(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RS xA(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
244 |
O uDiss,vDiss, |
O uDiss,vDiss, |
245 |
& myThid) |
& myThid) |
246 |
ENDIF |
ENDIF |
247 |
|
C or in terms of tension and strain |
248 |
|
IF (viscAstrain.NE.0. .OR. viscAtension.NE.0.) THEN |
249 |
|
CALL MOM_CALC_TENSION(bi,bj,k,uFld,vFld, |
250 |
|
O tension, |
251 |
|
I myThid) |
252 |
|
CALL MOM_CALC_STRAIN(bi,bj,k,uFld,vFld,hFacZ, |
253 |
|
O strain, |
254 |
|
I myThid) |
255 |
|
CALL MOM_HDISSIP(bi,bj,k, |
256 |
|
I tension,strain,hFacZ,viscAtension,viscAstrain, |
257 |
|
O uDiss,vDiss, |
258 |
|
I myThid) |
259 |
|
ENDIF |
260 |
ENDIF |
ENDIF |
261 |
|
|
262 |
C---- Zonal momentum equation starts here |
C---- Zonal momentum equation starts here |
493 |
& myCurrentTime-deltaTClock) |
& myCurrentTime-deltaTClock) |
494 |
& ) THEN |
& ) THEN |
495 |
CALL WRITE_LOCAL_RL('Ph','I10',Nr,phi_hyd,bi,bj,1,myIter,myThid) |
CALL WRITE_LOCAL_RL('Ph','I10',Nr,phi_hyd,bi,bj,1,myIter,myThid) |
496 |
|
CALL WRITE_LOCAL_RL('Ds','I10',1,strain,bi,bj,k,myIter,myThid) |
497 |
|
CALL WRITE_LOCAL_RL('Dt','I10',1,tension,bi,bj,k,myIter,myThid) |
498 |
CALL WRITE_LOCAL_RL('fV','I10',1,uCf,bi,bj,k,myIter,myThid) |
CALL WRITE_LOCAL_RL('fV','I10',1,uCf,bi,bj,k,myIter,myThid) |
499 |
CALL WRITE_LOCAL_RL('fU','I10',1,vCf,bi,bj,k,myIter,myThid) |
CALL WRITE_LOCAL_RL('fU','I10',1,vCf,bi,bj,k,myIter,myThid) |
500 |
CALL WRITE_LOCAL_RL('Du','I10',1,uDiss,bi,bj,k,myIter,myThid) |
CALL WRITE_LOCAL_RL('Du','I10',1,uDiss,bi,bj,k,myIter,myThid) |
501 |
CALL WRITE_LOCAL_RL('Dv','I10',1,vDiss,bi,bj,k,myIter,myThid) |
CALL WRITE_LOCAL_RL('Dv','I10',1,vDiss,bi,bj,k,myIter,myThid) |
502 |
|
CALL WRITE_LOCAL_RL('Z3','I10',1,vort3,bi,bj,k,myIter,myThid) |
503 |
|
CALL WRITE_LOCAL_RL('W3','I10',1,omega3,bi,bj,k,myIter,myThid) |
504 |
|
CALL WRITE_LOCAL_RL('KE','I10',1,KE,bi,bj,k,myIter,myThid) |
505 |
|
CALL WRITE_LOCAL_RL('D','I10',1,hdiv,bi,bj,k,myIter,myThid) |
506 |
ENDIF |
ENDIF |
507 |
|
|
508 |
RETURN |
RETURN |