regrid/scrip_examples/split_CS32_faces.m

%
%  Ed Hill
%
%  Split a C32 grid along its faces and produce ASCII output.
%


nc = netcdf('rmp_CS32_to_LL128x64_conserv.nc','nowrite');

isrc = nc{'src_address'}(:);
idst = nc{'dst_address'}(:);
wts = nc{'remap_matrix'}(:);

nc = close(nc);

icull = find(wts(:,1) < 1.0e-18);
ikeep = find(wts(:,1) > 1.0e-18);
isrck = isrc(ikeep);
idstk = idst(ikeep);
wtsk = wts(ikeep,:);

pface = {};
fi = 1;
for fi = 1:6
  fname = sprintf('CS32_LL128x64_%03d.regrid.ascii',fi);
  fid = fopen(fname,'wt');
  imin = 32^2 * (fi - 1) + 1;
  imax = 32^2 * fi;
  pface(fi).inds = find(imin<=isrck & isrck<=imax);
  for j = 1:length(pface(fi).inds)
    ii = pface(fi).inds(j);
    fprintf(fid,' %10d %10d %28.22E\n', ...
            (isrck(ii) - imin + 1),idstk(ii),wtsk(ii,1));
  end
  fclose(fid);
end

1	edhill	1.1	%
2			% Ed Hill
3			%
4			% Split a C32 grid along its faces and produce ASCII output.
5			%
6
7
8			nc = netcdf('rmp_CS32_to_LL128x64_conserv.nc','nowrite');
9
10			isrc = nc{'src_address'}(:);
11			idst = nc{'dst_address'}(:);
12			wts = nc{'remap_matrix'}(:);
13
14			nc = close(nc);
15
16			icull = find(wts(:,1) < 1.0e-18);
17			ikeep = find(wts(:,1) > 1.0e-18);
18			isrck = isrc(ikeep);
19			idstk = idst(ikeep);
20			wtsk = wts(ikeep,:);
21
22			pface = {};
23			fi = 1;
24			for fi = 1:6
25			fname = sprintf('CS32_LL128x64_%03d.regrid.ascii',fi);
26			fid = fopen(fname,'wt');
27			imin = 32^2 * (fi - 1) + 1;
28			imax = 32^2 * fi;
29			pface(fi).inds = find(imin<=isrck & isrck<=imax);
30			for j = 1:length(pface(fi).inds)
31			ii = pface(fi).inds(j);
32			fprintf(fid,' %10d %10d %28.22E\n', ...
33			(isrck(ii) - imin + 1),idstk(ii),wtsk(ii,1));
34			end
35			fclose(fid);
36			end
37