PHOTON USE
  p;p5;5 1
 
  view x;up z; norm
  msg horizontal dimension has been multiplied by 5
  msg Press  after each pause
  gr x 1;pause;  gr off;red;gr ou x 1; set vec comp; u1 av1 cw1
  msg Gas velocity vectors at time = 5 second
  vec x 1 sh; pause; vec off;red
  msg Liquid velocity vectors
  set vec comp; u2 bv2 dw2; vec x 1 sh; pause; vec off;red
  msg pressure contours
  con p1 x 1 fi;0.01; pause; con off;red
  msg gas volume-fraction contours
  con gas x 1 fi;0.01; pause
 
  p;p10;5 1
 
  view x;up z; norm
  msg horizontal dimension has been multiplied by 5
  msg Press  after each pause
  gr x 1;pause;  gr off;red;gr ou x 1; set vec comp; u1 av1 cw1
  msg Gas velocity vectors at time = 10 second
  vec x 1 sh; pause; vec off;red
  msg Liquid velocity vectors
  set vec comp; u2 bv2 dw2; vec x 1 sh; pause; vec off;red
  msg pressure contours
  con p1 x 1 fi;0.01; pause; con off;red
  msg gas volume-fraction contours
  con gas x 1 fi;0.01; pause
 
  p;p15;5 1
 
  view x;up z; norm
  msg horizontal dimension has been multiplied by 5
  msg Press  after each pause
  gr x 1;pause;  gr off;red;gr ou x 1; set vec comp; u1 av1 cw1
  msg Gas velocity vectors at time = 15 second
  vec x 1 sh; pause; vec off;red
  msg Liquid velocity vectors
  set vec comp; u2 bv2 dw2; vec x 1 sh; pause; vec off;red
  msg pressure contours
  con p1 x 1 fi;0.01; pause; con off;red
  msg gas volume-fraction contours
  con gas x 1 fi;0.01; pause
 
  p;p20;5 1
 
  view x;up z; norm
  msg horizontal dimension has been multiplied by 5
  msg Press  after each pause
  gr x 1;pause;  gr off;red;gr ou x 1; set vec comp; u1 av1 cw1
  msg Gas velocity vectors at time = 20 second
  vec x 1 sh; pause; vec off;red
  msg Liquid velocity vectors
  set vec comp; u2 bv2 dw2; vec x 1 sh; pause; vec off;red
  msg pressure contours
  con p1 x 1 fi;0.01; pause; con off;red
  msg gas volume-fraction contours
  con gas x 1 fi;0.01; pause;con off; red
  msg contours of vertical slip velocity
  con slpw x 1 fi;0.01; pause;con off; red
  msg contours of horizontal slip velocity
  con slpv x 1 fi;0.01; pause;con off; red
 
  msg Press e to END
  enduse
 
    GROUP 1. Run title
TEXT(transient air-lift pump            :W430
TITLE
    GROUPs 2, 3, 4 & 5
STEADY=F; lstep=20; tlast=20.0
CARTES=F; NX=1; XULAST=0.1; NY=15; YVLAST=0.25; NZ=30; ZWLAST=2.0
l(unigrid
 
    GROUP 7. Variables stored, solved & named
ONEPHS=F;SOLVE(P1,V1,V2,W1,W2,R1,R2);SOLUTN(P1,Y,Y,Y,P,P,P)
NAME(V1)=AV1;NAME(V2)=BV2;NAME(W1)=CW1;NAME(W2)=DW2
NAME(R1)=GAS;NAME(R2)=SOL
  ** provide for the calculation and storage of slip velocities
STORE(SLPV,SLPW)
    GROUP 8. Terms (in differential equations) & devices
TERMS(AV1,Y,Y,N,Y,P,P);TERMS(BV2,Y,Y,N,Y,P,P)
TERMS(CW1,Y,Y,N,Y,P,P);TERMS(DW2,Y,Y,N,Y,P,P)
    GROUP 9. Properties of the medium (or media)
ENUL=0.0;RHO2=2500.0
     Compressible gas with density =1.0 at p = PRESS0
RHO1=COMPRESS;RHO1B=1.0/1.4;RHO1C=0.0;PRESS0=1.E5
RHO1A=PRESS0**(-RHO1B);DRH1DP=COMPRESS
    GROUP 10. Inter-phase-transfer processes and properties
  ** Set a constant inter-phase friction coefficient
cfips=5.e5
    GROUP 11. Initialization of variable or porosity fields
FIINIT(GAS)=0.1;FIINIT(SOL)=0.9
INIADD=T
   CONPOR(0.0,CELL,1,1,NY/2+1,NY,NZ/2+1,NZ)
patch(top,inival,1,1,1,ny,nz-2,nz-2,1,1)
init(top,gas,0.9,0.0);init(top,sol,-0.9,0.0)
    GROUP 13. Boundary conditions and special sources
  ** outlet1
PATCH(OUTLET1,HIGH,1,1,1,NY/2,NZ,NZ,1,LSTEP)
COVAL(OUTLET1,P1,FIXP,0.0);COVAL(OUTLET1,P2,FIXP*RHO2,0.0)
  ** inlet1
PATCH(INLET1,LOW,1,1,1,NY/5,NZ/4,NZ/4,1,LSTEP)
COVAL(INLET1,P1,FIXFLU,50.0);COVAL(INLET1,CW1,0.0,-50.0)
COVAL(INLET1,P1,FIXFLU,5.0);COVAL(INLET1,CW1,0.0,5.0)
  ** Gravity
PATCH(GRAVITY,PHASEM,1,NX,1,NY,1,NZ,1,LSTEP)
COVAL(GRAVITY,CW1,FIXFLU,-9.81);COVAL(GRAVITY,DW2,FIXFLU,-9.81)
  ** dividing wall
PATCH(PANEL,NORTH,1,1,NY/2,NY/2,2,3*NZ/4,1,LSTEP)
COVAL(PANEL,AV1,FIXVAL,0.0);COVAL(PANEL,BV2,FIXVAL,0.0)
    GROUP 15. Termination of sweeps
lsweep=20
    GROUP 17. Under-relaxation devices
RELAX(GAS,LINRLX,0.1);RELAX(GAS,LINRLX,0.1);RELAX(AV1,FALSDT,0.01)
RELAX(CW1,FALSDT,0.01);RELAX(BV2,FALSDT,0.01);RELAX(DW2,FALSDT,0.01)
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
    GROUP 18. Limits on variables or increments to them
VARMIN(GAS)=1.E-5;VARMIN(SOL)=1.E-5
    GROUP 22. Spot-value print-out
TSTSWP=-1;IXMON=NX/2;IYMON=NY/2;NPRMNT=1
    GROUP 23. Field print-out and plot control
IPLTL=LSWEEP;NTPRIN=LSTEP/8;NXPRIN=NX/20;NYPRIN=NY/20
NTPRIN=LSTEP/5;NZPRIN=NZ/5;NZPRIN=NZ/5
PATCH(4,PROFIL,NX/2,NX/2,1,NY,1,1,1,LSTEP)
PLOT(4,GAS,0.0,0.);PATCH(5,PROFIL,NX/2,NX/2,NY/2,NY/2,1,1,1,LSTEP)
PLOT(5,CW1,0.,0.);PLOT(5,DW2,0.,0.)
PATCH(6,CONTUR,1,NX,1,NY,1,1,1,LSTEP)
PLOT(6,P1,0.0,0.);PLOT(6,GAS,0.0,0.)
    GROUP 24. Dumps for restarts
  ** provision for dumping grid and field files for photon
CSG1=PHI;CSG2=XYZ;IDISPA=5