high_res_cube/matlab-grid-generator/write_tile.m

function [] = write_tile(fn,A,varargin)
% Write tiled field to single file "fn.bin".
%
% Data is in IEEE-bigendian real*8 format by default.
%
% Dimensions of A must be (n,n)
%
% Usage:
% >>  write_tiles('LATC',latC);
% >>  write_tiles('LATC',latC,'real*4');

% Default precision
prec='real*8';
ieee='b';

% Check optional arguments
args=char(varargin);
while (size(args,1) > 0)
 if deblank(args(1,:)) == 'real*4' | deblank(args(1,:)) == 'float32'
  prec='real*4';
 elseif deblank(args(1,:)) == 'real*8' | deblank(args(1,:)) == 'float64'
  prec='real*8';
 else   
  sprintf(['Optional argument ' args(1,:) ' is unknown'])
  return
 end    
 args=args(2:end,:);
end     

fid=fopen(sprintf('%s.bin',fn),'w',ieee);
if fid<0
  error('Unable to create file')
end
fwrite(fid,A,prec);
fclose(fid);
1	function [] = write_tile(fn,A,varargin)
2	% Write tiled field to single file "fn.bin".
3	%
4	% Data is in IEEE-bigendian real*8 format by default.
5	%
6	% Dimensions of A must be (n,n)
7	%
8	% Usage:
9	% >> write_tiles('LATC',latC);
10	% >> write_tiles('LATC',latC,'real*4');
11
12	% Default precision
13	prec='real*8';
14	ieee='b';
15
16	% Check optional arguments
17	args=char(varargin);
18	while (size(args,1) > 0)
19	if deblank(args(1,:)) == 'real*4' \| deblank(args(1,:)) == 'float32'
20	prec='real*4';
21	elseif deblank(args(1,:)) == 'real*8' \| deblank(args(1,:)) == 'float64'
22	prec='real*8';
23	else
24	sprintf(['Optional argument ' args(1,:) ' is unknown'])
25	return
26	end
27	args=args(2:end,:);
28	end
29
30	fid=fopen(sprintf('%s.bin',fn),'w',ieee);
31	if fid<0
32	error('Unable to create file')
33	end
34	fwrite(fid,A,prec);
35	fclose(fid);