CHLNDDEV / OceanMesh2D

A two-dimensional triangular mesh generator with pre- and post-processing utilities written in pure MATLAB (no toolboxes required) designed specifically to build models that solve shallow-water equations or wave equations in a coastal environment (ADCIRC, FVCOM, WaveWatch3, SWAN, SCHISM, Telemac, etc.).
https://github.com/sponsors/krober10nd
GNU General Public License v3.0
182 stars 65 forks source link

How to set up parameters of hybrid bottom friction in OM #313

Closed Jiangchao3 closed 4 months ago

Jiangchao3 commented 5 months ago

Hi @krober10nd and @WPringle,

I am wondering how to set up parameters for hybrid bottom friction in OM. The following is the description in the manual:

If NOLIBF = 2, the hybrid bottom friction formulation is used FFACTOR = CF*[1+(HBREAK/H)FTHETA](FGAMMA/FTHETA and CF is as specified in this expression. Note, that FFACTOR approaches CF in deep water (H > HBREAK) and the hybrid friction formulation reverts to a standard quadratic formulation. This option is not available for friction coefficients specified using nodal attributes.

So four parameters is required, CF, HBREAK, FTHETA and FGAMMA.

If I want to use this function, should I set one line CF = CF, HBREAK, FTHETA and FGAMMA to assign this four parameters when using Make_f15.m, or in other ways ?

Looking forward to your reply. Thanks very much.

Jiangchao

Jiangchao3 commented 4 months ago

Just got it, it is easy to accomplish. Sorry for disturbing.

krober10nd commented 4 months ago

Maybe if you could post how you did it briefly here for others that would be helpful. Thanks!

Jiangchao3 commented 4 months ago

Yes,

When useing the hybrid bottom fraction,

m.f15.bolibf = 2; m.f15.taucf = [0.0025 2 10 1/3]; % [CF HBREAK FTHETA FGAMMA]