[MITgcm-support] Uniform diffkrfile
Geoffroy Aubry
gaubry at rip.ens-cachan.fr
Tue Aug 7 11:58:07 EDT 2007
Dimitris Menemenlis a écrit :
> Geoffroy, when you specify a diffkrfile with uniform 1.E-5 m2.s-1
> diffusivity, the situation would (should!) be equivalent to your case 3
> below, that is, diffKzT=diffKzS=1.E-5. Do your case 1 and case 3 below
> give equivalent results?
>
Yes you're right, the outputs are exactly the same (I did a diff -q for
all the files in the input directory, and the only differences are the
data files).
But I have a uniform salinity map (I don't give any hydrogSaltFile in
the data file). So if I put any diffusivity coefficient for salt, I
should have the same results isn't it? But it is not the case... I did
an experiment with diffKzS at 1.E-5 and 5.E-5 and I dont get the same
results (all the Eta.n, n>2, PH.n, n>3, T.n, n>48, U.n, n>2, W.n, n>2
differ).
For the moment, I try to understand why with a uniform vertical
diffusivity, with a temperature profile of θ = −1.10^−3 z + 5.5 °C (z>0)
and with a topography everywhere -1500m but at one place a wall (see
enclosed picture), without doing any forcing, I get everything mixed.
If the initial situation is stable (which I think it is), it should stay
the same over the time.
But after 50000s, I get some huge velocities (W~4m.s^-1).
If I remove the wall, then I get a stable situation again (test2 here :
http://perso.crans.org/gaubry/dossier/stage/index.html#bug ) and if I
put the vertical diffusivity for salt to 0, the situation is also stable
( http://perso.crans.org/gaubry/dossier/stage/bug.html )
Do you have any clues?
Thank you for your support,
--
Geoffroy
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image_50000.png
Type: image/png
Size: 32563 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20070807/e9b6027e/attachment.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image_00000.png
Type: image/png
Size: 13441 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20070807/e9b6027e/attachment-0001.png>
More information about the MITgcm-support
mailing list