`example run file to compute wave forcing statistics in a seaway
`requires SK.LIB
`requires GHS v17.32 or later

clear all | clear

`read geometry and establish loading condition
read fv
we 250 0 0 5 /gyr:gf,gf,gf
so

`specify seaway
wave (spe) br 10 3

`run seakeeping and generate wave and forcing data file
report 
sea /head:90 /data:wa,fo
report off

`use SK.LIB to parse data files
run sk.lib /call /quiet 
.sk.getforcingdata "sk-forcing.dat"
.sk.getwavedata "sk-waves.dat"

macro superimpose `%1 is mode
`superimpose force/moment "RAO" with wave spectrum
 macro superimpose0 `%1 is wave index
  .sk.makevari "wfs%1_%%91"
  sk.wfs%1_%%91:=sk.Ord_%%91*sk.FTAmp%1_%%91^2   
 //
 .superimpose0({sk.nwaves},1)1
/

macro integrate_spectrum `%1 is mode
`integrate excitation spectrum to compute variance
.sk.makevari "spwf"
.sk.makevari "r1" | .sk.makevari "r2"
.sk.makevari "dw"
 macro integrate_spectrum1 `%1 is moment 0,1,2
  macro integrate_spectrum0 `%1 is wave index
   sk.r1:=sk.EncFreq_{%%%991+1}^%%91*sk.wfs%1_{%%%991+1} 
   sk.r2:=sk.EncFreq_%%%991^%%91*sk.wfs%1_%%%991
   sk.spwf:=sk.r1+sk.r2
   sk.dw:=sk.Freq_{%%%991+1}-sk.Freq_{%%%991}
   sk.r1:=0.5*sk.dw*sk.spwf
   sk.m%%91_%1:=sk.m%%91_%1+sk.r1
  ///
  .sk.makevari "m%%91_%1"
  sk.m%%91_%1:=0
  .integrate_spectrum0({{sk.nwaves}-1},1)1
 //
 .integrate_spectrum1 0 `m0
/

macro compute_excitation_variance `%1 is mode
 .superimpose %1
 .integrate_spectrum %1
/

.compute_excitation_variance(6,1)1 `compute for all modes

vari factor=3.85 `for max1000 statistic, see SK Manual pg. 77
\ Max1000 Surge Force: {{factor}*({sk.m0_1}^0.5)} lbs \
\ Max1000 Sway Force: {{factor}*({sk.m0_2}^0.5)} lbs \
\ Max1000 Heave Force: {{factor}*({sk.m0_3}^0.5)} lbs \
\ Max1000 Roll Moment: {{factor}*({sk.m0_4}^0.5)} lb-ft \
\ Max1000 Pitch Moment: {{factor}*({sk.m0_5}^0.5)} lb-ft \
\ Max1000 Yaw Moment: {{factor}*({sk.m0_6}^0.5)} lb-ft \