```
** Macro used in some 2-phase combustion cases
GROUP 7. Variables stored, solved & named
SOLUTN(11,Y,Y,Y,P,P,P)
** provide storage for inter-phase mass transfer
STORE(MDOT,CFIP)
** Solve additionally for the mixture fraction, i.e. the quantity
of phase-2 material which has entered phase 1.
SOLUTN(C1,Y,Y,Y,P,P,P); NAME(C1)=MIXF
STORE(YCO,YO2,YCO2,YN2,YH2,YH2O,RHO1,RHO2)
SOLUTN(H1,Y,Y,Y,P,P,P); SOLUTN(H2,Y,Y,Y,P,P,P)
STORE(TMP1,TMP2)

GROUP 8. Terms (in differential equations) & devices
TERMS(H1,N,P,P,P,P,P); TERMS(H2,N,P,P,P,P,P)

GROUP 9. Properties of the medium (or media)
RHO1=7GASES;  RHO2=1.E3 ; PRESS0=1.E5 ; TEMP0=0.0
RHO1A=CINCL ; RHO1B=HINCL

GROUP 10. Inter-phase-transfer processes and properties
** Set constant interphase friction factor and activate
the calculation of the interphase mass transfer by:
CFIPS=GRND1; CFIPC=1.E5
CMDOT=GRND3; CMDTA=BURNRATE; CMDTC=FS

** Note that grnd3 is the mdot option which makes the mass-
transfer rate proportional to (cmdtc-mixf), where cmdtc
stands for the saturation value of mixf, i.e. the largest
value which can be attained as a result of mass transfer.

CINT(MIXF)=0.0; PHINT(MIXF)=1.0

GROUP 11. Initialization of variable or porosity fields
RSIN=FLOS/(RHO2*VELO);  RHOGIN=PRESS0/(287.41*TGIN)
RGIN=1.-RSIN;           VELG=FLOG/(RHOGIN*RGIN)