[modeller_usage] No PDB output

Hi,

I ran a test modeller script to generate 2 models - script attached to this
e-mail. From running this script, I only got the Dxxxxx files and not the PDB
output files. The log file didn't contain any error messages for 
generating PDB
files - might anyone be able to suggest a solution please?

Thanks,

Yalini

from modeller.automodel import *

log.verbose()    # for debugging

env = environ(rand_seed=-12312) # random seed

env.io.atom_files_directory = '.'  # Working directory

def defsym(mdl, aln, seg1, seg2):
    """Constrain segments |seg1| and |seg2| to be identical"""
    for (set,seg) in [(2,seg1), (3,seg2)]:
        mdl.pick_atoms(aln, pick_atoms_set=set, selection_segment=seg,
                       atom_types='ALL', selection_status='INITIALIZE',
                       selection_search='SEGMENT')
    mdl.symmetry.define(symmetry_weight=1.0, add_symmetry=(True, False))


def makealpha(mdl, aln, residue_ids):
    """Enforce alpha-helical structure on |residue_ids|"""
    mdl.restraints.make(aln, restraint_type='ALPHA', residue_ids=residue_ids,
                        spline_on_site=False)


class mymodel(automodel):
    def special_restraints(self, aln):
        segs = [('1', '99'), ('100', '198'), ('199', '297'), ('298', '396')]
        for n in range(len(segs)):
            defsym(self, aln, segs[n-1], segs[n])

        makealpha(self, aln, ('5', '31'))  	#M1 region (consensus)
        makealpha(self, aln, ('40', '53')) 	#P region (jalview p/out)
        makealpha(self, aln, ('64', '95')) 	#M2 region (consensus)
        makealpha(self, aln, ('104', '130'))
        makealpha(self, aln, ('139', '152'))
        makealpha(self, aln, ('163', '194'))
	makealpha(self, aln, ('203', '229'))
        makealpha(self, aln, ('238', '251'))
        makealpha(self, aln, ('262', '293'))
	makealpha(self, aln, ('302', '328'))
        makealpha(self, aln, ('337', '350'))
        makealpha(self, aln, ('361', '392'))


a = mymodel(env,
            alnfile  = 'test.ali',  # The alignment filename
            knowns   = ('target1', 'target2'),  # The structure filenames, the pdb
                                     # is called filename.atm
            sequence = 'template')      # The sequence filename
a.starting_model= 1                  # Starting and ending model indices
a.ending_model  = 2
a.deviation = 4.0                    # Deviation in models

# Call the routine model which in turn calls all the other routines required.
a.make()