from ase.io import read
from mbpol_ase import mbpol

ice = read('ice.vasp')
calc = mbpol(atoms=ice,X=2,Y=2,Z=2)
ice.set_calculator(calc)

print (ice.get_potential_energy())