action="store",
type="string",
default="CG_bilayer",
dest="filename")
parser.add_option("-t", action="store_true", dest="timtails", default=False)
(options, args) = parser.parse_args()
filename = options.filename
# load a structure, generated from mbuild
#hoomd.init.read_gsd(filename, restart= 'trajectory.gsd')
hoomd.init.read_gsd(filename, restart='restart.gsd')
atomtypes = martini_hoomd.get_atom_types()
hoomd.util.quiet_status()
martini_hoomd.set_bonds()
#martini_hoomd.set_harmonic_angles(timtails = True) # For using harmonic-fit potentials
martini_hoomd.set_cosine_angles(timtails=options.timtails)
nl = hoomd.md.nlist.cell()
# Set exclusions only if using martinitim
martini_hoomd.set_exclusions(nl=nl)
if options.timtails:
table_dir = "lambda0.9-martinitim"
else:
table_dir = "lambda0.9-martinifull"
lamb0 = martini_hoomd.set_pairs(table_dir=table_dir, nl=nl)
# Define some groups
all = hoomd.group.all()
hoomd.util.unquiet_status()
# Output parameters
action="store",
type="string",
default="CG_bilayer",
dest="filename")
(options, args) = parser.parse_args()
filename = options.filename
# load a structure, generated from mbuild
hoomd.init.read_gsd(filename, restart='restart.gsd')
#hoomd.init.read_gsd('trajectory.gsd')
atomtypes = martini_hoomd.get_atom_types()
ENERGY_TABLE, DISTANCE_TABLE, FORCE_TABLE = martini_hoomd.get_table_units()
hoomd.util.quiet_status()
martini_hoomd.set_bonds()
#martini_hoomd.set_harmonic_angles(timtails = True) # For using harmonic-fit potentials
martini_hoomd.set_cosine_angles(timtails=False)
nl = hoomd.md.nlist.cell()
# Set exclusions only if using martinitim
#martini_hoomd.set_exclusions(nl=nl)
#lamb0 = martini_hoomd.set_pairs(table_dir = "lambda0-martinitim",nl=nl)
lamb0 = martini_hoomd.set_pairs(table_dir="lambda0-martinifull", nl=nl)
# Set reaction field parameters
#martini_hoomd.set_reactionfield_electrostatics(nl=nl)
# Define some groups
all = hoomd.group.all()
hoomd.util.unquiet_status()
# Output parameters