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 */