windSpeed_kts = 29.93
windSpeed_mps = windSpeed_kts * 1852.0 / 3600

g = 9.8;
U195 = windSpeed_mps;
T1 = 10;

f =  0.02:.001:0.5;
w = 2 * pi * f;
alpha = 8.1e-3;
beta = 0.74;
w0 = g / U195;


PM_Sw = (2*pi)*((alpha * g^2) ./ (w.^5)) .* exp((-1*beta)*((w0./w).^4));
modPM_Sw = (0.11/(2*pi))*(((w.*T1)/(2*pi)).^-5).*exp(-.44*((w.*T1)./(2*pi)).^-4);

% figure
% plot(f,PM_Sw,f,modPM_Sw);

%figure
%plot(f(30:end),PM_Sw(30:end),f(30:end),modPM_Sw(30:end));