0

ERROR: STOP ITOUGH2 : Signal 1 recevied.

Hello everyone,

i'm a student, and i'm working on changing the primary variable in iTOUGH2, after some changing in file EOSX and T2F.F and IT2MAIN.F , i compile my executable and when i want to test my modification i got this 'ERROR:  STOP ITOUGH2 : Signal 1 recevied.' i know this is a message of subroutine signal 1, due to NBREAK=1.

1-i hope that someone will tell me where i made the error , 

2-i'm not familliar with the code but i run my code on mingw ( gcc compiler) i don't know the comande for just a forward solving and not inverse solving, so wich parameter i should fixe to get a forward solving like tOUGH2

 

thank you ! 

14 replies

null
    • Finsterle GeoConsulting
    • Stefan_Finsterle
    • 4 yrs ago
    • Reported - view

    Hi Amri,

    The error message most likely occurs because of an array-bound violation somewhere in the new code section you added. Check your modifications carefully. Here are some other suggestions:

    (1) Compile in debug mode, so you may get traceback information.

    (2) In case you use a newer version of iTOUGH2, add the following command to the iTOUGH2 input file:

    > COMPUTATION

      >> STOP

         >>> ignore SIGNAL

         <<<

    This disables the signal handler, so you should get the line number in the source code where the exception occurs.

    (3) I don't know whether you are running (i)TOUGH2 using the script files you find in ~/itough2/Scripts, which is the preferred way. If so, just type

    tough2 input_file 3

    to run iTOUGH2 in forward mode (add directory Scripts to your command search path).

    If you run the executable directly, provide a dummy iTOUGH2 input file (it actually can be an empty file; or see file ~/itough2/SpecialFiles/invdir) to run iTOUGH2 in forward mode.

    All of this is described in great detail in the iTOUGH2 read.me file - please read it!

     

    Best,

    Stefan

      • woemn_geoscience
      • 4 yrs ago
      • Reported - view

      Stefan Finsterle 

      Dear Finsterle,

      First I want to thank you for your answer,

      I can't find the file read.me in my folder, and i want to know if it's possible to solve with (PG,PL,T) and not (PG,Sg,T) in the EOS3 ?, There are places where I change other than EOS3.f and PC-SL in T2F.f ?

      Thanks ! 

    • Finsterle GeoConsulting
    • Stefan_Finsterle
    • 4 yrs ago
    • Reported - view

    Amri,

    I don't know which version of the code you have. The distribution license by LBNL includes the read.me file.

    Please study the primary variable choices for single-phase and two-phase conditions, which explains why (PG,PL,T) is not a suitable combination (specifically if capillary pressures are zero).

    It is not quite clear why you want to change primary variables. I hope you have a compelling reason, as it requires substantial changes in the EOS module and all subroutines that access the primary variable vector X.

    Best,

    Stefan

      • woemn_geoscience
      • 4 yrs ago
      • Reported - view

      Stefan Finsterle 

      Dear Finsterle,

      I succed to change the primary variable in the eos3  my new code solve with (P,X,T) in single phase and (Pg,Pl,T) in two-phase, i have some difference betwen the old calcul and the new ( in Sl  and not PCAP, Xair) (perehapes it's du to some call of incremented and last update variables) can you send me please the file where there's explication why (PG,PL,T) is not a suitable combination (specifically if capillary pressures are zero). i knew that but  when we suppose THERE IS a residual saturation of gas and we can do so ? 

       

      i want i's to take into accont of the entry air  pressure in the  Cox

      this is a capillary equilibrum -10 m =PC  in the bottom 

    • woemn_geoscience
    • 4 yrs ago
    • Reported - view
    • Finsterle GeoConsulting
    • Stefan_Finsterle
    • 4 yrs ago
    • Reported - view

    Amri,

     

    If you don't have capillary pressures, Pl=Pg for all saturations, so you end up with only two equations for three primary variables.

    I'm sorry, but I did not understand the other questions/requests you had.

    Stefan

    • woemn_geoscience
    • 4 yrs ago
    • Reported - view

    Dear Finsterle, 

    I would like to thank you for your answers, and I am taking into account your answers, but I can not understand when Pc = 0.0 the code changes the resolution in (P, X, T) so why there will be a problem ( number of equations and number of unknowns)?

    • Finsterle GeoConsulting
    • Stefan_Finsterle
    • 4 yrs ago
    • Reported - view

    Amri,

    I assume with "resolution" you mean "solution variables", i.e., the set of primary variables, which in TOUGH2 depends on the phase state, i.e., whether you have two-phase (gas and liquid) conditions or single-phase (gas or liquid) conditions. The choice does not depend on Pc, as you seem to imply.

    Your choice of the set (Pg, Pl, T) for two-phase conditions is acceptable, but collapses if Pc(0<Sl<1) = 0.0 or very small or very flat.

    Again, assuming I understand you correctly, you wanted this alternative set because you want to handle a finite air-entry pressure at Sg=0. Your choice is fine for that. There may be other alternatives, but I don't know about the relative stability of each of those options. Please let us know how it goes!

    Stefan

      • woemn_geoscience
      • 4 yrs ago
      • Reported - view

      Stefan Finsterle 

      Dear Stefan, 

      thank you again for your relevant answer, and i'm wondering to know wich alternative you think i can use?

      • Finsterle GeoConsulting
      • Stefan_Finsterle
      • 4 yrs ago
      • Reported - view

      Abdellah Amri 

      (PL, Sg, T) or no change in primary variables but adjustments of the Pc curve near full saturation or simply specifying a small irreducible gas saturation to mimic the effect of a finite air-entry pressure.

      Stefan

      • woemn_geoscience
      • 4 yrs ago
      • Reported - view

      Stefan Finsterle 

      in two-phase flow, my code is good, when there is an iteration where the code computes in liquid. I started to have bad results :(

      • woemn_geoscience
      • 4 yrs ago
      • Reported - view

      Stefan Finsterle thanks stefan

      • Finsterle GeoConsulting
      • Stefan_Finsterle
      • 4 yrs ago
      • Reported - view

      Abdellah Amri 

      Amri,

      Check proper variable initialization after a change from two-phase to single-phase liquid conditions.

      (I hope you understand that I cannot effectively "demote-debug" iTOUGH2 after changes have been made by a user.)

      Good luck!

      Stefan 

      • woemn_geoscience
      • 4 yrs ago
      • Reported - view

      Stefan Finsterle 

      Finsterle,

      Thank you for your suggestion, i try to do in other way thanks ! 

Content aside

  • 4 yrs agoLast active
  • 14Replies
  • 66Views
  • 2 Following