Welcome, Guest
Username: Password: Remember me

TOPIC: Empy Drogue File

Empy Drogue File 6 years 6 months ago #29945

  • Beya
  • Beya's Avatar
Hello,

I am trying to put drogues in a simulation and I am getting an empty DROGUE FILE. Can anyone help me with this?

See the CAS and flot.f files attached.

Also, I don't know where to locate the flot.f file. I looked at the example flotteurs and it has it in the sub-folder user_fortran-v2p0 so I created one with the version I am using (user_fortran-v7p3) but I am not sure if this is correct. I don't see where you specify the folder location for the flot.f file. For my last run I put it in the same folder as the CAS file.

I don't know very much of fortran and there are things I don't understand in the syntaxes of the float.f file such as the .D0 parts next which appear in the float.f file in the ...telemac2D/sources/ folder.

Ideally I would like to specify 100 Drogues to be released at a certain time and location but the example was for releasing 1 drogue at a time at the same position and I wouldn't have a clue how to programme that in fortran, but for now I would like to get some results.


Many thanks in advance!
Attachments:
The administrator has disabled public write access.

Empy Drogue File 6 years 6 months ago #29948

  • shenh
  • shenh's Avatar
  • OFFLINE
  • Senior Boarder
  • Posts: 148
  • Thank you received: 37
Yes in .cas file you need to specify the location of flot.f with keyword FORTRAN FILE. Before version v7.3, you must specify the name of one single user Fortran file. If you have multiple subroutines to customize, you need to combine them into one Fortran file. However with version 7.3, you do not have to combine subroutines in one file but can store all Fortran files in one folder and set keyword FORTRAN FILE to this folder.

So you just need to specify FORTRAN FILE = 'your-flot-file-folder-name' in your .cas file. Example flotteurs also includes this keyword.

Hailiang
The administrator has disabled public write access.

Empy Drogue File 6 years 6 months ago #29964

  • Beya
  • Beya's Avatar
Thank you Hailiang,

That worked!

It is now producing a non-empty drogue file using the float.f file from the flotteurs-example with new coordinates but when I want to change the time when the drogues are released it crashes at iteration 1010 with an error:

"NUMBER OF PARTICLES GREATER THAN 100"
PROGRAM STOPPED AFTER AN ERROR
RETURNING EXIT CODE : 2


The code I modified from the original float.f file from the flotteurs example and that now make telemac2d crash looks like


IF(LT.LE.2600.AND.(10*(LT/10).EQ.LT.OR.LT.EQ.2001)) THEN
CALL ADD_PARTICLE(1993783.D0,5668457.D0,0.D0,LT,NFLOT,
& NFLOT_MAX,XFLOT,YFLOT,YFLOT,TAGFLO,
& SHPFLO,SHPFLO,ELTFLO,ELTFLO,MESH,1,
& 0.D0,0.D0,0.D0,0.D0,0,0)
ENDIF

float.f file and cas file attached

File Attachment:

File Name: flot_2018-05-04.f
File Size: 3 KB

File Attachment:

File Name: Whakaki_50m_test_2018-05-04.cas
File Size: 3 KB



The code that was working before changing the times the drogues were released was:

IF(LT.LE.600.AND.(10*(LT/10).EQ.LT.OR.LT.EQ.1)) THEN
CALL ADD_PARTICLE(1993783.D0,5668457.D0,0.D0,LT,NFLOT,
& NFLOT_MAX,XFLOT,YFLOT,YFLOT,TAGFLO,
& SHPFLO,SHPFLO,ELTFLO,ELTFLO,MESH,1,
& 0.D0,0.D0,0.D0,0.D0,0,0)
ENDIF


Any clues on why it doesn't work?


Thanks again!
The administrator has disabled public write access.

Empy Drogue File 6 years 6 months ago #29967

  • c.coulet
  • c.coulet's Avatar
  • OFFLINE
  • Moderator
  • Posts: 3722
  • Thank you received: 1031
Hi

This is clearly indicate in the error message!
At iteration 1010, you try to add a particle which is the 101th particle.
By default, there is some harcoded limit in telemac to fix the maximum number of particle to 100.
So the program stop.
In your first program, you released 60 particles in the simulation.
In the latest version you plan to release 260 particles
You should adapt the code to change the limitation at least to allow 260 particles

Kind regards
Christophe
The administrator has disabled public write access.
The following user(s) said Thank You: Beya
Moderators: pham

The open TELEMAC-MASCARET template for Joomla!2.5, the HTML 4 version.