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

  Q1 uses the In-Form functions:
      SOURCE and STORED.

  The Q1 contains PHOTON USE commands
  ENDDIS  
  PHOTON USE
  p
 
 
 
  msg             Numerical solution  contours
  con C1 Z M fi;0.001
  con C1 Y M fi;0.001
  con C1 X M fi;0.001
  gr ou x m;gr ou y m;gr ou z m
  msg Press Enter to continue
  PAUSE
  CLEAR
  con EXAC Z M fi;0.001
  con EXAC Y M fi;0.001
  con EXAC X M fi;0.001
  gr ou x m;gr ou y m;gr ou z m
  msg               Exact solution  contours
  msg Press e to END
  ENDUSE

load(z118) 
text(In-Form equivalent of PLANT case Z118
libref=718

  INFORM13BEGIN
    ** The next line de-activates PLANT **
namsat=none
DIFWES=SKIP;DIFNOR=SKIP;DIFEAS=SKIP
DIFSOU=SKIP;DIFHIG=SKIP;DIFLOW=SKIP
  INFORM13END

     ** The following settings have been chosen as being
        preferable to those of z118, in respect of convergence
        and printout.                                         **            
LSWEEP=10
  RELAX(C1,LINRLX,1.0)
  
     ** In-Form property formulae which are equivalent to the
        PLANT formulae of z118                               **
  INFORM13BEGIN
nx=20;ny=20;nz=20
#unigrid  
PATCH(iDIFWES,WWALL,1,1,1,NY,1,NZ,1,LSTEP)
(SOURCE of C1 at iDIFWES is 2.*RHO1*ENUL*$
(-YG^2-ZG^2-C1)/DXG/:PRNDTL(C1): with line)
 
PATCH(iDIFNOR,NWALL,1,NX,NY,NY,1,NZ,1,LSTEP)
(SOURCE of C1 at iDIFNOR is 2.*RHO1*ENUL*$
(2.*XG^2-1.-ZG^2-C1)/DYG/:PRNDTL(C1): with line)
 
PATCH(iDIFEAS,EWALL,NX,NX,1,NY,1,NZ,1,LSTEP)
(SOURCE of C1 at iDIFEAS is 2.*RHO1*ENUL*$
(2.-YG^2-ZG^2-C1)/DXG/:PRNDTL(C1): with line)
 
PATCH(iDIFSOU,SWALL,1,NX,1,1,1,NZ,1,LSTEP)
(SOURCE of C1 at iDIFSOU is 2.*RHO1*ENUL*$
(2.*XG^2-ZG^2-C1)/DYG/:PRNDTL(c1): with line)
 
PATCH(iDIFHIG,HWALL,1,NX,1,NY,NZ,NZ,1,LSTEP)
(SOURCE of C1 at iDIFHIG is 2.*RHO1*ENUL*$
(2.*XG^2-YG^2-1.-C1)/DZG/:PRNDTL(c1): with line)
 
PATCH(iDIFLOW,LWALL,1,NX,1,NY,1,1,1,LSTEP)
(SOURCE of C1 at iDIFLOW is 2.*RHO1*ENUL*$
(2.*XG^2-YG^2-C1)/DZG/:PRNDTL(c1): with line)
  INFORM13END

  !!!! EXAC is the exact solution !!!!

  INFORM7BEGIN
(STORED of EXAC is 2.*XG^2 - YG^2 - ZG^2)
  INFORM7END
DISTIL=T
EX(C1)=5.964E-01; EX(EXAC)=5.906E-01
STOP