Issue with HISTORY file logical switch for Hbbl

[osean]

I am trying to make use of the logical switch for 'depth of bottom boundary layer' - i.e. Hbbl switch in the roms.in file.

I have set the following switches under the Logical switches (TRUE/FALSE) to activate writing of fields into HISTORY output file section to the following -

Hout(idDano) == F ! rho density anomaly
Hout(idVvis) == F ! AKv vertical viscosity
Hout(idTdif) == F ! AKt vertical T-diffusion
Hout(idSdif) == F ! AKs vertical Salinity diffusion
Hout(idHsbl) == T ! Hsbl depth of surface boundary layer
Hout(idHbbl) == T ! Hbbl depth of bottom boundary layer
Hout(idMtke) == T ! tke turbulent kinetic energy
Hout(idMtls) == T ! gls turbulent length scale

But after the run is complete, my roms_his.nc file doesn't have any of Hsbl, Hbbl, tke, or gls variables in it. Is there something else that I need to switch on / modify for the logical switch to work for them?

Similarly, I tried switching off salinity but setting the following line to

Hout(idTvar) == T F ! temp, salt temperature and salinity

But the output his.nc file still have salt. To me the logical switch doesn't seem to be working.

Any idea why so?

Additionally, is there somewhere I can find the definition of Hbbl that would be written into the output? Like how is Hbbl defined?

Cheers

[jcwarner]

gls and tke will only be output if you activate GLS_MIXING or maybe also MY25_MIXING.

for the hsbl, those go with the LMD options.

-j

[osean]

thanks - never done this before -

If I activate GLS_MIXING, do I need to activate all the 'related' ones as well that are described on this page? ==> https://www.myroms.org/wiki/GLS_MIXING

For LMD - do I simply activate it or is there anything I should be aware of? I can't find any instruction material on LMD. thanks

[osean]

Similarly, Ur and Vr (bottom velocity above bed) writing also don't appear to be activated

Qout(idUbur) == T ! Ur bottom U-velocity above bed
Qout(idVbvr) == T ! Vr bottom V-velocity above bed

Do they also need to activate some options? I can't find any info on this

[jcwarner]

for gls i typically use

# define GLS_MIXING
# ifdef GLS_MIXING
# define RI_SPLINES
# define KANTHA_CLAYSON
# define N2S2_HORAVG
# endif

it can get more options, for surface mixing etc, not sure what type of applications you are using.

for LMD i will let someone else answer.

[huangzhy36]

gls and tke will only be output if you activate GLS_MIXING or maybe also MY25_MIXING.

for the hsbl, those go with the LMD options.

-j

Is this means that Hsbl and Hbbl can only output when we use LMD mixing? It could not output when we use GLS or MY-2.5 mixing?

[jcwarner]

in ROMS/Utility/def_his.F I see

# ifdef LMD_SKPP
!
! Define depth of surface boundary layer.
!
IF (Hout(idHsbl,ng)) THEN
.........
status=def_var(ng, iNLM, HIS(ng)%ncid, HIS(ng)%Vid(idHsbl), &
& NF_FOUT, nvd3, t2dgrd, Aval, Vinfo, ncname)
IF (FoundError(exit_flag, NoError, __LINE__, MyFile)) RETURN
END IF
# endif
# ifdef LMD_BKPP
!
! Define depth of bottom boundary layer.
!
IF (Hout(idHbbl,ng)) THEN
...
status=def_var(ng, iNLM, HIS(ng)%ncid, HIS(ng)%Vid(idHbbl), &
& NF_FOUT, nvd3, t2dgrd, Aval, Vinfo, ncname)


so you need LMD_SKPP fo rthe Hsbl and LMD_BKPP for the Hbbl