Hello. My CROCO setup consists of a 300 m resolution grid with variable terrain, 72 vertical levels with theta_s=6, theta_b=6. The Haidvogel number is 0.2. My model setup is running fine in Hydrostatic Mode (undef NBQ), but crashing at the first time step in non-Hydrostatic Mode (define NBQ).

I think I am not using the correct cpp definitions for advection/diffusion while compiling in non-hydrostatic mode. Below, is a list of cpp definitions that I am using for the run. The only difference is that the model is running fine when NBQ is turned off, but crashes when NBQ is turned on.

I need to know which cpp definitions are the best suited for a non-hydrostatic (NBQ) run. Below, is the list of the cpp definitions that I am using. If anyone has used NBQ before, please let me know where I am doing wrong here. What are the correct cpp definitions to use with NBQ? Thanks.

                  /* Parallelization */

define OPENMP

                  /* Non-hydrostatic options */

define NBQ

ifdef NBQ

define W_HADV_TVD

define W_VADV_TVD

endif

                  /* Open Boundary Conditions */

define TIDES

define OBC_EAST

define OBC_WEST

define OBC_NORTH

define OBC_SOUTH

                  /* Grid configuration */

define CURVGRID

define SPHERICAL

define MASKING

define WET_DRY

define NEW_S_COORD

                  /* Model dynamics */

define SOLVE3D

define UV_COR

define UV_ADV

                  /* Equation of State */

define SALINITY

define NONLIN_EOS

                  /* Lateral Momentum Advection (default UP3) */

undef UV_HADV_UP3

undef UV_HADV_UP5

define UV_HADV_WENO5

undef UV_HADV_TVD

                  /* Lateral Explicit Momentum Mixing */

define UV_VIS2

ifdef UV_VIS2

define UV_VIS_SMAGO

endif

                  /* Vertical Momentum Advection  */

undef UV_VADV_SPLINES

define UV_VADV_WENO5

undef UV_VADV_TVD

                  /* Lateral Tracer Advection (default UP3) */

undef TS_HADV_UP3

undef TS_HADV_RSUP3

undef TS_HADV_UP5

define TS_HADV_WENO5

                  /* Lateral Explicit Tracer Mixing */

define TS_DIF2

undef TS_DIF4

undef TS_MIX_S

                  /* Vertical Tracer Advection  */

undef TS_VADV_SPLINES

define TS_VADV_AKIMA

undef TS_VADV_WENO5

                  /* Sponge layers for UV and TS */

define SPONGE

                  /* Vertical Mixing */

define LMD_MIXING

ifdef LMD_MIXING

define LMD_SKPP

define LMD_BKPP

define LMD_RIMIX

define LMD_CONVEC

undef LMD_DDMIX

define LMD_NONLOCAL

undef MLCONVEC

endif