Python script to calculate the electrostatic potential in a CG system