[MITgcm-devel] nonhydrostatic code

[Alistair Adcroft]

Martin Losch wrote:
> Hi,
> 
> probably a stupid question: As far as I can see, the 
> momentum_correction_step calls integr_continuity, which calls 
> integrate_for_w even wenn I use ALLOW_NONHYDROSTATIC and 
> nonhydrostatic=.TRUE.. That is, the vertical velocities, that have been 
> updated in calc_gw are overwritten with the integral of the horizontal 
> divergence? That can't be true. What am I missing?

If the solver had iterated to round-off tolerence then the two different 
estimates of W would be exactly the same. When they are only 
approximately the same, we find it more consistent to use the W that is 
exactly non-divergent.

A.
-- 
Dr Alistair Adcroft            http://www.mit.edu/~adcroft
MIT Climate Modeling Initiative        tel: (617) 253-5938
EAPS 54-1624,  77 Massachusetts Ave,  Cambridge,  MA,  USA