** Macro used in some 2-phase combustion cases
    GROUP 7. Variables stored, solved & named
SOLUTN(11,Y,Y,Y,P,P,P)
NAME(9)=GAS; NAME(R2)=FUEL; NAME(11)=SHAD
  ** 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)
FIINIT(GAS)=RGIN;       FIINIT(FUEL)=RSIN; FIINIT(SHAD)=RSIN
FIINIT(GVEL)=VELG;      FIINIT(CVEL)=VELO; FIINIT(MDOT)=0.01*FLOS
FIINIT(RHO1)=RHOGIN;    FIINIT(C4)=0.0;    FIINIT(C6)=0.0
FIINIT(H1)=HGIN;        FIINIT(H2)=HSIN;   FIINIT(MIXF)=0.
    GROUP 13. Boundary conditions and special sources
REAL(OUTCO1); OUTCO1=1.E3
REAL(MULT);MULT=2.0