Post a reply || Back to TOUGH2.com - TOUGH2 and TOUGH2-MP Forum

ECO2N

Submitted by Keni on 04/Mar/2009
128.3.13.140

Message:

In Compu_Eos.f,

there is a code fragment as follows:
c For small CO2 mass fractions, or for forced-isothermal calculation
c (NEQ=NK), jump over this section (1905051).
if(k.eq.1) x30=x3
if(x3.gt.1.d-6.and.(neq.eq.nk1.or.nb.ge.8)) then

In the last line of this code fragment, we should have "x30" instead
of "x3"; so the line should be changed to
if(x30.gt.1.d-6.and.(neq.eq.nk1.or.nb.ge.8)) then

Making this change will improve code performance for problems that are non-isothermal and/or include aqueous diffusion of
dissolved CO2.

Replies to this post

No replies yet