Global Index (short | long) | Local contents | Local Index (short | long)
shade_ccoef(dat);
shade_ccoef(data); This function pcolors the NaN's of data
This function calls | This function is called by |
---|---|
function shade_ccoef(dat); if iselement(size(dat), 1); dat = squeeze(dat); 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)]; 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, -1*ones(size(data2)), data2); colormap([0.9; 0.2]*[1 1 1]); shading flat box on hold off;