Hello
I have implemented a CROCO-WW3 coupled simulation with BBL activated (as per the default coupling cpp keys). I noticed that the bottom stress seen by the model goes to zero in the shallow regions where one would instead expect the waves to increase the bottom stress.
On inspecting bbl.F, it appears that when waves are significant (Ub > Ubc) (bbl.F:345) the wave-modulated bottom stress taucw is only updated inside the #ifdef Z0_BL block (bbl.F:470). Since Z0_BL is defined only when SEDIMENT is defined, and I am running BBL without SEDIMENT, taucw is never updated and is left at its initial value of zero (bbl.F:177), which is what I see in my outputs.
Proposed fix: move bbl.F:471 (endif /* Z0_BL */) up to bbl.F:453, so that the bottom-stress is updated whether Z0_BL is defined or not.
Thanks for taking a look.
Cheers
Giles