clear all
load('sim1')
ne = 2; % choose the number of experiment
[ndata, text, ~] = xlsread('BinMatEvol.xlsx','BinRange');
min_U = ndata(:,3);
step_U = ndata(:,4);
max_U = ndata(:,5);
param_V = ndata(:,1);
ParName = text(2:end,2);
BinAge = fliplr(0:5:595);
        
k=1;
for ii=1:length(ParName)
    M=EXP(ne).S{ii};
    BinV = min_U(ii):step_U(ii):max_U(ii);
    dpar =( BinV(2) - BinV(1) ); % Parameter increment               
    binCen = BinV + dpar/2; % Bin centers
    figure
    contourf(BinAge,binCen,M ,10,'LineColor','none'); 
    ylabel(ParName{ii})

end
