cs_adj_advdiff/results_ad/mitgcm_single.m

ypre  = './'

ydate = '0000000000'
%
yvar = 'adxx_theta'
%
level=1
%
A=rdmds([ ypre, '/', yvar, '.', ydate ]);
size(A)
%
yexp = 'cub'

if yexp == '1x1'
%
nlonw = 0.5
nlone = 359.5
nlonres = 1
nlats = -79.5
nlatn = 79.5
nlatres = 1
%
elseif yexp == '2x2'
%
nlonw = 1
nlone = 359
nlonres = 2
nlats = -79
nlatn = 79
nlatres = 2
%
elseif yexp == '4x4'
%
nlonw = 2
nlone = 358
nlonres = 4
nlats = -78
nlatn = 78
nlatres = 4
%
elseif yexp == 'cub'
%
nlonres = 1
nlonw = 0.5
nlone = nlonw + 191.*nlonres
nlatres = 1
nlats = 0.5
nlatn = nlats + 31.*nlatres
%
elseif yexp == 'med'
%
nlonres = 0.25
nlonw = -8.875
nlone = nlonw + 183.*nlonres
nlatres = 0.25
nlats = 30.125
nlatn = nlats + 63.*nlatres
%
end
%
mlon  = (nlone-nlonw)/nlonres + 1
mlat  = (nlatn-nlats)/nlatres + 1
mvert = 15
%
%--------------------------------------------------------------
%
%%% set(gcf,'PaperOrientation','portrait')
%%% set(gcf,'PaperPosition',[0.5 0.5 7.5 10.])
set(gcf,'PaperOrientation','landscape')
set(gcf,'PaperPosition',[0.5 0.5 10. 7.5])
%
temp=squeeze(A(:,:,level));
%
if yvar == 'T'
  cmin = -16.
  cmax = 37.
  cstep = 4.
elseif yvar == 'S'
  cmin=34;
  cmax=38;
  cstep=0.5;
elseif yvar == 'U'
  cmin=-0.7;
  cmax=0.7;
  cstep=0.1;
elseif yvar == 'V'
  cmin=-0.2;
  cmax=0.2;
  cstep=0.02;
else
  cmin=min(min(temp(:,:)));
  cmax=max(max(temp(:,:)));
  cstep=(cmax-cmin)/10.
end
%
%cmin = -0.0
%cmax =  0.001
%cstep=(cmax-cmin)/10.
%
    missing = 0.;
    [xxx,yyy] = find(temp==missing);
    nnn = length(xxx);
    for iii = 1:nnn
      temp(xxx(iii),yyy(iii)) = NaN;
    end
