function [t,y]=cemrd2(fn)
% [t,y]=cemrdtd('filename')
% Read CEM time domain file
tdf=fopen(fn);
if tdf==-1 
disp('Can not open file!');

else
tdid=fscanf(tdf,'%*s %d %d %d',3);
td=fscanf(tdf,'%g %g',inf);
t=[td(1:2:size(td,1))];
y=[td(2:2:size(td,1))];
fclose(tdf);
end

end