enderton/Diagnostics/DiagLoadGridData.m

function [XC,XG,YC,YG,Ylat,ZC,ZF,RAC,drC,drF,HFacC,HFacW,HFacS,dxG,dyG,dxC,dyC] = ...
    DiagLoadGridData(LoadGridData,grd,gdf,flu,GridSuffix,ZcordFile)

% Load grid data.

% If LoadGridData is turned on, load grid data from the model output files.
% If it is turned off, load it from the data file where saved from a
% previous loading of the raw data.  The file name (and path) of this file
% is set in the DiagGenParam file.
if LoadGridData
    
    % Load y-axis for CS zonal average calculation.
    load('ZonAvgLat');

        if isequal(gdf,'MDS')
                XC  = rdmds([grd,'/','XC' ]);
                YC  = rdmds([grd,'/','YC' ]);
                XG  = rdmds([grd,'/','XG' ]);
                YG  = rdmds([grd,'/','YG' ]);
                dxC = rdmds([grd,'/','DXC']);
                dyC = rdmds([grd,'/','DYC']);
                dxG = rdmds([grd,'/','DXG']);
                dyG = rdmds([grd,'/','DYG']);
                RAC = rdmds([grd,'/','RAC']);
                HFacC = rdmds([grd,'/','hFacC']);
                HFacS = rdmds([grd,'/','hFacS']);
                HFacW = rdmds([grd,'/','hFacW']);
        
        % Information for z-axis.
        if isequal(ZcordFile,'')
                load(['Zcord',flu,GridSuffix,'.mat']);
        else
            load(ZcordFile);
        end     
        
        elseif isequal(gdf,'MNC')
        data = rdmnc_mod([grd,'grid.*']);
                XC = data.XC;                           % [192x32 double]
                YC = data.YC;                           % [192x32 double]
                XG = data.XG; XG = XG(1:end-1,1:end-1); % [193x33 double] -> [192x32 double]
                YG = data.YG; YG = YG(1:end-1,1:end-1); % [193x33 double] -> [192x32 double]
                dxC = data.dxC;                         % [192x32 double]
                dyC = data.dyC;                         % [192x32 double]
                dxG = data.dxG;                         % [192x32 double]
                dyG = data.dyG;                         % [192x32 double]
        RAC = data.rA;                          % [192x32 double]
                HFacC = data.HFacC;                     % [192x32x15 double]
                HFacW = data.HFacW;                     % [192x32x15 double]
                HFacS = data.HFacS;                     % [192x32x15 double]       
                ZC = data.RC;
                ZF = data.RF;
                drC = data.drC;
                drF = data.drF;
    end
    
    save(['GridData',flu,GridSuffix,'.mat'],...
         'XC','XG','YC','YG','Ylat','ZC','ZF','RAC','drC','drF',...
         'HFacC','HFacW','HFacS','dxG','dyG','dxC','dyC');
     
% If LoadGridData is turned off, load it from the data file where saved
% from a previous loading of the raw data.  The file name (and path) of
% this file is set in the DiagGenParam file.
else
    load(['GridData',flu,GridSuffix,'.mat'])
end
1	function [XC,XG,YC,YG,Ylat,ZC,ZF,RAC,drC,drF,HFacC,HFacW,HFacS,dxG,dyG,dxC,dyC] = ...
2	DiagLoadGridData(LoadGridData,grd,gdf,flu,GridSuffix,ZcordFile)
3
4	% Load grid data.
5
6	% If LoadGridData is turned on, load grid data from the model output files.
7	% If it is turned off, load it from the data file where saved from a
8	% previous loading of the raw data. The file name (and path) of this file
9	% is set in the DiagGenParam file.
10	if LoadGridData
11
12	% Load y-axis for CS zonal average calculation.
13	load('ZonAvgLat');
14
15	if isequal(gdf,'MDS')
16	XC = rdmds([grd,'/','XC' ]);
17	YC = rdmds([grd,'/','YC' ]);
18	XG = rdmds([grd,'/','XG' ]);
19	YG = rdmds([grd,'/','YG' ]);
20	dxC = rdmds([grd,'/','DXC']);
21	dyC = rdmds([grd,'/','DYC']);
22	dxG = rdmds([grd,'/','DXG']);
23	dyG = rdmds([grd,'/','DYG']);
24	RAC = rdmds([grd,'/','RAC']);
25	HFacC = rdmds([grd,'/','hFacC']);
26	HFacS = rdmds([grd,'/','hFacS']);
27	HFacW = rdmds([grd,'/','hFacW']);
28
29	% Information for z-axis.
30	if isequal(ZcordFile,'')
31	load(['Zcord',flu,GridSuffix,'.mat']);
32	else
33	load(ZcordFile);
34	end
35
36	elseif isequal(gdf,'MNC')
37	data = rdmnc_mod([grd,'grid.*']);
38	XC = data.XC; % [192x32 double]
39	YC = data.YC; % [192x32 double]
40	XG = data.XG; XG = XG(1:end-1,1:end-1); % [193x33 double] -> [192x32 double]
41	YG = data.YG; YG = YG(1:end-1,1:end-1); % [193x33 double] -> [192x32 double]
42	dxC = data.dxC; % [192x32 double]
43	dyC = data.dyC; % [192x32 double]
44	dxG = data.dxG; % [192x32 double]
45	dyG = data.dyG; % [192x32 double]
46	RAC = data.rA; % [192x32 double]
47	HFacC = data.HFacC; % [192x32x15 double]
48	HFacW = data.HFacW; % [192x32x15 double]
49	HFacS = data.HFacS; % [192x32x15 double]
50	ZC = data.RC;
51	ZF = data.RF;
52	drC = data.drC;
53	drF = data.drF;
54	end
55
56	save(['GridData',flu,GridSuffix,'.mat'],...
57	'XC','XG','YC','YG','Ylat','ZC','ZF','RAC','drC','drF',...
58	'HFacC','HFacW','HFacS','dxG','dyG','dxC','dyC');
59
60	% If LoadGridData is turned off, load it from the data file where saved
61	% from a previous loading of the raw data. The file name (and path) of
62	% this file is set in the DiagGenParam file.
63	else
64	load(['GridData',flu,GridSuffix,'.mat'])
65	end