matlab_class/gcmfaces_IO/interp2nctiles.m

function []=interp2nctiles(indFiles);
% INTERP2CNTILES creates netcdf files from interpolated 
%   fields that were created by process2interp. The 
%   input vector process2interp specifies the file subset
%   to be processed (1:length(listInterp) by default).
%
% note: this routine overrides mygrid with mygrid_latlon 
%   after storing the reference grid in mygrid_orig

gcmfaces_global; global mygrid_orig;

dirModel='./';
if isempty(mygrid_orig);
grid_load; mygrid_orig=mygrid;
cd(dirModel);
end;

lon=[-179.75:0.5:179.75]; lat=[-89.75:0.5:89.75];
[lat,lon] = meshgrid(lat,lon);

mygrid_latlon.nFaces=1;
mygrid_latlon.dirGrid='none';
mygrid_latlon.fileFormat='straight';
mygrid_latlon.ioSize=size(lon);
mygrid_latlon.XC=gcmfaces({lon});
mygrid_latlon.YC=gcmfaces({lat});
mygrid_latlon.RC=mygrid.RC;
mygrid_latlon.RF=mygrid.RF;
mygrid_latlon.DRC=mygrid.DRC;
mygrid_latlon.DRF=mygrid.DRF;
mygrid_latlon.mskC=1+0*repmat(mygrid_latlon.XC,[1 1  length(mygrid.RC)]);
mygrid_latlon.RAC=[];
mygrid_latlon.gcm2facesFast=0;
mygrid_latlon.facesExpand=[];

[listInterp,listNot]=process2interp;
if isempty(who('indFiles')); 
  indFiles=[1:length(listInterp)];
end;

mygrid=mygrid_latlon;
for ii=indFiles;
tic; process2nctiles(dirModel,listInterp{ii},[]);
fprintf(['DONE: ' listInterp{ii} ' (in ' num2str(toc) 's)\n']);
end;
mygrid=mygrid_orig;

1	function []=interp2nctiles(indFiles);
2	% INTERP2CNTILES creates netcdf files from interpolated
3	% fields that were created by process2interp. The
4	% input vector process2interp specifies the file subset
5	% to be processed (1:length(listInterp) by default).
6	%
7	% note: this routine overrides mygrid with mygrid_latlon
8	% after storing the reference grid in mygrid_orig
9
10	gcmfaces_global; global mygrid_orig;
11
12	dirModel='./';
13	if isempty(mygrid_orig);
14	grid_load; mygrid_orig=mygrid;
15	cd(dirModel);
16	end;
17
18	lon=[-179.75:0.5:179.75]; lat=[-89.75:0.5:89.75];
19	[lat,lon] = meshgrid(lat,lon);
20
21	mygrid_latlon.nFaces=1;
22	mygrid_latlon.dirGrid='none';
23	mygrid_latlon.fileFormat='straight';
24	mygrid_latlon.ioSize=size(lon);
25	mygrid_latlon.XC=gcmfaces({lon});
26	mygrid_latlon.YC=gcmfaces({lat});
27	mygrid_latlon.RC=mygrid.RC;
28	mygrid_latlon.RF=mygrid.RF;
29	mygrid_latlon.DRC=mygrid.DRC;
30	mygrid_latlon.DRF=mygrid.DRF;
31	mygrid_latlon.mskC=1+0*repmat(mygrid_latlon.XC,[1 1 length(mygrid.RC)]);
32	mygrid_latlon.RAC=[];
33	mygrid_latlon.gcm2facesFast=0;
34	mygrid_latlon.facesExpand=[];
35
36	[listInterp,listNot]=process2interp;
37	if isempty(who('indFiles'));
38	indFiles=[1:length(listInterp)];
39	end;
40
41	mygrid=mygrid_latlon;
42	for ii=indFiles;
43	tic; process2nctiles(dirModel,listInterp{ii},[]);
44	fprintf(['DONE: ' listInterp{ii} ' (in ' num2str(toc) 's)\n']);
45	end;
46	mygrid=mygrid_orig;
47