*HEADING
 Same as iw_1d_ac2d3_xpl_s_pa.inp except that now total wave
 formulation is used.
 1D transient wave propagation, acoustic medium, incident wave load, spherical
 and defined through pressure amplitude
**
** C = root (K/rho)= 343.05 m/s   K=bulk mod, rho=density.
**
** Analytical solution to prescribed constant acceleration of a0 at x=0
**
** p(x,t) = (K a0 / c^2) (ct - x)  for x < ct
**        = 0  			   for x > ct
**
*ACOUSTIC WAVE, TYPE=TOTAL
*PREPRINT,ECHO=NO,HISTORY=NO,MODEL=YES
*NODE, NSET=ODB_NSET
    1, 0,
  101, 1,
 1001, 0, 1.E-2
 1101, 1, 1.E-2
*NGEN,NSET=CHANNEL
    1,101
 1001,1101
*NSET,NSET=SAMPLE,GENERATE
 1,101
*ELEMENT,TYPE=AC2D3
 1,1,2,1002
 2,1002,1001,1
*ELGEN,ELSET=CHANNEL
 1,100,1,2
 2,100,1,2
**************************
*SURFACE,NAME=IW_LOAD
2,S2
**************************
** CROSS-SECTION IS 1.0E-4
*SOLID SECTION,MATERIAL=AIR,ELSET=CHANNEL
 1.E-2,
**
*MATERIAL,NAME=AIR
*DENSITY
 1.21,
*ACOUSTIC MEDIUM,BULK MODULUS
 0.1424E6,
**
**************************
*INCIDENT WAVE PROPERTY,NAME=IWPROP,TYPE=SPHERE
0.0,0.0,
-10.0E04,0.0,
*INCIDENT WAVE FLUID PROPERTY
0.1424E6,1.21,0.0
*AMPLITUDE,DEFINITION=TABULAR,NAME=PRESSUREVTIME
0.0,0.0,1.0,1.0
**************************
*RESTART,WRITE
**
*STEP
*DYNAMIC,EXPLICIT
 ,.0044,
** DEFINE NONREFLECTIVE BOUNDARY AT X=1
*IMPEDANCE
 199,I2
**
*INCIDENT WAVE,PRESSUREAMPLITUDE=PRESSUREVTIME,PROPERTY=IWPROP
IW_LOAD,415.09517
**
*NSET,NSET=QA_TEST
ODB_NSET,
*OUTPUT,FIELD,NUMBER INTERVAL=2
*NODE OUTPUT
 POR,
*NODE OUTPUT,NSET=QA_TEST
 POR,
*OUTPUT,HISTORY,TIME INTERVAL=0.1E-3
*NODE OUTPUT,NSET=SAMPLE
 POR,
*ENDSTEP