fid = fopen('/scratch/bliu17/Laurence/run_dir/polarMaps/03-07pol.okc', 'rt');
Data = textscan(fid, '%f %f %f %f %f %f', 'headerLines', 13, 'CollectOutput', true);
fclose(fid);
Data= cell2mat( Data);
X=Data(:,1);
Y=Data(:,2);
polx=Data(:,4);
poly=Data(:,5);
%quiver(X,Y,polx,poly);
quiver(Y,X,poly,polx);
xlabel('Z');
ylabel('Y');