# creates:  diffusion-path.png
import runpy

from ase.io import read, write

runpy.run_path('diffusion4.py')
images = [read('mep%d.traj' % i) for i in range(5)]
a = images[0] + images[1] + images[2] + images[3] + images[4]
del a.constraints
a *= (2, 1, 1)
a.set_cell(images[0].get_cell())
renderer = write('diffusion-path.pov', a,
                 rotation='-90x',
                 povray_settings=dict(transparent=False))
renderer.render()