| 484 | 
  | 
      } | 
| 485 | 
  | 
       | 
| 486 | 
  | 
 | 
| 487 | 
< | 
      //creat cutoff group for molecule | 
| 487 | 
> | 
      //create cutoff group for molecule | 
| 488 | 
  | 
      molInfo.myCutoffGroups.clear(); | 
| 489 | 
  | 
      for (j = 0; j < molInfo.nCutoffGroups; j++){ | 
| 490 | 
  | 
 | 
| 563 | 
  | 
  sprintf(checkPointMsg, "all molecules initialized succesfully"); | 
| 564 | 
  | 
  MPIcheckPoint(); | 
| 565 | 
  | 
#endif // is_mpi | 
| 566 | 
– | 
 | 
| 567 | 
– | 
  // clean up the forcefield | 
| 568 | 
– | 
 | 
| 569 | 
– | 
  if (!globals->haveRcut()){ | 
| 566 | 
  | 
 | 
| 571 | 
– | 
    the_ff->calcRcut(); | 
| 572 | 
– | 
 | 
| 573 | 
– | 
  } else { | 
| 574 | 
– | 
     | 
| 575 | 
– | 
    the_ff->setRcut( globals->getRcut() ); | 
| 576 | 
– | 
  } | 
| 577 | 
– | 
 | 
| 578 | 
– | 
  the_ff->cleanMe(); | 
| 567 | 
  | 
} | 
| 568 | 
  | 
 | 
| 569 | 
  | 
void SimSetup::initFromBass(void){ | 
| 1004 | 
  | 
#endif //is_mpi | 
| 1005 | 
  | 
 | 
| 1006 | 
  | 
    double theRcut, theRsw; | 
| 1007 | 
+ | 
 | 
| 1008 | 
+ | 
    if (globals->haveRcut()) { | 
| 1009 | 
+ | 
      theRcut = globals->getRcut(); | 
| 1010 | 
+ | 
 | 
| 1011 | 
+ | 
      if (globals->haveRsw())  | 
| 1012 | 
+ | 
        theRsw = globals->getRsw(); | 
| 1013 | 
+ | 
      else  | 
| 1014 | 
+ | 
        theRsw = theRcut; | 
| 1015 | 
+ | 
       | 
| 1016 | 
+ | 
      info[i].setDefaultRcut(theRcut, theRsw); | 
| 1017 | 
+ | 
 | 
| 1018 | 
+ | 
    } else { | 
| 1019 | 
+ | 
       | 
| 1020 | 
+ | 
      the_ff->calcRcut(); | 
| 1021 | 
+ | 
      theRcut = info[i].getRcut(); | 
| 1022 | 
+ | 
 | 
| 1023 | 
+ | 
      if (globals->haveRsw())  | 
| 1024 | 
+ | 
        theRsw = globals->getRsw(); | 
| 1025 | 
+ | 
      else  | 
| 1026 | 
+ | 
        theRsw = theRcut; | 
| 1027 | 
+ | 
       | 
| 1028 | 
+ | 
      info[i].setDefaultRcut(theRcut, theRsw); | 
| 1029 | 
+ | 
    } | 
| 1030 | 
  | 
 | 
| 1031 | 
  | 
    if (globals->getUseRF()){ | 
| 1032 | 
  | 
      info[i].useReactionField = 1; | 
| 1033 | 
< | 
 | 
| 1033 | 
> | 
       | 
| 1034 | 
  | 
      if (!globals->haveRcut()){ | 
| 1035 | 
  | 
        sprintf(painCave.errMsg, | 
| 1036 | 
  | 
                "SimSetup Warning: No value was set for the cutoffRadius.\n" | 
| 1107 | 
  | 
  strcpy(checkPointMsg, "post processing checks out"); | 
| 1108 | 
  | 
  MPIcheckPoint(); | 
| 1109 | 
  | 
#endif // is_mpi | 
| 1110 | 
+ | 
 | 
| 1111 | 
+ | 
  // clean up the forcefield | 
| 1112 | 
+ | 
  the_ff->cleanMe(); | 
| 1113 | 
  | 
} | 
| 1114 | 
  | 
   | 
| 1115 | 
  | 
void SimSetup::initSystemCoords(void){ |