[Mitgcm-support] advection code: calc_gt/gs

[mitgcm-support at dev.mitgcm.org]

Alistair,

short question: There are lines in calc_gs/gt like

#ifdef THIRD_ORDER
           SxP=(salt(i,j,k,bi,bj)-salt(i-1,j,k,bi,bj))*
     &          _maskW(i,j,k,bi,bj)
           SxM=(salt(i-1,j,k,bi,bj)-salt(i-2,j,k,bi,bj))*
     &          _maskW(i-1,j,k,bi,bj)
           af(i,j)= uTrans(i,j)*( -0.125*(SxP-SxM)
     &         0.5*( salt(i,j,k,bi,bj)+salt(i-1,j,k,bi,bj) ) )

The last line is wrong. Should it read "+0.5...", "*0.5..." or what?

Any new code on calc_mom_rhs.F? I suspect that momentum is the main part
that needs the high values for biharmonic option.

Thanks,
	Arne
-- 

Dr. Arne Biastoch
Physical Oceanography Research Group       phone: +1-858-822-3787
Scripps Institution of Oceanography        fax  : +1-858-534-9820
MS: 0230                                   email: abiastoch at ucsd.edu
8605 La Jolla Shores Dr.            
La Jolla, CA 92093-0230, U.S.A.   http://www.ecco.ucsd.edu/~biastoch