*heading
checking airy-wave with getwavevel utility subroutine    
*****************velocities from the wave theory
**
** h-d  -- horizontal distance from the crest
** v-d  -- vertical position
**
** velocity-vx
**
**h-d ->   0.0  20.7  41.5  62.2  83.0  103.7 124.5  145.2 166.0
**
**v-d
**200.00  16.12 15.88 15.16 13.98 12.37 10.40  8.10  5.57  2.85
**198.00  15.87 15.64 14.92 13.76 12.18 10.24  7.97  5.48  2.81
**192.00  15.15 14.92 14.24 13.13 11.62  9.77  7.61  5.23  2.68
**182.00  14.03 13.82 13.19 12.16 10.76  9.05  7.05  4.84  2.48
**168.00  12.62 12.43 11.87 10.94  9.68  8.14  6.34  4.36  2.23
**150.00  11.07 10.90 10.40  9.60  8.49  7.14  5.56  3.82  1.96
**128.00   9.51  9.37  8.94  8.24  7.30  6.13  4.78  3.28  1.68
**102.00   8.08  7.96  7.60  7.01  6.20  5.21  4.06  2.79  1.43
** 72.00   6.91  6.81  6.50  5.99  5.30  4.46  3.47  2.39  1.22
** 38.00   6.12  6.03  5.75  5.30  4.69  3.94  3.07  2.11  1.08
**  0.00   5.82  5.73  5.47  5.04  4.46  3.75  2.92  2.01  1.03
**
** velocity-vy
**
**h-d ->   0.0  20.7  41.5  62.2  83.0  103.7 124.5 145.2 166.0
**
**v-d
**200.00   0.00  2.60  5.13  7.50  9.65 11.50 13.00 14.11 14.80
**198.00   0.00  2.55  5.04  7.36  9.47 11.29 12.77 13.86 14.54
**192.00   0.00  2.42  4.77  6.97  8.97 10.69 12.10 13.13 13.77
**182.00   0.00  2.21  4.36  6.36  8.19  9.76 11.04 11.98 12.56
**168.00   0.00  1.94  3.82  5.58  7.19  8.56  9.69 10.51 11.02
**150.00   0.00  1.63  3.21  4.69  6.04  7.20  8.14  8.84  9.27
**128.00   0.00  1.30  2.57  3.75  4.83  5.75  6.50  7.06  7.40
**102.00   0.00  0.97  1.91  2.80  3.60  4.29  4.85  5.27  5.52
** 72.00   0.00  0.65  1.27  1.86  2.39  2.85  3.23  3.50  3.67
** 38.00   0.00  0.33  0.64  0.94  1.21  1.44  1.63  1.77  1.86
**  0.00   0.00  0.00  0.00  0.00  0.00  0.00  0.00  0.00  0.00
*****************
*aqua
0., 200., 32.2, 1.99
0., 0., 0., 0.
0., 0., 0., 200.
*wave, type=airy, wave period
30.0, 12.5, 180.0, 1.0
*node,nset=pr
101,
102,0.5
*ngen,nset=n2
101,102
*user element,nodes=2,type=u1,properties=4,coordinates=3,variables=396
1,2,3           
*element,type=u1,elset=utruss
1 , 101, 102
*uel property,elset=utruss
0.002, 2.1e11, 0.3, 7200.
*boundary
n2,1,3
**
*step
*dynamic
12.5,12.5
*elfile
sdv, 
*Output,field
*end step