clear all
clear

macro condition
  dr 8
  lcg 0
  tcg 0
  vcg 5
  so we
  we * /gyr:gf,gf,gf
/

macro scale_beam
  read fv.gf
  enter pm
    scale 1,%1,1
    write fv.gf1
  qu `pm
/

macro vary_beam
  b:=b+%2
  .scale_beam {b}
  .condition
  wave (spe) p2 3.5 `waves constant
  sea /speed:6 /head:135 /data:st /brief /noplot /samp:5
  .sk.getstatsdata "SK-STATS.DAT"
  vari RMS4%1={sk.RMS4_CG_P} `save RMS roll angle
  vari RMS5%1={sk.RMS5_CG_P} `save RMS roll angle
  vari B%1={woa}
/

run sk.lib /call /quiet

vari b=0.6 `starting beam scale
.vary_beam(10,1)1 0.1 `increment beam by 10%

macro plot
  me {B%1},{RMS4%1},{RMS5%1}
/

report RMSRollvsBeam.PF
  me plotstart "RMS Roll and Pitch with Variation in Beam" /norotate /page%:40
  me plotlabel "Beam (FT)|RMS (Deg)","Roll","Pitch"
  .plot(10,1)1
  me plotend
report close /thick:3,5 /preview