Documentation of dc_color

Global Index (all files) (short | long) | Local Index (files in subdir) (short | long)

Function Synopsis

hh = dc2(dat);

Help text

  dc2(data);

  This function pcolors the NaN's of data

Cross-Reference Information

Listing of function dc_color

function hh = dc2(dat);

dat = squeeze(dat);

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)];
  XAX2 = [XAX; (XAX(1)+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) & (FRAME(3) == -90);
    yticks = sort( [ 0:-gridy:(y0+gridy) gridy:gridy:(y1-gridy) ] );
  elseif FRAME(4) == 90;
    yticks = sort( [ 0:-gridy:y0 gridy:gridy:(y1 - gridy) ] );
  elseif FRAME(3) == -90;
    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;

hh = surface(XAX2, YAX2, ones(size(data2)), data2);
set(hh, 'FaceColor', [0.6 0.6 0.6]);
%colormap([0.6]*[1 1 1]);
shading flat
box on

hold off;