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.