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