&dimen -dimen meters m-tons
&set piledia = 72*25.4    $ pile diameter in millimeters
$ note: seg len should be approx greater then len/dia
&set fmul    = &number(real (3.1416*%piledia%)/1000**2)
&set pmul    = 9.8e-5     $ to convert kn/m to m.tons/mm
&set ymul    = 10         $ to convert cm to mm
&data soil s_china
depth 0
tz 0 0 0.11 0.7 0.14 1.5 0.21 3.7 0.26 7.3 0.28 11.0 0.30 18.3 0.30 228.6
py 0 0 30.5 0.34 39.4 1.02 44.4 1.70 48.1 2.37 51.0 3.05 66.1 6.86 66.1 500.0
depth 3.0
tz 0 0 0.11 0.7 0.14 1.5 0.21 3.7 0.26 7.3 0.28 11.0 0.30 18.3 0.30 228.6
py 0 0 30.5 0.34 39.4 1.02 44.4 1.70 48.1 2.37 51.0 3.05 66.1 6.86 66.1 500.0
depth 4.0
py 0 0 66.8 0.47 81.7 1.11 90.9 1.76 97.7 2.40 103.2 3.05 133.2 6.86 133.2 500.00
depth 5.50
tz 0 0 .37 .7 .50 1.5 .71 3.7 .88 7.3 .97 11.0 1.02 18.3 1.02 228.6
depth 6.0
py 0 0 138.8 0.55 166.9 1.18 185.2 1.80 199.1 2.42 210.6 3.05 274.8 6.86 +
   274.8 500.00
depth 8.00
tz 0 0 .51 .7 .69 1.5 .99 3.7 1.23 7.3 1.36 11.0 1.43 18.3 1.43 228.6
py 0 0 91.8 0.26 168.1 0.95 216.3 1.65 254.3 2.35 286.5 3.05 450.8 6.86 +
   450.8 500.00
depth 8.1
tz 0 0 .66 .7 .89 1.5 1.27 3.7 1.59 7.3 1.74 11.0 1.84 18.3 1.84 228.6
py 0.0 0.00 77.4 1.92 129.0 9.14 185.7 27.42 125.5 137.10 125.5 500.00
depth 10.00
tz 0 0 .73 .7 .99 1.5 1.42 3.7 1.76 7.3 1.94 11.0 2.04 18.3 2.04 228.6
py 0.0 0.00 90.8 1.92 151.3 9.14 217.9 27.42 191.7 137.10 191.7 500.00
depth 10.1
tz 0 0 .66 .7 .89 1.5 1.27 3.7 1.59 7.3 1.74 11.0 1.84 18.3 1.84 228.6
py 0 0 131.9 0.28 214.6 0.97 314.0 1.66 372.1 2.36 421.9 3.05 679.3 6.86 +
   679.3 500.00
depth 15.0
py 0 0 385.1 0.52 563.9 1.15 695.6 1.78 804.5 2.42 899.4 3.05 1439.1 6.86 +
   1439.1 500.00
depth 21.00
tz 0 0 1.28 .7 1.73 1.5 2.48 3.7 3.09 7.3 3.39 11.0 3.57 18.3 3.57 228.6
py 0 0 751.6 0.71 1004.6 1.29 1202.2 1.88 1369.5 2.46 1517.2 3.05 2427.6 6.86 +
   2427.6 500.00
depth 21.1
tz 0 0 1.21 .7 1.63 1.5 2.33 3.7 2.90 7.3 3.18 11.0 3.35 18.3 3.35 228.6
py 0.0 0.00 162.9 0.96 271.5 4.57 390.9 13.71 390.9 500.00
depth 30.00
tz 0 0 1.69 .7 2.28 1.5 3.26 3.7 4.06 7.3 4.45 11.0 4.69 18.3 4.69 228.6
py 0.0 0.00 227.0 0.96 378.4 4.57 544.9 13.71 544.9 500.00
depth 30.1
tz 0 0 1.76 .7 2.38 1.5 3.40 3.7 4.23 7.3 4.65 11.0 4.89 18.3 4.89 228.6
py 0 0 1082.2 0.70 1451.8 1.28 1739.6 1.87 1983.3 2.46 2198.1 3.05 3517.0 6.86 +
   3517.0 500.00
depth 38.00
tz 0 0 1.76 .7 2.38 1.5 3.40 3.7 4.23 7.3 4.65 11.0 4.89 18.3 4.89 228.6
py 0 0 1392.5 0.70 1865.0 1.29 2233.5 1.87 2545.5 2.46 2820.7 3.05 4513.1 6.86 +
   4513.1 500.00
depth 38.1
tz 0 0 2.13 .7 2.87 1.5 4.11 3.7 5.12 7.3 5.61 11.0 5.91 18.3 5.91 228.6
py 0.0 0.00 286.3 0.96 477.1 4.57 687.0 13.71 687.0 500.00
depth 55.00
tz 0 0 3.12 .7 4.21 1.5 6.02 3.7 7.50 7.3 8.23 11.0 8.67 18.3 8.67 228.6
py 0.0 0.00 419.5 0.96 699.2 4.57 699.2 500.00
qw -1,-10. 0 0 69.6 4.4 87.6 8.8 118.9 21.9 149.9 43.9 171.5 65.8 203.4 109.7 +
   203.4 365.6
