source: trunk/src/essai2.m @ 327

%+
%
% .. _essai2.m:
%
% ========
% essai2.m
% ========
%
% DESCRIPTION
% ===========
%
% EVOLUTIONS
% ==========
%
% $Id$
%
% $URL$
%
% - fplod 20110805T140251Z cratos (Linux)
%
%   * " is invalid in matlab
%
% - fplod 20110705T125628Z cratos.locean-ipsl.upmc.fr (Linux)
%
%   * minimal header
%
%-

mfilename
ifigure=1;
ifigue=ifigure+1;
figure(ifigure);
%nl=3;
%nc=3;
precip_journaliere=load('velingra_q_1961-2000.xsl.csv');
% initialise les annee
domaine_an=1961:2000;
nban=size(domaine_an,2);
% initialise le tableau des mois pour chaque annee
domaine_month=1:12;
nbmonth=size(domaine_month,2);
nbdecad=3;
% initialise le tableau contenant les sommes par decade
% pour chaques mois de  toutes les annee
sumdecad=zeros(nban,nbdecad*nbmonth);
xdata=zeros(nbdecad*nbmonth);
% boucle sur les annees
for ian=1:nban;
   % boucle sur les mois et sur les decads
   for imonth=domaine_month;
      for idecad=1:3;
          ndecad=(imonth-1)*nbdecad+idecad;
          ijour=(idecad-1)*10+1;
          xdata(ndecad)=datenum(ian+1961-1,imonth,ijour);
          a=(idecad-1)*10+3;
          b=(idecad-1)*10+12;
       if (idecad < 3);
            cmin=a;
            cmax=b;
       else
            cmin=a;
            cmax=b+1;
      % calcul de la ligne contenant les precipitations du mois
      % courant/annee courante
      iligne=imonth+((ian-1)*12);
            if (iligne <= size(precip_journaliere,1));
                sumdecad(ian,ndecad)=sum(precip_journaliere(iligne,cmin:cmax));
            end
       end
     end
   end
     start_date=datenum(ian+1961-1,1,1);
     end_date=datenum(ian+1961-1,12,31);
                   % for imonth=1:12;
                   %     for idecad=1:3;
                    %         ndecad=(imonth-1)*nbdecad+idecad;
                    %         sumdecad(ndecad)=sum(meandecad(precip_journaliere(1:36)));
                     %        ijour=(idecad-1)*10+1;
                     %       xdata(ndecad)=datenum(imonth,ijour);
                       % end
                    %end
 % ++todo+++ subplot(nl,nc,(ian+1)-1961);
   % dessine les precipitations des sommes partielles  decadaires de chaque mois   % de l'annee courante
   % chgt de dimension de sumdecad de 3 à 1 pour dessin
   % y3d=sumdecad(ian,:,:);
   % y1d=reshape(y3d,1,nbmonth*nbdecad);
    % whos y1d
    % x=1:nbmonth*nbdecad;
   % whos x
  % ifigure=ifigure+1;
   %figure(ifigure);
 %  plot(xdata(:,1),sumdecad(ian,:),'x');
  % set x values were we want tick:one tick for month
  % xticklabel=datestr(xdata,2);
  % set(gca,'xtick',xdata);
  % set(gca,'xticklabel',xticklabel);
  % xlabel(['domaine-month of the year' num2str(ian)]);
  % ylabel('sumdecad(precip-journaliere) (mm)');
 % title("sum decadaire of velingra precipitation");
end
% calcul du cycle annuel moyen
moydecad=zeros(nbmonth*nbdecad);
% boucle sur les decads
for idecad=1:nbdecad*nbmonth;
    moy(idecad)=mean(sumdecad(1:nban,idecad));
    ijour=(idecad-1)*10+1;
    xdata(idecad)=datenum(1961,imonth,ijour);
end

   % ifigue=ifigure+1;
    ifigure=1;
    plot(xdata(:,1),moy,'-');
   % set x values were we want tick:one tick for decad of month
    xticklabel=datestr(xdata,1);
    % tous les 10 jours de chaque mois de chaque année
    %xtick=start_date:10:end_date;
    set(gca,'xtick',xdata);
    set(gca,'xticklabel',xticklabel);
    xlabel(['decad of the year' num2str(ian)]);
    ylabel('mean(sumdecad(precip-journaliere)) (mm)');
    title('mean(sumdecad) of velingra precipitation');