[MITgcm-support] OBCS bug

Gianmaria Sannino gianmaria.sannino at casaccia.enea.it
Thu Jul 13 17:09:04 EDT 2006


Dear all,

there is a bug in obcs_calc.F, in particular when the 
nonlinear free surface option is defined on, the OBSeta 
and OBNeta are defined:

#ifdef NONLIN_FRSURF
               OBSeta(J,bi,bj)=0.
#endif

#ifdef NONLIN_FRSURF
               OBNeta(J,bi,bj)=0.
#endif


instead of:

#ifdef NONLIN_FRSURF
               OBSeta(I,bi,bj)=0.
#endif

#ifdef NONLIN_FRSURF
               OBNeta(I,bi,bj)=0.
#endif


ciao
gianmaria



More information about the MITgcm-support mailing list