open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

Classical error. As it is written "simulation en Francais" because the first keyword of the steering file is in French, but you have at least one keyword not known (at least WATER QUALITY PROCESS) because not in French, but English.

Do you know the dictionaries available in the sources folders (*.dico) or the reference manuals, available in the documentation folders or e.g. here : wiki.opentelemac.org/doku.php?id=documentation_v8p1r0

Moreover, if you use release 8.1, the numbering of the water quality processes has changed since a few versions. You should then read the WAQTEL user manual.

Hope this helps,

Chi-Tuan

Hi,
I have new error appear NaN
I think the problem can be in dynamic pressure, any idea?
Thanks a lot
No consistent outcomes

---

CALCULATION OF VERTICAL VELOCITY STEP

---

ADVECTION-DIFFUSION OF TRACERS
MURD3D OPTION: 4 1 ITERATIONS
DIFFUSION OF TRN1
GMRES (BIEF) : 1 ITERATIONS, RELATIVE PRECISION: NaN

---

ADVECTION-DIFFUSION OF VELOCITIES STEP
MURD3D OPTION: 4 1 ITERATIONS
MURD3D OPTION: 4 1 ITERATIONS

---

PROPAGATION AND DIFFUSION WITH WAVE EQUATION
GMRES (BIEF) : 1 ITERATIONS, RELATIVE PRECISION: NaN
POSITIVE DEPTHS OBTAINED IN 1 ITERATIONS

---

CALCULATION OF VERTICAL VELOCITY STEP

---

ADVECTION-DIFFUSION OF TRACERS
MURD3D OPTION: 4 1 ITERATIONS
DIFFUSION OF TRN1
GMRES (BIEF) : 1 ITERATIONS, RELATIVE PRECISION: NaN

---

ADVECTION-DIFFUSION OF VELOCITIES STEP
MURD3D OPTION: 4 1 ITERATIONS
MURD3D OPTION: 4 1 ITERATIONS

---

PROPAGATION AND DIFFUSION WITH WAVE EQUATION
GMRES (BIEF) : 1 ITERATIONS, RELATIVE PRECISION: NaN
POSITIVE DEPTHS OBTAINED IN 1 ITERATIONS

---

CALCULATION OF VERTICAL VELOCITY STEP

---

ADVECTION-DIFFUSION OF TRACERS
MURD3D OPTION: 4 1 ITERATIONS
DIFFUSION OF TRN1
GMRES (BIEF) : 1 ITERATIONS, RELATIVE PRECISION: NaN

---

ADVECTION-DIFFUSION OF VELOCITIES STEP
MURD3D OPTION: 4 1 ITERATIONS
MURD3D OPTION: 4 1 ITERATIONS

---

PROPAGATION AND DIFFUSION WITH WAVE EQUATION
GMRES (BIEF) : 1 ITERATIONS, RELATIVE PRECISION: NaN
POSITIVE DEPTHS OBTAINED IN 1 ITERATIONS

---

CALCULATION OF VERTICAL VELOCITY STEP

---

ADVECTION-DIFFUSION OF TRACERS
MURD3D OPTION: 4 1 ITERATIONS
DIFFUSION OF TRN1
GMRES (BIEF) : 1 ITERATIONS, RELATIVE PRECISION: NaN

---

ADVECTION-DIFFUSION OF VELOCITIES STEP
MURD3D OPTION: 4 1 ITERATIONS
MURD3D OPTION: 4 1 ITERATIONS

---

PROPAGATION AND DIFFUSION WITH WAVE EQUATION
GMRES (BIEF) : 1 ITERATIONS, RELATIVE PRECISION: NaN
POSITIVE DEPTHS OBTAINED IN 1 ITERATIONS

---

CALCULATION OF VERTICAL VELOCITY STEP

---

ADVECTION-DIFFUSION OF TRACERS
MURD3D OPTION: 4 1 ITERATIONS
DIFFUSION OF TRN1
GMRES (BIEF) : 1 ITERATIONS, RELATIVE PRECISION: NaN

---

ADVECTION-DIFFUSION OF VELOCITIES STEP
MURD3D OPTION: 4 1 ITERATIONS
MURD3D OPTION: 4 1 ITERATIONS

---

PROPAGATION AND DIFFUSION WITH WAVE EQUATION
GMRES (BIEF) : 1 ITERATIONS, RELATIVE PRECISION: NaN
POSITIVE DEPTHS OBTAINED IN 1 ITERATIONS

---

CALCULATION OF VERTICAL VELOCITY STEP

---

ADVECTION-DIFFUSION OF TRACERS
MURD3D OPTION: 4 1 ITERATIONS
DIFFUSION OF TRN1
GMRES (BIEF) : 1 ITERATIONS, RELATIVE PRECISION: NaN

