clear all | clear

vari sw=0

read fv.gf
we 250 0 0 5 /gyr:gf,gf,gf
so

`generate data first
macro gen_rao_data
  report skrep.$$$
    wave (spe) 1,20.944,15.708,12.566,10.472,8.976,7.854,6.981,6.283,5.712,5.236,4.833,4.488,4.189,3.927,3.696,3.491,3.307,3.142,2.992,2.856,2.732,2.618,2.513,2.417,2.327,2.244,2.167,2.094
    sea /brief /head:180 /speed:6 /data:re
  report off 
  erase skrep.$$$ /norecycle
/

.gen_rao_data

`seakeeping macro
macro seak
  sw:=sw+1
  sea data sk*-response.dat /brief /noplot /data:st
  rename skdata\SK-STATS.DAT skdata\SW{sw}-SK-STATS.DAT
  page
/

`use data to test different wave conditions
report waves.pf
  \\Pierson-Moskowitz using Wind Speed\
  \\
  wave (spe) PM 20 `knots
  .seak
  \\Pierson-Moskowitz using Significant Wave Height\
  \\
  wave (spe) P2 3
  .seak  
  \\Bretschneider using Peak Period and Significant Wave Height\
  \\
  wave (spe) BR 7,2
  .seak  
  \\JONSWAP using Wind Speed and Fetch\
  \\
  wave (spe) J1 30,500000 `fetch in current length units
  .seak  
  \\JONSWAP using Peak Period and Significant Wave Height\
  \\
  wave (spe) J2 9,3
  .seak  
  \\Ochi-Hubble Six-Parameter\
  \\
  wave (spe) O6 0.3,2,3.3,1.5,1.5,3.3
  .seak 
  \\Ochi-Hubble One-Parameter (Significant Wave Height)\
  \\
  wave (spe) O1 3
  .seak   
  
  `collect response statistics in a table (could do this efficiently in a macro)
  run sk.lib /call /quiet
  .sk.getstatsdata skdata/SW1-SK-STATS.DAT `P1
  vari SW1stat={sk.RMS3_CG_P}
  .sk.getstatsdata skdata/SW2-SK-STATS.DAT `P2
  vari SW2stat={sk.RMS3_CG_P}
  .sk.getstatsdata skdata/SW3-SK-STATS.DAT `BR
  vari SW3stat={sk.RMS3_CG_P}
  .sk.getstatsdata skdata/SW4-SK-STATS.DAT `J1
  vari SW4stat={sk.RMS3_CG_P}
  .sk.getstatsdata skdata/SW5-SK-STATS.DAT `J2
  vari SW5stat={sk.RMS3_CG_P}
  .sk.getstatsdata skdata/SW6-SK-STATS.DAT `O6
  vari SW6stat={sk.RMS3_CG_P}
  .sk.getstatsdata skdata/SW7-SK-STATS.DAT `O1
  vari SW7stat={sk.RMS3_CG_P}
  
  `make table
  me repstart 
  \RMS for Seaways\
  \Seaway-----------RMS\
  \P1        { SW1stat}\
  \P2        { SW2stat}\
  \BR        { SW3stat}\
  \J1        { SW4stat}\
  \J2        { SW5stat}\
  \O6        { SW6stat}\
  \O1        { SW7stat}\
  me repend
report off