matlab_class/gcmfaces_IO/write2meta.m

function []=write2meta(myFile,myDim,varargin);
%purpose:       write meta file to match binary file
%
%inputs:        myFile is the file name (must finish with '.data')
%               myDim is the size of a field (2D or 3D integer vector)
%(optional)     myPrec is the file myPrecision (32, by default, or 64)

%precision:
if nargin>2; myPrec=varargin{1}; else; myPrec=32; end;
%number of dimensions
nDim=length(myDim);
%number of records:
tmp1=dir(myFile);
nRec=tmp1.bytes/myDim(1)/myDim(2)/myPrec*8;
if nDim>2; nRec=nRec/myDim(3); end;
if nRec<1|floor(nRec)~=nRec; error('inconsistent dimensions'); end;
%check that creating a meta file makes sense w\r\t rdmds
if ~strcmp(myFile(end-4:end),'.data');
  error('file name must finish in ''.data''\n');
else;
  myFile=[myFile(1:end-5) '.meta'];
end;

%create the meta file:
fid=fopen(myFile,'wt'); 
%%
fprintf(fid,'nDims = [   %i ];\n',min(nDim,3));
fprintf(fid,' dimList = [\n');
fprintf(fid,' %5i, %5i, %5i,\n',myDim(1),1,myDim(1));
fprintf(fid,' %5i, %5i, %5i,\n',myDim(2),1,myDim(2));
if nDim>2; fprintf(fid,' %5i, %5i, %5i,\n',myDim(3),1,myDim(3)); end;
fprintf(fid,' ];\n');
fprintf(fid,' dataprec = [ ''float%2i'' ];\n',myPrec);
fprintf(fid,' nrecords = [ %5i ];\n',nRec);
%%
fclose(fid);

1	function []=write2meta(myFile,myDim,varargin);
2	%purpose: write meta file to match binary file
3	%
4	%inputs: myFile is the file name (must finish with '.data')
5	% myDim is the size of a field (2D or 3D integer vector)
6	%(optional) myPrec is the file myPrecision (32, by default, or 64)
7
8	%precision:
9	if nargin>2; myPrec=varargin{1}; else; myPrec=32; end;
10	%number of dimensions
11	nDim=length(myDim);
12	%number of records:
13	tmp1=dir(myFile);
14	nRec=tmp1.bytes/myDim(1)/myDim(2)/myPrec*8;
15	if nDim>2; nRec=nRec/myDim(3); end;
16	if nRec<1\|floor(nRec)~=nRec; error('inconsistent dimensions'); end;
17	%check that creating a meta file makes sense w\r\t rdmds
18	if ~strcmp(myFile(end-4:end),'.data');
19	error('file name must finish in ''.data''\n');
20	else;
21	myFile=[myFile(1:end-5) '.meta'];
22	end;
23
24	%create the meta file:
25	fid=fopen(myFile,'wt');
26	%%
27	fprintf(fid,'nDims = [ %i ];\n',min(nDim,3));
28	fprintf(fid,' dimList = [\n');
29	fprintf(fid,' %5i, %5i, %5i,\n',myDim(1),1,myDim(1));
30	fprintf(fid,' %5i, %5i, %5i,\n',myDim(2),1,myDim(2));
31	if nDim>2; fprintf(fid,' %5i, %5i, %5i,\n',myDim(3),1,myDim(3)); end;
32	fprintf(fid,' ];\n');
33	fprintf(fid,' dataprec = [ ''float%2i'' ];\n',myPrec);
34	fprintf(fid,' nrecords = [ %5i ];\n',nRec);
35	%%
36	fclose(fid);
37