*HEADING
 OBLIQUE SHOCK REFLECTIONS
 WEDGE HALF-ANGLE IS 15.13 DEGREE
*NODE
    1,   0.,      0.
   11,   8.,      0.    
 1651,   0.,     80.
 1661,   8.,     80. 
 2091,   4.056,  95.
 2101,   8.,     95.
 2751,   4.056, 150.
 2761,   8.,    150.
 5000,   0.,     75.
*NGEN,NSET=INFLOW
 1,11,1
*NGEN,NSET=WEDGEB
 1651,1661,1
*NGEN,NSET=WEDGET
 2091,2101,1
*NGEN,NSET=OUTFLOW
 2751,2761,1
*NFILL,NSET=ALL
 INFLOW,WEDGEB,150,11
 WEDGEB,WEDGET, 40,11
*NFILL,NSET=ALL,BIAS=0.975 
 WEDGET,OUTFLOW,60,11
*NSET,NSET=RIGHT,GENERATE
 11,2761,11
*NSET,NSET=BOT,GENERATE
 1,1661,1
*NSET,NSET=TOP,GENERATE
 1662,2761,1
**
*ELEMENT,TYPE=CPE4R
 1,1,2,13,12
*ELGEN,ELSET=GAS
 1,10,1,1,250,11,10
*ELSET,ELSET=BOT,GENERATE
 1,1500,1
*ELSET,ELSET=TOP,GENERATE
 1501,2500,1
*ELSET,ELSET=LEFT,GENERATE
 1,2491,10
*ELSET,ELSET=LOWER,GENERATE
 1,10,1
*ELSET,ELSET=UPPER,GENERATE
 2491,2500,1
*ELSET,ELSET=QA_TEST_UPPER,GENERATE
 2491,2500,1
**
*NSET,NSET=VIEWTOP,GENERATE
 2201,2211,1
*NSET,NSET=VIEWBOT,GENERATE
 1541,1551,1
*ELSET,ELSET=VIEW1,GENERATE
 1401,2000,1
*ELSET,ELSET=VIEW2,GENERATE
 1401,1900,1
**
*SOLID SECTION,ELSET=GAS,MATERIAL=GAS,CONTROLS=CONT
 1.,
*SECTION CONTROLS,NAME=CONT,HOURGLASS=VISCOUS
 1.,
*MATERIAL,NAME=GAS
*EOS,TYPE=IDEAL GAS
 0.2, 
*SPECIFIC HEAT
 0.5,
*DENSITY,DEP=1
 1., , 1.
 6., , 2. 
**
*INITIAL CONDITIONS,TYPE=STRESS
 TOP, -2.E-6, -2.E-6, -2.E-6
 BOT, -1.2, -1.2, -1.2
*INITIAL CONDITIONS,TYPE=FIELD,VAR=1
 TOP, 1.
 BOT, 2.
*INITIAL CONDITIONS,TYPE=VELOCITY
 BOT, 2, 1.
**
*SURFACE,TYPE=SEGMENT,NAME=WALL
 START, 4.056, 200.
 LINE,  4.056,  95.
 LINE,  0.,     80.
 LINE,  0.,    -50.
*SURFACE,TYPE=ELEMENT,NAME=IN,REGION TYPE=EULERIAN
 LOWER,S1 
*SURFACE,TYPE=ELEMENT,NAME=OUT,REGION TYPE=EULERIAN
 UPPER,S3 
*SURFACE,TYPE=ELEMENT,NAME=GAS
 LEFT,S4
*RIGID BODY,REF=5000, ANALYTICAL SURFACE =WALL
*STEP
*DYNAMIC,EXPLICIT
 ,150.
**
*SURFACE INTERACTION,NAME=GAS_WALL
*CONTACT PAIR,INTERACTION=GAS_WALL
 WALL,GAS
**
*BOUNDARY,OP=NEW,TYPE=VELOCITY,REGION TYPE=EULERIAN
 INFLOW,2,2,1.
*BOUNDARY,OP=NEW
 RIGHT,1,1,0.
 5000, 1,6,0.
**
*ADAPTIVE MESH,ELSET=GAS,MESHSWEEPS=5
*ADAPTIVE MESH CONSTRAINT
 VIEWBOT,2,2,0.0
 VIEWTOP,2,2,0.0
 INFLOW,1,2,0.0
 OUTFLOW,2,2,0.0
**
*FILE OUTPUT,TIMEMARKS=YES,NUMBER INTERVAL=4
*EL FILE,ELSET=UPPER
 TEMP,
 PRESS,
*OUTPUT, FIELD, VARIABLE=PRESELECT,NUMBER INTERVAL=2
*OUTPUT,FIELD,NUMBER INTERVAL=150
*ELEMENT OUTPUT,ELSET=GAS
PRESS,
*OUTPUT,FIELD,NUMBER INTERVAL=4
*ELEMENT OUTPUT,ELSET=QA_TEST_UPPER
 TEMP,
 PRESS,
**
*END STEP