--- MITgcm/pkg/mom_vecinv/mom_vecinv.F 2004/01/25 00:31:52 1.13 +++ MITgcm/pkg/mom_vecinv/mom_vecinv.F 2004/02/07 23:15:47 1.14 @@ -1,4 +1,4 @@ -C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/mom_vecinv/mom_vecinv.F,v 1.13 2004/01/25 00:31:52 dimitri Exp $ +C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/mom_vecinv/mom_vecinv.F,v 1.14 2004/02/07 23:15:47 dimitri Exp $ C $Name: $ #include "PACKAGES_CONFIG.h" @@ -240,7 +240,7 @@ IF (momViscosity) THEN C Calculate del^2 u and del^2 v for bi-harmonic term - IF (viscA4.NE.0.) THEN + IF (viscA4.NE.0. .OR. viscA4Grid.NE.0.) THEN CALL MOM_VI_DEL2UV(bi,bj,k,hDiv,vort3,hFacZ, O del2u,del2v, & myThid) @@ -250,7 +250,8 @@ ENDIF C Calculate dissipation terms for U and V equations C in terms of vorticity and divergence - IF (viscAh.NE.0. .OR. viscA4.NE.0.) THEN + IF (viscAh.NE.0. .OR. viscA4.NE.0. .OR. + & viscAhGrid.NE.0. .OR. viscA4Grid.NE.0. ) THEN CALL MOM_VI_HDISSIP(bi,bj,k,hDiv,vort3,hFacZ,dStar,zStar, O uDiss,vDiss, & myThid)