#!/usr/bin/env python

__author__ = "Magdalena A. Jonikas"
__version__ = "1.0"
__date__ = "2010-05-22"
__doc__ = """Example script to test Nast random unfold algorithm."""


seqInName  		    = "6TNA.seq"
pdbInName           = "pdb-start.pdb"
pdbOutName          = "6TNA-pdb"
helixFilename		= "6TNA_helix.txt"
numSteps		    = 20000 
stepsPerReport		= 1000
defaultTemperature	= 300.0
randomSeed		    = None
verbose			    = True
movieOut            = "6TNA-pdb"     # Can be set to None for no movie trace

import simtk.nast.molecule as molecule
molecule.runRandomUnfold(seqInName=seqInName,
             pdbOutName=pdbOutName,
             helixFilename=helixFilename,
             pdbInName=pdbInName,
             numSteps=numSteps,
             stepsPerReport=stepsPerReport,
             defaultTemperature=defaultTemperature,
             randomSeed=randomSeed,
             verbose=verbose,
             movieOut=movieOut)