problem using usini1 and CS_2.1.4

[sylvain_serra]

Hello everyone,

I'm trying to run a case on saturne 2,1,4.

This case use a particular physics « gas combustion ».
I don't use the GUI

so I add in the:
DATA directory, the dp_C3P
cs_user_scripts.py files

SRC direcory, the usini1.f90
usppmo.f90
usebu1.f90
usebui.f90
usebuc.f90 files

@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
@
@ @@ WARNING: STOP WHILE READING INPUT DATA
@ =========
@ IECAUX MUST BE AN INTEGER EQUAL 0 OR 1
@ IT HAS VALUE -512647296
@
@ The calculation could NOT run.
@
@ Check the input data given via User Interface or in usini1.
@
@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
@

@
@
@
@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
@
@ @@ WARNING: ABORT IN THE DATA SPECIFICATION
@ ========
@ THE CALCULATION PARAMETERS ARE INCOHERENT OR INCOMPLET
@
@ The calculation will not be run ( 1 errors).
@
@ See previous impressions for more informations.
@ Verify the provided data in the interface, usini1 or
@ the other initialization subroutines.
@
@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
@


Pile d'appels :
1: 0x2b9e436787ac <initi1_+0xec> (libsaturne.so.0)
2: 0x2b9e435d8117 <cs_run+0xc7> (libsaturne.so.0)
3: 0x2b9e435d8985 <main+0x1e5> (libsaturne.so.0)
4: 0x39aac1d994 <__libc_start_main+0xf4> (libc.so.6)
5: 0x4038c9 <main+0x131> (cs_solver)
Fin de la pile

-------------------------------------

In attachment, you could find my case
I checked that iecaux = 1 in usini1.f90 and it is ok

Can someone help me?
Thank you very much

Sylvain

[Yvan Fournier]

Hello,

Using a the code compiled in debug mode (with bounds checking enabled), I get an earlier error:

At line 105 of file /home/yvan/Code_Saturne/saturne-2.1/src/cogz/ebuver.f90
Fortran runtime error: Index '0' of dimension 1 of array 'visls0' below lower bound of 1.
Looking at the source code, it seems that the variable "ihm", used for enthalpy, is not set (it is set at 0, and should be >= 4, as we usually have at least pressure and velocity, + turbulence variables).

I would imagine that with combusting activated, Enthalpy would automatically be activated, but you may try to force this in usini1.f90.

I suspect a bug, as the update of specific physics from version 2.0 to 2.1 or 2.2 have not received exhaustive testing yet, but I'll try to check within a few days with the Gas combustion people, and suggest to them to take a look at the Forum.

Best regards,

Yvan

[Yvan Fournier]

I just tried running with bounds-checking enabled, and do not reproduce your error, so the issue I had in the previous post is probably only a minor bug.

I am not running with version 2.1.4, but the 2.1 branch head, which you could obtain using:
svn checkout http://code-saturne.org/svn/saturne/branches/Version2_1 saturne-2.1
then running ./sbin/bootstrap in saturne-2.1 (but which requires developper tools such as the GNU autotools, gettext, bison, flex, ... to be present)

As there have been very few changes since 2.1.5, I would simply recommend installing 2.1.5, which you can download from this site.

Best regards,

Yvan

[sylvain_serra]

Thank you very much for your answers Yvan.