I am trying to model a salt intrusion in a confined aquifer using TOUGH3 and EO7, following Henry's problem.

I am defining initial conditions (with multi 2 2 2 8: only water and brine) that allow simulating full saturation in water. When I add the diffusion parameters, I obtain the following error:

MESH: NEL=         800  NCON=        1540
 TOUGH STATUS: Preprocessing step completed.
 TOUGH STATUS: Total preprocessing time =   2.999612502753735E-002
forrtl: severe (408): fort: (2): Subscript #2 of the array FDDIAG has value 3 which is greater than the upper bound of 2

Image              PC                Routine            Line        Source
tough3-eos7        000000011045004D  for_emit_diagnost     Unknown  Unknown
tough3-eos7        000000010F2D9815  eos_                  Unknown  Unknown
tough3-eos7        000000010F357B23  cycit_                Unknown  Unknown
tough3-eos7        000000010F2E9C03  MAIN__                Unknown  Unknown
tough3-eos7        000000010F2D78EE  main                  Unknown  Unknown

However, if I initiate the simulation using  multi 3 3 2 8, TOUGH3 does recognize the diffusion parameters I assigned.

Could you explain why it does not work when using multi 2 2 2 8?

Please find in attachment my INFILE.


Thank you very much for taking the time to answer my question, I really appreciate it.

Best Regards,



Reply Oldest first
  • Oldest first
  • Newest first
  • Active threads
  • Popular
Like Follow
  • 7 mths agoLast active
  • 17Views
  • 1 Following