L(327)
    GROUP 1. Run title and other preliminaries
TEXT(1D Piston-With-Bowl Study         
TITLE
  In this case, 10 cells are used,and half of them have a volume
  porosity of 0.5, to represent the bowl in the piston.
 
    GROUP 5. Z-direction grid specification
GRDPWR(Z,10,0.2,1.0)
 
    GROUP 7. Variables stored, solved & named
SOLVE(P1,W1);SOLUTN(P1,Y,Y,Y,N,N,N);STORE(VPOR,HPOR)
 
    GROUP 11. Initialization of variable or porosity fields
   *** Initialize the porosities
FIINIT(VPOR)=1.0; FIINIT(HPOR)=1.0
 
   *** Set the porosities for the bowl region
INIADD=T
PATCH(BOWL1,INIVAL,1,1,1,1,NZ/2+1,NZ,1,1)
INIT (BOWL1,VPOR,0.0,-0.5)
PATCH(BOWL2,INIVAL,1,1,1,1,NZ/2,NZ,1,1)
INIT (BOWL2,HPOR,0.0,-0.5)
 
    GROUP 15. Termination of sweeps
LSWEEP=20
   *** The false-time-step mode of under-relaxation is used
       for W1
RELAX(W1,FALSDT,0.01)
 
    GROUP 19. Data communicated by satellite to GROUND
   *** The iz location of the last slab in part one of an n-part
       moving grid is changed from 1 to 5.
IZW1=5
 
    GROUP 22. Spot-value print-out
NTPRIN=5
TSTSWP=LSWEEP/4
 
    GROUP 23. Field print-out and plot control
PATCH(TIMEPLOT,PROFIL,1,1,1,1,NZ,NZ,1,LSTEP)
PLOT (TIMEPLOT,P1,0.0,0.0)
PATCH(TIM10,PROFIL,1,1,1,1,1,NZ,LSTEP/2,LSTEP/2)
PLOT (TIM10,P1,0.0,0.0)
PATCH(TIM20,PROFIL,1,1,1,1,1,NZ,LSTEP,LSTEP)
PLOT (TIM20,P1,0.0,0.0)
 
    GROUP 24. Dumps for restarts