--- MITgcm/pkg/mom_vecinv/mom_vecinv.F 2005/07/30 22:05:36 1.43 +++ MITgcm/pkg/mom_vecinv/mom_vecinv.F 2005/08/24 23:16:05 1.45 @@ -1,4 +1,4 @@ -C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/mom_vecinv/mom_vecinv.F,v 1.43 2005/07/30 22:05:36 jmc Exp $ +C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/mom_vecinv/mom_vecinv.F,v 1.45 2005/08/24 23:16:05 jmc Exp $ C $Name: $ #include "MOM_VECINV_OPTIONS.h" @@ -243,9 +243,7 @@ C Eddy component of vertical flux (interior component only) -> vrF IF (momViscosity.AND..NOT.implicitViscosity) THEN - IF ( k.LT.Nr ) THEN - CALL MOM_U_RVISCFLUX(bi,bj,k+1,uVel,KappaRU,vrF,myThid) - ENDIF + CALL MOM_U_RVISCFLUX(bi,bj,k+1,uVel,KappaRU,vrF,myThid) C Combine fluxes DO j=jMin,jMax @@ -305,9 +303,7 @@ C Eddy component of vertical flux (interior component only) -> vrF IF (momViscosity.AND..NOT.implicitViscosity) THEN - IF ( k.LT.Nr ) THEN - CALL MOM_V_RVISCFLUX(bi,bj,k+1,vVel,KappaRV,vrF,myThid) - ENDIF + CALL MOM_V_RVISCFLUX(bi,bj,k+1,vVel,KappaRV,vrF,myThid) C Combine fluxes -> fVerV DO j=jMin,jMax @@ -460,14 +456,12 @@ ENDIF #ifdef ALLOW_TIMEAVE -#ifndef MINIMAL_TAVE_OUTPUT IF (taveFreq.GT.0.) THEN CALL TIMEAVE_CUMUL_1K1T(uZetatave,vCf,deltaTClock, & Nr, k, bi, bj, myThid) CALL TIMEAVE_CUMUL_1K1T(vZetatave,uCf,deltaTClock, & Nr, k, bi, bj, myThid) ENDIF -#endif /* ndef MINIMAL_TAVE_OUTPUT */ #endif /* ALLOW_TIMEAVE */ C-- Vertical shear terms (-w*du/dr & -w*dv/dr)