| 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.2 |
%list0={'XC','XG','YC','YG','RC','RF','RAC','DRC','DRF',... |
| 13 |
|
|
% 'DXC','DXG','DYC','DYG','hFacC','hFacS','hFacW','Depth'}; |
| 14 |
gforget |
1.1 |
%list0={'AngleCS','AngleSN','Depth','DRC','DRF','DXC','DXG','DYC','DYG',... |
| 15 |
|
|
% 'hFacC','hFacS','hFacW','maskCtrlC','maskCtrlS','maskCtrlW',... |
| 16 |
|
|
% 'PHrefC','PHrefF','RAC','RAS','RAW','RAZ','RC','RF','XC','XG','YC','YG'}; |
| 17 |
|
|
|
| 18 |
|
|
list0={'XC','XG','YC','YG','RAC','RAZ','DXC','DXG','DYC','DYG',... |
| 19 |
|
|
'hFacC','hFacS','hFacW','Depth','AngleCS','AngleSN'}; |
| 20 |
|
|
for iFld=1:length(list0); |
| 21 |
|
|
eval(['mygrid.' list0{iFld} '=rdmds2gcmfaces([dirGrid ''' list0{iFld} '*''],[],nFaces);']); |
| 22 |
|
|
end; |
| 23 |
|
|
|
| 24 |
|
|
list0={'RC','RF','DRC','DRF'}; |
| 25 |
|
|
for iFld=1:length(list0); |
| 26 |
gforget |
1.2 |
eval(['mygrid.' list0{iFld} '=squeeze(rdmds([dirGrid ''' list0{iFld} '*'']));']); |
| 27 |
gforget |
1.1 |
end; |
| 28 |
|
|
|
| 29 |
|
|
mygrid.hFacCsurf=mygrid.hFacC; |
| 30 |
|
|
for ff=1:mygrid.hFacC.nFaces; mygrid.hFacCsurf{ff}=mygrid.hFacC{ff}(:,:,1); end; |
| 31 |
|
|
|
| 32 |
gforget |
1.2 |
mskC=mygrid.hFacC; mskC(mskC==0)=NaN; mskC(mskC>0)=1; mygrid.mskC=mskC; |
| 33 |
|
|
mskW=mygrid.hFacW; mskW(mskW==0)=NaN; mskW(mskW>0)=1; mygrid.mskW=mskW; |
| 34 |
|
|
mskS=mygrid.hFacS; mskS(mskS==0)=NaN; mskS(mskS>0)=1; mygrid.mskS=mskS; |
| 35 |
|
|
|
| 36 |
gforget |
1.1 |
end; |
| 37 |
|
|
|
Parent Directory
| 
Revision Log
| 
Revision Graph