Take the following scenario. I have two cube files, the first one holds some values of a function (electron localization function),
while in the second cube for each value in the first cube file i can have either one of the numbers from 0 to 6.
After loading the first cube file in Pymol, I create an isosurface of a given value for the contour level, lets say 0.8 a.u.
Next I want to color this isosurface based on the values in the second cube file, e.g. all points that have a value of 3 in the second cube
should be colored red etc.
Any ideas on how to accomplish this? If needed, I could provide the two cube files for a simple test case