$
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$
$
$   ----------    Finding the Maximum for Many Cases  --------
$
$  The test shows how to compute the maximum G load on cargo located
$  at several points for several periods and headings.
$
$
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$
$
$*********************************************      SET BASIC PARAMETERS
$
&dimen -dimen feet kips
&device   -g_default device
&TITLE Sample of Strip-Theory Frequency Domain
$
$*********************************************      Read Model
$
inmodel
$
$*********************************************      Set Condition
$
&instate -con 10 0 0000
&weight -compute 27  32 97 97
&status
$
$*********************************************      Frequency Domain
$
hydro
   g_press
end
$
$*********************************************      Define Report Pts.
$
medit
   *bowp    00 -50 25
   *bows    00  50 25
   *center 200   0 25
   *cwl    200   0 10
   *sternp 200 -50 25
   *sterns 200  50 25
end
&set points = *bowp *bows *center *cwl *sternp *sterns
$
$*********************************************      initialize
$
&set env      =
&set max_x    =
&set max_y    =
&set max_z    =
$
$*********************************************      Define env
$
&set Weight   = 1000
&set hs       = 10
&set periods  = 07 08 09  10  11  12
&set headings = 000 045 090 135 180 225 270 315
&data env
   &loop p ( %periods )
      &loop h ( %headings )
          &set e   = e|%p|%h
          &set env = %env %e
          env %e -sea issc %h %hs %p
      &endloop
   &endloop
end
$
$*********************************************      Frequency Post
$
freq_resp
$
$*********************************************      response
$
   rao
$
$*********************************************      loop points
$
   &loop p ( %points )
      fr_point *bowp
         end
      fr_fcargo %weight
         end
$
$*********************************************      loop envs
$
      &loop e ( %env )
$
$*********************************************      get statistics
$
         st_fcargo %e
            report
$           vlist
$
$*********************************************      get max
$
            set_var x  -column 2
            set_var y  -column 3
            set_var z  -column 4
            &set max_x = &number(max %max_x %x)
            &set max_y = &number(max %max_y %y)
            &set max_z = &number(max %max_z %z)
$
$*********************************************      end dispose
$
            end
$
$*********************************************      end loops
$
      &endloop
   &endloop
$
$*********************************************      report
$
&type
&type Maximum Cargo Force Statistics
&type     Longitudinal &format(f8.4 %max_x/%weight )
&type     Transverse   &format(f8.4 %max_y/%weight )
&type     Vertical     &format(f8.4 %max_z/%weight )
$
$*********************************************      All done
$
&finish