AUTOPLOT USE
  file
  phi 5
 
  da 1
  1v
  screen
  msg 1st-phase velocity at t=10s
  pl 1
  msg press  to continue
  pause
  cl
  da 1
  2v
  screen
  msg 2nd-phase velocity at t=10s
  pl 1
  msg press  to continue
  pause
  cl
  da 1
  r1
  screen
  msg 1st-phase volume fraction at t=10s
  pl 1
  msg press  to continue
  pause
  cl
  da 1
  r2
  screen
  msg 2nd-phase volume fraction at t=10s
  pl 1
  msg press  to continue
  pause
  cl
  da 1
  p1
  screen
  msg pressure at t=10s
  pl 1
  msg press e to END
  enduse
 
    GROUP 1. Run title and other preliminaries
TEXT(AS FOR RUN 771 BUT IN Y DIRECTION : W772
TITLE
    GROUP 2. Transience; time-step specification
STEADY=F;GRDPWR(T,10,10.0,1.0)
    GROUP 5. Z-direction grid specification
GRDPWR(Y,10,1.0,1.0)
    GROUP 7. Variables stored, solved & named
ONEPHS=F
SOLVE(P1,V1,V2,R1,R2);NAME(V1)=1V;NAME(V2)=2V
    GROUP 9. Properties of the medium (or media)
RHO2=1000.0
    GROUP 10. Inter-phase-transfer processes and properties
CFIPS=1.E4
    GROUP 11. Initialization of variable or porosity fields
FIINIT(R1)=0.1;FIINIT(R2)=0.9;FIINIT(P1)=1.0E2
    GROUP 13. Boundary conditions and special sources
PATCH(GRAVITY,PHASEM,1,1,1,NY,1,1,1,LSTEP)
COVAL(GRAVITY,1V,FIXFLU,-9.81)
COVAL(GRAVITY,2V,FIXFLU,-9.81)
PATCH(TOP,CELL,1,1,NY,NY,1,1,1,LSTEP)
COVAL(TOP,P1,1.0,0.0);COVAL(TOP,P2,RHO2,0.0)
    GROUP 15. Termination of sweeps
LSWEEP=30
RESREF(P1)=1.E-6;RESREF(R1)=1.E-6;RESREF(R2)=1.E-6
RESREF(1V)=1.E-6;RESREF(2V)=1.E-6
    GROUP 17. Under-relaxation devices
RELAX(R2,LINRLX,0.5)
RELAX(1V,FALSDT,0.1);RELAX(2V,FALSDT,0.1)
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
    GROUP 21. Print-out of variables
OUTPUT(R1,N,N,N,N,N,N)
OUTPUT(P1,Y,Y,Y,Y,Y,Y);OUTPUT(R2,Y,Y,Y,Y,Y,Y)
OUTPUT(1V,Y,Y,Y,Y,Y,Y);OUTPUT(2V,Y,Y,Y,Y,Y,Y)
    GROUP 22. Spot-value print-out
TSTSWP=-1;IYMON=2;NTPRIN=5
    GROUP 23. Field print-out and plot control
IPLTF=2;IPLTL=LSWEEP;YZPR=T;IPROF=3;ITABL=3
PATCH(PROFILE,PROFIL,1,1,1,NY,1,1,1,LSTEP)
PLOT(PROFILE,P1,0.0,0.0);PLOT(PROFILE,R1,0.0,0.0)
PLOT(PROFILE,1V,0.,0.);PLOT(PROFILE,2V,0.,0.)
PATCH(TIMEPLOT,PROFIL,1,1,NY,NY,1,1,1,LSTEP)
PLOT(TIMEPLOT,R1,0.0,0.0)
    GROUP 24. Dumps for restarts