clear
proj COW073

vari (string) x1,x2,x3,x4,p,f,r,row,lbls,lind
vari i,j,imax,jmax `i indexes frequency, j indexes heading

macro rdln
  read p,f,r /skip 
/

macro parseRAO
  .rdln
  if {r}={lind} then .rdln(3,1)
  `if {p}={ind} then .rdln(3,1)
  if {p}="RAO" then set imax={i} | exit
  set i={i} plus 1
  vari rao{i}{j}, freq{i}
  set rao{i}{j}={r}, freq{i}={f}
  exit parseRAO
/

macro parsePF
  read x1,x2,x3,x4 /skip
  if {x4}="DEGREES" then set j={j} plus 1 | vari head{j} | set head{j}={x3}
  if {x1}="HEAVE" then if {x2}="RAO" then .rdln | set i=0 | .parseRAO
  if eof then set jmax={j} | exit
  exit parsePF
/

macro columns
  set row="{row} {rao{i}{j}}"
  set j={j} plus 1
/

macro plot
  set row="{freq{i}}", j=1
  .columns({jmax},1)
  me {row}
  set i={i} plus 1
/

macro plot_labels
  set lbls="{lbls}, ""Heading={head{j}} Deg"""
  set j={j} plus 1
/

`--- Parse Report File ---
read (data) HEADINGS.PF
read x1,x2,x3 /skip
set lind={x3} | set j=0 | .parsePF
read (data) off

`--- Create Custom Plot ---
report
me plotstart "Variable Heading Heave RAO"
set lbls="""Wave Frequency (Rad/Sec)|RAO (Ft/Ft)"""
set j=1 | .plot_labels({jmax},1)
me plotlabel {lbls}
set i=1 | .plot({imax},1)
me plotend
report off