diff --git a/NonLinearSolver/nlsolver/nonLinearMechSolver.cpp b/NonLinearSolver/nlsolver/nonLinearMechSolver.cpp
index 811ad2ed6efc132357fdf22289568709aa529f01..b921f14e3580bdb929021321210773a80e3cc63f 100644
--- a/NonLinearSolver/nlsolver/nonLinearMechSolver.cpp
+++ b/NonLinearSolver/nlsolver/nonLinearMechSolver.cpp
@@ -1150,7 +1150,7 @@ void nonLinearMechSolver::oneStepPostSolve(const double curtime, const int numst
       int fieldNum = it->first;
       double val;
       this->getTotalIntegralByVolumeIntegral(fieldNum,val);
-      fprintf(it->second,"%e ; %e; \n",curtime,val);
+      fprintf(it->second,"%e;%e \n",curtime,val);
       fflush(it->second);
     }
   }
@@ -1160,7 +1160,7 @@ void nonLinearMechSolver::oneStepPostSolve(const double curtime, const int numst
       int fieldNum = it->first;
       double val;
       this->getHomogenizedPropertyByVolumeIntegral(fieldNum,val);
-      fprintf(it->second,"%e ; %e; \n",curtime,val);
+      fprintf(it->second,"%e;%e \n",curtime,val);
       fflush(it->second);
     }
   }
@@ -1172,7 +1172,7 @@ void nonLinearMechSolver::oneStepPostSolve(const double curtime, const int numst
       this->getHomogenizedPropertyByVolumeIntegralOnActiveDissipationDomain(fieldNum,valCur,IPStateBase::current);
       this->getHomogenizedPropertyByVolumeIntegralOnActiveDissipationDomain(fieldNum,valPrev,IPStateBase::previous);
       double curVal = lastVal + valCur - valPrev;
-      fprintf(it->second,"%e ; %e;\n",curtime,curVal);
+      fprintf(it->second,"%e;%e\n",curtime,curVal);
       fflush(it->second);
       lastVal = curVal;
     }