$
$*@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$*@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$*@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$
$ ------- Strip Theory for a Barge ---------------------
$ --------    405 X 90 X 20 FT BARGE    --------------
$ ---- To be used in compartison to Noble Denton curves
$
$*@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$*@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$*@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$
$******************************* SET BASIC PARAMETERS
$
&device -g_default device  -oecho no
&dimen -dimen feet kips
&title Comparison of MOSES Strip Theory w/Model Tests
&set subt = 405x90x20 ft Barge, 9.0ft draft
$
$*************************************      READ IN MODEL
$
inmodel
$
$*************************************      set condition
$
&instate -cond 9
&weight -total 19946 212.8 0 32 51.9 123 123
&status
$
$***********************************      plot the model
$
&channel  gra -file ndai.ans/ndai_iso.eps
&picture iso
$
$***********************************      frequency domain
$
hydro
    g_press -period %period -heading 180 90
end
$
$**********************************     frequency
$
freq_resp
   rao -steep 20
   fr_point 212 0 32
      report
      &channel  gra -file ndai.ans/ndai_surge.eps
      plot 2 15 -no -t_left 'Surge, ft/ft' -t_x 'Period, Sec' \
                -t_sub  '%subt%, Head Seas' -legend 1 Moses \
                -add 4               \
                -points "Model Test" \
                       7.0 .05       \
                       9.8 .17       \
                      11.6 .46       \
                      16.0 .83       \
                      19.5 .92
      &channel  gra -file ndai.ans/ndai_heave_180.eps
      plot 2 19 -no -t_left 'Heave, ft/ft' -t_x 'Period, Sec' \
                -t_sub  '%subt%, Head Seas' -legend 1 Moses \
                -add 4               \
                -points "Model Test" \
                       6.0 .07       \
                       8.0 .13       \
                       9.7 .21       \
                      11.8 .52       \
                      16.0 .86       \
                      19.8 .95
      &channel  gra -file ndai.ans/ndai_pitch.eps
      plot 2 23 -no -t_left 'Pitch deg/ft' -t_x 'Period, Sec' \
                -t_sub  '%subt%, Head Seas' -legend 1 Moses \
                -add 4               \
                -points "Model Test" \
                       6.0  .30/3.28  \
                       8.0  .52/3.28  \
                       9.8 1.25/3.28  \
                      11.7 1.29/3.28  \
                      16.0  .90/3.28  \
                      19.5  .65/3.28
      &channel  gra -file ndai.ans/ndai_sway.eps
      Range -left 0 1.5 -x 5 20
      plot 2 5  -no -t_left 'Sway, ft/ft' -t_x 'Period, Sec' \
                -t_sub  '%subt%, Beam Seas' -legend 1 Moses \
                -add 4               \
                -points "Model Test" \
                       6.0 .32       \
                       8.0 .56       \
                       9.6 .95       \
                      11.6 .875      \
                      17.0 .85       \
		      19.2 .93
      &channel  gra -file ndai.ans/ndai_heave_90.eps
      plot 2 7  -no -t_left 'Heave, ft/ft' -t_x 'Period, Sec' \
                -t_sub  '%subt%, Beam Seas' -legend 1 Moses \
                -add 4               \
                -points "Model Test" \
                       6.0 .675      \
                       8.0 .92       \
                       9.8 .78       \
                      11.8 .98       \
                      16.0 1.        \
                      20.0 1.
   end
end
$
$**********************************     different steepness
$
&subtitle
&loop I (33 52 75 133)
   &set roll%I% =
   freq_resp
      rao -steep %I% -heading 180 90
      fr_point 212 0 32
         &loop J (%period%)
           set_var roll -select 9 2 %J%
           &set    roll%I% = %(roll%(I)) %roll%
         &endloop
         set_var maxp -max 9 2
         &set roll%I% = %(roll%(I)) %maxp%
         report
      end
   end
&endloop
$
$**********************************     build graph
$
&channel  gra -file ndai.ans/ndai_roll.eps
&buildg -brief
   Period
   Moses Roll S33
   Moses Roll S52
   Moses Roll S75
   Moses Roll S133
 
&loop I 1 &token(n %period%)
   &token(%I% %period%) &token(%I% %roll33%) &token(%I% %roll52%)  \
   &token(%I% %roll75%)  &token(%I% %roll133%)
&endloop
 
 
      Range -left 0 8 -x 5 20
      plot 1 2 3 4 5  -no -t_left 'Roll deg/ft' -t_x 'Period, Sec' \
                -t_sub  '%subt%, Beam Seas' \
                -add 4 -clean               \
                -points "Model Test" \
                       6.0  .70/3.28  \
                       8.0 2.50/3.28  \
                      10.0 7.10/3.28  \
                      11.7 4.89/3.28  \
                      16.0 1.30/3.28
end
$
$**********************************     all done
$
&fini