[x, y, zone ] =  deg2utm(tp_lat, tp_long);

fid = fopen(['tethys_pings_bound.log'], 'wt');
for i=1:length(x)
  fprintf(fid, '%f,%f,%f\n', x(i), y(i), tp_sr(i));
end
fclose(fid)