Documentation of mshade

Global Index (short | long) | Local contents | Local Index (short | long)

Function Synopsis

mshade(aray,levs,mstyle,origin)

Help text

  mshade(data, levs, mstyle, origin)
  
  You may omit mstyle and origin and mcont will contour
  using polar stereographic with the North Pole as the
  origin.

Cross-Reference Information

Listing of function mshade

function mshade(aray,levs,mstyle,origin)

if nargin == 2
  mstyle = 'stereo';
  origin = [90 270];
elseif nargin == 3
  if isstr(mstyle);
    if strcmp(mstyle, 'stereo')
      origin = [90 270];
    else 
      origin = [0 180];
    end
  elseif (length(mstyle) == 2);
    origin = mstyle; mstyle = 'stereo'; lab = 0;
  end
end

global XAX YAX FRAME x_coasts y_coasts;

axesm('MapProjection',mstyle,'origin',origin,'grid','on')

[xk, yk] = keep_var(FRAME, XAX, YAX);
if (FRAME(2)-FRAME(1) == 360); 
  aray2 = [aray(yk,:) aray(yk,1)];
  XAX2 = [XAX;XAX(1)+360];
else
  aray2 = aray(yk, xk);
  XAX2 = XAX(xk);
end
YAX2 = YAX(yk);
XAX3 = ones(size(YAX2)) * XAX2';
YAX2 = YAX2 * ones(size(XAX2'));

%if length(levs) == 2;
%  colormap([1 1 1; .9 .9 .9]);
%  aray2(find(aray2 < levs(1))) = zeros(size(find(aray2 < levs(1))));
%  aray2(find(aray2 >= levs(1))) = levs(1) * ones(size(find(aray2 >= levs(1))));
%else
%  shade_vect = (levs - min(levs)) / (max(levs) - min(levs));
%  colormap([shade_vect' shade_vect' shade_vect']);
%end

%pcolorm(YAX2,XAX3,aray2);
%surfacem(YAX2,XAX3,aray2,levs);
greyshade(XAX3,YAX2,aray2,levs);

%  [b] = plotm(y_coasts,x_coasts,'-');
%  set(b,'linewidth',0.5);

%if ~strcmp(lower(mstyle), 'stereo')
  setm(gca, 'maplonlimit', FRAME(1:2), 'maplatlimit', FRAME(3:4), ...
            'flonlimit',   FRAME(1:2), 'flatlimit',   FRAME(3:4))
%end