From b7c477d34e4a70e8e3dcbc893548ddcadcf984dc Mon Sep 17 00:00:00 2001
From: Gabriel Studer <gabriel.studer@unibas.ch>
Date: Wed, 14 Jun 2017 14:29:43 +0200
Subject: [PATCH] script to generate the charmm and amber forcefields in ost

---
 modules/mol/mm/pymod/generate_forcefields.py | 33 ++++++++++++++++++++
 1 file changed, 33 insertions(+)
 create mode 100644 modules/mol/mm/pymod/generate_forcefields.py

diff --git a/modules/mol/mm/pymod/generate_forcefields.py b/modules/mol/mm/pymod/generate_forcefields.py
new file mode 100644
index 000000000..4dfe517b1
--- /dev/null
+++ b/modules/mol/mm/pymod/generate_forcefields.py
@@ -0,0 +1,33 @@
+import os
+from ost.mol import mm
+
+path_to_gromacs = "xyz"
+charmm_ff = "charmm27.ff"
+amber_ff = "amber03.ff"
+
+charmm_path = os.path.join(path_to_gromacs, "share", "top", charmm_ff)
+amber_path = os.path.join(path_to_gromacs, "share", "top", amber_ff)
+
+# Let's first do CHARMM
+reader = mm.FFReader(charmm_path)
+reader.ReadGromacsForcefield()
+reader.ReadResidueDatabase("aminoacids")
+reader.ReadResidueDatabase("rna")
+reader.ReadResidueDatabase("dna")
+reader.ReadResidueDatabase("lipids")
+reader.ReadITP("tip3p")
+reader.ReadITP("ions")
+ff = reader.GetForcefield()
+ff.Save("CHARMM27.dat")
+
+# Let's do amber
+reader = mm.FFReader(amber_path)
+reader.ReadGromacsForcefield()
+reader.ReadResidueDatabase("aminoacids")
+reader.ReadResidueDatabase("rna")
+reader.ReadResidueDatabase("dna")
+reader.ReadITP("tip3p")
+reader.ReadITP("ions")
+ff = reader.GetForcefield()
+ff.Save("AMBER03.dat")
+
-- 
GitLab