cs_adj_advdiff/results_ad/myplotcube.m


xcg=rdmds('XG');
ycg=rdmds('YG');
xcs=rdmds('XC');
ycs=rdmds('YC');

BBB=rdmds('adxx_theta.0000000000');

AAA=squeeze(BBB(:,:,1));

mymin=min(min(AAA));
mymax=max(max(AAA));

    missing = 0.;
    [xxx,yyy] = find(AAA==missing);
    nnn = length(xxx);
    for iii = 1:nnn
      AAA(xxx(iii),yyy(iii)) = NaN;
    end

grph_CS(AAA,xcs,ycs,xcg,ycg,-0.E-2,1.E-2,0)

title( ['cubed-spere adv.-diff. adjoint: adxx\_theta', ...
        ' min/max = ',num2str(mymin,3), ' / ', num2str(mymax,3)] )

filename= [ 'cs_advdiff_adjoint_adxx_theta' ]
%%%eval([ 'print -depsc ', filename ])
%%%eval([ 'print -djpeg ', filename ])
1	heimbach	1.1
2			xcg=rdmds('XG');
3			ycg=rdmds('YG');
4			xcs=rdmds('XC');
5			ycs=rdmds('YC');
6
7			BBB=rdmds('adxx_theta.0000000000');
8
9			AAA=squeeze(BBB(:,:,1));
10
11	heimbach	1.2	mymin=min(min(AAA));
12			mymax=max(max(AAA));
13
14	heimbach	1.1	missing = 0.;
15			[xxx,yyy] = find(AAA==missing);
16			nnn = length(xxx);
17			for iii = 1:nnn
18			AAA(xxx(iii),yyy(iii)) = NaN;
19			end
20
21	heimbach	1.2	grph_CS(AAA,xcs,ycs,xcg,ycg,-0.E-2,1.E-2,0)
22	heimbach	1.1
23	heimbach	1.2	title( ['cubed-spere adv.-diff. adjoint: adxx\_theta', ...
24			' min/max = ',num2str(mymin,3), ' / ', num2str(mymax,3)] )
25	heimbach	1.1
26			filename= [ 'cs_advdiff_adjoint_adxx_theta' ]
27			%%%eval([ 'print -depsc ', filename ])
28	heimbach	1.2	%%%eval([ 'print -djpeg ', filename ])