% Get the number of data sets

fname      = '../read_hdf_sd_demo/CISM_DX_1_day_data.hdf';
n_datasets = load_hdf(fname)
SDS        = load_hdf(fname,0);
SDS        = load_hdf(fname,0,1);

break

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

SDS