TELEMAC-MASCARET Forum: Telemac3d+Gaia Wet/Dry beach (1/1)

Telemac3d+Gaia Wet/Dry beach 4 years 3 days ago #37200

ekirinus
OFFLINE
Fresh Boarder
Posts: 21
Thank you received: 1

Hi all,

I am now trying to couple Gaia with Telemac3d (v8p1r2) into my simulation from my previous question [1], and I'm running into problems early in the parametrization stage. The sediment is one class 0,160mm in the entire domain with bedload only. The .cas files are attached here, but the entire setup is 90 MB, so I uploaded it on gdrive: files

Question 1: I was trying to run the attached case with SUSPENSION FOR ALL SANDS=NO (the Gaia equivalent of Sisyphe's SUSPENSION keyword), but it seems that with only non-cohesive sediment and SUSP=.false., Gaia will silently ignore it and later on there will be a Fortran runtime error from within telemac3d.F. Set SUSPENSION FOR ALL SANDS=NO in the attached gaia_cc_trial.cas to see that error. So the question is why can't I run non-cohesive without suspension?

Question 2: Yielding and using SUSPENSION FOR ALL SANDS=YES, I get three messages and the program stops:

WARNING IN DICTIONARY:
FOR KEYWORD: VARIABLES FOR GRAPHIC PRINTOUTS
THE NUMBER OF DEFAULT VALUES            1  IS DIFFERENT FROM THE DECLARED SIZE            7

WARNING IN DICTIONARY:
FOR KEYWORD: D90 SAND DIAMETER FOR ONLY ONE CLASS
THE NUMBER OF DEFAULT VALUES            1  IS DIFFERENT FROM THE DECLARED SIZE            2

GIVE THE KEY-WORD
INITIAL VALUES OF TRACERS
FOR ALL TRACERS

So what is wrong with the VARIABLES FOR GRAPHIC PRINTOUTS? There are 10 values to the keyword, so why is it saying there is only one? Same for D90 SAND DIAMETER FOR ONLY ONE CLASS?

About the INITIAL VALUES OF TRACERS, there is no keyword with that name in gaia.dico, so I suppose it has to be given in the telemac3d cas file? And what is this value, the height of the sediment layer, the mass, or the concentration? I couldn't find it in the Gaia manual.

[1]: opentelemac.org/index.php/kunena/21-tele...topobathymetric-data

The administrator has disabled public write access.

Telemac3d+Gaia Wet/Dry beach 4 years 3 days ago #37201

pavans
OFFLINE
Moderator
Posts: 46
Thank you received: 17

Hello,

at the moment I cannot open files that you uploaded on google drives so I cannot answer to your first question even if I can say that is possible to run simulation with non cohesive sediment and only bedload.

For your second question, the first two messages are only warning messages but the third message says that you have to set the keyword INITIAL VALUES OF TRACERS (in the TELEMAC-3D steering file). Indeed sediments are seen as tracers by TELEMAC-3D (in case of suspension) so you have to assign their initial values of concentration. For the next version it will be possible set them in the Gaia steering file using the keyword INITIAL SUSPENDED SEDIMENTS CONCENTRATION VALUES.

Hope it helps,

Sara

The administrator has disabled public write access.

Telemac3d+Gaia Wet/Dry beach 4 years 3 days ago #37202

ekirinus OFFLINE Fresh Boarder Posts: 21 Thank you received: 1	I update the link
	The administrator has disabled public write access.

Telemac3d+Gaia Wet/Dry beach 3 years 11 months ago #37212

pavans OFFLINE Moderator Posts: 46 Thank you received: 17	Hello, I've checked your steering files: you have 2 tracers (salinity and temperature) and 1 NCO sediment. This configuration has not been checked before and you've found a bug, this is why runtime error appears... I will come back asap with a fix. Best regards, Sara
	The administrator has disabled public write access.

Telemac3d+Gaia Wet/Dry beach 3 years 11 months ago #37213

pavans OFFLINE Moderator Posts: 46 Thank you received: 17	* I'm speaking about the case with only bedload
	The administrator has disabled public write access.

Telemac3d+Gaia Wet/Dry beach 3 years 10 months ago #37559

pavans
OFFLINE
Moderator
Posts: 46
Thank you received: 17

Hello again Eduardo,

I'm sorry for the late reply.
For your issue (question 1) you can try this quick fix : adding after line 1804 (else statement) of telemac3d.F these 2 lines :
ISUSP=1
NUM_ISUSP_ICLA(1)=1
(Then you should put the subroutine telemac3d.F in user_fortran folder, with the other user fortrans modified).
I hope it helps, this will be for sure fixed in the next release also, thanks for your reporting.

Best regards,
Sara

The administrator has disabled public write access.