[MITgcm-support] anelastic configuation for Jupiter

[Eli Galanti]

hello,

I'm running the deep anelastic case, but set to Jupiter parameters. The
goal is to revive the PhD work of Yohai Kaspi.

when running with a boussinesq approximation (like in the verification),
the model runs fine, but with the the full anelastic approximation (the
reference density changes dramatically with depth, and the prognostic
density is a function of temperature and pressure) the solution get out of
range within a few minutes.

The problems arises in the polar regions (the domain goes 80S to 80N),
where small scale velocities emerge, and grow exponentially.

I noticed a warning I get, one that did not appear in the code few years
ago:

(PID.TID 0001.0001) ** WARNING ** CONFIG_CHECK: Implicit viscosity applies
to provisional u,vVel
(PID.TID 0001.0001) ** WARNING ** => not consistent withfinal vertical
shear (after appling 3-D solver solution
(PID.TID 0001.0001) ** WARNING ** CONFIG_CHECK: Implicit viscosity not
implemented in CALC_GW
(PID.TID 0001.0001) ** WARNING ** CONFIG_CHECK: Explicit viscosity might
become unstable if too large

could that explains the problem, and if yes, what can I do?

thank you!

Eli Galanti
Weizmann Institute of Science
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20150316/a3aca999/attachment.htm>