Issue with HISTORY file logical switch for Hbbl

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
osean
Posts: 47
Joined: Sat Jan 23, 2021 3:46 pm
Location: MIT-WHOI

Issue with HISTORY file logical switch for Hbbl

#1 Unread post by 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
Posts: 1172
Joined: Wed Dec 31, 2003 6:16 pm
Location: USGS, USA

Re: Issue with HISTORY file logical switch for Hbbl

#2 Unread post by 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
Posts: 47
Joined: Sat Jan 23, 2021 3:46 pm
Location: MIT-WHOI

Re: Issue with HISTORY file logical switch for Hbbl

#3 Unread post by 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
Posts: 47
Joined: Sat Jan 23, 2021 3:46 pm
Location: MIT-WHOI

Re: Issue with HISTORY file logical switch for Hbbl

#4 Unread post by 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
Posts: 1172
Joined: Wed Dec 31, 2003 6:16 pm
Location: USGS, USA

Re: Issue with HISTORY file logical switch for Hbbl

#5 Unread post by 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
Posts: 5
Joined: Wed Feb 27, 2019 8:24 pm
Location: School of Marine Science, Sun Yat-sen University

Re: Issue with HISTORY file logical switch for Hbbl

#6 Unread post by huangzhy36 »

jcwarner wrote: Fri Jul 29, 2022 7:49 pm 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
Posts: 1172
Joined: Wed Dec 31, 2003 6:16 pm
Location: USGS, USA

Re: Issue with HISTORY file logical switch for Hbbl

#7 Unread post by 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

Post Reply