*heading checking stokes-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 21.9 43.8 65.7 87.6 109.6 131.5 153.4 175.3 ** **v-d **23.516 21.50 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 **232.81 21.08 20.68 0.00 0.00 0.00 0.00 0.00 0.00 0.00 **225.75 19.88 19.50 18.38 16.59 0.00 0.00 0.00 0.00 0.00 **214.00 18.05 17.72 16.72 15.11 12.98 10.41 0.00 0.00 0.00 **197.53 15.82 15.53 14.67 13.29 11.46 9.25 6.78 4.15 1.47 **176.37 13.42 13.18 12.48 11.34 9.81 7.97 5.90 3.69 1.42 **150.50 11.10 10.91 10.34 9.42 8.19 6.69 5.01 3.19 1.32 **119.93 9.05 8.90 8.45 7.72 6.74 5.54 4.19 2.72 1.19 ** 84.66 7.42 7.30 6.94 6.35 5.57 4.60 3.51 2.31 1.06 ** 44.68 6.34 6.24 5.93 5.44 4.78 3.97 3.04 2.03 0.96 ** 0.00 5.94 5.84 5.56 5.11 4.49 3.74 2.87 1.92 0.92 ** ** velocity-vy ** **h-d -> 0.0 21.9 43.8 65.7 87.6 109.6 131.5 153.4 175.3 ** **v-d **23.516 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 **232.81 0.00 3.83 0.00 0.00 0.00 0.00 0.00 0.00 0.00 **225.75 0.00 3.57 6.99 10.09 0.00 0.00 0.00 0.00 0.00 **214.00 0.00 3.19 6.24 9.02 11.42 13.35 0.00 0.00 0.00 **197.53 0.00 2.72 5.33 7.71 9.77 11.45 12.66 13.41 13.67 **176.37 0.00 2.21 4.34 6.28 7.98 9.37 10.39 11.03 11.29 **150.50 0.00 1.71 3.35 4.86 6.18 7.27 8.09 8.62 8.85 **119.93 0.00 1.23 2.42 3.51 4.47 5.27 5.88 6.28 6.48 ** 84.66 0.00 0.80 1.56 2.27 2.90 3.42 3.83 4.10 4.24 ** 44.68 0.00 0.39 0.77 1.13 1.44 1.70 1.90 2.04 2.11 ** 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 60.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,nlgeom *dynamic 12.5,12.5 *elfile sdv, *Output, field *end step