192 |
C use the same maskZ (and hFacZ) => needs 2 call(s) |
C use the same maskZ (and hFacZ) => needs 2 call(s) |
193 |
c CALL MOM_VI_HFACZ_DISS(bi,bj,k,hFacZ,r_hFacZ,myThid) |
c CALL MOM_VI_HFACZ_DISS(bi,bj,k,hFacZ,r_hFacZ,myThid) |
194 |
|
|
195 |
CALL MOM_CALC_KE(bi,bj,k,0,uFld,vFld,KE,myThid) |
CALL MOM_CALC_KE(bi,bj,k,selectKEscheme,uFld,vFld,KE,myThid) |
196 |
|
|
197 |
CALL MOM_CALC_HDIV(bi,bj,k,2,uFld,vFld,hDiv,myThid) |
CALL MOM_CALC_HDIV(bi,bj,k,2,uFld,vFld,hDiv,myThid) |
198 |
|
|
199 |
CALL MOM_CALC_RELVORT3(bi,bj,k,uFld,vFld,hFacZ,vort3,myThid) |
CALL MOM_CALC_RELVORT3(bi,bj,k,uFld,vFld,hFacZ,vort3,myThid) |
200 |
|
|
|
CALL MOM_CALC_TENSION(bi,bj,k,uFld,vFld,tension,myThid) |
|
|
|
|
|
CALL MOM_CALC_STRAIN(bi,bj,k,uFld,vFld,hFacZ,strain,myThid) |
|
|
|
|
201 |
IF (useAbsVorticity) |
IF (useAbsVorticity) |
202 |
& CALL MOM_CALC_ABSVORT3(bi,bj,k,vort3,omega3,myThid) |
& CALL MOM_CALC_ABSVORT3(bi,bj,k,vort3,omega3,myThid) |
203 |
|
|
204 |
IF (momViscosity) THEN |
IF (momViscosity) THEN |
205 |
C Calculate Viscosities |
C Calculate Viscosities |
206 |
CALL MOM_CALC_VISC( |
CALL MOM_CALC_TENSION(bi,bj,k,uFld,vFld,tension,myThid) |
207 |
|
|
208 |
|
CALL MOM_CALC_STRAIN(bi,bj,k,uFld,vFld,hFacZ,strain,myThid) |
209 |
|
|
210 |
|
CALL MOM_CALC_VISC( |
211 |
I bi,bj,k, |
I bi,bj,k, |
212 |
O viscAh_Z,viscAh_D,viscA4_Z,viscA4_D, |
O viscAh_Z,viscAh_D,viscA4_Z,viscA4_D, |
213 |
O harmonic,biharmonic,useVariableViscosity, |
O harmonic,biharmonic,useVariableViscosity, |
611 |
CALL DIAGNOSTICS_FILL(KE, 'momKE ',k,1,2,bi,bj,myThid) |
CALL DIAGNOSTICS_FILL(KE, 'momKE ',k,1,2,bi,bj,myThid) |
612 |
CALL DIAGNOSTICS_FILL(hDiv, 'momHDiv ',k,1,2,bi,bj,myThid) |
CALL DIAGNOSTICS_FILL(hDiv, 'momHDiv ',k,1,2,bi,bj,myThid) |
613 |
CALL DIAGNOSTICS_FILL(vort3, 'momVort3',k,1,2,bi,bj,myThid) |
CALL DIAGNOSTICS_FILL(vort3, 'momVort3',k,1,2,bi,bj,myThid) |
|
CALL DIAGNOSTICS_FILL(strain,'Strain ',k,1,2,bi,bj,myThid) |
|
|
CALL DIAGNOSTICS_FILL(tension,'Tension ',k,1,2,bi,bj,myThid) |
|
614 |
CALL DIAGNOSTICS_FILL(gU(1-Olx,1-Oly,k,bi,bj), |
CALL DIAGNOSTICS_FILL(gU(1-Olx,1-Oly,k,bi,bj), |
615 |
& 'Um_Advec',k,1,2,bi,bj,myThid) |
& 'Um_Advec',k,1,2,bi,bj,myThid) |
616 |
CALL DIAGNOSTICS_FILL(gV(1-Olx,1-Oly,k,bi,bj), |
CALL DIAGNOSTICS_FILL(gV(1-Olx,1-Oly,k,bi,bj), |
617 |
& 'Vm_Advec',k,1,2,bi,bj,myThid) |
& 'Vm_Advec',k,1,2,bi,bj,myThid) |
618 |
IF (momViscosity) THEN |
IF (momViscosity) THEN |
619 |
|
CALL DIAGNOSTICS_FILL(strain,'Strain ',k,1,2,bi,bj,myThid) |
620 |
|
CALL DIAGNOSTICS_FILL(tension,'Tension ',k,1,2,bi,bj,myThid) |
621 |
CALL DIAGNOSTICS_FILL(guDiss,'Um_Diss ',k,1,2,bi,bj,myThid) |
CALL DIAGNOSTICS_FILL(guDiss,'Um_Diss ',k,1,2,bi,bj,myThid) |
622 |
CALL DIAGNOSTICS_FILL(gvDiss,'Vm_Diss ',k,1,2,bi,bj,myThid) |
CALL DIAGNOSTICS_FILL(gvDiss,'Vm_Diss ',k,1,2,bi,bj,myThid) |
623 |
ENDIF |
ENDIF |