function [edata1,edata2]=extractRegion(data1,data2,region);
% Extract data for a given geographical region
% Input: [numeric data] [textual data] [region code]
%
% [region code] is compared to column 2 in the textual data matrix
%
% Note: rows in each data matrix must be equal

j=1;
% numeric data
edata1=[];
% textual data
edata2=cell(0);
% iterate over data
for i=1:size(data1,1)
    % compare region codes
    if (strcmp(data2(i,2),region) == 1)
        % extract row of data and save in numeric and textual output matrix
        edata1(j,:) = data1(i,:);
        edata2(j,:) = data2(i,:);
        j = j + 1;
    end
end