Documentation of mquiv

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

Function Synopsis

[h, a, c] = mquiv(aray1,aray2,len,mstyle,origin,lab)

Help text

  mquiv(xdat, ydat, len, 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 mquiv

function [h, a, c] = mquiv(aray1,aray2,len,mstyle,origin,lab)

cla

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

global XAX YAX FRAME x_coasts y_coasts;

[xk, yk] = keep_var(FRAME, XAX, YAX);

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

if (FRAME(2)-FRAME(1) == 360); 
  aray1 = [aray1(yk,:), aray1(yk,1)];
  aray2 = [aray2(yk,:), aray2(yk,1)];
  XAX2 = [XAX;XAX(1)+360];
else
  XAX2 = XAX(xk);
  aray1 = aray1(yk, xk);
  aray2 = aray2(yk, xk);
end

YAX2 = YAX(yk) * ones(1, length(XAX2));
XAX2 = ones(size(YAX(yk), 1), 1) * XAX2';
[h,a,b,c] = quiverm2(YAX2, XAX2, aray2, aray1, 'k', len);

plotm(a(:,[2 4]), a(:,[1 3]), '-k');
textm(a(1,2), a(1,1)+4, 'Max');

yc = y_coasts; xc = x_coasts;
yc(find(yc <= FRAME(3))) = NaN * ones(size(find(yc <= FRAME(3))));
xc(find(yc <= FRAME(3))) = NaN * ones(size(find(yc <= FRAME(3))));

plotm(yc, xc, '-');

if ~strcmp(lower(mstyle), 'stereo')
  setm(gca, 'maplonlimit', FRAME(1:2), 'maplatlimit', [b(3) FRAME(4)])
  setm(gca, 'grid', 'off')
end