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