talk=t;run(1,1)
    PHOTON USE
  autoplot
  upause 1
  file; phi 5

  cl; da 1; mpro; da 1; moxi; da 1; fmix; da 1; mfue; da 1; oxix;
  da 1; fex; da 1; fuex; da 1; prox; scale;
  col7 1; col2 2; colf 3; colC 4; col1 5; col1 6; col1 7; col1 8; bi
  msg
  msg           Mass fractions of oxidant - blue, fuel - yellow,
  msg           products - green and mixture fraction - red.
  msg           Exact solution - white ......................... mor
  pause
  cl; da 1; temp; da 1; temx; col7 1; col1 2; big;
  msg
  msg           Temperature, K
  msg           Exact solution - white.
  pause
  enduse
     GROUP 1. Run title and other preliminaries

TEXT(Convection-only steady combustion of CO
     ---------------------------------

  DISPLAY

    Convection-only steady combustion of CO:

    One-phase, SCRS, 1D, one space, constant properties,
    convection only, comparison with analitical solution.

  ENDDIS

    GROUP 4. Y-direction grid specification
NY=1;GRDPWR(Y,NY,1.,1.0)
    GROUP 5. Z-direction grid specification
NZ=30;GRDPWR(Z,NZ,10.,1.0)
    GROUP 7. Variables stored, solved & named
SOLVE(P1,W1,FMIX,H1)
REAL(FSTOI,CP,STOIC,RIMF,HFU,HFUEL,TFUEL,TOX,HOX)
TFUEL = 289.0;TOX   =289.0
CP=  1100
HFU=1.e7
HOX  = CP*TOX;HFUEL= CP*TFUEL + HFU
STOIC=0.5174/0.232
FSTOI=1/(1+STOIC)
RIMF=1./(1-FSTOI)
STORE(TEMP,MPRO,MFUE,MOXI)
SOLUTN(P1,Y,Y,Y,N,N,N)
    GROUP 8. Terms (in differential equations) & devices
TERMS(FMIX,N,Y,Y,Y,Y,Y)
TERMS(H1,N,Y,Y,N,Y,N)
    GROUP 9. Properties of the medium (or media)
REAL(RHOIN1,RHOIN2,WAIR,WFU,WPR)
PRESS0=1.e5
WAIR=29.;WFU=28.;WPR=29.
RHOIN1=PRESS0*WFU/(8314.*TFUEL)
RHOIN2=PRESS0*WAIR/(8314.*TOX)
RHO1=RHOIN1
    GROUP 13. Boundary conditions and special sources
  ** Inlet Boundaries
INLET(INLET,LOW,1,NX,1,NY,1,1,1,1)
VALUE(INLET,P1 , RHOIN1*1.)
VALUE(INLET,W1 , 1.)
VALUE(INLET,FMIX , 1.0)
VALUE(INLET,H1,HFUEL)
  ** Volumetric mass source
PATCH(MASS,VOLUME,1,NX,1,NY,1,NZ,1,1)
COVAL(MASS,P1,FIXFLU,.4)
COVAL(MASS,FMIX,ONLYMS,0.0)
COVAL(MASS,H1,ONLYMS,HOX)
  **Outlet boundary
PATCH(OUTLET,HIGH,1,NX,1,NY,NZ,NZ,1,1)
COVAL(OUTLET,P1,fixp,0.0)
    GROUP 15. Termination of sweeps
LSWEEP=200
RESFAC=0.0001
    GROUP 16. Termination of iterations
LITHYD=10
    GROUP 22. Monitor print-out
IZMON=NZ-1;IYMON=NY-1;UWATCH=T
    GROUP 23. Field print-out and plot control
NPLT=1;NYPRIN=1;NZPRIN=1
NYPRIN=1;IYPRF=1;IYPRL=30
TSTSWP=-1

namsat=mosg

   FMIX=AMAX1(0.,AMIN1(1.,FMIX))
   MFUE=AMAX1(0.,:RIMF:*FMIX-:RIMF:*:FSTOI:)
   MOXI=AMAX1(0.,1.-FMIX/:FSTOI:)
   MPRO=1.-MOXI-MFUE
   TEMP=(H1-:HFU:*MFUE)/:CP:
STORE(FEX,FUEX,OXIX,PROX,H1X,TEMX)
  FEX=:RHOIN1:*1/(:RHOIN1:*1+.4*ZGNZ)
  OXIX=AMAX1(0.,1.-FEX/:FSTOI:)
  FUEX=AMAX1(0.,:RIMF:*FEX-:RIMF:*:FSTOI:)
  PROX=1.-OXIX-FUEX
  H1X=:HFUEL:*FEX+:HOX:*(1.-FEX)
  TEMX=(H1X-:HFU:*FUEX)/:CP:
STOP