0

"EOS CANNOT FIND PARAMETERS AT ELEMENT" error message

Hi, dear all

I'm working on geothermal model by using TOUGH2 with PetraSim.

Now I had generated input file from PetraSim (0114.dat) and recompile the eos1.exe file.

I try to output the temperature restlt at last time setp at  3.1536E13 (s).

But the model seems to not complete executed and the out file showed error message like below:

++++++++++ EOS CANNOT FIND PARAMETERS AT ELEMENT *410 9*      XX(M) = -0.70049E+07  0.31783E+02                                            
+++++++++ REDUCE TIME STEP AT (   1,  2) ++++++++++++++++   NEW DELT = 0.250000E+03                                                        
++++++++++ EOS CANNOT FIND PARAMETERS AT ELEMENT *412 5*      XX(M) = -0.74834E+06  0.23701E+02                                            
+++++++++ REDUCE TIME STEP AT (   1,  2) ++++++++++++++++   NEW DELT = 0.625000E+02                                                        
++++++++++ EOS CANNOT FIND PARAMETERS AT ELEMENT *41185*      XX(M) = -0.17859E+10  0.23710E+02                                            
+++++++++ REDUCE TIME STEP AT (   1,  2) ++++++++++++++++   NEW DELT = 0.156250E+02                                                        
++++++++++ EOS CANNOT FIND PARAMETERS AT ELEMENT *41185*      XX(M) = -0.87538E+08  0.24417E+02                                            
+++++++++ REDUCE TIME STEP AT (   1,  2) ++++++++++++++++   NEW DELT = 0.390625E+01                                                        
++++++++++ EOS CANNOT FIND PARAMETERS AT ELEMENT *41185*      XX(M) = -0.28586E+11  0.23710E+02                                            
+++++++++ REDUCE TIME STEP AT (   1,  2) ++++++++++++++++   NEW DELT = 0.976562E+00                                                        
39221(   1,  1) ST = 0.976562E+00 DT = 0.976562E+00 DX1= 0.000000E+00 DX2= 0.000000E+00 T =  31.834 P =  2335912. S = 0.000000E+00 ICURV=0 
39221(   2,  1) ST = 0.292969E+01 DT = 0.195312E+01 DX1= 0.000000E+00 DX2= 0.000000E+00 T =  31.834 P =  2335912. S = 0.000000E+00 ICURV=0 
++++++++++ EOS CANNOT FIND PARAMETERS AT ELEMENT *41185*      XX(M) = -0.28586E+11  0.23710E+02                                                                                        
+++++++++ CONVERGENCE FAILURE ON TIME STEP #   3 WITH DT = 0.390625E+01 SECONDS, FOLLOWING TWO STEPS THAT CONVERGED ON ITER = 1            
          STOP EXECUTION AFTER NEXT TIME STEP                                                                                              
+++++++++ REDUCE TIME STEP AT (   3,  2) ++++++++++++++++   NEW DELT = 0.976562E+00                                                        
39221(   3,  1) ST = 0.390625E+01 DT = 0.976562E+00 DX1= 0.000000E+00 DX2= 0.000000E+00 T =  31.834 P =  2335912. S = 0.000000E+00 ICURV=0 
 

As mention before, the input file was generated from PetraSim.  In PetraSim, the model can execute complete. so I don't know what's the input file error could caused these error message. 

 

In addition,when I examined my output file, I saw the WARNING message like below:

 WWWWWWWWWW WARNING WWWWWWWWWW: NUMBER OF SIGNIFICANT DIGITS IS INSUFFI
 WWWWWWWWWWWWWWWWWWWWWWWWWWWWW: CODE SHOULD BE RUN IN DOUBLE PRECISION!
And I also don't know whether these warning message will caused the error message  as mentioned before or not.

 I would be appreciated if someone could let me know where I am doing wrong.

Many Thanks

Yane

3replies Oldest first
  • Oldest first
  • Newest first
  • Active threads
  • Popular
    • Andri Arnaldsson
    • Senior Geothermal Scientist - Vatnaskil Consulting Engineers
    • Andri_Arnaldsson
    • 6 yrs ago
    • Reported - view

    What compiler do you use? Make sure to use the appropriate flags or options to generate a double precision executable.

    Like
  • Hi,

    I used the g95 for compile the .exe file, and it's my first time to compile so I don't know how to solve the double precision executable problem. Did that the reason caused  the model didn't execute complete?

    Like
    • Andri Arnaldsson
    • Senior Geothermal Scientist - Vatnaskil Consulting Engineers
    • Andri_Arnaldsson
    • 6 yrs ago
    • Reported - view

    I would recommend using gfortran (via MinGW or cygwin on Windows) rather than g95 since I think the g95 project is now semi-dormant. But you can try to add the -r8 option.

    http://www.g95.org/docs.shtml

    Like
Like Follow
  • 6 yrs agoLast active
  • 3Replies
  • 1176Views
  • 2 Following