PHOTON USE
  p;parphi;1 1 1000
 
  do kk=1,100,2
  msg pressure contours
  con p1 z kk fi;1
  enddo
  pause;con off;red
 
  do kk=1,100,2
  msg marker concentrations
  con c3 z kk fi;1
  enddo
  pause;con off;red;view z
  enduse
 
l($332
    GROUP 1. Run title and other preliminaries
TEXT(Shock Tube With 2D Obstacle 
TITLE
  DISPLAY
  This case is similar to case 332; but only half of the diaphragm
  breaks, with the result that a two-dimensional flow ensues.
 
  Case 332 is loaded first, so that the Q1 contains only those
  settings which the two-dimensionality necessitates.
 
  ENDDIS
 
    GROUP 2. Y-direction grid specification
GRDPWR(Y,20,1.0,1.0)
 
    GROUP 7. Variables stored, solved & named
SOLVE(V1)
 
    GROUP 8. Terms (in differential equations) & devices
  ** De-activate the diffusion term in the V1 equation
TERMS(V1,Y,Y,N,Y,Y,Y)
 
 
    GROUP 11. Initialization of variable or porosity fields
 
FIINIT(V1)=0.0
PATCH(START,INIVAL,1,NX/2,1,NY,1,1,1,1);INIADD=F
 
    GROUP 13. Boundary conditions and special sources
PATCH(OBSTACLE,CELL,NX/2,NX/2,1,NY/2,1,1,1,LSTEP)
COVAL(OBSTACLE,U1,FIXVAL,0.0)
 
LSWEEP=5;IXMON=NX/2;IYMON=NY/2+1
 
IDISPA=5;IDISPB=1;IDISPC=LSTEP