Ocean Modeling Discussion

Hi. This might not be a satisfying answer but I would solve this problem in a simple way. I would not modify the code. That is taking the risk of including bugs and it makes it difficult to move to newer code version in the future. I would instead make a bash or python script that runs the model for...

The problem is solved! It was indeed NDTFAST, I sought it was the barotropic time step in seconds so I was decreasing it to remove the instability which made things worse.
Thank you for your help!

Using GLS vertical mixing parameterization with standard k-epsilon implementation and N2S2_HORAVG does not remove the instability.

Dear jcwarner,

For the moment I am using LMD_MIXING so from what I understand N2S2_HORAVG is not used in this case.
But I am going to try GLS_MIXING with N2S2_HORAVG maybe it solves my issue.

Thanks for the tip,
Dewi

Dear all, I should be doing something wrong but I can't find out so any help is welcome. I develop a 3D ROMS (version 3.7) configuration of the european coast forced with ERA-interim fields and bulk formula with horizontal resolution of 20km and 30 vertical layers. I get what looks like a checkerboa...