function null = mccgramplot3dmesh(correlogram)
% function null = mccgramplot3dmesh(correlogram)
%
%--------------------------------------------------------------------
% Plots a 3-dimensional binaural correlogram uising 'mesh'
%--------------------------------------------------------------------
%
% Input:
%    correlogram = 'correlogram' structure as defined in mccgramcreate.m
%
% Output parameters:
%    none
%
% Examples:
% to plot a previously-made correlogram cc1, type:
% >> mccgramplot3dmesh(cc1);
%
% See mccgramplot4panel.m for another example.
%
%
% MAA Winter 2001 2i01
%--------------------------------------------------------

% ****************************************************************
% This MATLAB software was developed by Michael A Akeroyd for 
% supporting research at the University of Connecticut
% and the University of Sussex.  It is made available
% in the hope that it may prove useful. Any support is extremely
% limited; in particular, there is no guarantee that any help 
% can be offered in its use nor that any bugs will be fixed.
% 
% Any for-profit use or redistribution is prohibited. No warranty
% is expressed or implied. If you modify this software,
% please include a statement indicating what has been changed.
% 
%    Contact address:
%      Dr Michael A Akeroyd,
%      Laboratory of Experimental Psychology, 
%      University of Sussex, 
%      Falmer, 
%      Brighton, BN1 9QG, 
%      United Kingdom.
%    email:   maa@biols.susx.ac.uk 
%    webpage: http://www.biols.susx.ac.uk/Home/Michael_Akeroyd/
%    
% ******************************************************************


hold off;


% plot!
mesh(correlogram.delayaxis, correlogram.freqaxiserb, correlogram.data);
shading interp;
view(-10,70);
zlabel('Crossproduct');

switch correlogram.type
case 'binauralcorrelogram' ; 
   title('Binaural correlogram');
case 'autocorrelogram' ; 
   title('Autocorrelogram');
end;

colormap(jet);

% reset freq axes to filterbank limits
oldaxes = axis;
[filter_q, filter_bwmin] = mstandarderbparameters;
xmin = correlogram.mindelay;  
xmax = correlogram.maxdelay;  
ymin = mhztoerb(correlogram.mincf,filter_q, filter_bwmin, 0);
ymax = mhztoerb(correlogram.maxcf,filter_q, filter_bwmin, 0);
axis([xmin xmax ymin ymax oldaxes(5) oldaxes(6)]);
maxislabels_freq('ytick', 'yticklabel');
ylabel('Frequency (Hz)');
maxislabels_delay(correlogram.type);

% the end!
%---------------------------------------------------------------------------