--- trunk/OOPSE/libmdtools/SimSetup.cpp	2003/09/21 05:09:16	779
+++ trunk/OOPSE/libmdtools/SimSetup.cpp	2003/09/22 21:23:25	780
@@ -22,8 +22,6 @@
 #define NVT_ENS        1
 #define NPTi_ENS       2
 #define NPTf_ENS       3
-#define NPTim_ENS      4
-#define NPTfm_ENS      5
 
 #define FF_DUFF 0
 #define FF_LJ   1
@@ -603,12 +601,6 @@ void SimSetup::gatherInfo(void){
   }
   else if (!strcasecmp(ensemble, "NPTf")){
     ensembleCase = NPTf_ENS;
-  }
-  else if (!strcasecmp(ensemble, "NPTim")){
-    ensembleCase = NPTim_ENS;
-  }
-  else if (!strcasecmp(ensemble, "NPTfm")){
-    ensembleCase = NPTfm_ENS;
   }
   else{
     sprintf(painCave.errMsg,
@@ -1350,9 +1342,7 @@ void SimSetup::makeIntegrator(void){
 
   NVT<RealIntegrator>* myNVT = NULL;
   NPTi<NPT<RealIntegrator> >* myNPTi = NULL;
-  NPTf<RealIntegrator>* myNPTf = NULL;
-  NPTim<RealIntegrator>* myNPTim = NULL;
-  NPTfm<RealIntegrator>* myNPTfm = NULL;
+  NPTf<NPT<RealIntegrator> >* myNPTf = NULL;
   
   for (k = 0; k < nInfo; k++){
     switch (ensembleCase){
@@ -1430,10 +1420,10 @@ void SimSetup::makeIntegrator(void){
       case NPTf_ENS:
         if (globals->haveZconstraints()){
           setupZConstraint(info[k]);
-          myNPTf = new ZConstraint<NPTf<RealIntegrator> >(&(info[k]), the_ff);
+          myNPTf = new ZConstraint<NPTf<NPT <RealIntegrator> > >(&(info[k]), the_ff);
         }
         else
-          myNPTf = new NPTf<RealIntegrator>(&(info[k]), the_ff);
+          myNPTf = new NPTf<NPT <RealIntegrator> >(&(info[k]), the_ff);
 
         myNPTf->setTargetTemp(globals->getTargetTemp());
 
@@ -1459,47 +1449,6 @@ void SimSetup::makeIntegrator(void){
 
         if (globals->haveTauBarostat())
           myNPTf->setTauBarostat(globals->getTauBarostat());
-        else{
-          sprintf(painCave.errMsg,
-                  "SimSetup error: If you use an NPT\n"
-                  "    ensemble, you must set tauBarostat.\n");
-          painCave.isFatal = 1;
-          simError();
-        }
-        break;
-
-      case NPTim_ENS:
-        if (globals->haveZconstraints()){
-          setupZConstraint(info[k]);
-          myNPTim = new ZConstraint<NPTim<RealIntegrator> >(&(info[k]), the_ff);
-        }
-        else
-          myNPTim = new NPTim<RealIntegrator>(&(info[k]), the_ff);
-
-        myNPTim->setTargetTemp(globals->getTargetTemp());
-
-        if (globals->haveTargetPressure())
-          myNPTim->setTargetPressure(globals->getTargetPressure());
-        else{
-          sprintf(painCave.errMsg,
-                  "SimSetup error: If you use a constant pressure\n"
-                  "    ensemble, you must set targetPressure in the BASS file.\n");
-          painCave.isFatal = 1;
-          simError();
-        }
-
-        if (globals->haveTauThermostat())
-          myNPTim->setTauThermostat(globals->getTauThermostat());
-        else{
-          sprintf(painCave.errMsg,
-                  "SimSetup error: If you use an NPT\n"
-                  "    ensemble, you must set tauThermostat.\n");
-          painCave.isFatal = 1;
-          simError();
-        }
-
-        if (globals->haveTauBarostat())
-          myNPTim->setTauBarostat(globals->getTauBarostat());
         else{
           sprintf(painCave.errMsg,
                   "SimSetup error: If you use an NPT\n"
@@ -1509,47 +1458,6 @@ void SimSetup::makeIntegrator(void){
         }
         break;
 
-      case NPTfm_ENS:
-        if (globals->haveZconstraints()){
-          setupZConstraint(info[k]);
-          myNPTfm = new ZConstraint<NPTfm<RealIntegrator> >(&(info[k]), the_ff);
-        }
-        else
-          myNPTfm = new NPTfm<RealIntegrator>(&(info[k]), the_ff);
-
-        myNPTfm->setTargetTemp(globals->getTargetTemp());
-
-        if (globals->haveTargetPressure())
-          myNPTfm->setTargetPressure(globals->getTargetPressure());
-        else{
-          sprintf(painCave.errMsg,
-                  "SimSetup error: If you use a constant pressure\n"
-                  "    ensemble, you must set targetPressure in the BASS file.\n");
-          painCave.isFatal = 1;
-          simError();
-        }
-
-        if (globals->haveTauThermostat())
-          myNPTfm->setTauThermostat(globals->getTauThermostat());
-        else{
-          sprintf(painCave.errMsg,
-                  "SimSetup error: If you use an NPT\n"
-                  "    ensemble, you must set tauThermostat.\n");
-          painCave.isFatal = 1;
-          simError();
-        }
-
-        if (globals->haveTauBarostat())
-          myNPTfm->setTauBarostat(globals->getTauBarostat());
-        else{
-          sprintf(painCave.errMsg,
-                  "SimSetup error: If you use an NPT\n"
-                  "    ensemble, you must set tauBarostat.\n");
-          painCave.isFatal = 1;
-          simError();
-        }
-        break;
-
       default:
         sprintf(painCave.errMsg,
                 "SimSetup Error. Unrecognized ensemble in case statement.\n");