[MITgcm-support] DST3 advection scheme

[Jean-Michel Campin]

Hi Samar,

The DST3 could be written sligtly differently with no division by
Rj, but for now, I suggest to un-comment those 2 lines and comment 
the 3rd one:

c       thetaP=0.
c       IF (Rj.NE.0.) thetaP=Rjm/Rj
        thetaP=Rjm/(1.D-20+Rj)

And the same with thetaM, 
and probably the same in adv_y.F & adv_z.F

See you,

Jean-Michel