From 0d17ebd453ac8025b1d7b6b9a6eb2b7e931e0f8b Mon Sep 17 00:00:00 2001 From: Jean-Denis Vauguet Date: Mon, 8 Jun 2009 17:10:40 +0200 Subject: [PATCH] yeah --- .../MATLAB/vectorial/3D/data_processing.m | 6 ++--- .../MATLAB/vectorial/3D/wind_effect.m | 26 ++++++++++++---------- .../IPGP Saint Maur/tex/rapport/presentation.tex | 2 +- .../stage/IPGP Saint Maur/tex/rapport/rapport.tex | 4 ++-- 4 files changed, 20 insertions(+), 18 deletions(-) diff --git a/docs/fac/stage/IPGP Saint Maur/MATLAB/vectorial/3D/data_processing.m b/docs/fac/stage/IPGP Saint Maur/MATLAB/vectorial/3D/data_processing.m index e8c35d2..77598da 100644 --- a/docs/fac/stage/IPGP Saint Maur/MATLAB/vectorial/3D/data_processing.m +++ b/docs/fac/stage/IPGP Saint Maur/MATLAB/vectorial/3D/data_processing.m @@ -5,11 +5,11 @@ function [w,u,p,r] = data_processing(thefile, varargin) % - inferior limit % - resolution % - superior limit +% - 3D flag optargin = size(varargin, 2); stdargin = nargin - optargin; -%pause data = load(thefile, '-ascii'); @@ -22,7 +22,7 @@ if optargin > 4 if varargin{5} == '3D' threeD = 1; else - error('bad value (try 3D or no value)'); + error('last argument has bad value (try 3D or let empty)'); end end @@ -78,7 +78,7 @@ if optargin > 0 bsup = varargin{4}; end else - % on demand… + % on demand for k = 1:optargin % TODO disp 'ERROR: not implemented yet!' diff --git a/docs/fac/stage/IPGP Saint Maur/MATLAB/vectorial/3D/wind_effect.m b/docs/fac/stage/IPGP Saint Maur/MATLAB/vectorial/3D/wind_effect.m index 5c05fe0..973adcf 100644 --- a/docs/fac/stage/IPGP Saint Maur/MATLAB/vectorial/3D/wind_effect.m +++ b/docs/fac/stage/IPGP Saint Maur/MATLAB/vectorial/3D/wind_effect.m @@ -5,8 +5,9 @@ function wind_effect(varargin) % TODO s'assurer de l'existence des dossiers, les creer sinon -% TODO indiquer la res dans les noms de fichiers, ça eviterait de les ouvrir :) regexp ^^ -% TODO supprimer les output graphiques à chaque nouveau run ? +% TODO demander les profils a utiliser (uigetfile) +% TODO indiquer la res dans les noms de fichiers, ca eviterait de les ouvrir :) regexp ^^ +% TODO supprimer les output graphiques a chaque nouveau run ? % {{{ parameters handling @@ -139,7 +140,7 @@ else if exist(name, 'file') fprintf('%s %s %s %s', '"', name, '"', ' exists'); if check_res - fprintf('%s', ', checking its resolution…'); + fprintf('%s', ', checking its resolution...'); do_check = load(name); % force computation if the provided file has wrong resolution if size(do_check(:,1),1) ~= (bsup-binf/res+1) @@ -226,7 +227,7 @@ else if exist(name, 'file') fprintf('%s %s %s %s', '"', name, '"', ' exists'); if check_res - fprintf('%s', ', checking its resolution…'); + fprintf('%s', ', checking its resolution...'); do_check = load(name); % force computation if the provided file has wrong resolution if size(do_check(:,1),1) ~= (bsup-binf/res+1) @@ -301,7 +302,7 @@ fprintf('\n') % {{{ amplitudes processing -% data dumb sample: +% dumb data sample: %data(:,:,1,1) = [1 2; 3 4; 5 6]; %data(:,:,1,2) = [11 22; 33 44; 55 66]; %data(:,:,1,3) = [-111 -222; 333 444; 555 -666]; @@ -314,14 +315,14 @@ disp '> processing amplitudes' amp_max = []; amp_ref = []; if ~auto_winds - n_param = 4 + n_param = 4; else - n_param = 5 + n_param = 5; end for x = 1:n_param - amp_tmp = squeeze(data(:,x,:,:)); - amp_tmp2 = squeeze(wind_ref(:,x,:,:)); + amp_tmp = squeeze(data(:,x,:,:)); + amp_tmp2 = squeeze(wind_ref(:,x,:,:)); % 1. get max with respect to the signs %test = squeeze(min(amp_tmp) +max(amp_tmp)) < 0; % look for the sign @@ -348,8 +349,8 @@ for x = 1:n_param %end % 2. get abs max - amp_max(:,:,x) = squeeze(max(abs(amp_tmp))) - amp_ref(:,:,x) = squeeze(max(abs(amp_tmp2))) + amp_max(:,:,x) = squeeze(max(abs(amp_tmp))); + amp_ref(:,:,x) = squeeze(max(abs(amp_tmp2))); % perform some divergence checks n_global_div_max = find(abs(amp_max(:,:,x)) > threshold3); @@ -378,6 +379,7 @@ for x = 1:n_param amp_max(amp_ref_bad) = 0; amp_ref(amp_max_bad) = 0; + % feedback outputs amp_ref amp_max amp_max-amp_ref @@ -475,7 +477,7 @@ if plot_pcolors set(gca, 'fontsize', font_sizeN); % horizontal speed disgression - % FIXME 3D -> 2D mean… + % FIXME 3D -> 2D mean... subplot(2, 2, 2); if flag_plot == 0 pcolor(wn, kx, (amp_max_u - amp_ref_u)./amp_ref_u.*100); diff --git a/docs/fac/stage/IPGP Saint Maur/tex/rapport/presentation.tex b/docs/fac/stage/IPGP Saint Maur/tex/rapport/presentation.tex index 4cd1cf6..0a7d954 100644 --- a/docs/fac/stage/IPGP Saint Maur/tex/rapport/presentation.tex +++ b/docs/fac/stage/IPGP Saint Maur/tex/rapport/presentation.tex @@ -40,7 +40,7 @@ A & = i \lp \Omega + \frac{g}{\Omega} \dz{\ln{\rho_0}} \rp & -\frac{1}{2} \dz{\ln{\rho_0}} \notag \end{bmatrix} \end{align} -Les $*$ dénotent une renormalisation des deux paramètres $\tuz$ et $\tp$ : $\tuz^* = \sqrt{\rho_0} \tuz$ et $\tp^* = \frac{\tp}{\sqrt{\rho_0}}$ (d'où les termes en $\pm \frac{1}{2} \dz{\ln{\rho_0}}$ de $A$). Cette renormalisation assure que les perturbations de la vitesse verticale et de la pression sont indépendantes de l'amplification exponentielle avec $z$. Au final, le système ainsi construit permet une étude modale atemporelle pour ces deux inconnues, en fonction de l'état moyen de l'atmosphère (profils des vents, de la densité, …) et du choix des paramètres spectraux, \ie périodes et longueurs d'onde. +Les $*$ dénotent une renormalisation des deux paramètres $\tuz$ et $\tp$ : $\tuz^* = \sqrt{\rho_0} \tuz$ et $\tp^* = \frac{\tp}{\sqrt{\rho_0}}$ (d'où les termes en $\pm \frac{1}{2} \dz{\ln{\rho_0}}$). Cette renormalisation assure que les perturbations de la vitesse verticale et de la pression sont indépendantes de l'amplification exponentielle avec $z$. Au final, le système ainsi construit permet une étude modale atemporelle pour ces deux inconnues, en fonction de l'état moyen de l'atmosphère (profils des vents, de la densité, …) et du choix des paramètres spectraux, \ie périodes et longueurs d'onde. Dans ces équations n'apparaissent pas explicitement de termes dissipatifs, lesquels seraient liés par exemple à des transferts thermiques, des interactions d'ondes, du déferlement, ou bien de la viscosité. La théorie linéaire dérivée dans ces conditions se prête bien à une étude modale mono-paramétrique, mais il faut garder à l'esprit que les ordres de grandeur absolus déduits de simulations opérées sous ces conditions n'ont qu'une valeur limitée, car les mécanismes dissipatifs et les interactions multi-échelles sont vraisemblablement primordiaux dans les cas réalistes (\cite{Fritts2003}). Par ailleurs, en ne considérant ni les transferts au flux moyen (conservation du pseudo-moment), ni la turbulence induite, ni les échanges de chaleur, on ne s'autorise aucun \emph{feedback} lors de la propagation pseudo-spectrale. C'est donc avec un modèle « découplé » que le travail numérique a été effectué. diff --git a/docs/fac/stage/IPGP Saint Maur/tex/rapport/rapport.tex b/docs/fac/stage/IPGP Saint Maur/tex/rapport/rapport.tex index 9b75159..b783eab 100755 --- a/docs/fac/stage/IPGP Saint Maur/tex/rapport/rapport.tex +++ b/docs/fac/stage/IPGP Saint Maur/tex/rapport/rapport.tex @@ -22,9 +22,9 @@ rapport de stage… \include{visco-elastique} -\clearpage +\include{conclusion} -\include{annexes} +%\include{annexes} \bibliographystyle{agu} \bibliography{/home/jd/biblio/jabref/jd.bib} -- 2.11.4.GIT