LoadData = 1;

if LoadData 
D=load('9641C_201301_F52.max');

z = [];
for i = 1:size(D,1);
if ~isempty(find(D(i,:) < 0));
z = [z i];
disp(i);
end
end
Z = setdiff(1:size(D,1),z);
E = D(Z,:);

end

k = 0;


k = 2;
j = 1;
Ix = 1;
clear Station
Station{1}(1) = 1;

while k < size(E,1)-1;

	a = num2str(E(k,1));
	b = num2str(E(k-1,1));
	if ~strncmp(a,b,6);
		Ix = Ix + 1;
		M(j,:) = [k Ix];
		j = j+1;
	end
	k = k+1;
end

for i = 1:size(M,1)-1
A = [M(i,1):M(i+1,1)-1];
s = num2str(E(A,1));
s = s(:,end-3:end);
F(1:length(A),1) = str2num(s);
F(1:length(A),2:14) = round(E(A,2:14)/10);
size(M,1)-i

%save(['~/Teaching/MA211/ClimateData/MaxTempsStation_',num2str(i),'.mat'],'F','-mat');
end


% To Load Paste this:
% 
% FF = load(['~/Teaching/MA211/ClimateData/MaxTempsStation_',num2str(k),'.mat'],'F','-mat');


