clear
fname = '../read_hdf_sd_demo/CISM_DX_1_day_data.hdf';
fname = ['/mnt/hda7/tmp/tests/CISM_DX_DATA/measurements/' ,...
	 'uniform_time_series/CISM_DX_1_day_data.hdf'];

[X,T] = uniform_time_series_slice('Bx from OMNI',[2000 1 1],[2000 1 1],'min');

break
Info  = hdfinfo(fname);

tic
temp1 = hdfread(Info.SDS(end));
temp1(100)
toc

tic
temp2 = double(hdfread(Info.SDS(end),'Index',{[100 1],[1 1],[10 1]}));
temp2(1)
toc

if (max((temp2-temp1(100:109))) > 0)
  error('hdfread error');
else
  fprintf('test hdfread: ok\n');fflush(stdout);
end

break

SDS        = load_hdf(fname,0);

n_datasets = load_hdf(fname);

% Load the information from the first data set

% Load the data from the first data set
[SDS,data] = load_hdf(fname,0);

SDS