%
pcolor([nlonw:nlonres:nlone],[nlats:nlatres:nlatn],temp');
shading flat
%
caxis([cmin cmax]);
colorbar
%
%%%hold on
%%%[c,h]=contour([nlonw:nlonres:nlone],[nlats:nlatres:nlatn], ...
%%%                 temp',[cmin:cstep:cmax],'w');
%%%caxis([cmin cmax]);
%%%clabel(c,[cmin:cstep:cmax]);
%
tempmin=min(min(temp(:,:)));
tempmax=max(max(temp(:,:)));
tempmin
tempmax
%
ytime = eval('date');
title( [ ypre, '.', yvar, '.', ydate, '.', yexp, '.lev', ...
   int2str(level) , ...
   '  min/max=', num2str(tempmin,3), ' / ', num2str(tempmax,3) ] )
%

filename= [ ypre, '.', yvar, '.', ydate, '.', yexp, '.lev', int2str(level), '.eps' ]
%%%eval([ 'print -depsc ', filename ])

1	heimbach	1.1	ypre = './'
2
3			ydate = '0000000000'
4			%
5			yvar = 'adxx_theta'
6			%
7			level=1
8			%
9			A=rdmds([ ypre, '/', yvar, '.', ydate ]);
10			size(A)
11			%
12			yexp = 'cub'
13
14			if yexp == '1x1'
15			%
16			nlonw = 0.5
17			nlone = 359.5
18			nlonres = 1
19			nlats = -79.5
20			nlatn = 79.5
21			nlatres = 1
22			%
23			elseif yexp == '2x2'
24			%
25			nlonw = 1
26			nlone = 359
27			nlonres = 2
28			nlats = -79
29			nlatn = 79
30			nlatres = 2
31			%
32			elseif yexp == '4x4'
33			%
34			nlonw = 2
35			nlone = 358
36			nlonres = 4
37			nlats = -78
38			nlatn = 78
39			nlatres = 4
40			%
41			elseif yexp == 'cub'
42			%
43			nlonres = 1
44			nlonw = 0.5
45			nlone = nlonw + 191.*nlonres
46			nlatres = 1
47			nlats = 0.5
48			nlatn = nlats + 31.*nlatres
49			%
50			elseif yexp == 'med'
51			%
52			nlonres = 0.25
53			nlonw = -8.875
54			nlone = nlonw + 183.*nlonres
55			nlatres = 0.25
56			nlats = 30.125
57			nlatn = nlats + 63.*nlatres
58			%
59			end
60			%
61			mlon = (nlone-nlonw)/nlonres + 1
62			mlat = (nlatn-nlats)/nlatres + 1
63			mvert = 15
64			%
65			%--------------------------------------------------------------
66			%
67			%%% set(gcf,'PaperOrientation','portrait')
68			%%% set(gcf,'PaperPosition',[0.5 0.5 7.5 10.])
69			set(gcf,'PaperOrientation','landscape')
70			set(gcf,'PaperPosition',[0.5 0.5 10. 7.5])
71			%
72			temp=squeeze(A(:,:,level));
73			%
74			if yvar == 'T'
75			cmin = -16.
76			cmax = 37.
77			cstep = 4.
78			elseif yvar == 'S'
79			cmin=34;
80			cmax=38;
81			cstep=0.5;
82			elseif yvar == 'U'
83			cmin=-0.7;
84			cmax=0.7;
85			cstep=0.1;
86			elseif yvar == 'V'
87			cmin=-0.2;
88			cmax=0.2;
89			cstep=0.02;
90			else
91			cmin=min(min(temp(:,:)));
92			cmax=max(max(temp(:,:)));
93			cstep=(cmax-cmin)/10.
94			end
95			%
96			%cmin = -0.0
97			%cmax = 0.001
98			%cstep=(cmax-cmin)/10.
99			%
100			missing = 0.;
101			[xxx,yyy] = find(temp==missing);
102			nnn = length(xxx);
103			for iii = 1:nnn
104			temp(xxx(iii),yyy(iii)) = NaN;
105			end
106			%
107			pcolor([nlonw:nlonres:nlone],[nlats:nlatres:nlatn],temp');
108			shading flat
109			%
110			caxis([cmin cmax]);
111			colorbar
112			%
113			%%%hold on
114			%%%[c,h]=contour([nlonw:nlonres:nlone],[nlats:nlatres:nlatn], ...
115			%%% temp',[cmin:cstep:cmax],'w');
116			%%%caxis([cmin cmax]);
117			%%%clabel(c,[cmin:cstep:cmax]);
118			%
119			tempmin=min(min(temp(:,:)));
120			tempmax=max(max(temp(:,:)));
121			tempmin
122			tempmax
123			%
124			ytime = eval('date');
125			title( [ ypre, '.', yvar, '.', ydate, '.', yexp, '.lev', ...
126			int2str(level) , ...
127			' min/max=', num2str(tempmin,3), ' / ', num2str(tempmax,3) ] )
128			%
129
130			filename= [ ypre, '.', yvar, '.', ydate, '.', yexp, '.lev', int2str(level), '.eps' ]
131			%%%eval([ 'print -depsc ', filename ])
132