global_ocean.128x64x15/diags_matlab/mit_readfield.m

function h = mit_readfield(fname,dims,accuracy)
%function h = mit_readfield(fname,dims,accuracy)
  
  ieee='ieee-be';

  [fid message] = fopen(fname,'r',ieee);
  if fid <= 0
    error([message ', filename: ', [fname]])
  end

  v = fread(fid,accuracy);
  if length(v) == prod(dims)
    h = reshape(v,dims);
  elseif length(v) == prod(dims(1:end-1))
    h = reshape(v,dims(1:end-1));
  else
    error('dimensions do not match')
  end
  fclose(fid);
  
  return
1	cnh	1.1	function h = mit_readfield(fname,dims,accuracy)
2			%function h = mit_readfield(fname,dims,accuracy)
3
4			ieee='ieee-be';
5
6			[fid message] = fopen(fname,'r',ieee);
7			if fid <= 0
8			error([message ', filename: ', [fname]])
9			end
10
11			v = fread(fid,accuracy);
12			if length(v) == prod(dims)
13			h = reshape(v,dims);
14			elseif length(v) == prod(dims(1:end-1))
15			h = reshape(v,dims(1:end-1));
16			else
17			error('dimensions do not match')
18			end
19			fclose(fid);
20
21			return