[MITgcm-devel] NONLIN_FRSURF and SIDEDRAG

[Baylor Fox-Kemper]

Does anybody understand what this is about?

Also, why is this set to be nonnegative only when sidedragfactor.le.0 ?

#ifdef NONLIN_FRSURF
C-      this will not give any side-drag along thin wall.
C       (but this might just be what we want ...)
         hFacZClosedS = h0FacW(i,j,k,bi,bj) - h0FacW(i,j-1,k,bi,bj)
         hFacZClosedN = h0FacW(i,j,k,bi,bj) - h0FacW(i,j+1,k,bi,bj)
         hFacZClosedS = MAX( 0. _d 0, hFacZClosedS )
         hFacZClosedN = MAX( 0. _d 0, hFacZClosedN )
#else
         hFacZClosedW = _hFacS(i,j,k,bi,bj) - hFacZ(i,j)
         hFacZClosedE = _hFacS(i,j,k,bi,bj) - hFacZ(i+1,j)
#endif

   Cheers,
     -Baylor