| 916 |
|
painCave.isFatal = 1; |
| 917 |
|
simError(); |
| 918 |
|
} |
| 919 |
< |
|
| 920 |
< |
// get the ensemble |
| 919 |
> |
if (globals->haveForceFieldVariant()) { |
| 920 |
> |
strcpy(forcefield_variant, globals->getForceFieldVariant()); |
| 921 |
> |
has_forcefield_variant = 1; |
| 922 |
> |
} |
| 923 |
> |
|
| 924 |
> |
// get the ensemble |
| 925 |
|
|
| 926 |
|
strcpy(ensemble, globals->getEnsemble()); |
| 927 |
|
|
| 1518 |
|
void SimSetup::createFF(void){ |
| 1519 |
|
switch (ffCase){ |
| 1520 |
|
case FF_DUFF: |
| 1521 |
< |
the_ff = new DUFF(); |
| 1521 |
> |
the_ff = new DUFF(); |
| 1522 |
|
break; |
| 1523 |
|
|
| 1524 |
|
case FF_LJ: |
| 1526 |
|
break; |
| 1527 |
|
|
| 1528 |
|
case FF_EAM: |
| 1529 |
< |
the_ff = new EAM_FF(); |
| 1529 |
> |
if (has_forcefield_variant) |
| 1530 |
> |
the_ff = new EAM_FF(forcefield_variant); |
| 1531 |
> |
else |
| 1532 |
> |
the_ff = new EAM_FF(); |
| 1533 |
|
break; |
| 1534 |
|
|
| 1535 |
|
case FF_H2O: |
| 1543 |
|
simError(); |
| 1544 |
|
} |
| 1545 |
|
|
| 1546 |
+ |
|
| 1547 |
|
#ifdef IS_MPI |
| 1548 |
|
strcpy(checkPointMsg, "ForceField creation successful"); |
| 1549 |
|
MPIcheckPoint(); |