open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi all,

I'm trying to implement a hydrodynamic model coupled with a wave model for a lake system, i.e., Chittaway Bay in southern lake Tuggerah in NSW, Australia. From the telemac-tomawac model manuals descriptions, it appears that I can use this tomawac to model the waves in this lake (no harbour like topography and no reflecting structures).

The lake has a lot of vegetation along its shoreline (and not much vegetation deeper in the lake). From looking at QVEG in the tomawac V7p0, the vegetation seems to be applied as a domain uniform feature, i.e., a constant dissipation value, R0.

If I want to do a 2D grid of vegetation, can I modify R0=constant to R0=some_analytic_function(x,y,z). This would allow me to specify 0.D0 at deeper ends and 1.D0 at shallower ends.

The change seems fairly straightforward and simple, so I'm wondering if I'm breaking the model in some way by doing this (that is not obvious).

Cheers,
Shiv

hello

Of course it is possible, you will find an example of such a dependance in the subroutine qveg of the test case dean.

that will not break anything

Hi everybody,
In version v8p4r0 of Tomawac. The example case "veget" don't use the QVEG subroutine for calculate wave dissipation on a wave flume. I want to re-use this example for a coastal purpose, but I just can't figure out if it's possible to define a vegetation zone within my domain. I have the impression that it's applied to the whole estate.

Can you help me with this?

Thank you in advance,
Mathieu S.

Hi Mathieu,
Attached is what I did. You can set polygon vertices in qveg.f file and call that fortran file in the tomawac cas file.
Cheers,
Shiv

Thank you so much, have a good day
Mathieu

Do you think writing vegetation parameters in the tomawac.cas file as in the example "Veget" overrides the parameters specified in the qveg subroutines?

/

---

/ VEGETATION
/

---

VEGETATION TAKEN INTO ACCOUNT = YES
VEGETATION DISSIPATION COEFFICIENT = 1.
NUMBER OF PLANTS M2= 20.
STEM DIAMETER=0.25
BULK DRAG COEFFICIENT=0.2
VEGETATION HEIGHT=2.

thanks

to be honest, I'm not sure. I hard coded the fortran file and it was simple enough since we were not changing those values.

Anything hard coded would likely not change but you can try to set CAS file to veg diss coeff=0 & hard code the same variable in the fortran file to 1 and see if it listens to CAS or fortran.

OR not fill in the NV, BV,CD, DV values in the fortran file and only specific them in CAS and see if it works.

I've tried it, so qveg.f takes over. My lines in tomawac.cas are useless

yeah, hard coded lines always override.

I assume you still run your simulations by entering the parameters into qveg.f?

yes, that's how I do it

Thank you for your help