two different color axis for two surfaces on the same plot?

Hi everyone,

I´ve been running some simple wavepacket simulations, which are running nice, but I am having some issues for plotting them. I am plotting the potential energy with meshc, so I get some nice contour lines at the bottom of the plot, and mesh for the wavepacket itself on the same plot using hold on. What I would like to do is to plot the potential energy with a preseted color axis and the wavepacket with a different color axis that adapts to the size of the wavepacket as its develop. But so far I have failed miserably at just trying to set two distinctive color schemes for each one.

The code I am using for this is

levels =  5.8e-20:0.2e-20:8.6e-20;
figure;
hold on
sc=meshc(r_ab,r_bc,V_e);
axis('square')
set(sc(2),'LevelList',levels)
zlim([1e-21 1e-19])
sc(2).ContourZLevel = 1e-21;
caxis(1.0e-20*[5 10])
global Plothandle2, Plothandle2 = mesh(r_ab,r_bc,abs(psi_ex));
view([120 30])

where r_ab and r_bc are coordinates generated with meshgrid, V_e is the potential surface and psi_ex is the wavepacket. The following is a snapshot, clearly the features of the wavepacket cant be seen when both surfaces share the same color axis.

Thanks!