[MITgcm-support] non-hydrostatic mode

[nkimura at cae.wisc.edu]

    Dear MITgcm users
    Hello.
    I have a fundamental question about non-hydrostatic mode in the ocean code.
    According to the manual (page 209), in order to set non-hydrostatic mode,
“nonHydrostatic=.TRUE.” in input/data and “#define ALLOW_NONHYDROSTATIC” in
code/CPP_OPTIONS.h should be done. I simply followed it, and tested wave
propagation.
    However, my result is not similar to the analytical solution and looks like
the result of quasi-hydrostatic mode (Causlli 1999), so I am wondering if or
not my setups of non-hydrostatic mode are correct.
    I would appreciate any suggestion.

    Nobuaki Kimura, Graduate student,
    CEE, UW-Madison

In my computational conditions of the oscillating basin (Causlli 1999,
International Journal for Numerical Methods in Fluid 30: 425-440)
- square basin with 10[m] in x and 10[m] in z
- initial free-surface of constant slope = 0.02 x – 0.1 [m]
- no viscosity, no forcing without initial free-surface, homogenous density
- no Coriolis force
 (Note)  My data, CPP_OPTIONS.h are attached.


-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: CPP_OPTIONS.h
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20060202/5ef102de/attachment.h>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: data
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20060202/5ef102de/attachment.el>