下载源码后解压，直接编译即可。

运行时需要输入这个几个参数：

第一个是包含电荷密度或者静电势原始数据的文件名称。第二个是输出文件名称。第三个是针对原始数据的缩放因子，因为cube格式使用原子单位制，需要转换为国际单位制（立方埃或eV）则指定该因子，不需转换时则设置为1。最后一个是积分算法，1为梯形法，2为简单求和，取1即可。

输出文件中有三列数据：第一列是z轴坐标，单位为bohr；第二列是积分后的格点数据；第三列为平均后的格点数据。做面内积分电荷密度采用第一、第二列数据，做平均电势取第一、第三列数据。运行效果如下(以黑磷加电场后，面内积分电荷密度为例)：

程序目前只支持cube格式的数据，且只能在xOy平面内积分。VASP CHGCAR可以用VTST附带的chg2cube.pl转换成cube格式。

potslice.tar.bz2