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 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69
| close all;
clear all;
clc;
%Le load correspond à la matrice désiré
%bien faire attention à la taille de la matrice
load ('C:\Users\benjamin\Documents\CDD_inge_climato\donnees_Hobo_matlab\matrices\hobo_loc_dates\comparaison\terrain_TCmean.mat');
Station={'Station 01';'Station 02';'Station 03';'Station 04';'Station 05';'Station 06';'Station 07';'Station 08';'Station 09';'Station 10';'Station 11';'Station 12';'Station 13';'Station 14';'Station 15';'Station 16';'Station 17';'Station 18';'Station 19';'Station 20';'Station 21';'Station 22';'Station 23';'Station 24';'Station 25';'Station 26';'Station 27';'Station 28';'Station 29';'Station 30';'Station 31';'Station 32';'Station 33';'Station 34';'Station 35';'Station 36';'Station 37';'Station 38';'Station 39';'Station 40';'Station 41';'Station 42';'Station 43';'Station 44';'Station 45';'Station 46';'Station 47';'Station 48';'Station 49';'Station 50';'Station 51'};
Stations={'Station 01 - Meteo France Longvic';'Station 02 - Meteo France Longvic';'Station 03 - Meteo France Longvic';'Station 04 - Meteo France Longvic';'Station 05 - Meteo France Longvic';'Station 06 - Meteo France Longvic';'Station 07 - Meteo France Longvic';'Station 08 - Meteo France Longvic';'Station 09 - Meteo France Longvic';'Station 10 - Meteo France Longvic';'Station 11 - Meteo France Longvic';'Station 12 - Meteo France Longvic';'Station 13 - Meteo France Longvic';'Station 14 - Meteo France Longvic';'Station 15 - Meteo France Longvic';'Station 16 - Meteo France Longvic';'Station 17 - Meteo France Longvic';'Station 18 - Meteo France Longvic';'Station 19 - Meteo France Longvic';'Station 20 - Meteo France Longvic';'Station 21 - Meteo France Longvic';'Station 22 - Meteo France Longvic';'Station 23 - Meteo France Longvic';'Station 24 - Meteo France Longvic';'Station 25 - Meteo France Longvic';'Station 26 - Meteo France Longvic';'Station 27 - Meteo France Longvic';'Station 28 - Meteo France Longvic';'Station 29 - Meteo France Longvic';'Station 30 - Meteo France Longvic';'Station 31 - Meteo France Longvic';'Station 32 - Meteo France Longvic';'Station 33 - Meteo France Longvic';'Station 34 - Meteo France Longvic';'Station 35 - Meteo France Longvic';'Station 36 - Meteo France Longvic';'Station 37 - Meteo France Longvic';'Station 38 - Meteo France Longvic';'Station 39 - Meteo France Longvic';'Station 40 - Meteo France Longvic';'Station 41 - Meteo France Longvic';'Station 42 - Meteo France Longvic';'Station 43 - Meteo France Longvic';'Station 44 - Meteo France Longvic';'Station 45 - Meteo France Longvic';'Station 46 - Meteo France Longvic';'Station 47 - Meteo France Longvic';'Station 48 - Meteo France Longvic';'Station 49 - Meteo France Longvic';'Station 50 - Meteo France Longvic';'Station 51 - Meteo France Longvic'};
mois={'15/06';'1/07';'15/07';'1/08';'15/08';'1/09';'15/09'};
C=nanmean(terrain_TCmean,2);
for i=1:51
Tmean_ete2014=terrain_TCmean(:,i);
S51=C;
min(S51)
max(S51)
S51=reshape(S51,[24,115]);
colormap('jet(16)')
figure(i); clf
subplot (3,1,1);
contourf(1:115,1:24,S51,[6:2:38]);
c=colorbar('southoutside');
set(gca,'xtick',[9;24;39;55;70;86;101],'xticklabel',mois,'xgrid','on')
ylabel('Heures');
title 'Moyenne globale des stations';
axes(c);
xlabel('Degrés Celsius');
caxis([0 40]);
xlim([0 40]);
S01=Tmean_ete2014;
min(S01)
max(S01)
S01=reshape(S01,[24,115]);
figure(i); hold on;
subplot (3,1,2);
colormap('jet(16)')
contourf(1:115,1:24,S01,[6:2:38]);
c=colorbar('southoutside');
set(gca,'xtick',[9;24;39;55;70;86;101],'xticklabel',mois,'xgrid','on')
ylabel('Heures');
title (Station{i});
axes(c);
xlabel('Degrés Celsius');
caxis([0 40]);
xlim([0 40]);
print -depsc2 -tiff s51&s01;
DIFF=S01-S51;
min(DIFF)
max(DIFF)
figure(i);
subplot (3,1,3);
contourf(1:115,1:24,DIFF,[-9:1:6]);
c=colorbar('southoutside');
colormap('jet(16)');
set(gca,'xtick',[9;24;39;55;70;86;101],'xticklabel',mois,'xgrid','on')
ylabel('Heures');
title (Stations{i});
axes(c);
xlabel('Différences de températures en degrés Celsius');
caxis([-9 6]);
xlim([-9 6]);
clear DIFF
clear S01
clear S51
clear c
end |
Partager