clear vari axis_inc=5 `degrees vari axis_angle=0 `degrees vari n set n=360 DIVIDE {axis_inc} PLUS 1 macro vary_axis `1=windf, 2=hmmt, 3=heel, 4=gmtra set axis_angle={axis_angle} PLUS {axis_inc} axis {axis_angle} hmmt wind /band /quiet if %1=1 then me {axis_angle} {windf} if %1=2 then me {axis_angle} {hmmt} if %1=3 then so | me {axis_angle} {heel} if %1=4 then so | me {axis_angle} {gmtra} / macro polar_plot `1=windf, 2=hmmt, 3=heel, 4=gmtra if %1=1 then vari plot_title="Force, Wind Speed: {wind:1} knots" if %1=2 then vari plot_title="Heeling Moment, Wind Speed: {wind:1} knots" if %1=3 then vari plot_title="Heel Angle, Wind Speed: {wind:1} knots" if %1=4 then vari plot_title="Metacentric Height, GMa" me plotstart "{plot_title}" /polar:-1 me plotstyle "L" .vary_axis({n}) %1 me plotend / read rig.gf tr 0 he 0 dr 24 vcg 15 so we lcg tcg wind (SP) 30 report RIG_POLAR.PF `.polar_plot 1 `wind force plot .polar_plot 2 `wind hmmt plot `.polar_plot 3 `heel angle plot .polar_plot 4 `GM plot report /preview /thick:5,5