Mercurial > repos > agpetit > calculate_diameter
comparison cut_trajectory.py @ 7:f1dd5d99ea2d draft
"planemo upload for repository https://github.com/mesocentre-clermont-auvergne/aubi_piaf commit b6488400d4478d46697019485e912c38ea2202a5-dirty"
author | agpetit |
---|---|
date | Mon, 30 May 2022 15:35:48 +0000 |
parents | e504457035e5 |
children | d45cecb7cfd0 |
comparison
equal
deleted
inserted
replaced
6:a7128cd8f0e1 | 7:f1dd5d99ea2d |
---|---|
257 out_traj = (str(out_dir) + prefix_name_file + "_traj_" + str(k) + ".xtc") | 257 out_traj = (str(out_dir) + prefix_name_file + "_traj_" + str(k) + ".xtc") |
258 bash_command = ("echo " + n_group + " | gmx trjconv -f " + xtc_file | 258 bash_command = ("echo " + n_group + " | gmx trjconv -f " + xtc_file |
259 + " -s " + gro_file + " -b " + str(v.split(",")[0]) | 259 + " -s " + gro_file + " -b " + str(v.split(",")[0]) |
260 + " -e " + str(v.split(",")[1]) + " -o " + out_traj) | 260 + " -e " + str(v.split(",")[1]) + " -o " + out_traj) |
261 if arguments.verbose: | 261 if arguments.verbose: |
262 log_directory = "log/tmp/" | 262 log_directory = "log/tmp/log/" |
263 logging_file = (log_directory | 263 logging_file = (log_directory |
264 + logging_file.rsplit("/", 1)[1].split(".")[0] | 264 + logging_file.rsplit("/", 1)[1].split(".")[0] |
265 + "_" + str(k) + ".log") | 265 + "_" + str(k) + ".log") |
266 if not os.path.exists(log_directory): | 266 if not os.path.exists(log_directory): |
267 os.makedirs(log_directory) | 267 os.makedirs(log_directory) |