TALK=T;RUN( 1, 1)
 ** LOAD(334) from the PHOENICS Input Library
  PHOTON USE
  p; parphi; 1 1 500
 
  gr ou x 1
  gr off;red;gr ou y 1
  msg pressure; distance is vertical, time horizontal
  con p1 x 1 fi;1
  msg Press  to continue
  pause
  con off
  red
  msg velocity; distance is vertical, time horizontal
  con v1 x 1 fi;1
  msg trajectory of interface
  con c3 x 1 col 1; 1.1  1.11  2
  enduse
 
    GROUP 1. Run title and other preliminaries
TEXT(Spherical Blast Wave
TITLE
  DISPLAY
  High-pressure air is separated from low-pressure air by a
  spherical diaphragm concentric with a spherical outer vessel.
 
  The diaphragm suddenly breaks.
 
  The subsequent motion is predicted by a one-dimensional analysis.
 
  The spherical geometry is simulated by the use of VPOR and NPOR,
  which are made proportional to iy by the setting iporib=3.
 
  ENDDIS
 
 
STEADY=F; tlast=0.001; lstep=200; yvlast=1.0; ny=200
  TEXT(Uniform Grid
  #pause
mesg(unigrid macro has been loaded
grdpwr(t,lstep,tlast, 1.0)
grdpwr(x,nx   ,xulast,1.0)
grdpwr(y,ny   ,yvlast,1.0)
if(parab) then
 parab
 zfrac(1)= -nz;zfrac(2)=1/nz
else
 grdpwr(z,nz   ,zwlast,1.0)
endif
if(cartes) then
 mesg(uniform cartesian grid
else
 mesg(uniform polar grid
endif
  end of unigrid
 ** LOAD( 50) from the PHOENICS Input Library
 ** LOAD( 50) from the PHOENICS Input Library
IPORIB=3
 
    Group 7. Variables stored, solved & named
SOLVE(P1,V1,C3); STORE(RHO1,NPOR,VPOR)
 
    GROUP 8. Terms (in differential equations) & devices
TERMS(v1,Y,Y,N,Y,Y,Y)
TERMS(C3,Y,Y,N,Y,Y,Y)
    GROUP 9. Properties of the medium (or media)
  TEXT(Choose Fluid Materials
    71 start of ....fluidmat
store(prps)
    71 end of ....fluidmat
 ** LOAD( 71) from the PHOENICS Input Library
 ** LOAD( 71) from the PHOENICS Input Library
SETPRPS(1,AIRISENT)
REAL(REFPRE,REFRHO);REFPRE=1.E5;REFRHO=1.0
MESG(Isentropic air, Refrho=:REFRHO:, Refpre=:REFPRE:
RHO1A=REFRHO*(1./REFPRE)**RHO1B
 
    GROUP 11. Initialization of variable or porosity fields
FIINIT(P1)=REFPRE; FIINIT(RHO1)=REFRHO
FIINIT(C3)=1.0
PATCH(START,INIVAL,1,1,1,ny/20,1,1,1,1);INIADD=F
COVAL(START,P1,0.0,REFPRE*20.0)
COVAL(START,RHO1,0.0,REFRHO*20.0)
COVAL(START,C3,0.0,20.0/(20.0**RHO1B))
 
    GROUP 15. Termination of sweeps
LSWEEP=5
    GROUP 16. Relaxation
CONWIZ=T
 
    GROUP 22. Spot-value print-out
IXMON=Ny/2; TSTSWP=-1
SPEDAT(SET,GXMONI,TRANSIENT,L,F)  
    GROUP 23. Field print-out and plot control
NTPRIN=LSTEP/5;NXPRIN=NX/5
 
PATCH(ALLY,PROFIL,1,1,1,NY,1,1,1,LSTEP)
PLOT(ALLY,P1,0.0,0.0); PLOT(ALLY,V1,0.0,0.0)
 
PATCH(TIMEPLOT,PROFIL,1,1,1,1,1,1,1,LSTEP)
PLOT(TIMEPLOT,P1,0.0,0.0); PLOT(TIMEPLOT,V1,0.0,0.0)
 
  create parphi file for display via PHOTON
IDISPA=1;IDISPB=1;IDISPC=LSTEP
 LIBREF  =       334
STOP