open TELEMAC-MASCARET The mathematically superior suite of solvers

hello all

if we want to implement instationary sources into the modell the dos-window tells them siphon-sources no. x (siphon is commented out)

the sources are steered with a fortran-file. and in the steering file of computation we defined the sources as discharge with default values.

with v5p9 the fortran files of sources overwrite the default values of discharge in the steering file, so instationary sources could be handled.
With v6p1 in parallel mode the default values in the steering file overwrite the instationary values of fortran file.

In scalar mode the instationary sources (including siphon) are working properly. the steering files are the same except the outcommented siphon and the number of processors.

cas file of parallel-mode is attached

so my question: how to implement instationary sources into v6p1-parallel mode.

thank you in advance

edit: attaching seems not to work. so in text
/ FILES USED FOR SIMULATION
/

---

/
FORTRAN FILE = stetten_100.f /f-file: this is the source-file with sources Re08 and St01
GEOMETRY FILE = Modell_stetten.net
BOUNDARY CONDITIONS FILE = Modell_stetten.lim
STEERING FILE = stetten_100_1.ste
LIQUID BOUNDARIES FILE = stetten_100.prn
RESULTS FILE = stetten_100_1.sel
/FORMATTED RESULTS FILE = schaenis_dorf_vorlauf.txt
/FORMATTED DATA FILE 1 = durchlaesse_stetten.txt


.
.
.
.
TYPE OF SOURCES = 2 /Dirac delta function should improve mass conservationmass conservation
/NUMBER OF CULVERTS = 1
/
ABSCISSAE OF SOURCES =
/665734.2500 ; 665707.9375 ; / siphon
666378.81 ; / Re08
666291.88 / St01
ORDINATES OF SOURCES =
/250514.6406 ; 250480.9688 ; / siphon
251111.95 ; / Re08
250794.66 / St01
/
WATER DISCHARGE OF SOURCES =
/0; /siphon
/0; / FORMATTED DATA FILE 1
1.0;
0.3
/
/

---

/
&ETA

Hello,

I had a look at what is done and you might be right :

in Telemac-2D DSCE is the discharge of sources as given in the steering file
and DSCE2 is the discharge varying in time. DSCE is kept because it may be used to build DSCE2 when calling function DEBSCE. Then when calling subroutine siphon.f in telemac-2d we should probably have DSCE2 instead of DSCE in the arguments. You can try to add subroutine telemac2d.f in your Fortran file and do this correction and see what happens. If it works we will correct it. However I do not catch really how parallelism can interfere with this.

With best regards,

Jean-Michel Hervouet

Dear Jean-Michel

I was wrong.... forgot to change numbering of sources after commenting the siphon out.

As well in v6p1-parallel mode sources work the same, if they are properly defined.

Thank you for your inputs. To analyse the source.f helped a lot.

best regards
daniel

Hello,

OK, but you helped me pointing out a problem, it is really DSCE2 that should be in the call to siphon, though the error is only a delay of one time step, which remained overlooked given the relaxation coefficient that also delays a lot.

With best regards,

Jean-Michel Hervouet