diff --git a/src/fix_gyration.cpp b/src/fix_gyration.cpp
index 0b8ed1668a..aee1221130 100644
--- a/src/fix_gyration.cpp
+++ b/src/fix_gyration.cpp
@@ -42,8 +42,8 @@ FixGyration::FixGyration(LAMMPS *lmp, int narg, char **arg) :
   }
 
   if (me == 0) {
-    fprintf(fp,"Radius-of-gyration for group %s\n",group->names[igroup]);
-    fprintf(fp,"TimeStep Rg\n");
+    fprintf(fp,"# Radius-of-gyration for group %s\n",group->names[igroup]);
+    fprintf(fp,"# TimeStep Rg\n");
   }
 }
 
diff --git a/src/fix_msd.cpp b/src/fix_msd.cpp
index f593c25cc8..916fb57242 100644
--- a/src/fix_msd.cpp
+++ b/src/fix_msd.cpp
@@ -46,9 +46,9 @@ FixMSD::FixMSD(LAMMPS *lmp, int narg, char **arg) :
   }
 
   if (me == 0) {
-    fprintf(fp,"Mean-squared Displacement for group %s\n",
+    fprintf(fp,"# Mean-squared Displacement for group %s\n",
 	    group->names[igroup]);
-    fprintf(fp,"TimeStep x y z total\n");
+    fprintf(fp,"# TimeStep x y z total\n");
   }
 
   // perform initial allocation of atom-based array
diff --git a/src/fix_print.cpp b/src/fix_print.cpp
index 1ac89724ba..b4848d75aa 100644
--- a/src/fix_print.cpp
+++ b/src/fix_print.cpp
@@ -68,7 +68,7 @@ FixPrint::FixPrint(LAMMPS *lmp, int narg, char **arg) :
 
   // print header into file
 
-  if (fp && me == 0) fprintf(fp,"Fix print output for fix %s\n",id);
+  if (fp && me == 0) fprintf(fp,"# Fix print output for fix %s\n",id);
 
   copy = new char[MAXLINE];
   work = new char[MAXLINE];
diff --git a/src/fix_rdf.cpp b/src/fix_rdf.cpp
index 1e94b53768..dc78212f08 100644
--- a/src/fix_rdf.cpp
+++ b/src/fix_rdf.cpp
@@ -260,7 +260,7 @@ void FixRDF::end_of_step()
     double constant = 4.0*PI / (3.0*domain->xprd*domain->yprd*domain->zprd);
     int irdf;
 
-    fprintf(fp,"%s %d \n","TIMESTEP", update->ntimestep);
+    fprintf(fp,"# Timestep %d\n",update->ntimestep);
     fprintf(fp,"%s","r");
 
     for (int i = 1; i <= atom->ntypes; i++)
@@ -300,7 +300,7 @@ void FixRDF::end_of_step()
 
     if (update->ntimestep + nevery > update->endstep) {
 
-      fprintf(fp,"%s \n","RUN AVERAGE");
+      fprintf(fp,"# Run Average\n");
       fprintf(fp,"%s","r");
 
       for (int i = 1; i <= atom->ntypes; i++)