GROUP 1. Run title and other preliminaries
TEXT(1D Venturi, Y-Direction           
TITLE
  This run performs the same analysis as above using the
  y-directed coordinate
 
    GROUP 4. Y-direction grid specification
GRDPWR(Y,5,1.0,1.0)
 
    GROUP 7. Variables stored, solved & named
SOLVE(P1,V1)
  ** Provide storage for porosities on the north faces.
STORE(NPOR)
 
    GROUP 11. Initialization of variable or porosity fields
FIINIT(V1)=1.0
  ** Set varying-cross section by means of non-uniform
     porosities (varying linearly with IY).
INIADD=T     
PATCH(SECT1,LINVLY,1,1,1,NY,1,1,1,1);INIT(SECT1,NPOR,-0.1,1.0)
  ** Modify porosity field from IY=2 to IY=5
PATCH(SECT2,LINVLY,1,1,2,NY,1,1,1,1);INIT(SECT2,NPOR,0.2,0.0)
 
    GROUP 13. Boundary conditions and special sources
  ** Inlet
INLET(IN,CELL,1,1,1,1,1,1,1,1)
VALUE(IN,P1,1.0);VALUE(IN,V1,0.9804)
  ** Outlet
PATCH(OUTLET,CELL,1,1,NY,NY,1,1,1,1);COVAL(OUTLET,P1,FIXP,0.0)
COVAL(OUTLET,V1,ONLYMS,0.0)
 
    GROUP 15. Termination of sweeps
LSWEEP=20
 
    GROUP 17. Under-relaxation devices
RELAX(V1,FALSDT,1.0)
 
    GROUP 23. Field print-out and plot control
NPRINT=20;ITABL=1
PATCH(LONGPLOT,PROFIL,1,1,1,NY,1,1,1,1)
PLOT(LONGPLOT,P1,0.0,0.0);PLOT(LONGPLOT,V1,0.0,0.0)