%% Paramètres physiques
g= 9.81; gamma= 0.074; rho= 1000; nu= 1.31*10^-6;

%% Nombre d'onde
lambda= 40; k= 2*pi/lambda

%% Domination du terme de gravité / terme de tension sup dans rdd
R= (g/k)/(gamma*k/rho)

%% Onde de gravité parce que k est petit !
c= sqrt(g/k)
omega= c*k
T= 2*pi/omega

%% Domination du terme de gravité / terme visqueux dans la cond. d'interf.
R= nu*k*omega/g

%% Interface : définition et tracé
% Valeurs tabulées de x puis z
x= linspace(0, 40, 30); A= 0.5; z= A * cos(k*x);
% Tracé
plot(x,z,'k','LineWidth',2)

%% Lignes de courant
% Grille = valeurs tabulées de x et z
[x,z]= meshgrid(linspace(0, 40, 30), linspace(-10, 0, 20));
% Valeurs tabulées de psi à un facteur près
psi= exp(k*z) .* cos(k*x);
% Lignes de courant : isocontours de la fonction courant psi
hold on; contour( x, z, psi, 'k')
% Imposer les intervalles voulus et un rapport d'aspect physique
axis([0 40 -10 0.5]); axis equal;

%% Champ de vitesse
% Nouvelle grille moins fine
[x,z]= meshgrid(linspace(0, 40, 18),linspace(-10, -2, 4));
% Valeurs tabulées de la vitesse à un facteur près
vx= exp(k*z) .* cos(k*x); vz= exp(k*z) .* sin(k*x);
% Tracé de ce champ de vecteurs
quiver(x,z,vx,vz)
% On n'oublie pas de labelliser les axes
xlabel('x en m'); ylabel('z en m');
