tutorial_global_oce_latlon/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
% $Header:  $
% $Name:  $

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