1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
|
#! /usr/bin/env -S bash -e
#equilibration time in Gromacs units (ps)
equi=200
echo equi = $equi
if [[ ! -f ../atomistic/topol.tpr || ! -f ../atomistic/traj.trr ]]; then
echo "Run atomistic simulation in ../atomistic first"
exit 1
fi
echo "Running force matching"
csg_fmatch --top ../atomistic/topol.tpr --trj ../atomistic/traj.trr --begin $equi --options fmatch.xml --cg water.xml
#integrate force table to get potential
csg_call table integrate CG-CG.force CG-CG.pot
csg_call table linearop CG-CG.pot CG-CG.pot -1 0
#copy CG-CG.pot to new file to prevent deletion by clean command
cp CG-CG.pot input.pot
#convert to gromacs potential
csg_call --options fmatch.xml --ia-name CG-CG --ia-type non-bonded convert_potential gromacs --clean input.pot table_CG_CG.xvg
|
