0

Error in reading INCON - Problem 5, eco2n_1D Radial

Hello, 

In problem 5, eco2n_1D Radial, I am using the input files provided (flow, solute, chemical, CO2TAB, .dat, and executable) to run a q=0 kg/s flow in order to find the equilibrium conditions (pH specifically). I run this simulation for 3000 years and find that the pH is approximately around the 9.1 region. 

I would like to use the 3000th year data as my initial conditions to run another simulation that involves injecting CO2 at a rate of 90 kg/s for 10 years and letting it equilibrate for 1000 years. After the 3000 year simulation is run,  I copy the input files along with SAVE and savechem and change their names to INCON and inchem respectively, and change the q to 90 in GENER and run the simulation to see how much pH changes with the CO2 injection, I get an error that says INCON. 

Can somebody help me figure out what I might be doing that is wrong?

Thank you!

6replies Oldest first
  • Oldest first
  • Newest first
  • Active threads
  • Popular
  • Hi Majid,

    Make sure the last time in GENER is 3000 years, and is equal to or greater than the final time in flow.inp. If you are still having troubles,  can you send the modified input files, and I can look at them? 

    regards,

    Eric

    Like
      • Majid Trf
      • Majid_Trf
      • 4 days ago
      • Reported - view

      Eric Sonnenthal Thank you for your quick response Eric. 

      Here are the modified files.

      In PARAM, the length of time to run the simulation to should be 3000 years. 

      In GENER, I have 90 kg/s flowing for 10 years and equilibrating for 2990 years. I made sure GENER and PARAM have the same simulation length of time. 

      I have also deleted INCON in flow.inp and made SAVE become INCON. 

      Like
    • Hi Majid,

      At the bottom of SAVE, there are *** because the number of time steps and iterations exceeded the fortran format. To restart a file from time 0, delete these two lines, and enter at least 2 blank lines instead.

      +++  
      ************  1 0.00000000E+00 0.94608000E+11

      If you need to restart from a specific time, then reset the time steps and iterations to 1, so it doesn't crash on reading it. Make sure the format is the same (2i6), so it looks like:

      +++  
           1     1  1 0.00000000E+00 0.94608000E+11

      regards,

      Eric

      Like 1
      • Majid Trf
      • Majid_Trf
      • 4 days ago
      • Reported - view

      Eric Sonnenthal This worked perfectly. 

      I am a little bit confused about what it means to restart from a specific time. If I want the conditions of the 3000th year to be used in the next simulation I have to follow the second part (2i6) with the time steps and iterations, correct? 

      Thank you, 

      Majid

      Like
    • Hi Majid,

      If you remove the last two lines, then the starting time is back to 0, but the conditions are whatever is in the INCON. TOUGHREACT looks at the starting simulation time and then looks in GENER to see what the injection rate should be. If the time at the bottom of INCON is 3000 years, and the last time in GENER is 3000 years, then TOUGHREACT will think that simulation is over. If you want to restart a simulation from time 0 and run it another 3000 years, then you need to remove those 2 lines. If you want to run it for 100 years, then restart it and run it another 2900 years, then the last time on the bottom of INCON needs to be 100 years, so TOUGHREACT knows what to values to take from GENER.

      regards,

      Eric

      Like 1
      • Majid Trf
      • Majid_Trf
      • 3 days ago
      • Reported - view

      Eric Sonnenthal Thank you for your thorough answer. This was very helpful. 

      Kind regards,

      Majid

      Like
Like Follow
  • 3 days agoLast active
  • 6Replies
  • 17Views
  • 2 Following