Andrea, Without trying to answer your questions substantially, I would just like to point out that you have, in your flash.par, sim_pAtm = 0. sim_pTarg = 0. so it seems you are initializing various pressures to 0. That shouldn't really work, and I am a bit surprised your simulation got even as far as it did. Klaus