Extract last frame of LAMMPS xyz file

extract_last.py

#!/bin/env python3
# Extract last frame from LAMMPS xyz dump file
# inpath and write it to outpath

inpath = "AgNWSiO2.xyz"
outpath = "data.AgNWSiO2"
df = open(inpath, "r")
buf = []
i = 0
for l in df:
    # ITEM: TIME
    # ITEM: TIMESTEP
    if "ITEM: TIME" in l:
        buf.clear()
        print("Step:", i)
        i += 1
        
    buf.append(l)
    
df.close()
df = open(outpath, "w")

for l in buf:
    df.write(l)
    
df.close()