trajin ../33/prod_all_fixed.crd
# Generate H-bond info 
hbond hbonds out HBOND/HBONDS_t33.dat angle 100 dist 3.0 avgout HBOND/HBONDS-Avg_t33.dat
go
