open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi, Elodie,

Thanks for your advice and my FOUD.t3s represents the bathymetry data in geo.slf. I don't know which name (BOTTOM or FOUD) is best to use in TOMAWAC?

I try to reduce the minimal frequency=0.03 or vary time step, but wave height HM0 is still much low than Mike module.

SUBROUTINE ANAVEN in my Fortyran programm is modified as blow;
VITES=VX_CTE
DIR_1=270.D0
C
CNVR=3.141592654D0/180.D0
C
UCONST=VITES*SIN(DIR_1*CNVR)
VCONST=VITES*COS(DIR_1*CNVR)
!
DO IP=1,NPOIN2
!UV(IP)=VX_CTE
!VV(IP)=VY_CTE
UV(IP)=UCONST
VV(IP)=VCONST
ENDDO

Is right or not?

Thanks and regards,

Zheng

Hi, Elodie,

Could you kind to take a look at my model and find out why I get such results, please? of course, when you have a time.

The geometry file, boundary condition file and steering file all been sent before.

Thanks,

Zheng

Hi

I think I have an idea about the reason of the problem!
In your geometry file, the bottom describe in the variable FOND doesn't have any attribute.
Could you check in the result file if the bottom is well take in account by adding B to the list of graphic output variables...

Hope this helps

Hi, Christophe,

Thank for your help

I changed FOND to the bottom in the geometry file and add B to the list of VARIABLES FOR 2D GRAPHIC PRINTOUTS as below:

VARIABLES FOR 2D GRAPHIC PRINTOUTS ='HM0,TM01,TM02,TMOY,FX,POW,DMOY,B'

but there is no value with B. what is problem?

VARIABLES FOR 2D GRAPHIC PRINTOUTS ='HM0,TM01,TM02,TMOY,FX,POW,DMOY,ZF'

but Bottom values all are -70 with ZF. is it right or not?


The new geometry file attached.

Thanks,

Zheng

Hi Zheng

Sorry for the mistake, the avlue is ZF not B for the graphic printout

I had a look into your case without success for the moment but I had some questions and remarks

For me the fortran file is not necessary, you give a constant value for the wind and it's sufficient.
For your steering file, why different type for initial and boundary values?
type 1 and 6
fetch 0 and 10000
...
It's strange for me but i'm not a specialist of Tomawac

Regards

Hi, Christophe,

Thank you very much.

1) I don’t know whether my Fortran file is need or not. In my Fortran file (subroutine Avaven) I only add :
VITES=VX_CTE (VX_CTE=40m/s, wind speed)
DIR_1=270.D0 (wind direction)
CNVR=3.141592654D0/180.D0
UCONST=VITES*SIN(DIR_1*CNVR)
VCONST=VITES*COS(DIR_1*CNVR)
DO IP=1,NPOIN2
!UV(IP)=VX_CTE
!VV(IP)=VY_CTE
UV(IP)=UCONST
VV(IP)=VCONST
ENDDO
But in my steering file,
WIND VELOCITY ALONG X =40.0
WIND VELOCITY ALONG Y = 0.0
CONSIDERATION OF A WIND = true ( may be wrong?)
STATIONARY WIND = true ( may be wrong?)

But there is not wind direction defination?

2) TYPE OF BOUNDARY DIRECTIONAL SPECTRUM = 6
TYPE OF INITIAL DIRECTIONAL SPECTRUM =1

Which are based on 8.5.2. presscribing the boundary condition in the CAS file and 8.4. Presscribing the initial conditions in guide_tomawac_v6p1.pdf. I rally don’t know how do set up the initial and boundary value in steering file? This is my first case to test Tomawac and do the comrasion on Mike21 SW.

3) BOUNDARY MEAN FETCH VALUE=0.
INITIAL MEAN FETCH VALUE =10000.
(Mike defination initial max.FETCH length=100000.)


May be some wrong in my steering file and could you help me to run the case (about 1 hour CPU used) and help me to find the problem, please?

Thank you very much
Zheng

Hi,

I agree with Christophe, your Fortran file is not needed : the wind direction is given when you impose the two keywords WIND VELOCITY ALONG X and WIND VELOCITY ALONG Y.

CONSIDERATION OF A WIND = true and STATIONARY WIND = true are OK.

I will try to have a look at your case file or run your case and I will tell you if I see something wrong.

Regards,

Elodie

Thanks to Elodie.

Hello again,

I ran your case and here are the different problems I saw :

- First, in your boundary condition, you put BOUNDARY PEAK FREQUENCY=9.334 : I think 9.334 is the peak period (in seconds) and your peak frequency then has to be 0.107 Hz.

- Secondly, I don't understand very well why your boundary condition is an imposed spectrum on the west boundary, as your wind is blowing to the west : this boundary condition is not very useful, and on the contrary some condition would be needed on the east boundary. Or maybe you made a mistake about the direction convention in Tomawac and your wind is actually blowing to the east?

- As said before, the fortran file is not needed, you only have to put the keywords :
CONSIDERATION OF A WIND = true
STATIONARY WIND = true
WIND VELOCITY ALONG X=40.0 (blowing to the east) or -40 (blowing to the west)
WIND VELOCITY ALONG Y=0.0

I hope this helps.

Best regards,

Elodie

Hi, Elodie,

Thank you very much!

It is my mistake for the BOUNDARY PEAK FREQUENCY sholud be 0.107Hz.

Secondly,the west boundary is an imposed spectrum with significant wave height=6.5, spectral peak period=9.334s, mean wave direction=270 degree. And directional standard deviation=0. and Wind forcing speed=40 and direction=270 degree. The Wind is blowing to the east direction?

I only modified BOUNDARY PEAK FREQUENCY =0.107 and fortran code is not used and then run the case again, but it is not right. Some keyword are list as below:

BOUNDARY DIRECTIONAL SPREAD 1 = 25.
BOUNDARY MAIN DIRECTION 1 =40.

WIND VELOCITY ALONG X =40.0
WIND VELOCITY ALONG Y =0.0

INITIAL DIRECTIONAL SPREAD 1 =10.
INITIAL MAIN DIRECTION 1 =270.

could you help me to look at them, please? and my result are attached.

Thanks,

Zheng