[MITgcm-support] question...
Martin Losch
mlosch at awi-bremerhaven.de
Fri Sep 29 02:48:09 EDT 2006
Dear Rima,
this check doesn't have anything to do with Tref or Sref.
S/R (=subroutine) config_check.F tests the consistency of a few
things. In this case the message is quite clear, isn't it:
"nonhydrostatic NOT SAFE with non fully implicit barotopic solver".
As far as I can see, this message is related to these lines in data:
implicSurfPress=0.5,
implicDiv2DFlow=0.5,
which make the solver semi-implicit (I am guessing, I have never used
these flags myself), they make the model use a Crank-Nickelson time
stepping scheme (according to model/inc/PARAMS.h, where most
parameters are briefly described).
You have two options:
1. comment out these lines to go back to full implicit solver (which
i would do, maybe someone can comment on WHY the solver is not fully
implicit in this experiment)
2. Ignore the warning and comment out the corresponding STOP
statement in model/src/check_config.F (as suggested in the second
line of the warning, but it's probably not safe.)
Martin, Germany
PS. there is also a flag "implicitIntGravWave" (default=false), which
you may want to try in conjunction with internal gravity wave .
On Sep 29, 2006, at 6:09 AM, Riema Rachmayani wrote:
> dear all,
>
> i've already run verification of internal wave, i'm not change the
> input...i just run it!!....when i check the input, Tref is various
> to the vertical layer (20 layer), and Sref is uniform to the
> vertical layer (20 layer), nonhydrostatic=.false....(it mean the
> equation still hydrostatic right?)...and then i change
> nonhydrostatic=.true. , when i run it again...there was a message :
>
> config_check : nonhydrostatic NOT SAFE with non fully implicit
> barotopic solver
> config_check : to by-pass this STOP, comment this test and re-
> compile config_check
> STOP ABNORMAL END : S/R Confid_check statement executed
>
> i thought that it happen because Sref uniform, so i change Sref
> became various like Tref, but when i run it again, there was
> appears the same message...
>
> what should i do??
>
> thx mitgcm
>
> best regard,
> rima, indonesia
More information about the MITgcm-support
mailing list