SeaSwellFun

[Hm0,Hm0sea,Hm0swell,Tp,Tpsea,Tpswell,fp,fseparation,f,Syy]=SeaSwellFun(input,fs,duration,nfft,h,fmin,fmax,ftailcorrection,tailpower,fminswell,fmaxswell,mincutoff,maxcutoff,tailcorrection,dispout)

DESCRIPTION

Separate sea wave from swell wave

INPUT

input=importdata(‘h.mat’);

Load water depth (h)/surface elevation (Eta) data and rename it “input” in (m)

fs=10;

Sampling frequency that data collected at in (Hz)

duration=1024;

Duration time that data collected in input in each burst in second

nfft=2^10;

NFFT for Fast Fourier Transform

h=1;

Mean water depth in (m)

fmin=0.04;

Minimum frequency for cut off the lower part of spectra

fmax=1;

Maximum frequency for cut off the upper part of spectra

ftailcorrection=1;

Frequency that diagnostic tail apply after that (typically set at 2.5fm, fm=1/Tm01)

tailpower=-5;

Power that diagnostic tail apply based on that (-3 for shallow water to -5 for deep water)

fminswell=0.1;

Minimum frequency that is used for Tpswell calculation

fmaxswell=0.25;

Maximum frequency that swell can have, It is about 0.2 in Gulf of Mexico

mincutoff=’off’;

Define if to cut off the spectra below fmin

mincutoff=’off’: Cutoff off

mincutoff=’on’: Cutoff on

maxcutoff=’off’;

Define if to cut off the spectra beyond fmax

maxcutoff=’off’: Cutoff off

maxcutoff=’on’: Cutoff on

tailcorrection=’off’;

Define if to apply diagnostic tail correction or not

tailcorrection=’off’: Not apply

tailcorrection=’jonswap’: JONSWAP Spectrum tail

tailcorrection=’tma’: TMA Spectrum tail

dispout=’on’;

Define to display outputs or not (‘off’: not display, ‘on’: display)

OUTPUT

Hm0

Zero-Moment Wave Height (m)

Hm0sea

Sea Zero-Moment Wave Height (m)

Hm0swell

Swell Zero-Moment Wave Height (m)

Tp

Peak wave period (second)

Tpsea

Peak Sea period (second)

Tpswell

Peak Swell Period (second)

fp

Peak Wave Frequency (Hz)

f

Frequency (Hz)

fseparation

Sea and Swell Separation Frequency (Hz)

Syy

Wave Surface Elevation Power Spectrum (m^2s)

EXAMPLE

[Hm0,Hm0sea,Hm0swell,Tp,Tpsea,Tpswell,fp,fseparation,f,Syy]=SeaSwellFun(water_pressure/(1000*9.81),10,1024,256,1.07,0.05,5,1,-5,0.1,0.25,'on','on','off','on')