timour_matlab/mscripts/Readmodel.m

%
%  Load M.I.T. GCM model data output files.
%
function array = Readmodel( fileName, dims, format )
if nargin  < 2
  fprintf(2,'**ERROR** Wrong number of arguments.\n')
  fprintf(2,'Usage: Readmodel(''pathname'', [dim1 dim2 dim3 ...])\n')
  return
end
if nargin == 2 
  format = 'float32';
end
if size(dims,1) ~= 1
  fprintf(2,'**ERROR** Array dimensions mis-specified.\n')
  fprintf(2,'Usage: Readmodel(''pathname'', [dim1 dim2 dim3 ...])\n')
  return
end
fid=fopen(fileName,'r','b');
if fid < 0 
  fprintf(2,'**ERROR** Unable to open input file "%s".\n',fileName)
  fprintf(2,'Usage: Readmodel(''pathname'', [dim1 dim2 dim3 ...])\n')
  return
end
arraySize=1+size(dims,2)+prod(dims);
array=zeros(1,arraySize);
array(1)=size(dims,2);
for i = 1 : size(dims,2);
  array(i+1)=dims(1,i);
end
[dummy,cnt]=fread(fid,1,'int');
if cnt ~= 1 
  fprintf(2,'**ERROR** Unable to read first byte of file "%s".\n',fileName)
  fprintf(2,'Usage: Readmodel(''pathname'', [dim1 dim2 dim3 ...])\n')
  fclose(fid)
  return
end
[array(1+size(dims,2)+1:size(array,2)),cnt]=fread(fid,prod(dims),format);
if cnt ~= prod(dims)
  fprintf(2,'**ERROR** Unable to read field data from file "%s".\n',fileName)
  fprintf(2,'Usage: Readmodel(''pathname'', [dim1 dim2 dim3 ...])\n')
  fclose(fid)
  return
end
return
1	edhill	1.1	%
2			% Load M.I.T. GCM model data output files.
3			%
4			function array = Readmodel( fileName, dims, format )
5			if nargin < 2
6			fprintf(2,'ERROR Wrong number of arguments.\n')
7			fprintf(2,'Usage: Readmodel(''pathname'', [dim1 dim2 dim3 ...])\n')
8			return
9			end
10			if nargin == 2
11			format = 'float32';
12			end
13			if size(dims,1) ~= 1
14			fprintf(2,'ERROR Array dimensions mis-specified.\n')
15			fprintf(2,'Usage: Readmodel(''pathname'', [dim1 dim2 dim3 ...])\n')
16			return
17			end
18			fid=fopen(fileName,'r','b');
19			if fid < 0
20			fprintf(2,'ERROR Unable to open input file "%s".\n',fileName)
21			fprintf(2,'Usage: Readmodel(''pathname'', [dim1 dim2 dim3 ...])\n')
22			return
23			end
24			arraySize=1+size(dims,2)+prod(dims);
25			array=zeros(1,arraySize);
26			array(1)=size(dims,2);
27			for i = 1 : size(dims,2);
28			array(i+1)=dims(1,i);
29			end
30			[dummy,cnt]=fread(fid,1,'int');
31			if cnt ~= 1
32			fprintf(2,'ERROR Unable to read first byte of file "%s".\n',fileName)
33			fprintf(2,'Usage: Readmodel(''pathname'', [dim1 dim2 dim3 ...])\n')
34			fclose(fid)
35			return
36			end
37			[array(1+size(dims,2)+1:size(array,2)),cnt]=fread(fid,prod(dims),format);
38			if cnt ~= prod(dims)
39			fprintf(2,'ERROR Unable to read field data from file "%s".\n',fileName)
40			fprintf(2,'Usage: Readmodel(''pathname'', [dim1 dim2 dim3 ...])\n')
41			fclose(fid)
42			return
43			end
44			return