% opendap url
opdfile = 'c1440_NR.inst01hr_3d_T_Cv.20060918_0900z.nc4'

% bounding box
%   lons = -130:0.5:-65
%   lats =   25:0.5:50
imin = round((-130. + 180.)/0.5);
imax = round(( -65. + 180.)/0.5);
jmin = round((  25 +  90.)/0.5);
jmax = round((  50 +  90.)/0.5);

% corresponding array sizes
im = imax-imin+1;
jm = jmax-jmin+1;
lm = 72; % read all 72 levels

% start/count for ncread
start = [imin, jmin, 1, 1];
count = [im, jm, lm, 1];

% now read
fprintf('Reading T (subset)...');
T = ncread(opdfile, 'T', start, count);
fprintf('done.\n')

% compute max/min of T
fprintf('max(T): %f\n', max(T(:)));
fprintf('min(T): %f\n', min(T(:)));