*HEADING
 Same as iwr_1d_1_xpl_s_pa.inp except that the total wave formulation is used
 and also the incident wave load sign is reversed.  There should not be any 
 change in the solution. 
**
 1D transient wave propagation, acoustic medium, incident wave 
 load and incident wave reflection, spherical and defined through 
 pressure amplitude
**
*ACOUSTIC WAVE, TYPE=TOTAL
*PREPRINT,ECHO=NO,HISTORY=NO,MODEL=YES
*NODE,NSET=PLATE
 1,0.,0.,0.
 2,1.,0.,0.
 3,1.,1.,0.
 4,0.,1.,0.
 5,0.,0.,-1.
 6,1.,0.,-1.
 7,1.,1.,-1.
 8,0.,1.,-1.
*ELEMENT,TYPE=AC3D8R,ELSET=ACOUSTIC
 1,1,4,3,2,5,8,7,6
**************************
*SURFACE,NAME=IW_LOAD
1,S1
**************************
*SOLID SECTION,MATERIAL=AIR,ELSET=ACOUSTIC
**
*MATERIAL,NAME=AIR
*DENSITY
 1.21,
*ACOUSTIC MEDIUM,BULK MODULUS
 0.1424E6,
**
**************************
*INCIDENT WAVE PROPERTY,NAME=IWPROP,TYPE=SPHERE
0.0,0.0,0.
0.0,0.0,10.
*INCIDENTWAVE FLUID PROPERTY
1.E8,1.E2
*AMPLITUDE,DEFINITION=TABULAR,NAME=PRESSUREVTIME
0.,1.,1.,1.
**
**************************
**************************
*STEP
*DYNAMIC,EXPLICIT,DIRECT USER CONTROL
 0.00001,0.02
**
*IMPEDANCE,non=pl
1,I2
*INCIDENTWAVE,PRESSUREAMPLITUDE=PRESSUREVTIME,PROPERTY=IWPROP
IW_LOAD,415.09517
**
*INCIDENT WAVE REFLECTION
5.,0.,1.,0.,1.0
**
*OUTPUT,FIELD,NUMBER INTERVAL=2
*NODE OUTPUT
POR,
*OUTPUT,HISTORY,TIME INTERVAL=0.1E-3
*NODE OUTPUT,NSET=PLATE
POR,
**
*END STEP