---

ADVECTION-DIFFUSION OF VELOCITIES STEP
MURD3D OPTION: 4 1 ITERATIONS
MURD3D OPTION: 4 1 ITERATIONS

---

PROPAGATION AND DIFFUSION WITH WAVE EQUATION
GMRES (BIEF) : 1 ITERATIONS, RELATIVE PRECISION: NaN
POSITIVE DEPTHS OBTAINED IN 1 ITERATIONS

---

CALCULATION OF VERTICAL VELOCITY STEP

---

ADVECTION-DIFFUSION OF TRACERS
MURD3D OPTION: 4 1 ITERATIONS
DIFFUSION OF TRN1
GMRES (BIEF) : 1 ITERATIONS, RELATIVE PRECISION: NaN

---

ADVECTION-DIFFUSION OF VELOCITIES STEP
MURD3D OPTION: 4 1 ITERATIONS
MURD3D OPTION: 4 1 ITERATIONS

---

PROPAGATION AND DIFFUSION WITH WAVE EQUATION
GMRES (BIEF) : 1 ITERATIONS, RELATIVE PRECISION: NaN
POSITIVE DEPTHS OBTAINED IN 1 ITERATIONS

---

CALCULATION OF VERTICAL VELOCITY STEP

---

ADVECTION-DIFFUSION OF TRACERS
MURD3D OPTION: 4 1 ITERATIONS
DIFFUSION OF TRN1
GMRES (BIEF) : 1 ITERATIONS, RELATIVE PRECISION: NaN

---

MASS BALANCE

WATER
MASS AT THE PREVIOUS TIME STEP : 21029.69
MASS AT THE PRESENT TIME STEP : 21029.69
MASS LEAVING THE DOMAIN DURING THIS TIME STEP : NaN
ERROR ON THE MASS DURING THIS TIME STEP : NaN
BALANCE RAIN-EVAPORATION : 0.000000
FLUX BOUNDARY 1: NaN M3/S ( >0 : ENTERING <0 : EXITING )
FLUX BOUNDARY 2: NaN M3/S ( >0 : ENTERING <0 : EXITING )

TRACER 1: TEMPERATURE C , UNIT : * M3)

ADVECTIVE FLUX THROUGH BOUNDARIES OR SOURCES : NaN
DIFFUSIVE FLUX THROUGH THE BOUNDARIES : NaN
MASS AT THE PREVIOUS TIME STEP : NaN
MASS AT THE PRESENT TIME STEP : NaN
MASS EXITING (BOUNDARIES OR SOURCE) : NaN
ERROR ON THE MASS DURING THIS TIME STEP : NaN

---

FINAL MASS BALANCE
T = 0.0300

--- WATER ---
INITIAL MASS : 21029.61
FINAL MASS : 21029.69
MASS LEAVING THE DOMAIN (OR SOURCE) : NaN
MASS LOSS : NaN

--- TRACER 1: TEMPERATURE C , UNIT : * M3)
INITIAL MASS : 358222.3
FINAL MASS : NaN
MASS EXITING (BOUNDARIES OR SOURCE) : NaN
MASS LOSS : NaN

EXITING MPI

CORRECT END OF RUN

ELAPSE TIME :
4 SECONDS


~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

... merging separated result files

+> cas_rejet_exchange.cas
collecting: T3DRES
STOP 0


~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

... handling result files

+> cas_rejet_exchange.cas
moving: res3D_stratif_21.slf


My work is done

I tried already:

1- When ran hydrodynamics is working.
2- When ran hydrodynamics and temperature without waqtel is working.
3- Differents time step
4- Refine the mesh
NOTHING WORK

Hello,

With dry zones, the ATMOSPHERE-WATER EXCHANGE MODEL = 2 currently does not work.

You could try ATMOSPHERE-WATER EXCHANGE MODEL = 1 either, it may better work with dry zones.

Chi-Tuan

Hi Ch-Tuan

Does ATMOSPHERE-WATER EXCHANGE MODEL = 2 work for the latest version on the trunk ?

If not, do you know what it is that is causing the problem ? Is your team actively looking at it ? I can't see anything on CUE.

Elfed

Hello Elfed,

No fix for ATMOSPHERE-WATER EXCHANGE MODEL = 2 and tidal flats for the trunk.
It is on our todo list of bugs to solve. But any help is welcome.

Chi-Tuan

Does not work either with ATMOSPHERE-WATER EXCHANGE MODEL = 1
I would like to ask something, What is the difference between run exchange without coupling with water quality, and run with waqtel? The way to board boundary conditions? why you are suggesting tidal flat?
For me is not clear which exactly it is the numerical problem, I am working in my thesis and I want to understand and help.
I am doing a real case in a reservoir and first of all, I was trying with a small example.
Thank you very much

