--- trunk/OOPSE/libmdtools/SimInfo.cpp	2003/07/16 21:30:56	626
+++ trunk/OOPSE/libmdtools/SimInfo.cpp	2003/08/07 00:47:33	669
@@ -33,6 +33,7 @@ SimInfo::SimInfo(){
   the_integrator = NULL;
   setTemp = 0;
   thermalTime = 0.0;
+  currentTime = 0.0;
   rCut = 0.0;
   ecr = 0.0;
   est = 0.0;
@@ -56,6 +57,16 @@ SimInfo::SimInfo(){
   wrapMeSimInfo( this );
 }
 
+SimInfo::~SimInfo(){
+
+  map<string, GenericData*>::iterator i;
+  
+  for(i = properties.begin(); i != properties.end(); i++)
+    delete (*i).second;
+
+  
+}
+
 void SimInfo::setBox(double newBox[3]) {
   
   int i, j;
@@ -291,6 +302,8 @@ void SimInfo::calcBoxL( void ){
   dsq = dx*dx + dy*dy + dz*dz;
   boxL[2] = sqrt( dsq );
   if( (0.5 * boxL[2]) < maxCutoff ) maxCutoff = 0.5 * boxL[2];
+  
+  checkCutOffs();
 
 }
 
@@ -457,11 +470,13 @@ void SimInfo::checkCutOffs( void ){
 
   int cutChanged = 0;
 
+
+
   if( boxIsInit ){
     
     //we need to check cutOffs against the box
-    
-    if( maxCutoff > rCut ){
+   
+    if(( maxCutoff > rCut )&&(usePBC)){
       if( rCut < origRcut ){
 	rCut = origRcut;
 	if (rCut > maxCutoff) rCut = maxCutoff;
@@ -490,7 +505,7 @@ void SimInfo::checkCutOffs( void ){
     }
 
 
-    if (rCut > maxCutoff) {
+    if ((rCut > maxCutoff)&&(usePBC)) {
       sprintf( painCave.errMsg,
                "New Box size is setting the long range cutoff radius "
 	       "to %lf\n",
@@ -527,4 +542,53 @@ void SimInfo::checkCutOffs( void ){
 
   oldEcr = ecr;
   oldRcut = rCut;
+}
+
+void SimInfo::addProperty(GenericData* prop){
+
+  map<string, GenericData*>::iterator result;
+  result = properties.find(prop->getID());
+  
+  //we can't simply use  properties[prop->getID()] = prop,
+  //it will cause memory leak if we already contain a propery which has the same name of prop
+  
+  if(result != properties.end()){
+    
+    delete (*result).second;
+    (*result).second = prop;
+      
+  }
+  else{
+
+    properties[prop->getID()] = prop;
+
+  }
+    
 }
+
+GenericData* SimInfo::getProperty(const string& propName){
+ 
+  map<string, GenericData*>::iterator result;
+  
+  //string lowerCaseName = ();
+  
+  result = properties.find(propName);
+  
+  if(result != properties.end()) 
+    return (*result).second;  
+  else   
+    return NULL;  
+}
+
+vector<GenericData*> SimInfo::getProperties(){
+
+  vector<GenericData*> result;
+  map<string, GenericData*>::iterator i;
+  
+  for(i = properties.begin(); i != properties.end(); i++)
+    result.push_back((*i).second);
+    
+  return result;
+}
+
+