def setup(): if tm.mod_seed.nff == 1: #Forward tm.mod_time.intstart = tm.mod_time.intmin tm.mod_time.intend = tm.mod_time.intmax else: #Backward tm.mod_time.intstart = tm.mod_time.intmin tm.mod_time.intend = tm.mod_time.intmax tm.mod_grid.kmt[:] = tm.mod_param.km from os.path import abspath, join, curdir filename = abspath(join(curdir, 'results-new', 'data')) if tm.mod_seed.nqua == 1: # number of trajectories (per time resolution) # num=NTRACMAX tm.mod_seed.num = tm.mod_param.partquant elif tm.mod_seed.nqua == 2: tm.mod_param.voltr = tm.mod_param.partquant elif tm.mod_seed.nqua == 3: tm.mod_param.voltr = tm.mod_param.partquant tm.fortran_file(56, filename + '_run.asc') # trajectory path tm.fortran_file(57, filename + '_out.asc') # exit position tm.fortran_file(58, filename + '_in.asc') # entrance position tm.fortran_file(59, filename + '_err.asc') # Error position
def __init__(self, start, end, delta): Project.__init__(self, start, end, delta) from os.path import abspath, join, curdir filename = abspath(join(curdir, 'results-new', 'data')) tm.fortran_file(56, filename + '_run.asc') # trajectory path tm.fortran_file(57, filename + '_out.asc') # exit position tm.fortran_file(58, filename + '_in.asc') # entrance position tm.fortran_file(59, filename + '_err.asc') # Error position
def main(): tm.init_params() tm.coordinat() writesetup() print 'tm.time.intmin', tm.mod_time.intmin print 'tm.time.intmax', tm.mod_time.intmax if tm.mod_seed.nff == 1: #Forward print "Forward" tm.mod_time.intstart = tm.mod_time.intmin tm.mod_time.intend = tm.mod_time.intmax else: #Backward print "Backward" tm.mod_time.intstart = tm.mod_time.intmin tm.mod_time.intend = tm.mod_time.intmax print "tm.time.intstart", tm.mod_time.intstart print "tm.time.intend", tm.mod_time.intend tm.init_seed() filename = abspath(join(curdir, 'results-old', 'data')) print "filename", filename print "tm.mod_seed.nqua", tm.mod_seed.nqua if tm.mod_seed.nqua == 1: # number of trajectories (per time resolution) # num=NTRACMAX print "tm.params.part_quant", tm.mod_param.partquant, tm.mod_seed.num tm.mod_seed.num = tm.mod_param.partquant elif tm.mod_seed.nqua == 2: tm.mod_param.voltr = tm.mod_param.partquant elif tm.mod_seed.nqua == 3: tm.mod_param.voltr = tm.mod_param.partquant tm.fortran_file(56, filename + '_run.asc') # trajectory path tm.fortran_file(57, filename + '_out.asc') # exit position tm.fortran_file(58, filename + '_in.asc') # entrance position tm.fortran_file(59, filename + '_err.asc') # Error position tm.loop.readfields = readfields print_state() # tm.loop() print "Done!"