[MITgcm-devel] twice-yearly TAMC reminder
Patrick Heimbach
heimbach at MIT.EDU
Fri Dec 10 15:22:57 EST 2004
Hi there,
here's an annual Christmas-reminder,
that comments interlaced within ongoing lines
break TAMC compatibility.
Examples of what I mean are in
pkg/mom_vecinv/mom_vi_u_coriolis_c4.F
In the example below, the first #ifdef block is BAD,
the second one is GOOD.
#undef THIS_IS_BAD_FOR_TAMC
#ifdef THIS_IS_BAD_FOR_TAMC
uCoriolisTerm(i,j)=
C high order vorticity advection term
& +vort3u*vBarXY*recip_dxc(i,j,bi,bj)
C linear Coriolis term
c & +0.5*(fCoriG(I,J,bi,bj)+fCoriG(I,J+1,bi,bj))*vBarXY
C full nonlinear Coriolis term
c & +0.5*(omega3(I,J)+omega3(I,J+1))*vBarXY
C correct energy conserving form of Coriolis term
c & +0.5*( fCori(I ,J,bi,bj)*vBarY(I ,J,K,bi,bj) +
c & fCori(I-1,J,bi,bj)*vBarY(I-1,J,K,bi,bj) )
C original form of Coriolis term (copied from calc_mom_rhs)
c & +0.5*(fCori(i,j,bi,bj)+fCori(i-1,j,bi,bj))*vBarXY
& *_maskW(I,J,K,bi,bj)
#else
uCoriolisTerm(i,j)=
& +vort3u*vBarXY*recip_dxc(i,j,bi,bj)
& *_maskW(I,J,K,bi,bj)
#endif
Cheers
-Patrick
--
----------------------------------------------------------
Patrick Heimbach Massachusetts Institute of Technology
FON: +1/617/253-5259 EAPS, Room 54-1518
FAX: +1/617/253-4464 77 Massachusetts Avenue
mailto:heimbach at mit.edu Cambridge MA 02139
http://www.mit.edu/~heimbach/ USA
More information about the MITgcm-devel
mailing list