% dump HDF SRF as ascii data and check results one more time
cmd = sprintf('!hdp dumpsds -dx -n srfval %s > srfval.asc', hfile);
eval(cmd);
fid = fopen('srfval.asc', 'r');
% [srf1, count] = fread(fid, inf, 'double');
[srf1, count] = fread(fid, inf, 'single');
srf1 = reshape(srf1, 471, 2378);
isequal(srfval, srf1)