vis/sf.m

   1 function sf(f,u,v,l)\r
   2     [m,n]=size(f);\r
   3     %f(f==0)=NaN;\r
   4     dx=6;dy=6;\r
   5     clf,hold off\r
   6     s=20;\r
   7     ix=1:ceil(m/s):m;\r
   8     iy=1:ceil(n/s):n;\r
   9     x=([0:m-1])*dx;\r
  10     y=([0:n-1])*dy;\r
  11     xh=([0:m-1]+0.5)*dx;\r
  12     yh=([0:n-1]+0.5)*dy;\r
  13     [xx,yy]=ndgrid(x,y);\r
  14     sc=2*min(m,n)/s;\r
  15     hsc=1e5; % scale for flux labels\r
  16     quiver(xx(ix,iy),yy(ix,iy),sc*u(ix,iy),sc*v(ix,iy),0,'k')\r
  17     hold on\r
  18     %contour(y,x,l',[0 0],'y')\r
  19     h=pcolor(xh,yh,(1/hsc)*f');\r
  20     %h=surf(xh,yh,f');\r
  21     c=hot;\r
  22     cc=c;\r
  23     c(2:31,:)=cc(22:51,:);\r
  24     for i=32:64\r
  25         t=(i-32)/(64-32);\r
  26         c(i,:)=t*cc(64,:)+(1-t)*cc(52,:);\r
  27     end\r
  28     c(1,:)=0.45;\r
  29     colormap(c);\r
  30     if(0),\r
  31         hc=colorbar\r
  32         % title(hc,sprintf('%e J/m^2/s',hsc))\r
  33         title(hc,'10^5 J/m^2/s')\r
  34     end\r
  35     mx=4e5; % max color displayed\r
  36     caxis([0,mx/hsc]);\r
  37     hold on\r
  38     set(h,'edgecolor','none')\r
  39     %alpha(h,0.5)\r
  40     %colorbar\r
  41     xlabel('dist (m)')\r
  42     ylabel('dist (m)')\r
  43     zlabel('heat flux (J/m^2/s)')\r
  44     axis equal\r
  45     axis([0,m*dx,0,n*dy,0,mx])\r
  46     %axis tight\r
  47 end\r