DELETE TYPES;
READ FILE "flash2.a4c";
COMPILE t OF do_a_flash;
BROWSE {t};
RUN {t.seqmod_P_Calculate_T};
SOLVE {t} WITH QRSlv;
BROWSE {t};

STUDY {Flash.T} {Flash.x['benzene']} {Flash.y['benzene']} {gammabenzene} \
{gammatoluene} \
IN t \
VARYING {Feed.y['benzene'] {0.01} {0.1} {0.2} \
{0.3} {0.4} {0.5} {0.6} {0.7} {0.8} {0.9} {0.99}} \
USING QRSlv \
OUTFILE /tmp/flash2PT.dat \
ERROR STOP; 

DELETE TYPES;
READ FILE "flash2.a4c";
COMPILE t OF do_a_flash;
BROWSE {t};
RUN {t.seqmod_T_Calculate_P};
SOLVE {t} WITH QRSlv;
BROWSE {t};

STUDY {Flash.P} {Flash.x['benzene']} {Flash.y['benzene']} {gammabenzene} \
{gammatoluene} \
IN t \
VARYING {Feed.y['benzene'] {0.01} {0.1} {0.2} \
{0.3} {0.4} {0.5} {0.6} {0.7} {0.8} {0.9} {0.99}} \
USING QRSlv \
OUTFILE /tmp/flash2TP.dat \
ERROR STOP; 

DELETE TYPES;
READ FILE "flash2.a4c";
COMPILE t OF do_a_flash;
BROWSE {t};
RUN {t.seqmod_PT_Flash};
SOLVE {t} WITH QRSlv;
BROWSE {t};

STUDY {Flash.state.phase_fraction['vapor']} \
{Flash.x['benzene']} {Flash.y['benzene']} {gammabenzene} \
{gammatoluene} \
IN t \
VARYING {Feed.y['benzene'] {0.50} {0.55} {0.60}} \
USING QRSlv \
OUTFILE /tmp/flash2f.dat \
ERROR STOP;