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!"
