[MITgcm-devel] mom_vi_hdissip

[Patrick Heimbach]

Hi Baylor,

I am wondering whether in mom_vi_hdissip
for biharmonic friction, both alternatives

#ifdef MOM_VI_ORIGINAL_VISCA4
          Dim=dyF( i ,j-1,bi,bj)*dStar( i ,j-1)
          Dij=dyF( i , j ,bi,bj)*dStar( i , j )
          Dmj=dyF(i-1, j ,bi,bj)*dStar(i-1, j )
          Zip=dxV( i ,j+1,bi,bj)*hFacZ( i ,j+1)*zStar( i ,j+1)
          Zij=dxV( i , j ,bi,bj)*hFacZ( i , j )*zStar( i , j )
          Zpj=dxV(i+1, j ,bi,bj)*hFacZ(i+1, j )*zStar(i+1, j )
#else
          Dim=dStar( i ,j-1)
          Dij=dStar( i , j )
          Dmj=dStar(i-1, j )
          Zip=hFacZ( i ,j+1)*zStar( i ,j+1)
          Zij=hFacZ( i , j )*zStar( i , j )
          Zpj=hFacZ(i+1, j )*zStar(i+1, j )
#endif

are still in use, or whether one has become obsolete now
(maybe there just for historical/backward compatibility reasons?).

Cheers
-Patrick

---
Dr Patrick Heimbach | heimbach at mit.edu | http://www.mit.edu/~heimbach
MIT | EAPS, 54-1518 | 77 Massachusetts Ave | Cambridge, MA 02139, USA
FON: +1-617-253-5259 | FAX: +1-617-253-4464 | SKYPE: patrick.heimbach