[MITgcm-devel] vermix test experiment

[Martin Losch]

Hi Jean-Michel,

great, go ahead with eosType='MDJWF'. The reason, why the wrong k- 
level gives the same result is probably because the tref and sref are  
constant with k.

sigmaR: when I wrote these pkgs, sigmaR was not around. Maybe it's  
really worth replacing the buoyancy computation with this field. I  
won't have time for this in the next 2 weeks, so if you want to go  
ahead, feel free.

Martin

On 12 Aug 2008, at 14:34, Jean-Michel Campin wrote:

> Hi Martin,
>
> I was thinking of changing verification/vermix/input/data
> to switch from Linear EOS to eosType='MDJWF', (it's already
> there, but commented out):
> The reason is that when I was changing all those "call find_rho"
> in all those vertical mixing scheme pkg (by the way, most of
> those calls could be removed just by using sigmaR),
> I intentionally added a bug (wrong k reference level) to check
> how wrong it was, and with linear EOS, the results were identical
> (no change, with or without the bug).
> Since those tests (vermix.*) are just there to check
> all those pkg, it would be better to make it more selective
> (and indeed, with MDJWF, the bug shows up clearly in the output).
> Plus, I don't see any experiment which is testing this EOS,
> is it right ? this could be an additional good reason.
>
> Cheers,
> Jean-Michel
>
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel