fafany0503
Messages postés1Date d'inscriptionmercredi 16 décembre 2015StatutMembreDernière intervention16 décembre 2015
-
Modifié par Strumpfette le 16/12/2015 à 15:44
bonjour
j'ai un code compliqué matlab mais qui tourne pour une seule ligne de donnée alors que moi j'ai des données toutes les demi heure qu'est ce que je peux faire sachant que j'ai droit à modifier seulement le main:
% PARAMETRES D'ENTREE
input= xlsread('datainput_EC');
%n = size(input,1);
Xtr = 3484110.393;
Ytr = 1220635.907;
zEC = 9.52;
for i=1:n
ustar(i,1)= input(i,1);
dir_wind(i,1) = input(i,2);
L (i,1)= input(i,3);
sigmav(i,1)= input(i,4);
h = 5;
resol = 0.5;
save var
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% lecture des données d'entrée du fichier data_EC
%[Ytr,Xtr,zEC,ustar,dir_wind,L,sigmav,h,resol]=donnees_EC(data);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
d=0.66*h; %d en m
z0=0.1*h; %z0 en m
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Ty=200; %Ty temps lagrangien (à voir si on peut le calculer avec EDIRE)
zm=zEC-d;
c=0.56; %c constante HW 92
% cette constante varie d'après Horst&Weil 1994 : c=0.66 pour r=2; c=0.63
% pour r=1.5 et 0.56 pour r=1
p=1.55;%p constante HW 92
kappa=0.41;%kappa constante
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%Debut des calculs%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% appel de la fonction plumeheight
[x,z]=plumheight(p,L,kappa,z0,resol,dir_wind);
% appel de la fonction plumeadvection
[u,U,dzdx,r]=plumadvection(ustar,zm,L,z0,z,kappa,c,p);
% appel de la fonction Fetch
Fy=FetchEC(zm,z,dzdx,u,U,r);
% appel de la fonction crosswind
Dy=crosswind(x,resol,U,Ty,sigmav);
% calcul du footprint
Fp=FprintEC(Fy,Dy);
% calcul de l'image du footprint
% on donne comme coordonnées de référence le milieu du bord de l'image et
% le milieu de l'image
[out,coord]=rotate_image(-90+dir_wind, Fp, [1;75] );
end
c'est moi qui a essayé de faire les boucles for et voila l'erreur
??? Error using ==> rdivide
Matrix dimensions must agree.
merci d'avance