TALK=T;RUN(1,1)
  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
  PHOTON USE
  p; parphi;  1 1 500

  gr ou y 1
  gr off;red;gr ou y 1
  msg pressure; distance is vertical, time horizontal
  con p1 y 1 fi;1
  msg Press  to continue
  pause
  con off
  red
  msg velocity; distance is vertical, time horizontal
  con u1 y 1 fi;1
  msg trajectory of interface
  con c3 y 1 col 1; 1.1  1.11  2
  enduse
  PHOTON USE
  p; parphi;  1 1 500

  gr ou y 1
  msg pressure; distance is vertical, time horizontal
  con p1 y 1 fi;1
  msg Press  to continue
  pause;con off;red
  msg velocity; distance is vertical, time horizontal
  con u1 y 1 fi;1
  enduse
  DISPLAY
  Air is moving steadily along a pipe. Valves are suddenly closed
  at each end, bringing the air to rest there. The subsequent
  distributions of velocity and pressure in the pipe are to be
  predicted.

  The analysis is one-dimensional, and the process is supposd to
  be isentropic, ie without either friction or heat transfer.

  ENDDIS
  DISPLAY
  High-pressure air is separated from low-pressure air by a
  diaphragm at the mid-section of a tube with closed ends.

  The diaphragm suddenly breaks.

  The subsequent motion is predicted by a one-dimensional analysis,
  similar to that of case 321, which is first loaded.

  The differences from case 321 are:-

  (1) the initial conditions of pressure and velocity; and
  (2) the solution of the "marker" variable C3, which serves as a
      multiplier in the density formula.

  ENDDIS
  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 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(Shock Tube With 2D Obstacle             )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 332
 ************************************************************
  Group 2. Time dependence
 STEADY = F
    * Set overall time and no. of steps
 TFIRST =0. ;TLAST =1.0E-02
 FSTEP = 1 ;LSTEP = 100
 TFRAC(1)=1.0E-02 ;TFRAC(21)=0.21
 TFRAC(41)=0.41 ;TFRAC(61)=0.61
 TFRAC(81)=0.81
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 100
 XULAST =2.
 XFRAC(1)=1.0E-02 ;XFRAC(21)=0.21
 XFRAC(41)=0.41 ;XFRAC(61)=0.61
 XFRAC(81)=0.81
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 20
 YVLAST =1.
 YFRAC(1)=0.05 ;YFRAC(2)=0.1
 YFRAC(3)=0.15 ;YFRAC(4)=0.2
 YFRAC(5)=0.25 ;YFRAC(6)=0.3
 YFRAC(7)=0.35 ;YFRAC(8)=0.4
 YFRAC(9)=0.45 ;YFRAC(10)=0.5
 YFRAC(11)=0.55 ;YFRAC(12)=0.6
 YFRAC(13)=0.65 ;YFRAC(14)=0.7
 YFRAC(15)=0.75 ;YFRAC(16)=0.8
 YFRAC(17)=0.85 ;YFRAC(18)=0.9
 YFRAC(19)=0.95 ;YFRAC(20)=1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 1
 ZWLAST =1.
 ZFRAC(1)=1.
 ************************************************************
  Group 6. Body-Fitted Coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = T
 NAME(1)=P1 ;NAME(3)=U1
 NAME(5)=V1 ;NAME(18)=C3
 NAME(149)=PRPS ;NAME(150)=RHO1
    * Y in SOLUTN argument list denotes:
    * 1-stored 2-solved 3-whole-field
    * 4-point-by-point 5-explicit 6-harmonic averaging 
 SOLUTN(P1,Y,Y,N,N,N,Y)
 SOLUTN(U1,Y,Y,N,N,N,Y)
 SOLUTN(V1,Y,Y,N,N,N,Y)
 SOLUTN(C3,Y,Y,N,N,N,Y)
 SOLUTN(PRPS,Y,N,N,N,N,Y)
 SOLUTN(RHO1,Y,N,N,N,N,Y)
 DEN1 = 150
 PRPS = 149
 ************************************************************
  Group 8. Terms & Devices
    * Y in TERMS argument list denotes:
    * 1-built-in source 2-convection 3-diffusion 4-transient
    * 5-first phase variable 6-interphase transport         
 TERMS(P1,Y,Y,Y,N,Y,Y)
 TERMS(U1,Y,Y,N,Y,Y,Y)
 TERMS(V1,Y,Y,N,Y,Y,Y)
 TERMS(C3,Y,Y,N,Y,Y,Y)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 NEWRH1 = T
 ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
 ************************************************************
  Group 9. Properties used if PRPS is not
  stored, and where PRPS = -1.0 if it is!
 RHO1 = GRND3 ;TMP1 =0. ;EL1 =0.
 TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
 DVO1DT =3.33E-03 ;DRH1DP = GRND3
 RHO1A =2.682563E-04 ;RHO1B =0.71429 ;RHO1C =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 ENUL =1.569E-05 ;ENUT =0.
 CP1 =1007. ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(P1)=1.0E+05 ;FIINIT(U1)=0.
 FIINIT(V1)=0. ;FIINIT(C3)=1.
 FIINIT(PRPS)=-1. ;FIINIT(RHO1)=1.
 
 PATCH(START ,INIVAL, 1, 50, 1, 20, 1, 1, 1, 1)
 INIT(START ,P1 ,0. ,2.0E+05 )
 INIT(START ,C3 ,0. ,1.21901 )
 INIT(START ,RHO1,0. ,2. )
 INIADD = F
 FSWEEP = 1
 NAMFI =CHAM
 ************************************************************
  Group 12. Patchwise adjustment of terms
  Patches for this group are printed with those
  for Group 13.
  Their names begin either with GP12 or &
 ************************************************************
  Group 13. Boundary & Special Sources
 
 PATCH(OBSTACLE,CELL , 50, 50, 1, 10, 1, 1, 1, 100)
 COVAL(OBSTACLE,U1 , FIXVAL ,0. )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 5 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-05
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=20 ;LITER(U1)=10
 LITER(V1)=10 ;LITER(C3)=20
 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
 ENDIT(V1)=1.0E-03 ;ENDIT(C3)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,1.)
 RELAX(U1,FALSDT,1.)
 RELAX(V1,FALSDT,1.)
 RELAX(C3,FALSDT,1.0E+09)
 RELAX(PRPS,LINRLX,1.)
 RELAX(RHO1,LINRLX,1.)
 OVRRLX =0.
 EXPERT = F ;NNORSL = F
 ************************************************************
  Group 18. Limits
 VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
 VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06
 VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
 VARMAX(C3)=1.0E+10 ;VARMIN(C3)=-1.0E+10
 VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
 VARMAX(RHO1)=1.0E+10 ;VARMIN(RHO1)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 PARSOL = F
 ISG62 = 1
 SPEDAT(SET,DOMAIN,PHASE_1_MAT,I,1)
 SPEDAT(SET,GXMONI,TRANSIENT,L,F)
 SPEDAT(SET,GXMONI,PLOTALL,L,T)
 ************************************************************
  Group 20. Preliminary Printout
 ************************************************************
  Group 21. Print-out of Variables
 INIFLD = F ;SUBWGR = F
    * Y in OUTPUT argument list denotes:
    * 1-field 2-correction-eq. monitor 3-selective dumping      
    * 4-whole-field residual 5-spot-value table 6-residual table
 OUTPUT(P1,Y,N,Y,Y,Y,Y)
 OUTPUT(U1,Y,N,Y,Y,Y,Y)
 OUTPUT(V1,Y,N,Y,Y,Y,Y)
 OUTPUT(C3,Y,N,Y,Y,Y,Y)
 OUTPUT(PRPS,Y,N,Y,N,N,N)
 OUTPUT(RHO1,Y,N,Y,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 50 ;IYMON = 11 ;IZMON = 1
 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000 ;NUMCLS = 5
 NTPRIN = 20 ;ISTPRF = 1 ;ISTPRL = 100000
 NXPRIN = 20 ;IXPRF = 1 ;IXPRL = 10000
 NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000
 IPLTF = 1 ;IPLTL = -1 ;NPLT = -1
 ISWPRF = 1 ;ISWPRL = 100000
 ITABL = 3 ;IPROF = 1
 ABSIZ =0.5 ;ORSIZ =0.4
 NTZPRF = 1 ;NCOLPF = 50
 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
 
 PATCH(ALLX ,PROFIL, 1, 100, 1, 1, 1, 1, 1, 100)
 PLOT(ALLX ,P1 ,0. ,0. )
 PLOT(ALLX ,U1 ,0. ,0. )
 
 PATCH(TIMEPLOT,PROFIL, 50, 50, 1, 1, 1, 1, 1, 100)
 PLOT(TIMEPLOT,P1 ,0. ,0. )
 PLOT(TIMEPLOT,U1 ,0. ,0. )
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
 IDISPA = 5 ;IDISPB = 1 ;IDISPC = 100
STOP