function make_AA_POCA_h5(in_dir, out_dir, out_file)

if ~exist('out_file','var') && ~isempty(in_dir);

disp('make_AA_POCA_h5: Generating bin centers and queue');

dx0=2e5;

[~, d]=unix(['ls ', in_dir,' | grep -P "CS.*00..mat"']);

d=strsplit(deblank(d));

for k=1:length(d);

d{k}=strrep(d{k}, in_dir,'');

temp=load([in_dir, '/', d{k}],'X0');

if ~isempty(temp);

X0{k}=temp.X0;

end

end

if ~exist(out_dir,'dir');

mkdir(out_dir);

end

uBins=unique(round_to(cat(1, X0{:}), dx0));

fid=fopen([out_dir,'/bin_list.txt'],'w');

for k=1:length(uBins);

fprintf(fid,'POCA_%dE_%dN.h5\n', real(uBins(k))/1000, imag(uBins(k))/1000);

end

fclose(fid);

save([out_dir,'/bins_for_CS_files.mat'], 'X0', 'd', 'dx0', 'in_dir');

%fid_list=fopen('POCA_h5_C/south/AA_bin_list.txt','r');

fid_list=fopen([out_dir,'/bin_list.txt'],'r');

fid_queue=fopen('queue_for_AA_POCA.txt','w');

while ~feof(fid_list);

fname=deblank(fgetl(fid_list));

fprintf(fid_queue, 'make_AA_POCA_h5(''%s'',''%s'', ''%s'');\n', in_dir, out_dir, fname);

end

fclose(fid_list);

fclose(fid_queue);

disp(sprintf('now run the queue in queue_for_AA_POCA.txt\n'));

return

if ~isempty(in_dir) && exist('out_file','var');

fprintf(1,'make_AA_POCA_h5: working on %s\n', out_file);

load([out_dir,'/bins_for_CS_files.mat'])

temp=regexp(out_file,'([-+]*\d+)E_([-+]*\d+)N','tokens');

x0=str2double(temp{1}{1})*1000;

y0=str2double(temp{1}{2})*1000;

% out_file=sprintf('%s/POCA_%dE_%dN.h5', out_dir, x0(kx)/1000, y0(ky)/1000)

% if exist(out_file,'file');

% continue

% end

% fid=fopen(out_file,'w'); fprintf(fid,'patience...\n'); fclose(fid);

clear D0

for k=1:length(X0);

these=real(X0{k})>x0-dx0/2 & real(X0{k}) <= x0+dx0/2 & imag(X0{k})>y0-dx0/2 & imag(X0{k}) <= y0+dx0/2;

if ~any(these);

continue

end

LL=load([in_dir,d{k}],'D');

if isfield(LL,'D') && ~isempty(LL.D);

D0(k)=LL.D;

end

end

if ~exist('D0','var'); return; end

f=fieldnames(D0);

for kf=1:length(f);

D.(f{kf})=cat(1, D0.(f{kf}));

end

these=real(D.x)>x0-dx0/2 & real(D.x) <= x0+dx0/2 & imag(D.x)>y0-dx0/2 & imag(D.x) <= y0+dx0/2 & isfinite(D.x);

if ~any(these); return; end

D=index_struct(D, these);

% check this

[M, INDEX]=index_point_data_h5('build_index_and_sort', D, [], out_dir, 1e4, 10.2);

index_point_data_h5('write_h5',[out_dir,'/', out_file], M, INDEX, true);

save([out_dir,'/', strrep(out_file,'.h5','_index.mat')], 'INDEX');

if isempty(in_dir) && exist(out_dir,'dir');

[~, h5_files]=unix(['ls ', out_dir,'/*.h5']);

h5_files=strsplit(deblank(h5_files));

for k=1:length(h5_files);

try

II=h5info(h5_files{k},'/INDEX');

disp(II.Filename);

catch

fprintf(1,'failure for %s\n', h5_files{k});

delete(h5_files{k})

end

end

INDEX_list=index_point_data_h5('collect_index', {out_dir}, '*.h5', {'./'});

master_index=index_point_data_h5('make_master_index', INDEX_list, [out_dir,'/master_index_h5.mat'], out_dir)