Documentation of drawmap2

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

Function Synopsis

hh = dc2(dat, colr, pn);

Help text

  dc2(data);

  This function pcolors the NaN's of data

Cross-Reference Information

Listing of function drawmap2

function hh = dc2(dat, colr, pn);

dat = squeeze(dat);

if nargin <= 1; colr = 0.35; end;
if nargin <= 2; pn = 1; end;

data = repmat(NaN, size(dat));
data(find(isnan(dat))) = 1;

global XAX YAX FRAME

if (size(YAX,1)==1); YAX = YAX'; end;
if (size(XAX,1)==1); XAX = XAX'; end;

[xk, yk] = keep_var(FRAME, XAX, YAX);
%xk = 1:length(XAX);
%yk = 1:length(YAX);

if (FRAME(2)-FRAME(1) == 360); 
  data2 = [data(yk,:) data(yk,1:2)];
  XAX2 = [XAX; (XAX(1:2)+360)];
  YAX2 = YAX(yk);
else
  data2 = data(yk,xk);
  XAX2 = XAX(xk);
  YAX2 = YAX(yk);
end
XAX2 = XAX2 - mean(diff(XAX2))/2;
YAX2 = YAX2 - mean(diff(YAX2))/2;

hold on;
cax = newplot;
hold_state = ishold;

global FRAME COASTWIDTH COASTTYPE COASTGRID GRDX_SPACING GRDY_SPACING
global USE_LANDMASK

%gridx = GRDX_SPACING;
%gridy = GRDY_SPACING;
%lgrd=':k';
%  xl=get(gca,'xlim'); yl=get(gca,'ylim');
%  x0=xl(1); x1=xl(2); y0=yl(1); y1=yl(2);
%  if ~USE_LANDMASK
%  xticks = sort([ 0:-gridx:x0 gridx:gridx:x1 ]);
%  if (FRAME(4) == 90.1) & (FRAME(3) == -90.1);
%    yticks = sort( [ 0:-gridy:(y0+gridy) gridy:gridy:(y1-gridy) ] );
%  elseif FRAME(4) == 90.1;
%    yticks = sort( [ 0:-gridy:y0 gridy:gridy:(y1 - gridy) ] );
%  elseif FRAME(3) == -90.1;
%    yticks = sort( [ 0:-gridy:(y0 + gridy) gridy:gridy:y1 ] );
%  else
%    yticks = sort( [ 0:-gridy:y0 gridy:gridy:y1 ] );
%  end
%  manual_ticks('x',xticks,label_eastwest(xticks));
%  manual_ticks('y',yticks,label_northsouth(yticks));
%  end
%  if lgrd
%    if FRAME(4) == 90 & FRAME(3) == -90;
%      vline(xticks,lgrd, [min(yticks)-gridy/2 max(yticks)+gridy/2]); 
%    elseif FRAME(4) == 90;
%      vline(xticks,lgrd, [min(yticks) max(yticks)+girdy/2]);
%    elseif FRAME(3) == -90;
%      vline(xticks,lgrd, [min(yticks)-gridy/2 max(yticks)]);
%    else
%      vline(xticks,lgrd,[min(yticks) max(yticks)]);
%    end
%    hline(yticks,lgrd);
%  end;

XAX2 = ones(length(YAX2), 1)*XAX2';
YAX2 = YAX2*ones(1, size(XAX2, 2));

hh = surfacem(YAX2, XAX2, data2, 100*pn*ones(size(data2)));
colormap(colr*[1 1 1]);
%shading flat
%box on

hold off;