Hello Sanaz,
What are the units of your y-axes in your figures, meters? For free surface, with respect to what the bottom elevation is related, Mean Sea Level?
Using OPTION FOR LIQUID BOUNDARIES = 2 (Thompson boundary conditions) does not prescribe exactly the water levels at the boundary conditions you give as input. If you want to prescribe exactly water levels at your boundary conditions, you should use the classical boundary conditions (OPTION FOR LIQUID BOUNDARIES = 1) and codes 5 4 4 in the BOUNDARY CONDITIONS FILE (LIHBOR = 5 prescribed water depth, LIUBOR=LIVBOR = 4, free velocity).
Anyway, can you show a figure of the extent of your domain with scale (e.g. longitude/latitude or a projection with scale)?
If using Thompson boundary, it is recommended to have boundaries almost perpendicular to the velocity currents.
The NUMBER OF BOTTOM SMOOTHINGS may change the water levels in your domain.
Moreover, you should clean your TELEMAC-3D steering file:
- some options are inconsistent (e.g. NON-HYDROSTATIC VERSION = NO and DYNAMIC PRESSURE IN WAVE EQUATION = YES),
- if modelling temperature and salinity and you want to see stratification or buoyancy effects, it is strongly recommended to use the NON-HYDROSTATIC VERSION (which is the default value),
- some default values for some keywords should be used first before tuning your TELEMAC-3D case. I recommend you to have a look at what my colleague Agnes has already written :
www.opentelemac.org/index.php/kunena/21-...se-with-some-figures
(IMPLICITATION FOR DEPTH = 0.55, IMPLICITATION FOR VELOCITIES = 0.55, SCHEME FOR ADVECTION OF... = 5 + SCHEME OPTION FOR ADVECTION OF... = 4 but also MINOR CONSTITUENTS INFERENCE = NO, NON-HYDROSTATIC VERSION = YES). When you change the default values, you should wonder the signification of the keyword and why you want to change it.
If you want ot compare your T2D and T3D model only for water levels, I would advise you to only compare tidal and hydrodynamics options.
As already written in a previous post, tpxo 9 includes 2 extra tidal constituents (2N2+S1) compared to previous versions and currently, TELEMAC does not take into account the S1 constituent as written in a warning message that you showed.
Did you have a careful look at the output file in particular before the first time step and at the final step just to be sure there is no warning from TELEMAC?
Hope this helps,
Chi-Tuan