Hello,
I’m working on an ESTUARY configuration with MRL_WCI; MRL_CEW and WKB_WWAVE.
Forcing is a sinusoidal tide (amplitude 2m, period 12h), river input with 400m3/s and west boundary wave (amplitude 0.5m, period 8s).
When MRL_CEW is activated, I observe a north/south asymmetry in hrm values, with higher values in the northern part of the estuary.
Without MRL_CEW, the hrm field is symmetrical (as expected).
I thought that enabling the WKB_KZ_FILTER and WKB_TIME_FILTER keys would solve the problem (from description in 1.3. Wave-averaged Equations — Croco Doc 1.3 documentation), but my simulation explodes at the first time step when I enable them.
I searched the code and found this line, which seems strange in my case with tidal forcing as Dstp is the water depth : OCEAN/wkb_wwave.F · master · croco-ocean / croco · GitLab
Dstp(i,j)=max(h(i,j)+zwave(i,j,wmod),eps)
If I change this line 531 to :
Dstp(i,j)=max(h(i,j)+zeta(i,j,knew),eps)
my simulation is stable and symmetrical with keys WKB_KZ_FILTER and WKB_TIME_FILTER.
Could it be a bug on this particular line?
