1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
| function [] = altitude()
clc
q5 = 0:0.25:1;
q50 = 0.02:0.02:1;
pathname = '~/Matlab/Codes/Altitudes/BV';
listbv = dir(pathname);
for i = 3:length(listbv)
if listbv(i).isdir
liststation = dir(fullfile(pathname,listbv(i).name));
for j = 3:length(liststation)
fid = fopen(fullfile(pathname,listbv(i).name,liststation(j).name));
data = textscan(fid,'%s %f %f %f %f %s %s','HeaderLines',3);
fclose(fid);
lat = data{3};
alt = data{4};
lat = mean(lat);
if lat>10000
lat = lat./1000;
end
medalt = quantile(alt,0.5)';
altz5 = quantile(alt,q5)';
altz50 = quantile(alt,q50)';
beta = 0; % pour le Québec selon Valéry 2010
Vmin = 0.1; % selon Valéry 2010
altitudes = fullfile('~/Matlab/Codes/Altitudes/BV',sprintf('data_alt_lat_%s.mat',liststation(j).name));
save(altitudes, 'medalt', 'altz5', 'altz50','lat','beta','Vmin','-mat');
end
end
end |
Partager