open TELEMAC-MASCARET The mathematically superior suite of solvers

I need to simulate eg. 21000 passive particles transport released in same time in a defined area.
The Telemac-3D have a limited number of particles, no much than 10000 released in my area. So, a think that putting 0 to 7000 particles at surface (Z=0), 7001 to 14000 particles at Z= -5, and 14001 to 21000 particles at Z=-10, can solve my problem. The Total depth is about 18 m.
In my create_3dtracers.py I changed the ID_Start and Z, using the same coordinates. So my simulation don't run, only lock.
Can someone help me? How can I set up, and simulate passive particle putting in same time particles in different depth?
Is necessary to change the NPLAN or something?

Hello,

Which version of TELEMAC do you use?
What message is written in the listing file?

Could it be possible to have your STEERING FILE + FORTRAN FILE?

Chi-Tuan

Thanks Plam for your replay

I used Telemac V7P0

See attached the simulation documents I already ran.
I put particles in Z=-5.
But I also want to put particles in Z = -10 and Z = -15 at the same coordinates and the same simulation.
How can I do that?

Hello,

Wow, I am not sure I have ever read a FORTRAN FILE with more thant 30,600 lines. I wonder if there is a maximum number of lines for this file or not.

Anyway, why have not you written a loop to define your ADD_PARTICLE statements in FORTRAN as you did in your external Python script? Maybe it can solve the problem?

Have you tried to run a simple simulation with 3 particles at different depths,
e.g.:

CALL ADD_PARTICLE(393553.D0,6437736.D0,-5.D0,1,NFLOT,
& NFLOT_MAX,XFLOT,YFLOT,ZFLOT,TAGFLO,
& SHPFLO,SHZFLO,ELTFLO,ETAFLO,MESH3D,NPLAN,
& 0.D0,0.D0,0.D0,0.D0,0,0)
CALL ADD_PARTICLE(393553.D0,6437736.D0,-10.D0,1,NFLOT,
& NFLOT_MAX,XFLOT,YFLOT,ZFLOT,TAGFLO,
& SHPFLO,SHZFLO,ELTFLO,ETAFLO,MESH3D,NPLAN,
& 0.D0,0.D0,0.D0,0.D0,0,0)
CALL ADD_PARTICLE(393553.D0,6437736.D0,-15.D0,1,NFLOT,
& NFLOT_MAX,XFLOT,YFLOT,ZFLOT,TAGFLO,
& SHPFLO,SHZFLO,ELTFLO,ETAFLO,MESH3D,NPLAN,
& 0.D0,0.D0,0.D0,0.D0,0,0)

Moreover, I would strongly recommend you to use the latest v8p0r3 version (or the new branch v8p1) as many bugs have been solved since V7P0. We cannot assist users for too old versions.

Hope this helps,

Chi-Tuan