utils/matlab/readmat.m

% read solver matrix
%--------------------
% input  : snx,sny,npx,npy
% output : as2d, aw2d, cent
% files  : aw2d*, as2d*, cent* in MITGCMUV mdsio format

nx=npx*snx;
ny=npy*sny;

aw2d=zeros(nx,ny);
as2d=zeros(nx,ny);
cent=zeros(nx,ny);
for i = 1:npx,
  ii=(i-1)*snx;
  xc1=ii+1:ii+snx;

  for j = 1:npy,
    jj=(j-1)*sny;
    yc1=jj+1:jj+sny;

    fname = sprintf('aw2d.%03i.%03i.data', i, j );
    fid=fopen(fname,'r');
    aw2d(xc1,yc1) = fread( fid, [snx,sny], 'float32' );
    fclose(fid);

    fname = sprintf('as2d.%03i.%03i.data', i, j );
    fid=fopen(fname,'r');
    as2d(xc1,yc1) = fread( fid, [snx,sny], 'float32' );
    fclose(fid);

    fname = sprintf('cent.%03i.%03i.data', i, j );
    fid=fopen(fname,'r');
    cent(xc1,yc1) = fread( fid, [snx,sny], 'float64' );
    fclose(fid);
  end
end
1	% read solver matrix
2	%--------------------
3	% input : snx,sny,npx,npy
4	% output : as2d, aw2d, cent
5	% files : aw2d, as2d, cent* in MITGCMUV mdsio format
6
7	nx=npx*snx;
8	ny=npy*sny;
9
10	aw2d=zeros(nx,ny);
11	as2d=zeros(nx,ny);
12	cent=zeros(nx,ny);
13	for i = 1:npx,
14	ii=(i-1)*snx;
15	xc1=ii+1:ii+snx;
16
17	for j = 1:npy,
18	jj=(j-1)*sny;
19	yc1=jj+1:jj+sny;
20
21	fname = sprintf('aw2d.%03i.%03i.data', i, j );
22	fid=fopen(fname,'r');
23	aw2d(xc1,yc1) = fread( fid, [snx,sny], 'float32' );
24	fclose(fid);
25
26	fname = sprintf('as2d.%03i.%03i.data', i, j );
27	fid=fopen(fname,'r');
28	as2d(xc1,yc1) = fread( fid, [snx,sny], 'float32' );
29	fclose(fid);
30
31	fname = sprintf('cent.%03i.%03i.data', i, j );
32	fid=fopen(fname,'r');
33	cent(xc1,yc1) = fread( fid, [snx,sny], 'float64' );
34	fclose(fid);
35	end
36	end