Hello everyone,
I’m runnung a simulation just to test how the tides module is working with TPX07. For this purpose I have a real pressure dataset to compare with the model ‘zeta’ output. The problem is that I don’t know how to match the phase of my simulation outputs with my “real” data. In other words, I want the model outputs to coincide in day and hour with my SSH data (highs and low tides).
I set all the parameters in crocotools_param, thats the Ymin,Mmin,Dmin,Hmin.
It’s not clear for me what I have to put on Yorig …
Thanks,
Nicolas
Hi Nicolas,
I think Yorig is not relevant. It is just to avoid having a very large number in your time vector variable (seconds since Yorig)
How bad is you tidal phase problem? Does shifting the time series by a constant value corrects it? Amplitude is ok?
You could also shift between linear, quadratic, and Roughness length formulation for the bottom stress in the croco.in file
bottom_drag:
RDRG [m/s]: Drag coefficient for linear bottom stress formulation
RDRG2: Drag coefficient for constant quadratic bottom stress formulation
Zob [m]: Roughness length for Von-Karman quadratic bottom stress formulation
Cdb_min: Minimum value of drag coefficient for Von-Karman quadratic bottom stress formulation
Cdb_max: Maximum value of drag coefficient for Von-Karman quadratic bottom stress formulation.
or activate the Bottom Boundary Layer module (BBL)
See more info here
https://croco-ocean.gitlabpages.inria.fr/croco_doc/model/model.parametrizations.friction.html
https://croco-ocean.gitlabpages.inria.fr/croco_doc/model/model.modules.bbl.html
I agree with Andres.
Note that the model time is in UTC. Your data is probably given in local time. You need to do the conversion.
If there is still a phase shift after the conversion, then check the sensitivity to bottom friction and also compare with TPXO tides, which are forcing CROCO at the boundaries: do you also have a phase shift with TPXO? …
Hello Andres, yes I also have problem with the amplitud. If I shift the time series a couple of hours the phase fits.
Also, I think the Yorig is being used in the phase corrections done by crocotools (in make_tides.m):
% Add a phase correction to be consistent with the ‘Yorig’ time
%
t0=t0-24.*mjd(Yorig,1,1);
then :
correc_phase=-phase_mkB(it)-pu(it)+360.*t0./periods(it);
The amplitude is like 1m higher in the model, so I m going to take a look to those bottom stress coefficients.
Thanks!