function CARPETsurfZ(d,z)
%
% function CARPETsurfZ(d,z)
%
% plot a surface representation of the plan with z Z-coordinate
% for the quantity d, using all the ref. level and components
%
% EXAMPLE:
%  str = CARPETListDatasetHDF5('data/hdf5','phi',1);
%  d   = CARPETLoadDatasetHDF5(str,2);
%  CARPETsurfZ(d,7)
%
% Author: Gian Mario Manca, manca@fis.unipr.it, 
%         http://www.fis.unipr.it/~manca
%
% These scripts are distributed under the GNU General Public License (GPL)
% 12 March 2006

 
 
 %figure
 
 
  for a=d.rlmin:d.rl
      for b=1:d.c(a)
        %a=1; b=1;
        n=d.num2name(a,b);
        Zr=d.dt(n).Zr;  
        dz=abs(Zr(2)-Zr(1))/2;
        iZ=find(Zr>z-dz-dz*0.0001&Zr<z+dz-dz*0.0001);
        if length(iZ)==1 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% HERE GRAPHICAL COMMANDS (LOOK AT THE MATLAB HELP FOR INFO)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
            surf(d.dt(n).Xr,d.dt(n).Yr,double(d.dt(n).val(:,:,iZ))); 
            %shading faceted; 
            hold on
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
        end
      end
  end
  
  caxis([d.vmin d.vmax]);
  %axis equal; %axis square;
  axis([d.xmin d.xmax d.ymin d.ymax]);
  
  hold off

  xlabel(['Time=' num2str(d.timeMSEC) ' ms and Z= ' num2str(z)]);
  title(d.names(d.num2name(d.rlmin,1)));
  disp(['DISPLAYING the plain with Z = ' num2str(z)]);
        
end