I'm currently working on configuring ROMS to include the calculation of the surface boundary layer (Hsbl) and bottom boundary layer height (Hbbl) in the output files rom_his.nc and roms_avg.nc. However, despite enabling these options, I'm unable to see the results in the output files.
Here is how I've attempted to enable the output for Hsbl and Hbbl:
Code: Select all
Hout(idHsbl) == T ! Hsbl depth of surface boundary layer
Hout(idHbbl) == T ! Hbbl depth of bottom boundary layer
For context, here are the definitions included in my .h file:
Code: Select all
#define UV_ADV
#define UV_COR
#define UV_QDRAG
#define UV_VIS2
#define MIX_S_UV
#define DJ_GRADPS
#define SPLINES_VDIFF
#define SPLINES_VVISC
#define TS_DIF2
#define MIX_GEO_TS
#define TS_U3HADVECTION
#define TS_C4VADVECTION
#define SPHERICAL
#define MASKING
#define SALINITY
#define SOLVE3D
#define BODYFORCE
#define ANA_DIAG
#undef ANA_GRID
#undef ANA_INITIAL
#define ANA_SMFLUX
#define ANA_STFLUX
#define ANA_SSFLUX
#define ANA_BTFLUX
#define ANA_BSFLUX
/* Boundary conditions */
#define RADIATION_2D
/*
#define ANA_M2OBC
#define ANA_FSOBC */
/* Mixing definitions */
/*#define GLS_MIXING*/
#define MY25_MIXING
#if defined GLS_MIXING || defined MY25_MIXING
# undef KANTHA_CLAYSON
# define CANUTO_A
# undef N2S2_HORAVG
# undef RI_SPLINES
#else
# define ANA_VMIX
#endif
Best regards,
Hsin-Yi