function CARPETpcolorY(d,y)
%
% function CARPETcontourfY(d,y)
%
% plot a surface representation of the plan with y Y-coordinate
% for the quantity d, using all the ref. level and components
%
% EXAMPLE:
%  str = CARPETListDatasetHDF5('data/hdf5','phi',1);
%  d   = CARPETLoadDatasetHDF5(str,2);
%  CARPETpcolorY(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
 
 
  for a=d.rlmin:d.rl
      for b=1:d.c(a)
        %a=1; b=1;
        n=d.num2name(a,b);
        Yr=d.dt(n).Yr;  
        dy=abs(Yr(2)-Yr(1))/2;
        iY=find(Yr>y-dy-dy*0.0001&Yr<y+dy-dy*0.0001);
        if length(iY)==1 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% HERE GRAPHICAL COMMANDS (LOOK AT THE MATLAB HELP FOR INFO)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
        pcolor(d.dt(n).Xr,d.dt(n).Zr,double(squeeze(d.dt(n).val(iY,:,:)))');             
        shading faceted; 
        hold on
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
        end
      end
  end
  
  caxis([d.vmin d.vmax]);
  axis equal; %axis square;
  axis([d.xmin d.xmax d.zmin d.zmax]);
  
  hold off

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

end