Hi there,
I have set up a 1km scale model and am having difficulty with overheating in the surface layer and throughout the water column. This heating up occurs within the first day and remains. I am using bulk fluxes. I have tried with both hourly ERA5 (from cds store using ERA5 scripts) and and the default bulk fluxes built into crocotools, but I get the same effect with both. I have attached my cpp file for inspection cppdefs.h;
I am currently working through a few possibilities related to vertical mixing parametrisation and OBC options, but nonetheless I would appreciate any suggestions based on your experience or knowledge.
Thanks in advance,
Joe
UPDATE: I am still asking for your suggestions and experience, but it seems (provisionally) that reverting to v1.1 solves this issue. I had switched everything to v1.2 (then v1.2.1) and it seems that there is some major difference, which either requires another CPP flag that isn’t contained in the default cppdefs.h, or there is something buried in the code which causes issues in my model domain, but not for other models used for benchmarking.
Here is the cppdefs.h file for v1.1 cppdefs.h