from modeller import *
from modeller.automodel import * # Load the automodel class
#from modeller import soap_protein_od

log.verbose()
env = environ()
a = automodel(env, alnfile='rpn11_human_yeast.aln', # alignment filename
              knowns='rpn11_yeast_23-306', # codes of the templates
              sequence='rpn11_human_27-310', # code of the target
              assess_methods=(	assess.DOPE,
                            	assess.DOPEHR,
                            	#soap_protein_od.Scorer(),
                            	assess.GA341))
a.starting_model = 1 	# index of the first model
a.ending_model = 10		# index of the last model
a.make()				# do homology modeling