trajin ../4/prod_all_fixed.crd
# Generate H-bond info 
hbond hbonds out NEW_HBOND/HBONDS_t4.dat angle 90 dist 3.9 avgout NEW_HBOND/HBONDS-Avg_t4.dat
go
