


fn = 'belem_2017.csv';

fh = fopen(fn);

% Read header line

line1 = fgetl(fh);

commas = strfind(line1,',');
commas = [0 commas numel(line1)+1];
nf = numel(commas)-1;

for ii=1:nf
  tmpstr = line1(commas(ii)+1:commas(ii+1)-1);
  tmpstr(isspace(tmpstr))=[];
  field{ii} = tmpstr;
end

data.mtime = [];
data.(field{1})='';
data.(field{2})='';
data.(field{3})='';
for ii=4:nf
  data.(field{ii}) = [];
end

% Start reading data file
done=false;

while (~done)

  [cc pos] = textscan(fh, '%s',nf, 'Delimiter',','); %,'TreatAsEmpty',{'////'});

  if(pos==-1 | numel(cc{1})==0)
    done=true;
  else 
    data.(field{1})(:,end+1) = cc{1}{1};
    data.(field{2})(1:10,end+1) = cc{1}{2}; %datenum(cc{1}{2},'dd.mm.yyyy'); 
    data.(field{3})(1:2,end+1) = cc{1}{3}; %str2num(cc{1}{3});
    data.mtime(end+1) = datenum([cc{1}{2} ' ' cc{1}{3}],'dd/mm/yyyy hh');
    for ii=4:nf
      str = cc{1}{ii};
      if(strcmp(str,'////'))
	val=NaN;
      else
	val = str2num(str);
      end
      data.(field{ii})(:,end+1) = val;
    end
  end

end

% FIX TIME: IT comes as GMT. Subtract 3hrs to get local time
data.mtime = data.mtime-3/24;




% Check wind direction as a function of hour and month
%
% April 2017 - March 2018
for im=1:12
  sdate = datenum([2017 04+(im-1) 01 0 0 0]);
  edate = datenum([2017 04+(im) 01 0 0 0]);
  idate{im} = find(data.mtime>=sdate & data.mtime<edate);

  dhour = mod(data.mtime,1)*24;
  % Early hours
  id1{im} = find(data.mtime>=sdate & data.mtime<edate & dhour>=0 & dhour<6);
  % Morning
  id2{im} = find(data.mtime>=sdate & data.mtime<edate & dhour>=6 & dhour<12);
  % Afternoon
  id3{im} = find(data.mtime>=sdate & data.mtime<edate & dhour>=12 & dhour<18);
  % Night
  id4{im} = find(data.mtime>=sdate & data.mtime<edate & dhour>=18 & dhour<24);
end

opts={'anglenorth',0,'angleeast',90,'vwinds',[0:.5:4],'freqlabelangle',210}; %,'maxfrequency',15};
% Four Day Roses
%
figure;
subplot(2,2,1);
WindRose(data.vento_vel([id1{:}]), data.vento_direcao([id1{:}]),{opts{:},'axes',gca,'titlestring','Madrugada'});
subplot(2,2,2);
WindRose(data.vento_vel([id2{:}]), data.vento_direcao([id2{:}]),{opts{:},'axes',gca,'titlestring','Manhã'});
subplot(2,2,3);
WindRose(data.vento_vel([id3{:}]), data.vento_direcao([id3{:}]),{opts{:},'axes',gca,'titlestring','Tarde'});
subplot(2,2,4);
WindRose(data.vento_vel([id4{:}]), data.vento_direcao([id4{:}]),{opts{:},'axes',gca,'titlestring','Noite'});

% Inverno - Nov Dec Jan Fev Mar Apr  - Verao: May Jun Jul Aug Sep Oct
%            8   9   10  11 12   1             2   3   4   5  6   7 
figure;
subplot(2,2,1);
WindRose(data.vento_vel([id1{[8 9 10 11 12 1]}]), data.vento_direcao([id1{[8 9 10 11 12 1]}]),{opts{:},'axes',gca,'titlestring','Madrugada (I)'});
subplot(2,2,2);
WindRose(data.vento_vel([id2{[8 9 10 11 12 1]}]), data.vento_direcao([id2{[8 9 10 11 12 1]}]),{opts{:},'axes',gca,'titlestring','Manhã (I)'});
subplot(2,2,3);
WindRose(data.vento_vel([id3{[8 9 10 11 12 1]}]), data.vento_direcao([id3{[8 9 10 11 12 1]}]),{opts{:},'axes',gca,'titlestring','Tarde (I)'});
subplot(2,2,4);
WindRose(data.vento_vel([id4{[8 9 10 11 12 1]}]), data.vento_direcao([id4{[8 9 10 11 12 1]}]),{opts{:},'axes',gca,'titlestring','Noite (I)'});


% Verao 
figure;
subplot(2,2,1);
WindRose(data.vento_vel([id1{[2 3 4 5 6 7]}]), data.vento_direcao([id1{[2 3 4 5 6 7]}]),{opts{:},'axes',gca,'titlestring','Madrugada (V)'});
subplot(2,2,2);
WindRose(data.vento_vel([id2{[2 3 4 5 6 7]}]), data.vento_direcao([id2{[2 3 4 5 6 7]}]),{opts{:},'axes',gca,'titlestring','Manhã (V)'});
subplot(2,2,3);
WindRose(data.vento_vel([id3{[2 3 4 5 6 7]}]), data.vento_direcao([id3{[2 3 4 5 6 7]}]),{opts{:},'axes',gca,'titlestring','Tarde (V)'});
subplot(2,2,4);
WindRose(data.vento_vel([id4{[2 3 4 5 6 7]}]), data.vento_direcao([id4{[2 3 4 5 6 7]}]),{opts{:},'axes',gca,'titlestring','Noite (V)'});

% Make 12 wind roses:
%

[fg1] = WindRose(data.vento_vel(idate{1}), data.vento_direcao(idate{1}),opts);
[fg2] = WindRose(data.vento_vel(idate{2}), data.vento_direcao(idate{2}),opts);



