Fortran Simulation Output¶
Bead positions and identity¶
Periodically while running the simulation will take a save point and save the configuration of the polymer(s). These are saved in the data/ folder in files such as data/r100 where the number 100 denotes what savepoint it is. When using MPI to parallel temper, the data will be saved in the format data/r110v6 where 100 refers to the save point and 6 refers to the replica. This data will have three columns for the x, y, and z positions of each bead. If there are multiple polymers they will be printed one after the next with no delineation.
If WLC_P__SAVEAB is specified as true than data/r100v6 will contain a fourth column specifying the bead type. For two tail methylation profiles a value of 0 in this column refers to either tail bound by HP1, a value of 1 or 2 refers to one or the other tail bound by HP1, while a value of 3 refers to both tails being bound by HP1. If WLC_P__CHANGINGCHEMICALIDENTITY is true than a fifth column will specify the number of methylated tails (0, 1 or 2).
For example
21.702 39.628 14.415 0 0
21.715 39.617 14.892 3 1
21.205 39.449 14.758 3 2
20.718 39.472 14.672 3 2
20.463 39.763 14.803 0 0
20.209 39.667 15.418 1 0
19.888 39.643 14.826 2 0
...
Other files¶
If WLC_P__SAVEU is true then the orientation unit vectors for each bead will be stored as in files of a similar format to the positions. For example data/u100 or data/u100v6. If WLC_P__LOCAL_TWIST is also true then there will be three additional columns in the data/u100.
The file data/energiesv6 contains the various energies and corresponding coefficients. All energies are in kT’s. There is one line in this file per savepoint. The file data/adaptationsv6 contains the adapted window size, success rate, and so on the various move types. These files are useful for diagnosing problems and determining if the simulation is equilibrated.
A file data/repHistory contains various information about replica coupling if it is in use.
Visualization¶
We provide visualization code based on pymol. To run, navigate to ‘visualization/polymer_mole’, enter the savepoint and other settings into run_visualization.py and then type python3 run_visualization.py. This will require both python3 and pymol. The way it works is that the output described above is converted into a pdb file which is then passed to pymol. Various sets of settings that have been used in the past are included in run_visualization.py. To use one of thes set image to the name of the corrisponing set.
If closePymol is set to true in run_visualization than the output will be saved to the data/ directiory. By editing the line for rep in … you can save many snapshots.