macro concerned with solid-fuel combustion
  called by cases 480, 481 and 482
  solid phase
COVAL(INLETS,CVEL,ONLYMS,VELO); COVAL(INLETS,H2,ONLYMS,HSIN)
COVAL(INLETS,P2,FIXFLU,FLOS)  ; COVAL(INLETG,P1,FIXFLU,FLOG)
  gas phase
COVAL(INLETG,GVEL,ONLYMS,VELG); COVAL(INLETG,MIXF,ONLYMS,0.)
COVAL(INLETG,H1,ONLYMS,HGIN)
  ** multiply OUTCO1 by estimated value of exit density
COVAL(OUTLET,P1,OUTCO1*0.1,0.0); COVAL(OUTLET,P2,OUTCO1*RHO2,0.0)
    GROUP 15. Termination of sweeps
LSWEEP=2000; SELREF=T; RESFAC=0.001
 
    GROUP 17. Under-relaxation devices
CHATIM=MULT*LENGTH/(VELG*NX*NY*NZ)
 
RELAX(P1,LINRLX,0.7);   RELAX(SHAD,LINRLX,0.6)
RELAX(FUEL,LINRLX,0.6); RELAX(GAS,LINRLX,0.6)
 
RELAX(GVEL,FALSDT,CHATIM); RELAX(CVEL,FALSDT,CHATIM)
 
RELAX(MIXF,FALSDT,CHATIM*100/burnrate)
RELAX(H1,FALSDT,CHATIM); RELAX(H2,FALSDT,CHATIM)
 
RELAX(RHO1,LINRLX,0.5); RELAX(MDOT,LINRLX,0.3)
 
    GROUP 18. Limits on variables or increments to them
VARMAX(MIXF)=CMDTC; VARMIN(MIXF)=0.
VARMIN(FUEL)=1.E-9
 
    GROUP 21. Print-out of variables
OUTPUT(SHAD,Y,Y,Y,Y,Y,Y); OUTPUT(GAS ,Y,Y,Y,Y,Y,Y)
OUTPUT(MIXF,Y,Y,Y,Y,Y,Y); OUTPUT(FUEL,Y,Y,Y,Y,Y,Y)
OUTPUT(H1  ,Y,Y,Y,Y,Y,Y); OUTPUT(H2  ,Y,Y,Y,Y,Y,Y)
    GROUP 22. Spot-value print-out
IYMON=NY-2
TSTSWP=-1
    GROUP 23. Field print-out and plot control
ITABL=1; ORSIZ=0.2; Nplt=1
 
PATCH(PROFIL1,PROFIL,1,NX,1,NY,1,NZ,1,1)
COVAL(PROFIL1,YCO,0,0)
COVAL(PROFIL1,YCO2,0,0)
COVAL(PROFIL1,MIXF,0,0)
COVAL(PROFIL1,FUEL,0,0)
 
PATCH(PROFIL2,PROFIL,1,NX,1,1,NY,NZ,1,1)
COVAL(PROFIL2,P1,0,0)
COVAL(PROFIL2,GVEL,0,0)
COVAL(PROFIL2,5,0,0)
COVAL(PROFIL2,RHO1,0,0)
 
PATCH(PROFIL3,PROFIL,1,NX,1,NY,1,NZ,1,1)
COVAL(PROFIL3,TMP1,0,0)
COVAL(PROFIL3,TMP2,0,0)