ECO2N, initial condition

Here's a quick comment and a few tips:

(1) The system clearly has not reached steady state after 1000 time steps (I can tell from ST, DT and DX1) - I would therefore be impossible for me to guess what pressure and temperature profile I would expect at this random time (i.e., I don't think anything is wrong - the inconsistency is just that your expectations are for steady-state conditions, whereas you stopped the simulation after a short time, T and even P cannot equilibrate over a distance of 2 km after a single day of simulated time!). 

(2) Your convergence criterion RE1=1E-9 is unnecessarily tight for a steady-state run; 1E-5 would suffice; this change should allow you to get to steady state much faster. (By the way: your SOLVR parameters are a little inconsistent with RE1=1E-9).

(3) Increase MCYC (I would try 2500, which should be enough to reach your specified endtime of 100 years; you still need to check whether you have actually reached steady state, though!).

(4) It would also be useful to actually observe how the system evolves from your initial conditions to steady state; so please set MCYPR to a reasonable value (e.g., 500) and/or add a TIMES block, so you get some output you can analyze and a feeling for the system behavior. This is also helpful if you submit questions to the forum, so your output file contains some useful information other than the raw convergence numbers.

I did not run your model, but these simple changes will hopefully help.

Stefan

Sawyer,

NaN means "not-a-number" and it indicates that something went really wrong (divide by zero; square root of a negative number, overflow/underflow, etc.). I would need to see the input file (but am not sure I have the time to debug this right away - anybody from the TOUGH community willing to have a look at this? Thanks!)

Stefan