open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

I have a problem with particular initial condition for tracer.
I need to have a certain initial value of tracer (e.g. C=10) when x is higher than x*, while a zero initial tracer value otherwise.
I modified the file condin.f as the following:

IF(NTRAC.GT.0) THEN
DO ITRAC=1,NTRAC
CALL OS( 'X=C ' , X=T%ADR(ITRAC)%P , C=TRAC0(ITRAC) )
DO I=1,NPOIN
IF(X.GE.560157.D0) THEN
T%ADR(1)%P%R(I)=10.D0
ELSE
T%ADR(1)%P%R(I)=0.D0
ENDIF
ENDDO
ENDDO
ENDIF

and I get this error:

IF(X.GE.560157.D0) THEN
1
Error: IF clause at (1) requires a scalar LOGICAL expression

Can someone please help me? It is really urgent that I solve this problem!
I also attach the cas file and the condin.f file.

Thank you in advance

Best regards

Hello,

It is just that you must put X(I) instead of X in the test.

REgards,

JMH

Thank you a lot!

Now I got this error...


INBIEF (BIEF): NOT A VECTOR MACHINE (ACCORDING TO YOUR DATA)
STRCHE (BIEF): NO MODIFICATION OF FRICTION

NUMBER OF LIQUID BOUNDARIES: 2

CORFON (TELEMAC2D): NO MODIFICATION OF BOTTOM

CONDIN : WITH SPECIAL INITIAL CONDITIONS
YOU HAVE TO MODIFY CONDIN

PLANTE: PROGRAM STOPPED AFTER AN ERROR
_____________
runcode::main:
/projekte/ematta/w_res_time/tracer_condin:
|runCode: Fail to run
|/usr/bin/mpiexec -wdir /projekte/ematta/w_res_time/tracer_condin/cas2_T.txt_2015-07-06-16h59min56s -n 8 /projekte/ematta/w_res_time/tracer_condin/cas2_T.txt_2015-
07-06-16h59min56s/out_condin_e

|mpiexec has exited due to process rank 1 with PID 122634 on
|node wahyd03 exiting improperly. There are two reasons this could occur:
|
|1. this process did not call "init" before exiting, but others in
|the job did. This can cause a job to hang indefinitely while it waits
|for all processes to call "init". By rule, if one process calls "init",
|then ALL processes must call "init" prior to termination.
|
|2. this process called "init", but exited without calling "finalize".
|By rule, all processes that call "init" MUST call "finalize" prior to
|exiting or it will be considered an "abnormal termination"
|
|This may have caused other processes in the application to be
|terminated by signals sent by mpiexec (as reported here).


I don't know what could be wrong...

Thank you

BR

Hello,

INITIAL CONITIONS : SPECIAL

is meant for depth, so you should select another option, like 'CONSTANT DEPTH' or 'CONSTANT ELEVATION', or if you want to keep SPECIAL, implement your value of the depth in condin.f, just where there is the error message that stops you.

With best regards,

Jean-Michel Hervouet

Dear JMH,

thank you again.
I see... therefore, one cannot just modify condin for the tracer, am I right? One can modify condin for depth and tracer, or for neither of them. I will try...

Best regards

Hello,

If you have standard conditions for the depth, you can just change condin for the tracers, there is no need for choosing SPECIAL in that case, SPECIAL is just if you want to implement special conditions for the depth. So implement the tracer and keep the old choice for the depth.

Regards,

JMH

Ok, now it is working, thank you.

Before it was not working, because I was using a previous computation file (a bit different from the second file). If I avoid the previous computation, then it is calculating the tracer concentration as I want.

Thanks again.

Best regards

Hello,

I would like to ask, though, about the previous computation file. Is it possible to compute with condin.f modified conditions just for the tracer, starting from a previous computation where there is no tracer?
So in the first file I would just compute the hydrodynamics until steady state and in the second computation I would start from steady state conditions of the previous file, regarding the flow field, and simulating a tracer with particular conditions modified in the condin.f. Is it possible?

Thank you.

Best regards

Hello

yes this is possible. You have to introduce the intitial conditon in the bloc

IF(NTRAC.GT.0)THEN
...
ENDIF
of subroutine condin.f


with my best regards
Riadh ATA

Hello!

I am back on this case because I am applying it for a low water level mesh and a daily variable discharge imposed at the inflow (liquid boundary file).
Can you check please the files?
Probably is something minor but I cannot understand what is the problem.

I get always this error:

uncontroled error from python:: KeyError('

---

BOUNDARY CONDITIONS FILE',)
~~~~~~~~~~~~~~~~~~
Traceback (most recent call last):
File "/gfs1/work/beimatta/telemac/v6p3r2/scripts/python27/runcode.py", line 713, in runCAS
cases,lang,ncsize = processCAS(casNames,MODFiles[codeName],MODFiles[codeName])
File "/gfs1/work/beimatta/telemac/v6p3r2/scripts/python27/runcode.py", line 234, in processCAS
cas = readCAS(scanCAS(casFile),dico,frgb)
File "/gfs1/work/beimatta/telemac/v6p3r2/scripts/python27/parsers/parserKeywords.py", line 145, in readCAS
kw = frgb[kw]
KeyError: '

---

BOUNDARY CONDITIONS FILE'

Thanks. Best regards