open TELEMAC-MASCARET The mathematically superior suite of solvers

I have run a model for a variety of wave heights however I now want to alter the incoming wave direction. Currently the model runs for waves entering at 90 degrees to the boundary however I want to alter this so they are entering at 45 degrees to the boundary. I have tried altering the Direction of Wave Propagation value in the .cas file (shown below) however it made no changes to the output results. How do I vary the direction?

/

---

/ NUMERICAL PARAMETERS, SOLVER
/

---

SOLVER =8
BREAKING =YES
WAVE PERIOD =11.
DIRECTION OF WAVE PROPAGATION =135.

Thanks for your help

Niamh

Hello Niamh,

Are you working with regular or irregular waves?
Which version of the code do you use?
Did you pay attention to phase definition?
Can you send me your .par and your fortran file ?

regards,

Christophe

Hi Christophe

I'm using version 6.3 and working with regular waves.
My fortran file is attached. What is the .par file? I haven't come across this so far.

Thanks for your help.

Niamh

OK, it seems like your phase is not set correctly.

when your incident wave is perpendicular to the boundary, you can assume that your phase is the same on the boundary and fix it to 0.
But in your case, your incident wave is not perpendicular to the boundary, so you have to define a phase shift in your fortran file. Basically, the phase will be <k,x>, where x is the coordinates of the node and k the wave number.

you can use the example files, as btwi (btwi-par.f file) for instance. The phase shift isgiven in ALFAP variable.

If you still have a problem, don't hesitate to send your case (parameter file, fortran file, mesh) for me to have a look.

Christophe
(sorry : I call .par or .cas the parameter file ^^)

Thanks for your help Christophe.

Is the x node just a single point along my incident wave boundary? Is this phase value then converted to degrees?

I have attached my files if you wouldn't mind having a look to give me some guidance on what needs to be changed. I have looked at the example you mentioned but still not too clear on it

Thanks

Niamh

Hi Niamh,

x is the position of the current node, in the loop on the boundary points.
Thus, you have to write a loop looking like :
do I=1,Number of nodes on the incident wave boundary
Phase = <k,x>, meaning :
ALFAP(I)= K * cos (TETAB) * X(I) + K * sin (TETAB) * Y(I)
enddo
You can also choose a reference point Xo,Yo for your phase and write :
ALFAP(I)= K * cos (TETAB) * (X(I)-Xo) + K * sin (TETAB) * (Y(I)-Yo)

The phase ALPHAP should be given in degrees, so you need to write ALPHAP=ALPHAP*180/PI

I'll have a look at your case on monday if you can you give me your mesh file for me to test.

Hope this helps,

Christophe

Thanks very much Christophe, my mesh files are attached if you wouldn't mind having a look.

Niamh

OK, I'll have a look.
Just a question : your file ".cas" calls a mesh file named "Bathymetry_4.ser".
I couldn't find this file in your mesh files...

Christophe

This file is attached now. Sorry it didn't upload with the rest.

Thanks,

Niamh

Hi Niamh,

I add a phase definition in your fortran file using the method proposed in btwi-par.f file.

As your depth is not constant on the boundary, I had to choose a reference depth (9m here) to define a reference wave number. I also chose a reference point for the phase (your node 1309) but you can use any other values.

It looks good.

I attached the case (OK.zip file) and hope it will be OK.

Regards,

Christophe