$ $@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ $@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ $@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ $ $ $ ---------- MODEL FOR BARGE ------------ $ $ $@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ $@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ $@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ $ $ $********************************************* read macros $ &insert sfun.mac $ $********************************************* define body $ &describe body info $ $********************************************* hydrostatic model $ pgen -tanaka 1 plane 0 50 100 150 200 250 300 -rect 0 20 90 end pgen $ $********************************************* input wind spectra $ &data curves f_spectrum w_f \ 0.01, 7079.67 \ 0.02, 5011.39 \ 0.03, 3938.40 \ 0.04, 3254.06 \ 0.05, 2769.01 \ 0.07, 2381.27 \ 0.08, 2066.94 \ 0.10, 1819.84 \ 0.11, 1619.22 \ 0.13, 1446.77 \ 0.14, 1294.72 \ 0.16, 1165.95 \ 0.18, 1046.27 \ 0.20, 940.28 \ 0.22, 864.78 \ 0.24, 797.16 \ 0.26, 727.80 \ 0.29, 665.53 \ 0.31, 588.58 \ 0.36, 539.54 \ 0.37, 509.89 \ 0.40, 467.32 \ 0.44, 427.96 \ 0.47, 391.90 \ 0.51, 359.22 \ 0.55, 328.10 \ 0.60, 298.86 \ 0.65, 273.12 \ 0.70, 249.04 \ 0.76, 226.91 \ 0.82, 206.47 \ 0.88, 186.16 \ 0.97, 167.86 \ 1.05, 153.03 \ 1.12, 138.10 \ 1.23, 123.62 \ 1.34, 110.81 \ 1.46, 98.20 \ 1.61, 87.50 \ 1.75, 76.72 \ 1.96, 66.49 \ 2.17, 57.81 \ 2.42, 50.90 \ 2.62, 43.90 \ 2.99, 37.42 \ 3.31, 30.86 \ 3.93, 24.94 \ 4.49, 20.30 \ 5.24, 15.91 \ 6.28, 11.20 &data curves p_spectrum w_p \ 1061.00, 7079.67 \ 353.00, 5011.39 \ 212.00, 3938.40 \ 151.00, 3254.06 \ 117.00, 2769.01 \ 95.00, 2381.27 \ 78.00, 2066.94 \ 66.00, 1819.84 \ 57.00, 1619.22 \ 50.00, 1446.77 \ 44.00, 1294.72 \ 39.00, 1165.95 \ 35.00, 1046.27 \ 31.00, 940.28 \ 28.00, 864.78 \ 26.10, 797.16 \ 23.80, 727.80 \ 21.80, 665.53 \ 20.00, 588.58 \ 17.30, 539.54 \ 16.90, 509.89 \ 15.60, 467.32 \ 14.40, 427.96 \ 13.30, 391.90 \ 12.30, 359.22 \ 11.40, 328.10 \ 10.50, 298.86 \ 9.70, 273.12 \ 9.00, 249.04 \ 8.30, 226.91 \ 7.70, 206.47 \ 7.10, 186.16 \ 6.50, 167.86 \ 6.00, 153.03 \ 5.60, 138.10 \ 5.10, 123.62 \ 4.70, 110.81 \ 4.30, 98.20 \ 3.90, 87.50 \ 3.60, 76.72 \ 3.20, 66.49 \ 2.90, 57.81 \ 2.60, 50.90 \ 2.40, 43.90 \ 2.10, 37.42 \ 1.90, 30.86 \ 1.60, 24.94 \ 1.40, 20.30 \ 1.20, 15.91 \ 1.00, 11.20 $ $********************************************* wind periods $ &set pw = 1061 353 212 151 117 95 78 66 57 50 44 39 35 31 28 \ 26.1 23.8 21.8 20 17.3 16.9 15.6 14.4 13.3 12.3 11.4 \ 10.5 9.7 9.0 8.3 7.7 7.1 6.5 6.0 5.6 5.1 4.7 4.3 3.9 \ 3.6 3.2 2.9 2.6 2.4 2.1 1.9 1.6 1.4 1.2 1.0 $ $********************************************* define env $ &data env &env e90 -sea issc 90 20 10 \ -time 1200 .5 -t_rei 20000 \ -wind 78.2 90 -w_spec api -w_period %pw &env i90f -sea issc 90 20 10 \ -time 1200 .5 -t_rei 20000 \ -wind 78.2 90 -w_spec w_f -w_period %pw &env i90p -sea issc 90 20 10 \ -time 1200 .5 -t_rei 20000 \ -wind 78.2 90 -w_spec w_p -w_period %pw end $ $********************************************* exit &data $ $ $********************************************* marine growth $ &data curves m_growth addit 0 100e3 500 100e3 1000 0 $ $********************************************* current profile $ &data curves c_profile loopcur 0.0 6.10 \ 300.0 6.10 \ 500.0 5.30 \ 1000.0 2.00 \ 2000.0 0.50 \ 6000.0 0.50 $ $********************************************* swell $ &set hs = 10 &set tp = 8 &set c1 = &set c2 = &set string = &data curves p_spectrum swell &set twopi = &number(real 2*3.1415927) &set stwopi = &number(sqrt %twopi) &set wp = &number(real %twopi/%tp ) &loop t 30 3 -.5 &set w = &number(real %twopi/%t) &set e = &number(real ((%w/%wp)-1)**2 ) &set e = &number(exp -50*%e ) &set d = &number(real .1*%wp*%stwopi) &set c = &number(real (%hs/4.)**2) &set s = &number(real (%c*%e)/%d) &set string = %string %t %s &set c1 = %c1 %t &set c2 = %c2 %s &endloop %string $ $********************************************* p_input $ &data convolution period p_input \ 25.00 37.55 0.00 217.43 0.97 3640.28 861.65 \ 20.00 37.55 0.00 227.18 2.93 3084.19 967.34 \ 19.00 37.55 0.00 230.02 3.76 2965.82 989.42 \ 18.00 37.55 0.00 233.27 4.90 2845.29 1011.40 \ 17.00 37.55 0.00 236.95 6.47 2722.82 1032.95 \ 16.00 37.55 0.00 241.09 8.66 2598.72 1053.62 \ 15.00 37.55 0.00 245.65 11.78 2473.44 1072.74 \ 14.50 37.55 0.00 248.05 13.81 2410.55 1081.45 \ 14.00 37.55 0.00 250.49 16.25 2347.62 1089.40 \ 13.50 37.55 0.00 252.90 19.20 2284.79 1096.40 \ 13.00 37.55 0.00 255.21 22.77 2222.21 1102.23 \ 12.50 37.55 0.00 257.29 27.09 2160.05 1106.65 \ 12.00 37.55 0.00 258.99 32.32 2098.55 1109.31 \ 11.50 37.55 0.00 260.07 38.64 2037.96 1109.83 \ 11.00 37.55 0.00 260.25 46.27 1978.64 1107.74 \ 10.50 37.55 0.00 259.13 55.40 1921.00 1102.45 \ 10.00 37.55 0.00 256.24 66.24 1865.57 1093.24 \ 9.50 37.55 0.00 251.05 78.92 1813.01 1079.22 \ 9.00 37.55 0.00 242.98 93.47 1764.15 1059.32 \ 8.50 37.55 0.00 231.51 109.74 1720.09 1032.21 \ 8.00 37.55 0.00 216.27 127.34 1682.20 996.30 \ 7.50 37.55 0.00 197.21 145.60 1652.26 949.72 \ 7.00 37.55 0.00 174.67 163.58 1632.49 890.36 \ 6.50 37.55 0.00 149.42 180.15 1625.68 816.05 \ 6.00 37.55 0.00 122.62 194.07 1635.10 724.91 \ 5.50 37.55 0.00 95.67 204.15 1664.28 616.15 \ 5.00 37.55 0.00 70.02 209.19 1716.43 491.28 \ 4.50 37.55 0.00 47.12 208.01 1793.42 355.64 \ 4.00 37.55 0.00 28.37 199.38 1896.63 218.47 \ 3.00 37.55 0.00 7.10 153.64 2123.13 79.54 \ -ndof 3 $ $********************************************* t_input $ &data convolution time t_input \ 0.25 0.00 -77.36 -194.12 \ 0.50 0.00 -60.42 -173.50 \ 0.75 0.00 -38.45 -145.45 \ 1.00 0.00 -14.85 -112.94 \ 1.25 0.00 6.92 -79.16 \ 1.50 0.00 23.94 -47.10 \ 1.75 0.00 34.31 -19.15 \ 2.00 0.00 37.45 3.19 \ 2.25 0.00 34.04 19.34 \ 2.50 0.00 25.86 29.63 \ 2.75 0.00 15.28 35.02 \ 3.00 0.00 4.76 36.88 \ 3.25 0.00 -3.63 36.62 \ 3.50 0.00 -8.57 35.44 \ 3.75 0.00 -9.72 34.16 \ 4.00 0.00 -7.59 33.17 \ 4.25 0.00 -3.39 32.46 \ 4.50 0.00 1.38 31.77 \ 4.75 0.00 5.30 30.73 \ 5.00 0.00 7.36 29.05 \ 5.25 0.00 7.12 26.64 \ 5.50 0.00 4.82 23.63 \ 5.75 0.00 1.22 20.34 \ 6.00 0.00 -2.66 17.21 \ 6.25 0.00 -5.79 14.67 \ 6.50 0.00 -7.38 13.02 \ 6.75 0.00 -7.10 12.35 \ 7.00 0.00 -5.13 12.56 \ 7.25 0.00 -2.03 13.32 \ 7.50 0.00 1.35 14.21 \ 7.75 0.00 4.18 14.81 \ 8.00 0.00 5.80 14.83 \ 8.25 0.00 5.90 14.11 \ 8.50 0.00 4.58 12.74 \ 8.75 0.00 2.31 10.96 \ 9.00 0.00 -0.25 9.10 \ 9.25 0.00 -2.40 7.54 \ 9.50 0.00 -3.57 6.54 \ 9.75 0.00 -3.50 6.24 \ 10.00 0.00 -2.25 6.58 \ 10.25 0.00 -0.20 7.36 \ 10.50 0.00 2.07 8.27 \ 10.75 0.00 3.98 9.00 \ 11.00 0.00 5.03 9.29 \ 11.25 0.00 4.99 9.02 \ 11.50 0.00 3.92 8.21 \ 11.75 0.00 2.14 7.03 \ 12.00 0.00 0.13 5.74 \ 12.25 0.00 -1.56 4.62 \ 12.50 0.00 -2.51 3.90 \ 12.75 0.00 -2.51 3.68 \ 13.00 0.00 -1.60 3.95 \ 13.25 0.00 -0.05 4.56 \ 13.50 0.00 1.70 5.30 \ 13.75 0.00 3.18 5.92 \ 14.00 0.00 4.00 6.22 \ 14.25 0.00 3.98 6.09 \ 14.50 0.00 3.14 5.54 \ 14.75 0.00 1.73 4.70 \ 15.00 0.00 0.14 3.77 \ 15.25 0.00 -1.21 2.96 \ 15.50 0.00 -1.96 2.46 \ 15.75 0.00 -1.95 2.34 \ 16.00 0.00 -1.20 2.61 \ 16.25 0.00 0.08 3.16 \ 16.50 0.00 1.52 3.79 \ 16.75 0.00 2.74 4.31 \ 17.00 0.00 3.43 4.57 \ 17.25 0.00 3.42 4.47 \ 17.50 0.00 2.73 4.03 \ 17.75 0.00 1.57 3.35 \ 18.00 0.00 0.25 2.59 \ 18.25 0.00 -0.87 1.94 \ 18.50 0.00 -1.50 1.55 \ 18.75 0.00 -1.49 1.49 \ 19.00 0.00 -0.86 1.76 \ 19.25 0.00 0.22 2.26 \ 19.50 0.00 1.44 2.84 \ 19.75 0.00 2.48 3.32 \ 20.00 0.00 3.06 3.57 \ 20.25 0.00 3.06 3.52 \ 20.50 0.00 2.47 3.16 \ 20.75 0.00 1.47 2.59 \ 21.00 0.00 0.34 1.95 \ 21.25 0.00 -0.62 1.40 \ 21.50 0.00 -1.17 1.07 \ 21.75 0.00 -1.16 1.03 \ 22.00 0.00 -0.62 1.27 \ 22.25 0.00 0.31 1.71 \ 22.50 0.00 1.37 2.22 \ 22.75 0.00 2.27 2.66 \ 23.00 0.00 2.78 2.90 \ 23.25 0.00 2.77 2.87 \ 23.50 0.00 2.26 2.59 \ 23.75 0.00 1.39 2.11 \ 24.00 0.00 0.40 1.57 \ 24.25 0.00 -0.44 1.11 \ 24.50 0.00 -0.92 0.83 \ 24.75 0.00 -0.92 0.81 \ 25.00 0.00 -0.44 1.03 \ 25.25 0.00 0.38 1.43 \ 25.75 0.00 2.11 2.28 \ 26.00 0.00 2.56 2.50 \ 26.25 0.00 2.56 2.48 \ 26.50 0.00 2.10 2.23 \ 26.75 0.00 1.33 1.81 \ 27.00 0.00 0.45 1.34 \ 27.25 0.00 -0.30 0.94 \ 27.50 0.00 -0.73 0.70 \ 27.75 0.00 -0.72 0.69 \ 28.00 0.00 -0.30 0.90 \ 28.25 0.00 0.44 1.26 \ 28.50 0.00 1.27 1.69 \ 28.75 0.00 1.98 2.04 \ 29.00 0.00 2.39 2.25 \ 29.25 0.00 2.38 2.23 \ 29.50 0.00 1.98 2.02 \ 29.75 0.00 1.28 1.65 \ 30.00 0.00 0.49 1.23 \ 30.25 0.00 -0.19 0.87 \ 30.50 0.00 -0.57 0.66 \ 30.75 0.00 -0.57 0.65 \ 31.00 0.00 -0.18 0.84 \ 31.25 0.00 0.48 1.17 \ 31.50 0.00 1.23 1.56 \ 31.75 0.00 1.88 1.88 \ 32.00 0.00 2.25 2.07 \ 32.25 0.00 2.24 2.06 \ 32.50 0.00 1.88 1.87 \ 32.75 0.00 1.25 1.54 \ 33.00 0.00 0.53 1.16 \ 33.25 0.00 -0.09 0.83 \ 33.50 0.00 -0.44 0.64 \ 33.75 0.00 -0.44 0.63 \ 34.00 0.00 -0.09 0.81 \ 34.25 0.00 0.52 1.11 \ 34.50 0.00 1.20 1.47 \ 34.75 0.00 1.79 1.77 \ 35.00 0.00 2.13 1.94 \ 35.25 0.00 2.13 1.93 \ 35.50 0.00 1.79 1.75 \ 35.75 0.00 1.21 1.45 \ 36.00 0.00 0.55 1.10 \ 36.25 0.00 -0.01 0.80 \ 36.50 0.00 -0.33 0.63 \ 36.75 0.00 -0.33 0.62 \ 37.00 0.00 -0.01 0.78 \ 37.25 0.00 0.55 1.06 \ 37.50 0.00 1.18 1.38 \ 37.75 0.00 1.72 1.66 \ 38.00 0.00 2.03 1.82 \ 38.25 0.00 2.03 1.81 \ 38.50 0.00 1.72 1.65 \ 38.75 0.00 1.19 1.37 \ 39.00 0.00 0.58 1.05 \ 39.25 0.00 0.05 0.77 \ 39.50 0.00 -0.24 0.61 \ 39.75 0.00 -0.24 0.60 \ 40.00 0.00 0.85 0.76 \ -ndof 3