MITgcm_contrib/tides/writebin.m

function writebin(fnam,fld,typ,prec,skip,mform)

% Function writebin(fnam,fld,typ,prec,skip,mform)
% write N-D binary field
%
% INPUTS
% fnam  input path and file name
% fld   input array
% typ   0: sequential FORTRAN;  1: plain binary (the default)
% prec  numeric precision (default 'real*4')
% skip  records to skip before writing (default 0)
% mform machine format (see fopen; default: 'ieee-be')
%
% SEE ALSO
% readbin

if nargin < 6, mform='ieee-be'; end
if nargin < 5, skip=0; end
if nargin < 4, prec='real*4'; end
if nargin < 3, typ=1; end
if skip>0 & typ==0, error('feature not implemented yet'); end
if nargin < 2, error('please specify array and output file name'); end
reclength=0;

if exist(fnam)==2
  fid=fopen(fnam,'r+',mform);
else
  fid=fopen(fnam,'w',mform);
end

switch typ
  
  case 0

    write_record(fid,fld,prec);
  
  case 1

    if skip==0
      
      fwrite(fid,fld,prec);
    
    else

      switch prec
        case {'int8','integer*1'}
          reclength=prod(size(fld));
        case {'int16','integer*2','uint16','integer*2'}
          reclength=2*prod(size(fld));
        case {'int32','integer*4','uint32','single','real*4','float32'}
          reclength=4*prod(size(fld));
        case {'int64','integer*8','uint64','double','real*8','float64'}
          reclength=8*prod(size(fld));
      end
      
      if fseek(fid,skip*reclength,'bof') == 0
        
        fwrite(fid,fld,prec);
        
      else
        
        fseek(fid,0,'eof');
        file_length=ftell(fid);
        fwrite(fid,zeros((skip*reclength-file_length),1),'int8');
        fwrite(fid,fld,prec);
        
      end
      
    end
    
end

fid=fclose(fid);
1	function writebin(fnam,fld,typ,prec,skip,mform)
2
3	% Function writebin(fnam,fld,typ,prec,skip,mform)
4	% write N-D binary field
5	%
6	% INPUTS
7	% fnam input path and file name
8	% fld input array
9	% typ 0: sequential FORTRAN; 1: plain binary (the default)
10	% prec numeric precision (default 'real*4')
11	% skip records to skip before writing (default 0)
12	% mform machine format (see fopen; default: 'ieee-be')
13	%
14	% SEE ALSO
15	% readbin
16
17	if nargin < 6, mform='ieee-be'; end
18	if nargin < 5, skip=0; end
19	if nargin < 4, prec='real*4'; end
20	if nargin < 3, typ=1; end
21	if skip>0 & typ==0, error('feature not implemented yet'); end
22	if nargin < 2, error('please specify array and output file name'); end
23	reclength=0;
24
25	if exist(fnam)==2
26	fid=fopen(fnam,'r+',mform);
27	else
28	fid=fopen(fnam,'w',mform);
29	end
30
31	switch typ
32
33	case 0
34
35	write_record(fid,fld,prec);
36
37	case 1
38
39	if skip==0
40
41	fwrite(fid,fld,prec);
42
43	else
44
45	switch prec
46	case {'int8','integer*1'}
47	reclength=prod(size(fld));
48	case {'int16','integer2','uint16','integer2'}
49	reclength=2*prod(size(fld));
50	case {'int32','integer4','uint32','single','real4','float32'}
51	reclength=4*prod(size(fld));
52	case {'int64','integer8','uint64','double','real8','float64'}
53	reclength=8*prod(size(fld));
54	end
55
56	if fseek(fid,skip*reclength,'bof') == 0
57
58	fwrite(fid,fld,prec);
59
60	else
61
62	fseek(fid,0,'eof');
63	file_length=ftell(fid);
64	fwrite(fid,zeros((skip*reclength-file_length),1),'int8');
65	fwrite(fid,fld,prec);
66
67	end
68
69	end
70
71	end
72
73	fid=fclose(fid);