close all; clear all

Res=10:500:10010;

OMEGAISP=[];
OMEGAISML=[];
N= 180;

for ii=1:length(Res)
    ii
    Re=Res(ii);
    [omegaisP]=OrrSommerfeld_Poiseuille_temporal(N,Re);
    [omegaisML]=OrrSommerfeld_tanh_temporal(N,Re);
    OMEGAISP=[OMEGAISP ; omegaisP];
    OMEGAISML=[OMEGAISML ; omegaisML];
    clear omegaisP omegaisML Re
    
end

alphas=0.1:0.01:1.5;
v=-0.000005:0.000005:0.000005%::0.01;

figure(1)
contour(Res,alphas,OMEGAISP.',v)
figure(2)
contour(Res,alphas,OMEGAISML.',v)