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','RAC','RAZ','DXC','DXG','DYC','DYG',...
       'hFacC','hFacS','hFacW','Depth'};
for iFld=1:length(list0);
   eval(['mygrid.' list0{iFld} '=rdmds2gcmfaces([dirGrid ''' list0{iFld} '*''],[],nFaces);']);
end;

list0={'AngleCS','AngleSN'};
test0=~isempty(dir([dirGrid 'AngleCS*']));
if test0;
  for iFld=1:length(list0);
     eval(['mygrid.' list0{iFld} '=rdmds2gcmfaces([dirGrid ''' list0{iFld} '*''],[],nFaces);']);
  end;
else;
  warning('\n AngleCS/AngleSN not found; set to 1/0 assuming lat/lon grid.\n');
  mygrid.AngleCS=mygrid.XC; mygrid.AngleCS(:)=1; 
  mygrid.AngleSN=mygrid.XC; mygrid.AngleSN(:)=0;
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.1	list0={'XC','XG','YC','YG','RAC','RAZ','DXC','DXG','DYC','DYG',...
15	gforget	1.4	'hFacC','hFacS','hFacW','Depth'};
16	gforget	1.1	for iFld=1:length(list0);
17			eval(['mygrid.' list0{iFld} '=rdmds2gcmfaces([dirGrid ''' list0{iFld} '*''],[],nFaces);']);
18			end;
19
20	gforget	1.4	list0={'AngleCS','AngleSN'};
21			test0=~isempty(dir([dirGrid 'AngleCS*']));
22			if test0;
23			for iFld=1:length(list0);
24			eval(['mygrid.' list0{iFld} '=rdmds2gcmfaces([dirGrid ''' list0{iFld} '*''],[],nFaces);']);
25			end;
26			else;
27			warning('\n AngleCS/AngleSN not found; set to 1/0 assuming lat/lon grid.\n');
28			mygrid.AngleCS=mygrid.XC; mygrid.AngleCS(:)=1;
29			mygrid.AngleSN=mygrid.XC; mygrid.AngleSN(:)=0;
30			end;
31
32	gforget	1.1	list0={'RC','RF','DRC','DRF'};
33			for iFld=1:length(list0);
34	gforget	1.2	eval(['mygrid.' list0{iFld} '=squeeze(rdmds([dirGrid ''' list0{iFld} '*'']));']);
35	gforget	1.1	end;
36
37			mygrid.hFacCsurf=mygrid.hFacC;
38			for ff=1:mygrid.hFacC.nFaces; mygrid.hFacCsurf{ff}=mygrid.hFacC{ff}(:,:,1); end;
39
40	gforget	1.2	mskC=mygrid.hFacC; mskC(mskC==0)=NaN; mskC(mskC>0)=1; mygrid.mskC=mskC;
41			mskW=mygrid.hFacW; mskW(mskW==0)=NaN; mskW(mskW>0)=1; mygrid.mskW=mskW;
42			mskS=mygrid.hFacS; mskS(mskS==0)=NaN; mskS(mskS>0)=1; mygrid.mskS=mskS;
43
44	gforget	1.1	end;
45