[MITgcm-devel] DST3, Leith and adjoint

[Patrick Heimbach]

Hi Jean-Michel and Baylor,

I made a few small changes to
o mom_vi_hdissip.F
o gad_dst3_adv_?.F
You might want to check whether I screwed up anything.

As a general reminder,

d/dx SQRT(x) = 1/(2*SQRT(x))

which is why we always exclude x=0 from our
definition range when using SQRT().

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