*HEADING
 Test to check enhanced drag (Geers & Park 2005) in bubble model; NEW interface
*PREPRINT,ECHO=NO,HISTORY=yes,MODEL=YES
*NODE, NSET=ALL
    1, 4.95, -2.35, -2.45
    2, 5.05, -2.35, -2.45
    3, 5.05, -2.45, -2.35
    4, 4.95, -2.45, -2.35
    5,13.660369766932588, -8.4735, -8.4735
    6,13.660369766932588, -8.4735, -8.5735
    7,13.660369766932588, -8.5735, -8.5735
    8,13.660369766932588, -8.5735, -8.4735
*NODE,NSET=STANDOFF
9,5., -2.3999, -2.3999
*NODE,NSET=SOURCE            
10,5., -8.5235, -8.5235    
*NSET,NSET=QA_TEST,GENERATE
1,8
*ELEMENT,TYPE=s4,ELSET=CHANNEL
**NOTE THAT AREA OF ELEMENT 1 IS SQRT(2)/100
**                          2 IS 1/100
 1,  1,2,3,4
 2,  5,6,7,8
**************************
*SURFACE,NAME=IW_LOAD1
1,sneg
2,sneg
**************************
*shell SECTION,MATERIAL=SOLID,ELSET=CHANNEL
 1.E-2,5
*MATERIAL,NAME=SOLID
*DENSITY
 8.0E3,
*ELASTIC
 2.322585E9,0.3
**************************
*INCIDENT WAVE INTERACTION PROPERTY,NAME=IWPROP,TYPE=SPHERE
1500.00,1.025E3,0.0
*UNDEX
52.117E6, .0895690e-3, 0.18, 0.185, 1045.e6, 1.27,
1.6E3, 305.9, 9.81E0, 1.0E5, 1.0, 2.25, 
1.3,
36.019,0.0, 0.0, 1.0
**************************
*BOUNDARY
ALL,ENCASTRE
**************************
*STEP,inc=100000
*DYNAMIC,DIRECT
6.7e-4,1.3
*INCIDENT WAVEINTERACTION,UNDEX,PROPERTY=IWPROP
IW_LOAD1,source,standoff,1.0
**
*OUTPUT,history, frequency=1
*NODE OUTPUT,nset=qa_test
RF,
** TO RECOVER INCIDENT PRESSURES ON THE ELEMENT FACES FROM THE
** RF'S, USE:
** (ELEMENT 1) RF * 4 /100 = P_INCIDENT
** (ELEMENT 2) RF * 4 /100 = P_INCIDENT
*ENDSTEP