trajin ../13/prod_all_fixed.crd
# Generate H-bond info 
hbond hbonds out NEW_HBOND/HBONDS_t13.dat angle 90 dist 3.9 avgout NEW_HBOND/HBONDS-Avg_t13.dat
go
