matlab_class/gcmfaces_IO/grid_load.m

function []=grid_load(varargin);

global mygrid;

if isempty(mygrid); 

if nargin==1; dirGrid=varargin{1}; nFaces=5;
elseif nargin==2; dirGrid=varargin{1}; nFaces=varargin{2};
else; dirGrid='/net/altix3700/raid4/gforget/mysetups/ecco_v4/RUNS/GRIDmds_90x50/'; nFaces=5;
end;

mygrid.dirGrid=dirGrid;

%list0={'XC','XG','YC','YG','RC','RF','RAC','DRC','DRF',...
%        'DXC','DXG','DYC','DYG','hFacC','hFacS','hFacW','Depth'};
%list0={'AngleCS','AngleSN','Depth','DRC','DRF','DXC','DXG','DYC','DYG',...
%       'hFacC','hFacS','hFacW','maskCtrlC','maskCtrlS','maskCtrlW',...
%       'PHrefC','PHrefF','RAC','RAS','RAW','RAZ','RC','RF','XC','XG','YC','YG'};

list0={'XC','XG','YC','YG','RAC','RAZ','DXC','DXG','DYC','DYG',...
       'hFacC','hFacS','hFacW','Depth','AngleCS','AngleSN'};
for iFld=1:length(list0);
   eval(['mygrid.' list0{iFld} '=rdmds2gcmfaces([dirGrid ''' list0{iFld} '*''],[],nFaces);']);
end;

list0={'RC','RF','DRC','DRF'};
for iFld=1:length(list0);
   eval(['mygrid.' list0{iFld} '=squeeze(rdmds([dirGrid ''' list0{iFld} '*'']));']);
end;

mygrid.hFacCsurf=mygrid.hFacC;
for ff=1:mygrid.hFacC.nFaces; mygrid.hFacCsurf{ff}=mygrid.hFacC{ff}(:,:,1); end;

mskC=mygrid.hFacC; mskC(mskC==0)=NaN; mskC(mskC>0)=1; mygrid.mskC=mskC;
mskW=mygrid.hFacW; mskW(mskW==0)=NaN; mskW(mskW>0)=1; mygrid.mskW=mskW;
mskS=mygrid.hFacS; mskS(mskS==0)=NaN; mskS(mskS>0)=1; mygrid.mskS=mskS;

end;

1	gforget	1.1	function []=grid_load(varargin);
2
3			global mygrid;
4
5			if isempty(mygrid);
6
7			if nargin==1; dirGrid=varargin{1}; nFaces=5;
8			elseif nargin==2; dirGrid=varargin{1}; nFaces=varargin{2};
9			else; dirGrid='/net/altix3700/raid4/gforget/mysetups/ecco_v4/RUNS/GRIDmds_90x50/'; nFaces=5;
10			end;
11
12	gforget	1.3	mygrid.dirGrid=dirGrid;
13
14	gforget	1.2	%list0={'XC','XG','YC','YG','RC','RF','RAC','DRC','DRF',...
15			% 'DXC','DXG','DYC','DYG','hFacC','hFacS','hFacW','Depth'};
16	gforget	1.1	%list0={'AngleCS','AngleSN','Depth','DRC','DRF','DXC','DXG','DYC','DYG',...
17			% 'hFacC','hFacS','hFacW','maskCtrlC','maskCtrlS','maskCtrlW',...
18			% 'PHrefC','PHrefF','RAC','RAS','RAW','RAZ','RC','RF','XC','XG','YC','YG'};
19
20			list0={'XC','XG','YC','YG','RAC','RAZ','DXC','DXG','DYC','DYG',...
21			'hFacC','hFacS','hFacW','Depth','AngleCS','AngleSN'};
22			for iFld=1:length(list0);
23			eval(['mygrid.' list0{iFld} '=rdmds2gcmfaces([dirGrid ''' list0{iFld} '*''],[],nFaces);']);
24			end;
25
26			list0={'RC','RF','DRC','DRF'};
27			for iFld=1:length(list0);
28	gforget	1.2	eval(['mygrid.' list0{iFld} '=squeeze(rdmds([dirGrid ''' list0{iFld} '*'']));']);
29	gforget	1.1	end;
30
31			mygrid.hFacCsurf=mygrid.hFacC;
32			for ff=1:mygrid.hFacC.nFaces; mygrid.hFacCsurf{ff}=mygrid.hFacC{ff}(:,:,1); end;
33
34	gforget	1.2	mskC=mygrid.hFacC; mskC(mskC==0)=NaN; mskC(mskC>0)=1; mygrid.mskC=mskC;
35			mskW=mygrid.hFacW; mskW(mskW==0)=NaN; mskW(mskW>0)=1; mygrid.mskW=mskW;
36			mskS=mygrid.hFacS; mskS(mskS==0)=NaN; mskS(mskS>0)=1; mygrid.mskS=mskS;
37
38	gforget	1.1	end;
39