From 466d13f302e1b6a84e8ad096586b1c97fc1c3708 Mon Sep 17 00:00:00 2001 From: Jean-Denis Vauguet Date: Mon, 1 Jun 2009 00:03:02 +0200 Subject: [PATCH] plop --- .../IPGP Saint Maur/MATLAB/vectorial/modesT_vect.m | 3 ++- .../IPGP Saint Maur/MATLAB/vectorial/wind_effect.m | 23 +++++++++++++++------- 2 files changed, 18 insertions(+), 8 deletions(-) diff --git a/docs/fac/stage/IPGP Saint Maur/MATLAB/vectorial/modesT_vect.m b/docs/fac/stage/IPGP Saint Maur/MATLAB/vectorial/modesT_vect.m index 13caef6..330952f 100644 --- a/docs/fac/stage/IPGP Saint Maur/MATLAB/vectorial/modesT_vect.m +++ b/docs/fac/stage/IPGP Saint Maur/MATLAB/vectorial/modesT_vect.m @@ -6,6 +6,7 @@ % c'est du 3D planaire disons % TODO ajouter un flag pour savoir si on calcule avec des L choisis ou bien des L dérivés des T via la vitesse de phase sismique +% TODO éventuellement, fichier ou vecteur de conf pour les arguments ? function modesT_vect(p_u0b, p_u0t, p_T, p_L, p_res, plotflags, exportflags) % arguments: @@ -74,7 +75,7 @@ do_export_append = exportflags(2); load modele1 % tests à résolution variable -res = 1 % resolution +res = p_res % resolution binf = 1 bsup = 9999 diff --git a/docs/fac/stage/IPGP Saint Maur/MATLAB/vectorial/wind_effect.m b/docs/fac/stage/IPGP Saint Maur/MATLAB/vectorial/wind_effect.m index 7f03785..eba8b71 100644 --- a/docs/fac/stage/IPGP Saint Maur/MATLAB/vectorial/wind_effect.m +++ b/docs/fac/stage/IPGP Saint Maur/MATLAB/vectorial/wind_effect.m @@ -95,14 +95,23 @@ figure(100); % mais il faut réfléchir à comment ploter ça… % idées : % faire les deux, ie. pour chaque L, le max en fonction de T, puis pour chaque T, le max en fonction de L, -% et essayer un pcolor pour avoir le max pour tous les couples L/T (mais ça doit plus ou moins revenir au même que le pcolor du script principal) +% et essayer un pcolor pour avoir le max pour tous les couples L/T (mais ça doit plus ou moins revenir au même que le pcolor du script principal : ses pcolor plottent l'amplification relative par rapport au sol, ici ce serait, intégré, par rapport à un cas de référence à vent nul, donc c'est pas exactement la même chose quand même) % FIXME ne pas faire de abs() en fait, sinon on ne voit pas si ça amplifie ou atténue ! -amp_data = squeeze(max(abs(data(:,:,1,:)))) -amp_ref = squeeze(max(abs(wind_ref(:,:,1,:)))) -plot(T, amp_data, 'r'); -plot(T, amp_ref, '--b'); -title('maximum w amplitude for each period'); -hold off; +% ne pas oublier de plotter le cas de référence à vent nul, d'ailleurs +% un truc qui serait chanmé : un pcolor en couche, avec une croix indiquant la position du cas de référence, une autre indiquant la position du cas étudié là, le vecteur vent moyen, en gros. À réfléchir… +% pour le moment, sur les L… +for j = 1:length(L) + for k = 1:length(T) + amp_data = data(:,:,j,k) % on récupère w, u, p et r pour les L et T courants + amp_data = amp_data(:,1); % on récupère le vecteur des w en particulier + amp_data = squeeze(max(amp_data)) + + plot(T, amp_data, 'r'); + plot(T, amp_ref, '--b'); + title('maximum w amplitude for each period'); + hold off; + end +end % }}} -- 2.11.4.GIT