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
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	mymin=min(min(AAA));
12	mymax=max(max(AAA));
13
14	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	grph_CS(AAA,xcs,ycs,xcg,ycg,-0.E-2,1.E-2,0)
22
23	title( ['cubed-spere adv.-diff. adjoint: adxx\_theta', ...
24	' min/max = ',num2str(mymin,3), ' / ', num2str(mymax,3)] )
25
26	filename= [ 'cs_advdiff_adjoint_adxx_theta' ]
27	%%%eval([ 'print -depsc ', filename ])
28	%%%eval([ 'print -djpeg ', filename ])