PHOTON USE
  p
 
 
 
  msg       ******************************************************
  msg       *  Note:                                             *
  msg       *               ========                             *
  msg       *  Inspect the result file of d_priv1 to see the     *
  msg       *  numerical data.                                   *
  msg       *  The contours illustrate the solution.             *
  msg       ******************************************************
  msg Press Enter to continue
  pause;cl
  gr ou z 1
  con c1 z 1 fil;0.01
  msg                     C1  CONTOURS
  msg Press Enter to continue
  pause
  con off;red
  con plnt z 1 fil;0.01
  msg                     PLANT CONTOURS
  MSG Press E to end
  ENDUSE
     GROUP 1. Run title and other preliminaries
TEXT(1D uniform flow dispersion:121
 
  DISPLAY
    This run  shows  the  PLANT  equivalents  of  standard PHOENICS
    settings  to  introduce  the  sources   of   linearized   form.
    One-dimensional convection-diffusion problem can also be tested
    against exact solution.  
 
     ==============================================
   Governing equation:
 
     d(RHO*U*C)/dX = d(GAMMA*dC/dX) -SIGMA*C + SOURCE
 
   Boundary conditions:
 
     X =  0   , C = Cfix
     X -> -oo , dC/dX = 0
     X -> +oo , dC/dX = 0
      =============================================
  ENDDIS

REAL(SIGMA,SOURC,RATE, UIN)
    West wind of 1 m/s
SIGMA=.3;SOURC=0.1;RATE=1.;UIN=1.
    GROUP 3. X-direction grid specification
GRDPWR(X,11,11.0,1.0)
    GROUP 4. Y-direction grid specification
GRDPWR(Y,2,2.0,1.0)
    GROUP 7. Variables stored, solved & named
SOLVE(C1,PLNT);STORE(U1,V1)
    GROUP 8. Terms (in differential equations) & devices
TERMS(C1,N,Y,Y,Y,Y,N);TERMS(PLNT,N,Y,Y,Y,Y,N)
    GROUP 9. Properties of the medium (or media)
ENUL=1.
    GROUP 11. Initialization of variable or porosity fields
FIINIT(C1)=1.0;FIINIT(PLNT)=1.0
FIINIT(U1)=UIN;FIINIT(V1)=0.0
    GROUP 13. Boundary conditions and special sources
 
     * PHOENICS settings and PLANT equivalents
 
RG(1)=SIGMA;RG(2)=SOURC;RG(3)=RATE

  PLANTBEGIN 
PATCH(INLET,WEST,1,1,1,NY,1,NZ,1,LSTEP)
COVAL(INLET,C1,1.*1.,SAME)
    VAL=PLNT
COVAL(INLET,PLNT,1.*1.,GRND)
 
PATCH(ABSORB,VOLUME,1,NX,1,NY,1,NZ,1,LSTEP)
COVAL(ABSORB,C1,RG(1),0.0)
    CO=RG(1)
COVAL(ABSORB,PLNT,GRND,0.0)
 
PATCH(SORCE,VOLUME,1,NX,1,NY,1,NZ,1,LSTEP)
COVAL(SORCE,C1,FIXFLU,RG(2))
    VAL=RG(2)
COVAL(SORCE,PLNT,FIXFLU,GRND)
 
PATCH(SPIL,CELL,(NX-1)/2+1,(NX-1)/2+1,1,NY,1,NZ,1,LSTEP)
COVAL(SPIL,C1,FIXVAL,RG(3))
    VAL=RG(3)
COVAL(SPIL,PLNT,FIXVAL,GRND)
  PLANTEND
 
    GROUP 15. Termination of sweeps
LSWEEP=5
    GROUP 19. Data communicated by satellite to GROUND
     * call to special SATELLITE
NAMSAT=MOSG
    GROUP 20. Preliminary print-out
tstswp=-1
dmpstk=t
DISTIL=T
EX(U1)=1.000E+00; EX(C1)=6.088E-01; EX(PLNT)=6.088E-01
 LIBREF=121
STOP