Documentation of wave_speed

Global Index (short | long) | Local contents | Local Index (short | long)

Help text

  Kelvin wave

Cross-Reference Information

Listing of script wave_speed

global RADUS
yr = 7:.2:33;
k = 2*pi/(RADUS*2*pi/3);
N = .1;
  beta = 2*7.292e-5/RADUS;
  omega = 2*pi./(3600*24*365*yr);
  kstar = k .* omega/beta;
n = 1;
m = beta*N./omega.^2 .* ( n+0.5 - ( (kstar+0.5).^2 + n*(n+1) ).^0.5 );
cp = sqrt((omega./m).^2 + (omega./k).^2);
cg = diff(omega)./diff(m);
cgom = omega(1:(length(omega)-1))+0.5*diff(omega);
plot(2*pi./(3600*24*365*cgom), 3600*24*365*cg, '-k')
hold on;
  plot(yr, 3600*24*365*omega./m, '--k');
hold off;
axis([5 35 -2 2]);
set(gca, 'YTick', -2:.5:2);
grid on;
ylabel('Vertical Wave Speed (m/yr)');
plot(yr, omega./m);
xlabel('Forcing period (yr)');
cg = 2e-3;
(1/1000)^3*((2*cg)/(1.5-sqrt(2)))^2/beta
1e-6 m/s
300/(2*3600*24*365)

clear

global RADUS
beta = 2*7.292e-5/RADUS;
N = 15;

%k = 2*pi/(RADUS*2*pi/2);
k = 2*pi/(30000000);
%yr = [.05:.05:1];
yr = 14/365;
omeg = 2*pi./(3600*24*365*yr);
kstar = k.*omeg./beta;
mstar = kstar;
%m = N*k./omeg;
m = mstar*beta*N./(omeg.^2);

cpz = omeg ./ m;
cgz = -1*N*k./(m.^2);
cpx = omeg ./ k;

figure(1); fl(1);
plot(yr, cpz, '-k', yr, cgz, '--k');