open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,
I'm very new to TELEMAC-3D and am trying to model with drogue/particle outputs.
I've copied the TELEMAC 3D particles example and edited the user_flot3d file to include UTM coordinates and depth (below). The job runs but my drogues output doesn't have any data. Each drogue section of the TELEMAC.out and sortie files show this error:
ERROR ON THE MASS DURING THIS TIME STEP

I'm guessing the issue may be due to either incoherence of time steps with the cas file? or do I need to change the numbers to do with LT?

Hope you can help, let me know if you need more :)

IF(LT.LE.600.AND.(10*(LT/10).EQ.LT.OR.LT.EQ.1)) THEN
CALL ADD_PARTICLE(346926.D0,5893245.D0+LT/3.D0,2.5D0+LT/100.D0,
& LT,1,NFLOT,
& NFLOT_MAX,XFLOT,YFLOT,ZFLOT,TAGFLO,CLSFLO,
& SHPFLO,SHZFLO,ELTFLO,ETAFLO,MESH3D,NPLAN,
& 0.D0,0.D0,0.D0,0.D0,0,0)
ENDIF

Hello Emma,

If you set MASS-BALANCE = YES, it is normal to have the message "ERROR ON THE MASS DURING THIS TIME STEP", this is information written periodically during the computation.

Are you sure the coordinates of particles you have given are inside the 3D volume of the domain? In particular for elevation (greater than 2.5 m)?

Which release do you use? Can you upload your steering file?

Hope this helps,

Chi-Tuan

Hi Chi-Tuan,
Thanks for your reply.
Ah, the coordinates should be within the domain but I am unsure of the depth at the moment!
Steering file attached .
Many Thanks,
Emma

Hello Emma,

There is no steering file attached. Be careful with the extension of the name of the file (.txt is OK but all do not work).

Chi-Tuan

Apologies, should be attached now.
Many Thanks,
Emma

Hello Emma,

Is your computation with particles working now?

Looking at your steering file:
- I wonder why you only use 2 planes (= 1 for the bottom and the other one for the free surface)? In that case, you may use TELEMAC-2D instead and it can be enough to track particles,
- if you use TREATMENT OF NEGATIVE DEPTHS = 1, you can have negative depths (h < 0 m), the comment just above is wrong,
- ACCURACY FOR PROPAGATION = 0.001 may be too coarse. Look at the listing file, the value may have to be decreased (1.E-4, 1.E-5 perhaps less).

Hope this helps,

Chi-Tuan

Hi Chi-Tuan,

Thanks for checking up and unfortunately not. This is something I'm working on in my spare time so taking me quite a while to get the hang of things.

Currently I am getting stuck with the correct wording in the user_flot file and I keep getting this error during runs:
PARTICLE 1 IN NONE OF THE SUB-DOMAINS
and no drogue output.
I believe my coordinates are definitely in the correct domain, so I think it must be the language I'm using. Would you be happy to have a quick look?

IF(LT.GE.1.AND.LT.LE.200) THEN
CALL ADD_PARTICLE(568280.2123D0,4926718.367D0,-1.5D0,
& LT,1,NFLOT,
& NFLOT_MAX,XFLOT,YFLOT,ZFLOT,TAGFLO,CLSFLO,
& SHPFLO,SHZFLO,ELTFLO,ETAFLO,MESH3D,NPLAN,
& 0.D0,0.D0,0.D0,0.D0,0,0)
ENDIF
With the above I was hoping to get *if timestep greater & equal to 1 & less than equal to 200 then call a particle on above coord* and expecting a particle/drogue to be released every timestep. Does this make sense?

Many Thanks again!

Hello Emma,

What you last sent means: for every time step between number 1 and 200, you add a particle at the same location: 568280.2123D0,4926718.367D0,-1.5D0

Can you run the following command (if you have already load a TELEMAC environment to run it) and copy/paste the result here please?
run_telfile.py scan name_geometry_file --data

If you can also do the same with one 3D result file:
run_telfile.py scan name_3D_result_file --data

Chi-Tuan