def main(nfGOtermgraph, taxid, randsimdir, outdir, windowsizes, outcols, dbname, dbuser, dbhost, dbpwd, aspects=aspects, metrics=metrics, blockscores=blockscores, silent=True): dbcon, dbcur = rec_to_db.dbconnect(dbhost=dbhost, dbuser=dbuser, dbname=dbname, dbpwd=dbpwd)
randomseed=int(dopt['-r'])
print "Protein order will be randomized on replicons. Random seed = %d"%randomseed
splitdiscordant = bool(int(dopt.get('-s', True)))
starttime = time.time()
print "Start job %d at %s"%(jobid, time.strftime("%a, %d %b %Y %H:%M:%S", time.localtime()))
if queue:
print "queue = %dh"%queue
tl = float(queue * 3600)
safetl = tl * 0.90
print "time limit = %fs"%tl
print "safe time limit = %fs"%safetl
# sets connection objects to AGROGENOMDB
db, cj = rec_to_db.dbconnect(dbname=dbname, dbuser=dbuser, dbhost=dbhost, dbpwd=dbpwd, dbclue=dbclue)
if buildancs:
dirout = "%s/getblockevent_out.a.%d"%(outdir,jobid)
else:
dirout = "%s/getblockevent_out.l.%d"%(outdir,jobid)
dirconcattrees = '%s/block_trees'%dirout
dirblockout = '%s/blocks'%dirout
dirsqlout = "%s/blocks_db_dump"%dirout
nfstat = '%s/transfered_blocks.stats'%dirout
nfprestat = '%s/before_split_blocks.stats'%dirout
nfancblocks = "%s/ancestral_blocks"%dirout
nfancstats = '%s/ancestral_blocks.stats'%dirout
nfmatrix = "%s/blocks_events.mat"%dirout
