global_ocean.128x64x15/diags_matlab/mit_loadglobaldata.m

% m-file: mit_loadglobaldata.m
% read all the data files for the 4x4 global run and take time averages
% of it
precision = mit_getparm('data','readBinaryPrec');
if isempty(precision); precision = 32.; end
if precision == 32
  acc = 'real*4';
elseif precision == 64;
  acc = 'real*8';
else
  error('readBinaryPrec contains unknown precision')
end

nx = grd.nx;
ny = grd.ny;
nz = grd.nz;

% load restoring fields and compute annual mean
thetaClimFile = mit_getparm('data','thetaClimFile');
if isempty(thetaClimFile) | thetaClimFile == ' ';
  tdmean = repmat(NaN,size(grd.hfacc(:,:,1)));
else
  tdmean = squeeze(mean(mit_readfield(thetaClimFile,[nx ny 12],acc),3));
end
saltClimFile = mit_getparm('data','saltClimFile');
if isempty(saltClimFile) | saltClimFile == ' ';
  sdmean = repmat(NaN,size(grd.hfacc(:,:,1)));
else
  sdmean = squeeze(mean(mit_readfield(saltClimFile,[nx ny 12],acc),3));
end
EmPmRFile = mit_getparm('data','EmPmRFile');
if isempty(EmPmRFile) | EmPmRFile == ' ';
  empr = repmat(NaN,nx,ny);
else
  empr = squeeze(mean(mit_readfield(EmPmRFile,[nx ny 12],acc),3));
end
empr = change(empr,'==',0,NaN);
surfQFile = mit_getparm('data','surfQFile');
if isempty(surfQFile) | surfQFile == ' ';
  qnet = repmat(NaN,nx,ny);
else
  qnet = squeeze(mean(mit_readfield(surfQFile,[nx ny 12],acc),3));
end
qnet = change(qnet,'==',0,NaN);
thetaFile = mit_getparm('data','hydrogThetaFile');
if isempty(thetaFile) | thetaFile == ' ';
  tdatamean = repmat(NaN,size(grd.hfacc));
else
  tdatamean = squeeze(mean(mit_readfield(thetaFile,[nx ny nz 12],acc),4));
end
msgbin = dir('sal.bin');
saltFile = mit_getparm('data','hydrogSaltFile');
if isempty(saltFile) | saltFile == ' ';
  sdatamean = repmat(NaN,size(grd.hfacc));
else
  sdatamean = squeeze(mean(mit_readfield(saltFile,[nx ny nz 12],acc),4));
end

if ~strcmp(grd.buoyancy,'OCEANIC');
  sdatamean = sdatamean(:,:,end:-1:1);
  tdatamean = tdatamean(:,:,end:-1:1);
end
1	cnh	1.1	% m-file: mit_loadglobaldata.m
2			% read all the data files for the 4x4 global run and take time averages
3			% of it
4			precision = mit_getparm('data','readBinaryPrec');
5			if isempty(precision); precision = 32.; end
6			if precision == 32
7			acc = 'real*4';
8			elseif precision == 64;
9			acc = 'real*8';
10			else
11			error('readBinaryPrec contains unknown precision')
12			end
13
14			nx = grd.nx;
15			ny = grd.ny;
16			nz = grd.nz;
17
18			% load restoring fields and compute annual mean
19			thetaClimFile = mit_getparm('data','thetaClimFile');
20			if isempty(thetaClimFile) \| thetaClimFile == ' ';
21			tdmean = repmat(NaN,size(grd.hfacc(:,:,1)));
22			else
23			tdmean = squeeze(mean(mit_readfield(thetaClimFile,[nx ny 12],acc),3));
24			end
25			saltClimFile = mit_getparm('data','saltClimFile');
26			if isempty(saltClimFile) \| saltClimFile == ' ';
27			sdmean = repmat(NaN,size(grd.hfacc(:,:,1)));
28			else
29			sdmean = squeeze(mean(mit_readfield(saltClimFile,[nx ny 12],acc),3));
30			end
31			EmPmRFile = mit_getparm('data','EmPmRFile');
32			if isempty(EmPmRFile) \| EmPmRFile == ' ';
33			empr = repmat(NaN,nx,ny);
34			else
35			empr = squeeze(mean(mit_readfield(EmPmRFile,[nx ny 12],acc),3));
36			end
37			empr = change(empr,'==',0,NaN);
38			surfQFile = mit_getparm('data','surfQFile');
39			if isempty(surfQFile) \| surfQFile == ' ';
40			qnet = repmat(NaN,nx,ny);
41			else
42			qnet = squeeze(mean(mit_readfield(surfQFile,[nx ny 12],acc),3));
43			end
44			qnet = change(qnet,'==',0,NaN);
45			thetaFile = mit_getparm('data','hydrogThetaFile');
46			if isempty(thetaFile) \| thetaFile == ' ';
47			tdatamean = repmat(NaN,size(grd.hfacc));
48			else
49			tdatamean = squeeze(mean(mit_readfield(thetaFile,[nx ny nz 12],acc),4));
50			end
51			msgbin = dir('sal.bin');
52			saltFile = mit_getparm('data','hydrogSaltFile');
53			if isempty(saltFile) \| saltFile == ' ';
54			sdatamean = repmat(NaN,size(grd.hfacc));
55			else
56			sdatamean = squeeze(mean(mit_readfield(saltFile,[nx ny nz 12],acc),4));
57			end
58
59			if ~strcmp(grd.buoyancy,'OCEANIC');
60			sdatamean = sdatamean(:,:,end:-1:1);
61			tdatamean = tdatamean(:,:,end:-1:1);
62			end