I tried to change the keyword FORMATTED DATA FILE 1 for ASCII ATMOSPHERIC DATA FILE And new error appear, I saw this suggestion in another post.




+> /PRODCOM/Ubuntu16.04/telemac/v7p2r3/builds/ubu16gnumpich/bin/partel < PARTEL.PAR >> partel_T3ATMB.log
Attempting to use an MPI routine before initializing MPICH
copyPartition:
|runPARTEL: Could not split your file T3ATMB (runcode=1) with the error as follows:
|
|... The following command failed for the reason above (or below)
|/PRODCOM/Ubuntu16.04/telemac/v7p2r3/builds/ubu16gnumpich/bin/partel < PARTEL.PAR >> partel_T3ATMB.log
|
| You may have forgotten to compile PARTEL with the appropriate compiler directive
| (add -DHAVE_MPI to your cmd_obj in your configuration file).
|
|Here is the log:
|
|
| +

---

+
|
| PARTEL/PARRES: TELEMAC METISOLOGIC PARTITIONER
|
|
|
| REBEKKA KOPMANN & JACEK A. JANKOWSKI (BAW)
|
| JEAN-MICHEL HERVOUET (LNHE)
|
| CHRISTOPHE DENIS (SINETICS)
|
| YOANN AUDOUIN (LNHE)
|
| PARTEL (C) COPYRIGHT 2000-2002
|
| BUNDESANSTALT FUER WASSERBAU, KARLSRUHE
|
|
|
| METIS 5.0.2 (C) COPYRIGHT 2012
|
| REGENTS OF THE UNIVERSITY OF MINNESOTA
|
|
|
| BIEF V7P2R0 (C) COPYRIGHT 2012 EDF
|
| +

---

+
|
|
|
|
|
| MAXIMUM NUMBER OF PARTITIONS: 100000
|
|
|
| +

---

+
|
|
|
| --INPUT FILE NAME <INPUT_NAME>:
|
| INPUT: T3ATMB
|
| --INPUT FILE FORMAT <INPFORMAT> [MED,SERAFIN,SERAFIND]:
|
| INPUT: SERAFIN
|
| --BOUNDARY CONDITIONS FILE NAME:
|
| INPUT: T3DCLI
|
|--NUMBER OF PARTITIONS <NPARTS> [2 -100000]:
|
| INPUT: 16
|
| PARTITIONING METHOD <PMETHOD> [1 (METIS) OR 2 (SCOTCH)]:
|
| --INPUT: 1
|
| --CONTROL SECTIONS FILE NAME (OR RETURN) :
|
| NO SECTIONS
|
| --CONTROL ZONES FILE NAME (OR RETURN) :
|
| NO ZONES
|
| --WEIR FILE NAME (OR RETURN) :
|
| NO WEIRS
|
| --GEOMETRY FILE NAME <INPUT_NAME>:
|
| INPUT: T3DGEO
|
| --GEOMETRY FILE FORMAT <GEOFORMAT> [MED,SERAFIN,SERAFIND]:
|
| INPUT: SERAFIN
|
| +---- PARRES: BEGINNING

---

+
|
| +---- USING GEOMETRY FILE: T3DGEO
|
| ERREUR -1 LORS DE L APPEL A IDENTIFY_ENDIAN_TYPE:READ:LITTLE_ENDIAN
|
| TEXTE DE L'ERROR : UNKONWN ERROR
|
|
|
|
|
|
|
| PLANTE: PROGRAM STOPPED AFTER AN ERROR
|
| RETURNING EXIT CODE:

Hello,

If you want to use the features described in the WAQTEL technical manual (or WAQTEL theory guide) available in wiki.opentelemac.org/doku.php?id=documentation_v8p1r1
In your case, if you may want to use the THERMIC module (see chapter 6), you will be able to compute boundary conditions and source terms for a complete balance of exchange fluxes, but you need to couple with WAQTEL. If not, you will only solve a classical advection-diffusion equation for temperature with classical boundary conditions at open and solid lateral boundaries and nothing special with the free surface (in particular the exchange with atmosphere).

I talked about tidal flats o know if you have some in your models or not (area that could be always or sometimes dry).
Because if you have some, the complete balance of heat fluxes does not currently work with dry areas.

Hope it is clearer now.

Chi-Tuan

Dear Juli,

I try to sumup:
- you use v8p1r1 version?
- your test case with hydro alone, and dry zones => OK
- your test case with hydro and temperature (as a passive or active tracer?) => OK
- your test case with hydro, temperature and waqtel (heat exchange at the atmosphere = 1) => NOK

Did you try a testcase with hydro+waqtel+heat exachange=1) without dry zones?

Best regards,