*heading
 shock tube (adaptive meshing with two domains)
 partially saturated sand
*parameter
 s = 1.93
 cs = 1480
 ce = 600
 ps = 6.5e6
 pe = 0.0
 alpha0 = 1.052364
 rhos0 = 2070
 gamma0 = 0.880
 el = 5e4
 pl = -4.23e8
 n0 = 1.0 - 1.0/alpha0
*restart,write,number=30
*node
   1, 0.0, 0.0
 101, 2.5, 0.0
 201, 5.0, 0.0
 202, 0.0, 0.1
 302, 2.5, 0.1
 402, 5.0, 0.1
*ngen,nset=botleft
 1,101,1
*ngen,nset=topleft
 202,302,1
*nset,nset=left
 botleft,topleft
*ngen, nset=qa_test_allnodes
 101,201,1
*ngen, nset=qa_test_allnodes
 302,402,1
*nset,nset=botright,generate
 102,201,1
*nset,nset=topright,generate
 303,402,1
*nset,nset=right
 botright,topright
*nset,nset=all
 left,right
*nset,nset=xsym
 1,201,202,402
*nset,nset=top
 topleft,topright
*nset,nset=bot
 botleft,botright
**
*element,type=cpe4r
 1,1,2,203,202
*elgen,elset=exp
 1,200,1,1
*elset,elset=expleft,generate
 1,100,1
*elset,elset=expright,generate
 101,200,1
*elset, elset=left, generate
 1,100,1
*elset, elset=right, generate
 101,200,1
**
*solid section,elset=exp,material=sand
 1.,
*material,name=sand
*eos,type=usup
 <cs>,<s>,<gamma0>
*eos compaction
 <ce>,<n0>,<pe>,<ps>
*density
 <rhos0>, 
*boundary
 top,ysymm
 bot,ysymm
 xsym,1
**
*initial conditions, type = specific energy
 left, <el>
*initial conditions, type = stress
 left, <pl>, <pl>, <pl>
*elset,elset=qa_test_allelems
 exp,
**
*nset, nset=qa_test_allnodes
 botleft, topleft, left, qa_test_allnodes
 botright, topright, right, all
 xsym, top, bot
**
*step
*dynamic,explicit
 ,0.001
*bulk viscosity
 0.3
**
*adaptive mesh,elset=expleft, frequency=1,initial mesh sweeps=1
*adaptive mesh,elset=expright,frequency=1,initial mesh sweeps=1
**
*output,field
*element output
 s, density, ener, palph, palphmin
*node output
 u
*output,history,variable=preselect
*end step