#cls
text(In-Form Equivalent Of PLANT Case Z112
title
libref=712
  
  DISPLAY
  
  This In-Form case does what Z112 does, but without fortran coding
  or re-compilation.
  
  First, case z112 is loaded. Then NAMSAT = NONE is set, to 
  deactivate.

  The Q1 uses the In-Form statements:
      STORED and SOURCE.
    
  The Q1 contains PHOTON USE commands
  ENDDIS  
  PHOTON USE
  p
 
 
 
  4. 1. 1.
 
 
  gr ou z 1
  vec z 1 sh
  msg                    Velocity distribution:
  msg
  msg Press Enter to continue
  pause;cl
  msg                  Concentration distribution :
  msg                     Numerical solution
  con c1 z 1 fil;0.01
  msg Press Enter to continue
  pause;con off;red
  msg                  Concentration distribution :
  msg                     Exact solution
  con exac z 1 fil;0.01
  msg Press e to END
  ENDUSE

load(z112) 
text(In-Form equivalent of PLANT case Z112
libref=712

  INFORM13BEGIN
    ** The next line de-activates PLANT **
namsat=none
SORC1=SKIP;SORC2=SKIP;SORC3=SKIP;EXACT=SKIP
  
     ** In-Form formulae which are equivalent to the
        PLANT formulae of z112                      **
PATCH(ISORC1,VOLUME,1,NX,1,NY,1,NZ,1,1)
(SOURCE of C1 at ISORC1 is 1.-2./PE+2.*XG)
 
PATCH(ISORC2,VOLUME,1,NX,1,NY,1,NZ,1,1)
(SOURCE of C1 at ISORC2 is PI*COS(PI*XG))
 
PATCH(ISORC3,VOLUME,1,NX,1,NY,1,NZ,1,1)
(SOURCE of C1 at ISORC3 is PI^2*SIN(PI*XG)/PE)
 
PATCH(IEXACT,CELL,1,NX,1,NY,1,NZ,1,1)
(SOURCE of AUX1 at IEXACT is 1.+XG+(XG)^2 with FIXV)
(SOURCE of AUX2 at IEXACT is EXP(-PE*(1.-XG)) with FIXV)
(SOURCE of AUX3 at IEXACT is SIN(PI*XG) with FIXV)
(SOURCE of EXCT at IEXACT is AUX1+AUX2+AUX3 with FIXV)
  INFORM13END
  INFORM7BEGIN
(STORED of EXAC is 1.+XG+(XG)^2+EXP(-PE*(1.-XG))+SIN(PI*XG))
  INFORM7END
DISTIL=T
EX(U1)=1.000E+00; EX(C1)=3.107E+00
EX(EXAC)=3.104E+00; EX(EXCT)=3.104E+00
STOP