open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello Hailiang,

I ran the "estu_gir_3D" example in Telemac3D and attempted using your script with the drogue result. However, I have encountered some error. Please find attached the error screen.

Regards,
Chijioke.

You just need to edit below
c = OilParticle2Pcl(r"C:\Users\Chijioke\Desktop\estu_gir_3D\test\t3d_particles.dat", 1)

I have the Nmass variable default value 4 which is for my case. In your example, you only have one mass. I should make the program smarter in next version, Nmass is actually the length of mass list in class OilDrogue.

If you have any other issues in using the script, just let me know.

Hailiang

Thank you Hailiang. You have been bery helpful.

Regards,
Chijioke.

That current_time 37800.0000 is the first time stamp in my displacement file. You could use another one time stamp in the file - not it must exist otherwise the program will fail (I did not catch this error in my code). BlueKenue supports reading xyz files. Function to_xyz_file(...) can output one that allows you to quickly visualize your drogues at a time stamp.

The other input values in function test() are also based on my displacement file. You need to edit them for your model.

Hailiang

Can you post your file t3d_particles.dat?

Hailiang

Hello Hialiang,

Find attached my t3d_particle file.
Using the script, I have successfully generated the xyz file. However, the pcl file produces an error when I try to read it with BlueKenue.

Regards,
Chijioke.

Chijioke, the two files are not what i wanted. I was referring to the .dat file which includes drogue locations, see attached example.





Hailiang

Sorry, about that. Didn't realise .dat will not be uploaded.
Please find attached the file.

Hello,

I am trying to run an oil spill simulation with TELEMAC3D, however, the oil spill (drogue) output file is empty. I have uploaded the oil_flot.f and steering file.

Please does anyone know what the problem might be.

Thank you.

Regards
Chijioke.

A few initial points:
- is particle release point (x=534303 and y=409994) inside your model domain?
- LT.EQ.7200 will never be hit and therefore no particles will be added, because your number of time steps in steering file is only 320. To quickly test the subroutine, change to a smaller number say LT.EQ.20.
- in subroutine ADD_PARTICLE, the third variable is z which should be the release elevation to identify vertical layer ID. Using COORD_Y is not correct.

Hope these help you get started,

Hailiang