[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