depth 55.1
tz 0 0 1.76 .7 2.38 1.5 3.40 3.7 4.23 7.3 4.65 11.0 4.89 18.3 4.89 228.6
py 0 0 2122.5 0.73 2810.6 1.31 3351.4 1.89 3811.0 2.47 4217.1 3.05 6747.3 6.86 +
   6747.3 500.00
depth 69.00
tz 0 0 1.76 .7 2.38 1.5 3.40 3.7 4.23 7.3 4.65 11.0 4.89 18.3 4.89 228.6
py 0 0 2694.8 0.73 3561.7 1.31 4244.0 1.89 4824.0 2.47 5336.8 3.05 8538.9 6.86 +
   8538.9 500.00
qw -1,-10. 0 0 173.9 4.4 219.1 8.8 297.4 21.9 374.7 43.9 428.9 65.8 508.5 109.7 +
   508.5 365.6
depth 69.1
tz 0 0 3.93 .7 5.30 1.5 7.58 3.7 9.44 7.3 10.36 11.0 10.91 18.3 10.91 228.6
py 0.0 0.00 528.1 0.96 880.2 4.57 880.2 500.00
qw -1,-10. 0 0 87.9 2.9 110.7 5.8 150.3 14.6 189.4 29.2 216.8 43.9 257.0 73.1 +
   257.0 365.6
depth 90.0
tz 0 0 5.18 .7 6.98 1.5 9.98 3.7 12.44 7.3 13.65 11.0 14.38 18.3 +
   14.38 228.6
py 0.0 0.00 695.9 0.96 1159.9 4.57 1159.9 500.00
qw -1,-10. 0 0 116.2 2.9 146.4 5.8 198.7 14.6 250.4 29.2 286.6 43.9 339.9 73.1 +
   339.9 365.6
depth 90.1
tz 0 0 1.76 .7 2.38 1.5 3.40 3.7 4.23 7.3 4.65 11.0 4.89 18.3 4.89 228.6
py 0 0 3642.4 0.76 4773.1 1.33 5668.3 1.90 6431.4 2.47 7107.1 3.05 11371.4 6.86 +
   11371.4 500.00
depth 100.00
tz 0 0 1.76 .7 2.38 1.5 3.40 3.7 4.23 7.3 4.65 11.0 4.89 18.3 4.89 228.6
py 0 0 3642.4 0.76 4773.1 1.33 5668.3 1.90 6431.4 2.47 7107.1 3.05 11371.4 6.86 +
   11371.4 500.00
qw -1,-10. 0 0 173.9 4.4 219.1 8.8 297.4 21.9 374.7 43.9 428.9 65.8 508.5 109.7 +
   508.5 365.6
depth 100.1
tz 0 0 5.80 .7 7.83 1.5 11.19 3.7 13.93 7.3 15.29 11.0 16.11 18.3 16.11 228.6
qw -1,-10. 0 0 130.0 2.9 163.8 5.8 222.3 14.6 280.0 29.2 320.6 43.9 380.1 73.1 +
   380.1 365.6
depth 117.00
tz 0 0 6.86 .7 9.26 1.5 13.24 3.7 16.49 7.3 18.10 11.0 19.07 18.3 19.07 228.6
qw -1,-10. 0 0 153.8 2.9 193.8 5.8 263.0 14.6 331.3 29.2 379.3 43.9 449.7 73.1 +
   449.7 365.6
depth 117.1
tz 0 0 1.76 .7 2.38 1.5 3.40 3.7 4.23 7.3 4.65 11.0 4.89 18.3 4.89 228.6
qw -1,-10. 0 0 173.9 4.4 219.1 8.8 297.4 21.9 374.7 43.9 428.9 65.8 508.5 109.7 +
   508.5 365.6
depth 120.00
tz 0 0 1.76 .7 2.38 1.5 3.40 3.7 4.23 7.3 4.65 11.0 4.89 18.3 4.89 228.6
qw -1,-10. 0 0 173.9 4.4 219.1 8.8 297.4 21.9 374.7 43.9 428.9 65.8 508.5 109.7 +
   508.5 365.6
&dimen -dimen feet kips
end_&data
 *p1104 0 0 -360
&describe body p
     ~pile tube 72.00 2.25 -soil s_china +
                           -refine 20        +
                           -pymult  %pmul%   %ymul% +
                           -tzmult  %fmul%  -len 100
     ~pile tube 72.00 2.25 -soil s_china +
                           -refine 20        +
                           -pymult  %pmul%   %ymul% +
                           -tzmult  %fmul%
     ~pile tube 72.00 2.00 -soil s_china +
                           -refine 52        +
                           -pymult  %pmul%   %ymul%  +
                           -tzmult  %fmul% -len 160
     ~30_100 tube 30 1.0
$
$ define the nodes on the piles
 *j1104 0 0 0
 *j1004 0 0 10
$
&set pf = 1.5
  beam 1 ~30_100 *j1004 *j1104