function CARPETcontourfX(d,x)
%
% function CARPETcontourfX(d,x)
%
% plot a surface representation of the plan with x X-coordinate
% for the quantity d, using all the ref. level and components
%
% EXAMPLE:
%  str = CARPETListDatasetHDF5('data/hdf5','phi',1);
%  d   = CARPETLoadDatasetHDF5(str,2);
%  CARPETcontourfX(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);
        Xr=d.dt(n).Xr;  
        dx=abs(Xr(2)-Xr(1))/2;
        iX=find(Xr>x-dx-dx*0.0001&Xr<x+dx-dx*0.0001);
        if length(iX)==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).Yr,d.dt(n).Zr,double(squeeze(d.dt(n).val(:,iX,:)))',isovalues);             
        %shading flat; 
        hold on    
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
        end
      end
  end
  
  caxis([minval maxval]);
  axis equal; %axis square;
  axis([d.ymin d.ymax d.zmin d.zmax]);
  
  hold off

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

end