GROUP 1. Run title and other preliminaries
TEXT(2D Channel Flow;Non Newtonian Fluid
TITLE
  DISPLAY
    This run illustrates the activation of a power-law
    laminar viscosity formula via the PROPS file
  ENDDIS
    GROUP 4. Y-direction grid specification
NREGY=2
IREGY=1;GRDPWR(Y,-5,0.5,1.5)
IREGY=2;GRDPWR(Y,-5,0.5,1.5)
    GROUP 5. Z-direction grid specification
NREGZ=3
IREGZ=1;GRDPWR(Z,5,1.0,-1.2)
IREGZ=2;GRDPWR(Z,-5,1.0,1.5)
IREGZ=3;GRDPWR(Z,5,2,2.0)
    GROUP 7. Variables stored, solved & named
SOLVE(P1,V1,W1)
#soltem
STORE(PRPS)
    GROUP 9. Properties of the medium (or media)
#use_props
GENK=T;STORE(ENUL,GEN1)
  
CSG10='q1'
    MATFLG=T;NMAT=1
    90 1000.0 GRND4 1000 .1428 0
    0.001 2
    GROUP 11. Initialization of variable or porosity fields
:fluid:=90
CONPOR(BLK,-1,CELL,1,1,-#2,#2,-#2,-#2);INIT(BLK,PRPS,0,101)
    GROUP 13. Boundary conditions and special sources
    Inlet
INLET(IN,LOW,1,1,1,NY,1,1,1,1)
VALUE(IN,P1,5.0*1000.0);VALUE(IN,W1,5.0)
VALUE(IN,V1,0.0);VALUE(IN,TEM1,9.0)
    Outlet
OUTLET(OUTLET,HIGH,1,1,1,NY,NZ,NZ,1,1)
    North wall
WALL (NORTH,NORTH,1,1,NY,NY,#1,#1,1,1)
COVAL(NORTH,TEM1,LOGLAW,0.0)
WALL (NORTH1,NORTH,1,1,NY,NY,#3,#3,1,1)
COVAL(NORTH1,TEM1,LOGLAW,0.0)
PATCH(SOLIDS,CELL,1,NX,1,NY,1,NZ,1,1)
    GROUP 15. Termination of sweeps
LSWEEP=50
RESREF(P1)=1.E-5;RESREF(V1)=1.E-5;RESREF(W1)=1.E-5
RESREF(TEM1)=1.E-5
    GROUP 17. Under-relaxation devices
RELAX(V1,FALSDT,ZWLAST/5/NZ/5);RELAX(W1,FALSDT,ZWLAST/5/NZ/5)
RELAX(P1,LINRLX,.7)
    GROUP 22. Spot-value print-out
IYMON=NY/2;IZMON=NZ-2
    GROUP 23. Field print-out and plot control