assistance to output Hsbl and Hbbl

[hyc006]

Hello everyone,

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

Could anyone advise if there are specific C-preprocessing options that need to be activated to successfully output Hsbl and Hbbl?

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

Any insights or suggestions would be greatly appreciated. Thank you in advance for your help!

Best regards,
Hsin-Yi

[wilkin]

I believe Hsbl and Hbbl are available for output only when using LMD_MIXING. You have MY25_MIXING.

Running with both BODYFORCE and a vertical turbulent mixing closure scheme might be inconsistent. Using MY25_MIXING (or GLS_MIXING or LMD_MIXING) assumes the stress is applied as the surface boundary condition to vertical momentum mixing.

[hyc006]

Running with both BODYFORCE and a vertical turbulent mixing closure scheme might be inconsistent. Using MY25_MIXING (or GLS_MIXING or LMD_MIXING) assumes the stress is applied as the surface boundary condition to vertical momentum mixing.

Thanks so much for your reply! For your comments on BODYFORCE and a vertical turbulent mixing closure scheme might be inconsistent, which kind of the mixing profile would you recommend me to used when running with BODYFORCE? Or this could be fixed if I'm not controlling the bodyforce by sustr and svstr but with bottom stress bustr and bvstr?

Thanks,
Hsin-Yi