Added check for domain decomposition and shells.

[alexxy/gromacs.git] / src / kernel / md.c

diff --git a/src/kernel/md.c b/src/kernel/md.c
index c9415eee962b465f76003d59b269d30f7b4ef0c7..4c4a88c174dbd8b47337a13fed1d23c3dadf4c76 100644 (file)
--- a/src/kernel/md.c
+++ b/src/kernel/md.c
@@ -351,6 +351,10 @@ double do_md(FILE *fplog, t_commrec *cr, int nfile, const t_filenm fnm[],
     {
         gmx_fatal(FARGS, "You have nstcalcenergy set to a value (%d) that is different from 1.\nThis is not supported in combinations with shell particles.\nPlease make a new tpr file.", ir->nstcalcenergy);
     }
+    if (shellfc && DOMAINDECOMP(cr))
+    {
+        gmx_fatal(FARGS, "In order to run parallel simulations with shells you need to use the -pd flag to mdrun.");
+    }
     if (shellfc && ir->eI == eiNM)
     {
         /* Currently shells don't work with Normal Modes */