Band Edge Alignment

[Chhatra]

Dear CRYSTAL developers and users,

I would like to perform a band edge alignment calculation by referencing the valence/conduction band edges to the vacuum level in the context of photocatalysis. For this purpose, I need the planar-averaged electrostatic potential along the slab normal direction.

I have already generated the electrostatic potential in Gaussian cube format using the POTCUBE keyword (file: POT_CUBE.DAT).

My questions are:

Does CRYSTAL provide any direct tool or utility to compute the planar average of the electrostatic potential from POT_CUBE.DAT (or other related output files)?

If not, is it recommended to post-process the cube file externally (e.g., with scripts, VESTA, or custom Python routines)?

Any guidance, example workflows, or scripts would be greatly appreciated.

Best regards,
Chhatra Subba