Intermolecular bonded interaction support added
[alexxy/gromacs.git] / src / gromacs / gmxpreprocess / grompp.c
index 30881cbf299ca5a8ce58338b5a48a69978c0975a..2c9ed586d79f77b70621aadcd9faf55990d2932f 100644 (file)
@@ -499,7 +499,8 @@ new_status(const char *topfile, const char *topppfile, const char *confin,
            t_gromppopts *opts, t_inputrec *ir, gmx_bool bZero,
            gmx_bool bGenVel, gmx_bool bVerbose, t_state *state,
            gpp_atomtype_t atype, gmx_mtop_t *sys,
-           int *nmi, t_molinfo **mi, t_params plist[],
+           int *nmi, t_molinfo **mi, t_molinfo **intermolecular_interactions,
+           t_params plist[],
            int *comb, double *reppow, real *fudgeQQ,
            gmx_bool bMorse,
            warninp_t wi)
@@ -524,7 +525,8 @@ new_status(const char *topfile, const char *topppfile, const char *confin,
     /* TOPOLOGY processing */
     sys->name = do_top(bVerbose, topfile, topppfile, opts, bZero, &(sys->symtab),
                        plist, comb, reppow, fudgeQQ,
-                       atype, &nrmols, &molinfo, ir,
+                       atype, &nrmols, &molinfo, intermolecular_interactions,
+                       ir,
                        &nmolblock, &molblock, bGB,
                        wi);
 
@@ -1510,7 +1512,7 @@ int gmx_grompp(int argc, char *argv[])
     t_gromppopts      *opts;
     gmx_mtop_t        *sys;
     int                nmi;
-    t_molinfo         *mi;
+    t_molinfo         *mi, *intermolecular_interactions;
     gpp_atomtype_t     atype;
     t_inputrec        *ir;
     int                natoms, nvsite, comb, mt;
@@ -1637,7 +1639,8 @@ int gmx_grompp(int argc, char *argv[])
     snew(state, 1);
     new_status(fn, opt2fn_null("-pp", NFILE, fnm), opt2fn("-c", NFILE, fnm),
                opts, ir, bZero, bGenVel, bVerbose, state,
-               atype, sys, &nmi, &mi, plist, &comb, &reppow, &fudgeQQ,
+               atype, sys, &nmi, &mi, &intermolecular_interactions,
+               plist, &comb, &reppow, &fudgeQQ,
                opts->bMorse,
                wi);
 
@@ -1806,7 +1809,8 @@ int gmx_grompp(int argc, char *argv[])
     }
 
     ntype = get_atomtype_ntypes(atype);
-    convert_params(ntype, plist, mi, comb, reppow, fudgeQQ, sys);
+    convert_params(ntype, plist, mi, intermolecular_interactions,
+                   comb, reppow, fudgeQQ, sys);
 
     if (debug)
     {