key old select all end
inquire char vari Str_arg text "ATOM COORDINATE FILE NAME::"
read file Str_arg atom pdb
key new select .not old end
calc bond resi sele new end
calc bond neigh sele new .a atom name C N end
calc bond disu sele new end
ima cent calc select new end
image select new end set bond 
return