% NOM: Rayonnement_Net.m
% 
% AUTEURS: Jean-Philippe PAQUIN & Camille GARNAUD
% 
% DATE: Mars 2011
% 
% REVISIONS: 
% 
% DESCRIPTION: Diagnostics de base pour les simulations >offline< de 
%              CLASS3.5 dans le contexte du Laboratoire pratique du cours
%              SCA7025
%              
% ENTREES: Fichiers ASCII comprenant le nom de la station (ORLANDO,
%          YUMA, FORT McMURRAY, FLINT) avec les extensions .OF(1-8)s
%          (voir le document Word pour la description detaillee des 
%          variables)
% 
% SORTIES: Figures
% 
%------------------------------------------------------------------------

clear

cd /skynet2_exec1/camille/SCA7025_CLASS/

%=================================================
%===  LECTURES DES FICHIERS ASCII NECESSAIRES  ===
%===  OFi sont des structures avec OFi.data    ===
%===  pour acceder aux donnees                 ===

OF1=importdata('YUMA.OF1',' ',3);
OF4=importdata('YUMA.OF4',' ',3);

% Figure montrant l'évolution sur 10 ans
hold on 
   plot(OF1.data(:,3),'b')                      % Net SW at sfc
   plot(OF1.data(:,4),'r')                      % Net LF at sfc
   title('Rayonnement Net')
   xlabel('Jours')
   ylabel('W/m^2')
   legend('Solaire Net','Infrarouge Net')
   box on ; grid on
hold off

% Figure montrant l'évolution sur 24 heures
figure(2)
hold on 
   plot(OF4.data(1:49,5),'b')                      % Net SW at sfc
   plot(OF4.data(1:49,6),'r')                      % Net LF at sfc
   title('Rayonnement Net')
   xlabel('Heures')
   ylabel('W/m^2')
   legend('Solaire Net','Infrarouge Net')
   box on ; grid on
hold off


%%% Exemple de cycle annuel moyen du rayonnement solaire

store=zeros(10,365);

for n = 1:10
    a = ((n-1)*365) + 1;  % Borne inférieure
    b = n*365;            % Borne supérieure
    store(n,:)=OF1.data(a:b,3);
end

annee_moyenne=mean(store);  % Moyenne de chaque jour
mois=(1:365)/30.41667;

figure(3)
hold on 
   plot(mois,annee_moyenne)                      % Net SW at sfc
   title('Rayonnement Solaire Net')
   xlabel('Mois')
   ylabel('W/m^2')
   box on ; grid on
hold off

