#cls
text(In-Form Equivalent Of PLANT Case Z244
title
libref=750
  
  DISPLAY
  
  This In-Form case does what Z244 does, but without fortran coding
  or re-compilation.
  
  This case use In-Form to create a flow resistance proportional
  to W1**1.8 in the central section of the pipe.
 
  First, case z244 is loaded. Then NAMSAT = NONE is set, to 
  deactivate.

  The Q1 uses SOURCE In-Form statement.

  The Q1 contains PHOTON USE commands
  ENDDIS  
  PHOTON USE
  p
  phi
  4 1
 
  view -x
  msg picture is enlarged 4 times in y direction
  gr ou x 1
  msg contours of concentration of material with Pr = 10.0
  con c x 1 fi;0.001;pause;con off; red
  msg contours of concentration of material with Pr = 0.1
  con b x 1 fi;0.001;pause;
  msg velocity vectors
  vec x 1
  enduse

load(z244)
text(In-Form equivalent of PLANT case Z244
libref=750

  INFORM13BEGIN
    ** The next line de-activates PLANT **
namsat=none
RESIST=SKIP
  INFORM13END
  
     ** The following settings have been chosen as being
        preferable to those of z244, in respect of convergence
        and printout.                                         **
ISWC1=LSWEEP-5
SELREF=F
RELAX(P1,LINRLX,0.5)
RELAX(V1,FALSDT,.001)
RELAX(W1,FALSDT,.001)

  INFORM13BEGIN
     ** In-Form formulae which are equivalent to the
        PLANT formulae of z244                      **

PATCH(iRESIST,VOLUME,1,1,1,NY,NZ/4+1,3*NZ/4,1,1)
(SOURCE of W1 at iRESIST is -RESCO*W1^(1+POWER) with LINE)
  INFORM13END
DISTIL=T
EX(P1)=7.609E-01; EX(W1)=1.057E-01; EX(H1)=5.241E-01
EX(A)=4.302E-01; EX(B)=9.120E-01; EX(C)=1.269E-01
EX(V1)=8.005E-04
STOP