#============================================================#
#          look at hydrogen bonds                            #
#============================================================#

system echo = off
system verbose = false

# set units
#----------
#units
units scale  name = mass   uname = reducedMass   value = 10000.0

# open pdb file
#--------------
database open name = db1  format = pdb  file = ../data/neck_xform.pdb

# read molecule
#--------------
molecule read  name = m1   database = db1

# create helix domain
#--------------------
domain create helix = C[783-791]
domain helix bonds atoms = true  atom_color = element  bond_color = element  
domain helix  backbone  color = [ 1 1 1 ]  planes = true  width = 3

domain helix hbonds color = [ 1 1 1 ]  width = 3 show