*heading
 shock tube  (lagrangian analysis)
 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
*ngen,nset=botleft
*ngen,nset=topleft
*nset,nset=left
 botleft,topleft
*ngen, nset=qa_test_allnodes
*ngen, nset=qa_test_allnodes
*nset,nset=botright,generate
*nset,nset=topright,generate
*nset,nset=right
 botright,topright
*nset,nset=all
 left,right
*nset,nset=xsym
*nset,nset=top
 topleft,topright
*nset,nset=bot
 botleft,botright
**
*element,type=cpe4r
*elgen,elset=exp
*elset, elset=left, generate
*elset, elset=right, generate
**
*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
**
*output,field
*element output
 s, density, ener, palph, palphmin
*node output
 u
*output,history,variable=preselect
*end step