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	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