#cls
text(In-Form chemical-reaction source
title
libref=751
  
  DISPLAY
  
  This In-Form case does without Fortran coding what Z251 does with
  it, namely  it creates a source which is dependent on both the 
  dependent variable (RCTD) and the independent variable (X), which is a 
  requirement for the method of simulation of laminar flame propagation 
  which employs
      x = exp(-const * distance)      as the independent variable.
 
  It also illustrates the use of the y dimension as a parametric
  multiplier, thus allowing the multiplying constant of the reaction
  rate to have ten (ie NY) values in a single run.
 
  This device is often useful when parametric studies have to be
  made; for it also permits PHOTON to be used for display purposes.
 
  ENDDIS  
  PHOTON USE
  p;;;;
 
  con rctd z 1 fi;0.02
  msg contours of reactedness; x=distorted distance; y=multiplier-1
  pause;con rate z 1 fi;0.02
  msg contours of reaction rate; x=distorted dist.; y=multiplier-1
  ENDUSE

load(z251) 
text(In-Form equivalent of PLANT case Z251
libref=751

  INFORM13BEGIN
    ** The next line de-activates PLANT **
namsat=none
REACRATE=SKIP
NAME(C1)= RCTD ! reactdeness, = 0 upstream and 1 downstream
NAME(C2)= RATE ! reaction rate,  
  INFORM13END
  
     ** In-Form formulae which are equivalent to the 
        PLANT formulae of z251                      **
  INFORM13BEGIN
PATCH(IREACRAT,VOLUME,1,NX,1,NY,1,1,1,1)
(SOURCE of RCTD at IREACRAT is RCONST*RCTD^EXPO*XG^(-2)$
*(1.0+YG)*(1-RCTD) with LINE)
(SOURCE of RATE at IREACRAT is RCONST*RCTD^EXPO$
*(1.0+YG)*(1-RCTD) with FIXV)
  INFORM13END
DISTIL=T
EX(RCTD)=6.872E-01; EX(RATE)=2.961E-01
STORE(NPOR); EX(NPOR)=1.000E-11
 STOP