open TELEMAC-MASCARET The mathematically superior suite of solvers

HELLO everybody
I am trying to run the overflow point of the pipe network as an input condition for TELEMAC.MY FORTRAN FILES WERE COPIED DIRECTLY FROM SOMEONE ELSE，I only changed the coordinates of the point source I needed inside。I would like to ask why the run has such an error。How can I change it to get it running?

Thank you so much

Hello,

As told many times in other posts, if a segmentation fault or a NaN (Not a Number) occurs, you should use a debug configuration with debug options to investigate.
See e.g. the S11.gfortran.debug configuration in the $HOMETEL/configs/systel.edf.cfg configuration file, in particular the flag fflags_debug_gfo for gfortran compiler.
fflags_debug_gfo: -g -Wall -fcheck=all -fbacktrace -fbounds-check -finit-integer=-1 -finit-real=nan -ffpe-trap=invalid,zero,overflow

It will show in which subroutine the issue/nan/segfault occurs, the first suspicious line and may help you to change something in your computation.

Anyway, I cannot try as there is no provided geometry file (be careful with the extensions of file allowed by this forum).

Whenever there are modifications in Fortran files, these modifications are suspicious as not standard TELEMAC implementation.
Looking quickly to your modified lecdon_telemac2d subroutine, ISCE is an array of integers and you give reals in SOURCEINDEX array, I cannot see how it can work at this step.

Hope this helps,

Chi-Tuan

this is my geometry

I will let you try to solve your issue because of many reasons:
- you do not use the latest official release (v8p5r0) or the main release and assistance is only done for these releases,
- when using Fortran files that are not official TELEMAC standard files, you have to understand every modification you do and cannot hope that taking modifications from someone else can be directly used by yourself. If you also do your own modifications, you have to understand exactly what you modify. You have to think a little bit and not expect other users can help you for basic issues,
- it will be a good experience for you to use debug options and discover how it is helpful to solve basic issues as the ones existing in your Fortran file. Segmentation faults are often due to memory allocations.

Chi-Tuan

Hi Chi-Tuan,
Thank you for your suggestion.I'M ashamed that I didn't really knowmuch about FORTRAN code and TELEMAC libraries involved in TELEMAC so far. I know this will be a waste of your time, but I would like to say that if you could please help me change this code to better understand the project I'm working on.
I'm so stupid, I don't really have any talent for it.

Hello,

To my point of view, learning Fortran is quite easier than learning TELEMAC or getting good results with TELEMAC.
You should make some efforts to learn it if you have to change some Fortran files.

Anyway, if you do not want, I can suggest you 2 possibilities:
- do not use Fortran files and try to use TELEMAC without them. It is also good advice to start from a simple model without so many things. If it works, you can try to model more complex features,
- you can contact services to assist you, you have a list (perhaps not up to date) here: opentelemac.org/index.php/support

Chi-Tuan