%----------------------------------------------------------------------- % [smin,smax,ibw]=colorband_fe(in,x,y,bnd,scalar,smin,smax,ibw) % This function creates a color banded plot using % Matlab5.1 commands and the jet color map. %----------------------------------------------------------------------- function [smin,smax,ibw]=colorband_fe(in,x,y,bnd,scalar,smin,smax,ibw) % % Echo scalar range to screen % fprintf(1,'Scalar Range: %f to %f\n',min(min(scalar)),max(max(scalar))) % % Set color banding settings (if 5 arguments are sent to function) % if nargin == 5 smin=input('Enter min contour level desired: '); smax=input('Enter max contour level desired: '); ibw =input('Enter the contour interval: '); end % % Set color bands % cval=smin:ibw:smax; nband=(2*ibw+smax-smin)/ibw; cmapjet=jet; size(cmapjet); njet=ans(1); clear cmap; cmap(1,:)=[0 0 0]; for i=2:nband-1 cmap(i,:)=cmapjet(1+round((njet-1)*(i-2)/(nband-3)),:); end cmap(nband,:)=[1 1 1]; % % Generate plot % scalara=max(scalar,smin-ibw); scalara=min(scalar,smax+ibw); hpb=plotbnd4(x,y,bnd);set(hpb,'color','k'); hl=lcontour4(in,x,y,scalara,cval); hc=colormeshm(in,x,y,scalara);colormap(cmap);caxis([smin-ibw smax+ibw]); hb=colorbar;set(hb,'ytick',cval);set(hb,'ticklength',[0.05 0.025]); set(hb,'FontSize',1.0); for l=1:nband-1 if cval(l)>=min(scalara) & cval(l)<=max(scalara) set(hl(l),'color','k'); end end clear cmap;