clear all

read sv.gf

he 0
tr 0
dr 4.2
vcg 4
so we lcg tcg
weight * /gyr: gf,gf,gf

crt (1) "bow" 61.949f,0,10.264

wave (spe) 6,8 `amplitude, period

report
 sea /crt:1 /rel /data:wa,re /resp:p
report off

`---- Use SEA.LIB ----
run sk.lib /call /quiet
set sk.nwaves=1 `number of wave components
set sk.ncrtpts=1 `number of Critical Points, not including CG
set sk.relative=1 `parse for relative motions
.sk.GetRespData "SK-RESPONSE.DAT"
.sk.GetWaveData "SK-WAVES.DAT"

`---- Make Time-Domain Visualization ----
vari nsteps,AbsPos,RelPos
vari dt=0.25 `time step in seconds
vari T=16 `total time in seconds

set nsteps={T} DIVIDE {dt}

macro TimeDomain 
 `%2 is wave component index 
 .sk.WaveElevation %2 -61.949 0 {t}
 .sk.Position %2 3 1 {t} 0 1 `%2 wave component, 3rd mode (heave), crt 1, additional phase shift, point location shift on
 set AbsPos={sk.Pos}
 .sk.Position %2 R 1 {t} 0 1 `%2 wave component, rel mode, crt 1, additional phase shift, point location shift on
 set RelPos={sk.Pos}
 me {t} {AbsPos} {sk.WELEV} {RelPos} 0 `send to plot
 set t={t} PLUS {dt}
/

report vmotion.pf
me plotstart "Vertical Bow Motion in a Sinusoidal Wave"
me plotlabel "Time (Sec)|Vertical Position (m.)","Bow Motion","Wave Elevation","Relative Bow Motion","Calm Water Surface"
vari t=0
.TimeDomain({nsteps},0)1 1
me plotend
report close /thick:5,5 /preview