oceanlyz.SeaSwellFun#
Hm0,Hm0sea,Hm0swell,Tp,Tpsea,Tpswell,fp,fseparation,f,Syy=oceanlyz.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')