Hello!
This is a Fortran90 code for the radial distribution function of a molecular trajectory averaged over time. This code uses molecular center of masses and periodic boundary conditions (3x3x3) for the calculation and uses unformatted files (gromacs: xtc). The part that reads the xtc format was written by James W. Barnett jbarnet4@tulane.edu.The code has comments all over, but feel free to reach out to me at leejunta@grinnell.edu if you have any more questions. Special thanks to Roger Rousseau (PNNL), David Cantu (PNNL), and Mal Soon Lee (PNNL) for assisting me in developing this code.
Here's how to run the code:
-
Install the xtc-library (http://www.gromacs.org/Developer_Zone/Programming_Guide/XTC_Library)
-
Copy the xtc-interface.f90 file into the appropriate folder
-
Edit the shell script to input the information you want inputted into the fortran code
-
gfortran -c xtc-interface.f90 -lxdrfile
-
gfortran xtc-interface.o rdf.f90 -lxdrfile -o rdf.x
-
chmod +x rdf.sh
-
./rdf.sh
Source code available at https://github.com/leejunta/Radial-Distribution-Function/