Nested 1/150° configuration crashes with neg dp
Alfredo Terrazas
Hello everyone,
I am working on a very high-resolution HYCOM configuration in the southern Gulf of Mexico with a horizontal resolution of 1/150°. This is the first time I have worked at such a high resolution, and I am encountering some difficulties with nesting.
My base experiment (expt_01.0), which runs without nesting, completes successfully (1-year run). However, when I run the nested configuration using boundary conditions from a parent simulation at 1/50° resolution, the model crashes after a few days with a "neg dp" error.
I have already experimented with reducing the time step, but the problem persists.
The model was compiled with:
setenv OCN_MISC "-DMASSLESS_1MM -DRDNEST_MASK -DLATBDT_NPLINE3"
I would like to ask specifically about the nesting mask. I am attaching the mask used for the nested domain and would appreciate any feedback on whether it looks reasonable or if there are obvious issues that could lead to negative layer thicknesses near the boundaries.
Has anyone experienced similar problems when nesting from 1/50° into 1/150°? Are there any particular checks you would recommend for the nesting mask or boundary treatment at this resolution?
I am also attaching the blkdat.input from the working non-nested experiment (expt_01.0) for reference and also from the nested experiment (expt_01.2).
Thank you for any suggestions.
