function CARPETcontourfZ(d,z)
%
% function CARPETcontourfZ(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);
%  CARPETcontourfZ(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 2007

 
 %figure

% Fix the interval of values used to generate the isodensity lines and the
% colormap: the default is the entire interval of values on the entire
% grid. Reduce it if you want to focus in a smaller interval of values.
minval=d.vmin;
maxval=d.vmax;
 
 
  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)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	isovalues = minval+(maxval-minval)*[0:0.1:1]; 
        % a simple linear scale with 10 point between minval and maxval        
        % isovalues = minval+(maxval-minval)*[0 2.^[-20:0]];
        % an example of logarithmic scale
        contourf(d.dt(n).Xr,d.dt(n).Yr,double(d.dt(n).val(:,:,iZ)),isovalues);             
        %shading flat; 
        hold on    
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
        end
      end
  end
  
  caxis([minval maxval]);
  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