diff --git a/doc/src/compute_displace_atom.txt b/doc/src/compute_displace_atom.txt
index 39c301cf911a159089189ee02cef15a5ffae4d1e..00e5f696c118029431bf5a36395bcd906fe432ee 100644
--- a/doc/src/compute_displace_atom.txt
+++ b/doc/src/compute_displace_atom.txt
@@ -15,7 +15,7 @@ compute ID group-ID displace/atom :pre
 ID, group-ID are documented in "compute"_compute.html command :ulb,l
 displace/atom = style name of this compute command :l
 zero or more keyword/arg pairs may be appended :l
-keyword = {refresh} :
+keyword = {refresh} :l
   {replace} arg = name of per-atom variable :pre
 
 :ule
diff --git a/doc/src/fix_reax_bonds.txt b/doc/src/fix_reax_bonds.txt
index 54aa7faef802f3061f795c54a7449d6d93265d0b..3f8f237de1c7d574bfc214b9dfec2892e11b3793 100644
--- a/doc/src/fix_reax_bonds.txt
+++ b/doc/src/fix_reax_bonds.txt
@@ -34,6 +34,8 @@ written to {filename} on timesteps that are multiples of {Nevery},
 including timestep 0.  For time-averaged chemical species analysis,
 please see the "fix reaxc/c/species"_fix_reaxc_species.html command.
 
+The specified group-ID is ignored by this fix.
+
 The format of the output file should be reasonably self-explanatory.
 The meaning of the column header abbreviations is as follows:
 
diff --git a/python/lammps.py b/python/lammps.py
index e798ef60718d7cd9763641f799b0dd8330da1f53..417427eb4b715bb64cdabeceb65b93ae036f492c 100644
--- a/python/lammps.py
+++ b/python/lammps.py
@@ -5,7 +5,7 @@
 #
 #   Copyright (2003) Sandia Corporation.  Under the terms of Contract
 #   DE-AC04-94AL85000 with Sandia Corporation, the U.S. Government retains
-#   certain rights in this software.  This software is distributed under 
+#   certain rights in this software.  This software is distributed under
 #   the GNU General Public License.
 #
 #   See the README file in the top-level LAMMPS directory.
@@ -37,7 +37,7 @@ def get_ctypes_int(size):
     return c_int32
   elif size == 8:
     return c_int64
-  return c_int 
+  return c_int
 
 class MPIAbortException(Exception):
   def __init__(self, message):
@@ -47,7 +47,7 @@ class MPIAbortException(Exception):
     return repr(self.message)
 
 class lammps(object):
-  
+
   # detect if Python is using version of mpi4py that can pass a communicator
 
   has_mpi4py = False
@@ -71,7 +71,7 @@ class lammps(object):
 
     # if a pointer to a LAMMPS object is handed in,
     # all symbols should already be available
-    
+
     try:
       if ptr: self.lib = CDLL("",RTLD_GLOBAL)
     except:
@@ -84,7 +84,7 @@ class lammps(object):
     #   so that LD_LIBRARY_PATH does not need to be set for regular install
     # fall back to loading with a relative path,
     #   typically requires LD_LIBRARY_PATH to be set appropriately
-      
+
     if not self.lib:
       try:
         if not name: self.lib = CDLL(join(modpath,"liblammps.so"),RTLD_GLOBAL)
@@ -110,15 +110,15 @@ class lammps(object):
     self.lib.lammps_gather_atoms.argtypes = \
       [c_void_p,c_char_p,c_int,c_int,c_void_p]
     self.lib.lammps_gather_atoms.restype = None
-    
+
     self.lib.lammps_gather_atoms_concat.argtypes = \
       [c_void_p,c_char_p,c_int,c_int,c_void_p]
     self.lib.lammps_gather_atoms_concat.restype = None
-    
+
     self.lib.lammps_gather_atoms_subset.argtypes = \
       [c_void_p,c_char_p,c_int,c_int,c_int,POINTER(c_int),c_void_p]
     self.lib.lammps_gather_atoms_subset.restype = None
-    
+
     self.lib.lammps_scatter_atoms.argtypes = \
       [c_void_p,c_char_p,c_int,c_int,c_void_p]
     self.lib.lammps_scatter_atoms.restype = None
@@ -137,7 +137,7 @@ class lammps(object):
     #   just convert it to ctypes ptr and store in self.lmp
 
     if not ptr:
-      
+
       # with mpi4py v2, can pass MPI communicator to LAMMPS
       # need to adjust for type of MPI communicator object
       # allow for int (like MPICH) or void* (like OpenMPI)
@@ -211,7 +211,7 @@ class lammps(object):
     self.c_imageint = get_ctypes_int(self.extract_setting("imageint"))
 
   # shut-down LAMMPS instance
-  
+
   def __del__(self):
     if self.lmp and self.opened:
       self.lib.lammps_close(self.lmp)
@@ -230,7 +230,7 @@ class lammps(object):
     self.lib.lammps_file(self.lmp,file)
 
   # send a single command
-    
+
   def command(self,cmd):
     if cmd: cmd = cmd.encode()
     self.lib.lammps_command(self.lmp,cmd)
@@ -250,13 +250,13 @@ class lammps(object):
     cmds = [x.encode() for x in cmdlist if type(x) is str]
     args = (c_char_p * len(cmdlist))(*cmds)
     self.lib.lammps_commands_list(self.lmp,len(cmdlist),args)
-    
+
   # send a string of commands
 
   def commands_string(self,multicmd):
     if type(multicmd) is str: multicmd = multicmd.encode()
     self.lib.lammps_commands_string(self.lmp,c_char_p(multicmd))
-    
+
   # extract lammps type byte sizes
 
   def extract_setting(self, name):
@@ -265,7 +265,7 @@ class lammps(object):
     return int(self.lib.lammps_extract_setting(self.lmp,name))
 
   # extract global info
-    
+
   def extract_global(self,name,type):
     if name: name = name.encode()
     if type == 0:
@@ -277,7 +277,7 @@ class lammps(object):
     return ptr[0]
 
   # extract global info
-    
+
   def extract_box(self):
     boxlo = (3*c_double)()
     boxhi = (3*c_double)()
@@ -286,11 +286,11 @@ class lammps(object):
     xz = c_double()
     periodicity = (3*c_int)()
     box_change = c_int()
-    
+
     self.lib.lammps_extract_box(self.lmp,boxlo,boxhi,
                                 byref(xy),byref(yz),byref(xz),
                                 periodicity,byref(box_change))
-    
+
     boxlo = boxlo[:3]
     boxhi = boxhi[:3]
     xy = xy.value
@@ -298,9 +298,9 @@ class lammps(object):
     xz = xz.value
     periodicity = periodicity[:3]
     box_change = box_change.value
-    
+
     return boxlo,boxhi,xy,yz,xz,periodicity,box_change
-    
+
   # extract per-atom info
   # NOTE: need to insure are converting to/from correct Python type
   #   e.g. for Python list or NumPy or ctypes
@@ -318,7 +318,7 @@ class lammps(object):
     else: return None
     ptr = self.lib.lammps_extract_atom(self.lmp,name)
     return ptr
-    
+
   @property
   def numpy(self):
     if not self._numpy:
@@ -371,7 +371,7 @@ class lammps(object):
     return self._numpy
 
   # extract compute info
-  
+
   def extract_compute(self,id,style,type):
     if id: id = id.encode()
     if type == 0:
@@ -384,9 +384,14 @@ class lammps(object):
       ptr = self.lib.lammps_extract_compute(self.lmp,id,style,type)
       return ptr
     if type == 2:
-      self.lib.lammps_extract_compute.restype = POINTER(POINTER(c_double))
-      ptr = self.lib.lammps_extract_compute(self.lmp,id,style,type)
-      return ptr
+      if style == 0:
+        self.lib.lammps_extract_compute.restype = POINTER(c_int)
+        ptr = self.lib.lammps_extract_compute(self.lmp,id,style,type)
+        return ptr[0]
+      else:
+        self.lib.lammps_extract_compute.restype = POINTER(POINTER(c_double))
+        ptr = self.lib.lammps_extract_compute(self.lmp,id,style,type)
+        return ptr
     return None
 
   # extract fix info
@@ -466,7 +471,7 @@ class lammps(object):
     cboxlo = (3*c_double)(*boxlo)
     cboxhi = (3*c_double)(*boxhi)
     self.lib.lammps_reset_box(self.lmp,cboxlo,cboxhi,xy,yz,xz)
-    
+
   # return vector of atom properties gathered across procs
   # 3 variants to match src/library.cpp
   # name = atom property recognized by LAMMPS in atom->extract()
@@ -475,7 +480,7 @@ class lammps(object):
   # returned data is a 1d vector - doc how it is ordered?
   # NOTE: need to insure are converting to/from correct Python type
   #   e.g. for Python list or NumPy or ctypes
-  
+
   def gather_atoms(self,name,type,count):
     if name: name = name.encode()
     natoms = self.lib.lammps_get_natoms(self.lmp)
@@ -487,7 +492,7 @@ class lammps(object):
       self.lib.lammps_gather_atoms(self.lmp,name,type,count,data)
     else: return None
     return data
-    
+
   def gather_atoms_concat(self,name,type,count):
     if name: name = name.encode()
     natoms = self.lib.lammps_get_natoms(self.lmp)
@@ -519,7 +524,7 @@ class lammps(object):
   # assume data is of correct type and length, as created by gather_atoms()
   # NOTE: need to insure are converting to/from correct Python type
   #   e.g. for Python list or NumPy or ctypes
-  
+
   def scatter_atoms(self,name,type,count,data):
     if name: name = name.encode()
     self.lib.lammps_scatter_atoms(self.lmp,name,type,count,data)
diff --git a/src/MOLECULE/angle_cosine_squared.cpp b/src/MOLECULE/angle_cosine_squared.cpp
index c261f3482d2a24917cebf64bbd21ed6df1329f62..c83ba90a60f35b1db16fde1c3a2cb53c31eafa82 100644
--- a/src/MOLECULE/angle_cosine_squared.cpp
+++ b/src/MOLECULE/angle_cosine_squared.cpp
@@ -34,7 +34,11 @@ using namespace MathConst;
 
 /* ---------------------------------------------------------------------- */
 
-AngleCosineSquared::AngleCosineSquared(LAMMPS *lmp) : Angle(lmp) {}
+AngleCosineSquared::AngleCosineSquared(LAMMPS *lmp) : Angle(lmp)
+{
+  k = NULL;
+  theta0 = NULL;
+}
 
 /* ---------------------------------------------------------------------- */
 
diff --git a/src/MOLECULE/angle_harmonic.cpp b/src/MOLECULE/angle_harmonic.cpp
index c7fe8d7a1fb47e6ec36744e90b88b7722be6c0af..d28afd76d688ca2a0a1a8d39d0724dc84a5210b9 100644
--- a/src/MOLECULE/angle_harmonic.cpp
+++ b/src/MOLECULE/angle_harmonic.cpp
@@ -30,7 +30,11 @@ using namespace MathConst;
 
 /* ---------------------------------------------------------------------- */
 
-AngleHarmonic::AngleHarmonic(LAMMPS *lmp) : Angle(lmp) {}
+AngleHarmonic::AngleHarmonic(LAMMPS *lmp) : Angle(lmp)
+{
+  k = NULL;
+  theta0 = NULL;
+}
 
 /* ---------------------------------------------------------------------- */
 
diff --git a/src/MOLECULE/dihedral_helix.cpp b/src/MOLECULE/dihedral_helix.cpp
index f50895dbee7b2addd92d27311256d6b4e35932f3..ae23b77951dda47600972f86df05636db8757001 100644
--- a/src/MOLECULE/dihedral_helix.cpp
+++ b/src/MOLECULE/dihedral_helix.cpp
@@ -332,3 +332,13 @@ void DihedralHelix::read_restart(FILE *fp)
 
   for (int i = 1; i <= atom->ndihedraltypes; i++) setflag[i] = 1;
 }
+
+/* ----------------------------------------------------------------------
+   proc 0 writes to data file
+------------------------------------------------------------------------- */
+
+void DihedralHelix::write_data(FILE *fp)
+{
+  for (int i = 1; i <= atom->ndihedraltypes; i++)
+    fprintf(fp,"%d %g %g %g\n",i,aphi[i],bphi[i],cphi[i]);
+}
diff --git a/src/MOLECULE/dihedral_helix.h b/src/MOLECULE/dihedral_helix.h
index 47756a26721d8129493d534103d72c9de801184d..745ab841156164e095c0ffbc55f93705d83f7662 100644
--- a/src/MOLECULE/dihedral_helix.h
+++ b/src/MOLECULE/dihedral_helix.h
@@ -33,6 +33,7 @@ class DihedralHelix : public Dihedral {
   void coeff(int, char **);
   void write_restart(FILE *);
   void read_restart(FILE *);
+  void write_data(FILE *);
 
  protected:
   double *aphi,*bphi,*cphi;
diff --git a/src/MOLECULE/dihedral_multi_harmonic.cpp b/src/MOLECULE/dihedral_multi_harmonic.cpp
index 97dab8b50173aee61c1e0958ba3f3799d459865e..f6461abb6e9fbf3ef84fbb82935755d73639c751 100644
--- a/src/MOLECULE/dihedral_multi_harmonic.cpp
+++ b/src/MOLECULE/dihedral_multi_harmonic.cpp
@@ -331,3 +331,13 @@ void DihedralMultiHarmonic::read_restart(FILE *fp)
 
   for (int i = 1; i <= atom->ndihedraltypes; i++) setflag[i] = 1;
 }
+
+/* ----------------------------------------------------------------------
+   proc 0 writes to data file
+------------------------------------------------------------------------- */
+
+void DihedralMultiHarmonic::write_data(FILE *fp)
+{
+  for (int i = 1; i <= atom->ndihedraltypes; i++)
+    fprintf(fp,"%d %g %g %g %g %g\n",i,a1[i],a2[i],a3[i],a4[i],a5[i]);
+}
diff --git a/src/MOLECULE/dihedral_multi_harmonic.h b/src/MOLECULE/dihedral_multi_harmonic.h
index 28a38e61ab4a0f2f52e29ffe7a88f171cc967090..32f17e75ae4a79571a6cdca27e199405d25069c8 100644
--- a/src/MOLECULE/dihedral_multi_harmonic.h
+++ b/src/MOLECULE/dihedral_multi_harmonic.h
@@ -33,6 +33,7 @@ class DihedralMultiHarmonic : public Dihedral {
   void coeff(int, char **);
   void write_restart(FILE *);
   void read_restart(FILE *);
+  void write_data(FILE *);
 
  protected:
   double *a1,*a2,*a3,*a4,*a5;
diff --git a/src/USER-MISC/angle_cosine_shift.cpp b/src/USER-MISC/angle_cosine_shift.cpp
index d782497d0f6876774930f598bea24824e1eafeb7..a361db497081441da026375482b8dcaef575bd4e 100644
--- a/src/USER-MISC/angle_cosine_shift.cpp
+++ b/src/USER-MISC/angle_cosine_shift.cpp
@@ -34,7 +34,10 @@ using namespace MathConst;
 
 /* ---------------------------------------------------------------------- */
 
-AngleCosineShift::AngleCosineShift(LAMMPS *lmp) : Angle(lmp) {}
+AngleCosineShift::AngleCosineShift(LAMMPS *lmp) : Angle(lmp)
+{
+  kcost = NULL;
+}
 
 /* ---------------------------------------------------------------------- */
 
diff --git a/src/USER-MISC/angle_cosine_shift_exp.cpp b/src/USER-MISC/angle_cosine_shift_exp.cpp
index 0312d3b36d5003c90ccf60052ebaac60dbe4974b..c87c73171a9509fc8895ab9ac1caaa58f0b64dde 100644
--- a/src/USER-MISC/angle_cosine_shift_exp.cpp
+++ b/src/USER-MISC/angle_cosine_shift_exp.cpp
@@ -34,7 +34,16 @@ using namespace MathConst;
 
 /* ---------------------------------------------------------------------- */
 
-AngleCosineShiftExp::AngleCosineShiftExp(LAMMPS *lmp) : Angle(lmp) {}
+AngleCosineShiftExp::AngleCosineShiftExp(LAMMPS *lmp) : Angle(lmp)
+{
+  doExpansion = NULL;
+  umin = NULL;
+  a = NULL;
+  opt1 = NULL;
+  theta0 = NULL;
+  sint = NULL;
+  cost = NULL;
+}
 
 /* ---------------------------------------------------------------------- */
 
diff --git a/src/USER-MISC/angle_dipole.cpp b/src/USER-MISC/angle_dipole.cpp
index 0fa8d3674dadcfb6e376f8bc1410ca140a2f19c1..bcb631c61f065353cbb60761a4e693c28800693a 100644
--- a/src/USER-MISC/angle_dipole.cpp
+++ b/src/USER-MISC/angle_dipole.cpp
@@ -32,7 +32,11 @@ using namespace MathConst;
 
 /* ---------------------------------------------------------------------- */
 
-AngleDipole::AngleDipole(LAMMPS *lmp) : Angle(lmp) {}
+AngleDipole::AngleDipole(LAMMPS *lmp) : Angle(lmp)
+{
+  k = NULL;
+  gamma0 = NULL;
+}
 
 /* ---------------------------------------------------------------------- */
 
diff --git a/src/USER-MISC/angle_fourier.cpp b/src/USER-MISC/angle_fourier.cpp
index 85e3e8f0367811162427a8b6e52fb331fbe5b46c..e6cc1f1a7e9945a55091002529051649f78186a9 100644
--- a/src/USER-MISC/angle_fourier.cpp
+++ b/src/USER-MISC/angle_fourier.cpp
@@ -35,7 +35,13 @@ using namespace MathConst;
 
 /* ---------------------------------------------------------------------- */
 
-AngleFourier::AngleFourier(LAMMPS *lmp) : Angle(lmp) {}
+AngleFourier::AngleFourier(LAMMPS *lmp) : Angle(lmp)
+{
+  k = NULL;
+  C0 = NULL;
+  C1 = NULL;
+  C2 = NULL;
+}
 
 /* ---------------------------------------------------------------------- */
 
diff --git a/src/USER-MISC/angle_fourier_simple.cpp b/src/USER-MISC/angle_fourier_simple.cpp
index 83fc00c1166f1bf73faefc86ea4f2dcdfe03f5bc..8464fe815c9db0e9948ec9a0c7185a4409edf6e3 100644
--- a/src/USER-MISC/angle_fourier_simple.cpp
+++ b/src/USER-MISC/angle_fourier_simple.cpp
@@ -35,7 +35,12 @@ using namespace MathConst;
 
 /* ---------------------------------------------------------------------- */
 
-AngleFourierSimple::AngleFourierSimple(LAMMPS *lmp) : Angle(lmp) {}
+AngleFourierSimple::AngleFourierSimple(LAMMPS *lmp) : Angle(lmp)
+{
+  k = NULL;
+  C = NULL;
+  N = NULL;
+}
 
 /* ---------------------------------------------------------------------- */
 
diff --git a/src/USER-MISC/dihedral_cosine_shift_exp.cpp b/src/USER-MISC/dihedral_cosine_shift_exp.cpp
index 253909530865fbf074459329d322b22f9065fb58..82da173f8e2c3d14bce6344f0f9d8080effd9594 100644
--- a/src/USER-MISC/dihedral_cosine_shift_exp.cpp
+++ b/src/USER-MISC/dihedral_cosine_shift_exp.cpp
@@ -26,6 +26,7 @@
 #include "force.h"
 #include "update.h"
 #include "memory.h"
+#include "math_const.h"
 #include "error.h"
 
 using namespace LAMMPS_NS;
@@ -285,9 +286,9 @@ void DihedralCosineShiftExp::coeff(int narg, char **arg)
     doExpansion[i]=(fabs(a_)<0.001);
     umin[i]  = umin_;
     a[i]     = a_;
-    cost[i]  = cos(theta0_*3.14159265/180);
-    sint[i]  = sin(theta0_*3.14159265/180);
-    theta[i] = theta0_*3.14159265/180;
+    cost[i]  = cos(theta0_*MathConst::MY_PI/180.0);
+    sint[i]  = sin(theta0_*MathConst::MY_PI/180.0);
+    theta[i] = theta0_*MathConst::MY_PI/180.0;
 
     if (!doExpansion[i]) opt1[i]=umin_/(exp(a_)-1);
 
@@ -338,3 +339,14 @@ void DihedralCosineShiftExp::read_restart(FILE *fp)
     if (!doExpansion[i]) opt1[i]=umin[i]/(exp(a[i])-1);
   }
 }
+
+/* ----------------------------------------------------------------------
+   proc 0 writes to data file
+------------------------------------------------------------------------- */
+
+void DihedralCosineShiftExp::write_data(FILE *fp)
+{
+  for (int i = 1; i <= atom->ndihedraltypes; i++)
+    fprintf(fp,"%d %g %g %g\n",i,umin[i],
+            theta[i]*180.0/MathConst::MY_PI,a[i]);
+}
diff --git a/src/USER-MISC/dihedral_cosine_shift_exp.h b/src/USER-MISC/dihedral_cosine_shift_exp.h
index ca4d08f1915a8a54e666fb765a4092ead605284c..fc4be5049b49c920d40d4fe83995a3b89f2b87db 100644
--- a/src/USER-MISC/dihedral_cosine_shift_exp.h
+++ b/src/USER-MISC/dihedral_cosine_shift_exp.h
@@ -33,6 +33,7 @@ class DihedralCosineShiftExp : public Dihedral {
   void coeff(int, char **);
   void write_restart(FILE *);
   void read_restart(FILE *);
+  void write_data(FILE *);
 
  protected:
   bool *doExpansion;
diff --git a/src/USER-MISC/dihedral_quadratic.cpp b/src/USER-MISC/dihedral_quadratic.cpp
index 9674f818f1d472b4dc93772362a53447d0c02a68..1b64b52faf34e50792105314277727edeb9efb1b 100644
--- a/src/USER-MISC/dihedral_quadratic.cpp
+++ b/src/USER-MISC/dihedral_quadratic.cpp
@@ -333,3 +333,12 @@ void DihedralQuadratic::read_restart(FILE *fp)
   for (int i = 1; i <= atom->ndihedraltypes; i++) setflag[i] = 1;
 }
 
+/* ----------------------------------------------------------------------
+   proc 0 writes to data file
+------------------------------------------------------------------------- */
+
+void DihedralQuadratic::write_data(FILE *fp)
+{
+  for (int i = 1; i <= atom->ndihedraltypes; i++)
+    fprintf(fp,"%d %g %g \n",i,k[i],phi0[i]*180.0/MY_PI);
+}
diff --git a/src/USER-MISC/dihedral_quadratic.h b/src/USER-MISC/dihedral_quadratic.h
index 45c3393e76ddd4ff57b5fd3215ead54d361db090..edc29c3cf8b5b87a84d291450dacdeb359d05ac1 100644
--- a/src/USER-MISC/dihedral_quadratic.h
+++ b/src/USER-MISC/dihedral_quadratic.h
@@ -33,6 +33,7 @@ class DihedralQuadratic : public Dihedral {
   void coeff(int, char **);
   void write_restart(FILE *);
   void read_restart(FILE *);
+  void write_data(FILE *);
 
  protected:
   double *k,*phi0;
diff --git a/src/USER-MISC/improper_cossq.cpp b/src/USER-MISC/improper_cossq.cpp
index 72cf0e79df352b55ef12361828ad32abc33e4be8..2483840e4234f37c14636937e4d1efe4d0498ced 100644
--- a/src/USER-MISC/improper_cossq.cpp
+++ b/src/USER-MISC/improper_cossq.cpp
@@ -312,3 +312,13 @@ void ImproperCossq::read_restart(FILE *fp)
 
   for (int i = 1; i <= atom->nimpropertypes; i++) setflag[i] = 1;
 }
+
+/* ----------------------------------------------------------------------
+   proc 0 writes to data file
+------------------------------------------------------------------------- */
+
+void ImproperCossq::write_data(FILE *fp)
+{
+  for (int i = 1; i <= atom->nimpropertypes; i++)
+    fprintf(fp,"%d %g %g\n",i,k[i],chi[i]/MY_PI*180.0);
+}
diff --git a/src/USER-MISC/improper_cossq.h b/src/USER-MISC/improper_cossq.h
index 46715433cd8432e0306c08df0dde7d677eff8f1e..ea880adfaf40845d301893abeb63849fb13a93e6 100644
--- a/src/USER-MISC/improper_cossq.h
+++ b/src/USER-MISC/improper_cossq.h
@@ -33,6 +33,7 @@ class ImproperCossq : public Improper {
   void coeff(int, char **);
   void write_restart(FILE *);
   void read_restart(FILE *);
+  void write_data(FILE *);
 
  protected:
   double *k, *chi;
diff --git a/src/USER-MISC/improper_distance.cpp b/src/USER-MISC/improper_distance.cpp
index 1c073cc791a4b792b4a1d46ddddeca7ffa354cf2..2edf37ec5c7b712c3c46a3da330ed4482edcd91d 100644
--- a/src/USER-MISC/improper_distance.cpp
+++ b/src/USER-MISC/improper_distance.cpp
@@ -223,7 +223,6 @@ void ImproperDistance::coeff(int narg, char **arg)
   int count = 0;
   for (int i = ilo; i <= ihi; i++) {
     k[i] = k_one;
-    //chi[i] = chi_one/180.0 * PI;
     chi[i] = chi_one;
     setflag[i] = 1;
     count++;
@@ -259,3 +258,13 @@ void ImproperDistance::read_restart(FILE *fp)
 
   for (int i = 1; i <= atom->nimpropertypes; i++) setflag[i] = 1;
 }
+
+/* ----------------------------------------------------------------------
+   proc 0 writes to data file
+------------------------------------------------------------------------- */
+
+void ImproperDistance::write_data(FILE *fp)
+{
+  for (int i = 1; i <= atom->nimpropertypes; i++)
+    fprintf(fp,"%d %g %g\n",i,k[i],chi[i]);
+}
diff --git a/src/USER-MISC/improper_distance.h b/src/USER-MISC/improper_distance.h
index 0f281803473ecd9624c02e2fbd8aad4a634ef857..57e4d671e94320503bd1497ed18849366cd983b3 100644
--- a/src/USER-MISC/improper_distance.h
+++ b/src/USER-MISC/improper_distance.h
@@ -33,6 +33,7 @@ class ImproperDistance : public Improper {
   void coeff(int, char **);
   void write_restart(FILE *);
   void read_restart(FILE *);
+  void write_data(FILE *);
 
  private:
   double *k,*chi;
diff --git a/src/USER-MISC/improper_fourier.cpp b/src/USER-MISC/improper_fourier.cpp
index 90596a1967feca597b9feb3041e99ef8b29692a9..11478d08eab135f70baa9ab6d69a762b4a302672 100644
--- a/src/USER-MISC/improper_fourier.cpp
+++ b/src/USER-MISC/improper_fourier.cpp
@@ -343,3 +343,13 @@ void ImproperFourier::read_restart(FILE *fp)
 
   for (int i = 1; i <= atom->nimpropertypes; i++) setflag[i] = 1;
 }
+
+/* ----------------------------------------------------------------------
+   proc 0 writes to data file
+------------------------------------------------------------------------- */
+
+void ImproperFourier::write_data(FILE *fp)
+{
+  for (int i = 1; i <= atom->nimpropertypes; i++)
+    fprintf(fp,"%d %g %g %g %g %d\n",i,k[i],C0[i],C1[i],C2[i],all[i]);
+}
diff --git a/src/USER-MISC/improper_fourier.h b/src/USER-MISC/improper_fourier.h
index 51421d6379667268fc486d632959db3a8ac26cfa..0525c45494b56741b100d2e545c21364f5e8955a 100644
--- a/src/USER-MISC/improper_fourier.h
+++ b/src/USER-MISC/improper_fourier.h
@@ -33,6 +33,7 @@ class ImproperFourier : public Improper {
   void coeff(int, char **);
   void write_restart(FILE *);
   void read_restart(FILE *);
+  void write_data(FILE *);
 
  protected:
   double *k, *C0, *C1, *C2;
diff --git a/src/USER-MISC/improper_ring.cpp b/src/USER-MISC/improper_ring.cpp
index 878c6dfc5e0d6fadb081187c93a606b5dea72df4..70124b2ed1d91de9d900aacdea1c0c828303c7fc 100644
--- a/src/USER-MISC/improper_ring.cpp
+++ b/src/USER-MISC/improper_ring.cpp
@@ -337,3 +337,13 @@ void ImproperRing::read_restart(FILE *fp)
 
   for (int i = 1; i <= atom->nimpropertypes; i++) setflag[i] = 1;
 }
+
+/* ----------------------------------------------------------------------
+   proc 0 writes to data file
+------------------------------------------------------------------------- */
+
+void ImproperRing::write_data(FILE *fp)
+{
+  for (int i = 1; i <= atom->nimpropertypes; i++)
+    fprintf(fp,"%d %g %g\n",i,k[i],acos(chi[i])/MY_PI*180.0);
+}
diff --git a/src/USER-MISC/improper_ring.h b/src/USER-MISC/improper_ring.h
index f721e8ff9eb71842f113ea7ee4f6120ed7569bd0..c31329816f806f573a044f77962eeab0c1cd418d 100644
--- a/src/USER-MISC/improper_ring.h
+++ b/src/USER-MISC/improper_ring.h
@@ -33,6 +33,7 @@ class ImproperRing : public Improper {
   void coeff(int, char **);
   void write_restart(FILE *);
   void read_restart(FILE *);
+  void write_data(FILE *);
 
  protected:
   double *k,*chi;
diff --git a/src/USER-MISC/pair_coul_shield.cpp b/src/USER-MISC/pair_coul_shield.cpp
index 5b321fe8ecd27447085d9f9c0f53874f30b4b0f7..deec7da81faaa4a2557f7e40bd21c1389c6aa0e2 100644
--- a/src/USER-MISC/pair_coul_shield.cpp
+++ b/src/USER-MISC/pair_coul_shield.cpp
@@ -231,6 +231,8 @@ void PairCoulShield::init_style()
 {
   if (!atom->q_flag)
     error->all(FLERR,"Pair style coul/shield requires atom attribute q");
+  if (!atom->molecule_flag)
+    error->all(FLERR,"Pair style coul/shield requires atom attribute molecule");
 
   neighbor->request(this,instance_me);
 }
diff --git a/src/USER-MISC/pair_ilp_graphene_hbn.cpp b/src/USER-MISC/pair_ilp_graphene_hbn.cpp
index 59dbc6142eca8424a1ef921e2111af13c253adcb..30ee2e7a5ae3c3aa42843f2ba07c1c8f36b7dbb9 100644
--- a/src/USER-MISC/pair_ilp_graphene_hbn.cpp
+++ b/src/USER-MISC/pair_ilp_graphene_hbn.cpp
@@ -636,7 +636,9 @@ void PairILPGrapheneHBN::calc_normal()
 void PairILPGrapheneHBN::init_style()
 {
   if (force->newton_pair == 0)
-    error->all(FLERR,"Pair style ILP requires newton pair on");
+    error->all(FLERR,"Pair style ilp/graphene/hbn requires newton pair on");
+  if (!atom->molecule_flag)
+    error->all(FLERR,"Pair style ilp/graphene/hbn requires atom attribute molecule");
 
   // need a full neighbor list, including neighbors of ghosts
 
diff --git a/src/USER-MISC/pair_kolmogorov_crespi_full.cpp b/src/USER-MISC/pair_kolmogorov_crespi_full.cpp
index b691150c46b7d5e30ac51ea243891d4d82bbf64c..be0e81d48d9e986e38247aae32bb60b0ac4fb377 100644
--- a/src/USER-MISC/pair_kolmogorov_crespi_full.cpp
+++ b/src/USER-MISC/pair_kolmogorov_crespi_full.cpp
@@ -640,7 +640,9 @@ void PairKolmogorovCrespiFull::calc_normal()
 void PairKolmogorovCrespiFull::init_style()
 {
   if (force->newton_pair == 0)
-    error->all(FLERR,"Pair style KC requires newton pair on");
+    error->all(FLERR,"Pair style kolmolgorov/crespi/full requires newton pair on");
+  if (!atom->molecule_flag)
+    error->all(FLERR,"Pair style kolmolgorov/crespi/full requires atom attribute molecule");
 
   // need a full neighbor list, including neighbors of ghosts
 
diff --git a/src/fix_restrain.cpp b/src/fix_restrain.cpp
index 2f75902d2cb1455aa0aa17975ca9f4eba50f9dd9..0f6c86ba5bc610f6e3d63e7921a47eca9432ec2e 100644
--- a/src/fix_restrain.cpp
+++ b/src/fix_restrain.cpp
@@ -261,7 +261,7 @@ void FixRestrain::restrain_bond(int m)
   if (r > 0.0) fbond = -2.0*rk/r;
   else fbond = 0.0;
 
-  energy = rk*dr;
+  energy += rk*dr;
 
   // apply force to each of 2 atoms
 
@@ -368,7 +368,7 @@ void FixRestrain::restrain_angle(int m)
   dtheta = acos(c) - target[m];
   tk = k * dtheta;
 
-  energy = tk*dtheta;
+  energy += tk*dtheta;
 
   a = -2.0 * tk * s;
   a11 = a*c / rsq1;
@@ -549,7 +549,7 @@ void FixRestrain::restrain_dihedral(int m)
   df1 *= -mult;
   p += 1.0;
 
-  energy = k * p;
+  energy += k * p;
 
   fg = vb1x*vb2xm + vb1y*vb2ym + vb1z*vb2zm;
   hg = vb3x*vb2xm + vb3y*vb2ym + vb3z*vb2zm;
diff --git a/src/library.cpp b/src/library.cpp
index 032b28d178c8a1911964bce545f06d1dad783b07..865bf2a0a5d41bda4e5f5c8ec349e1b2903b96af 100644
--- a/src/library.cpp
+++ b/src/library.cpp
@@ -480,10 +480,13 @@ void *lammps_extract_atom(void *ptr, char *name)
      compute's internal data structure for the entity
      caller should cast it to (double *) for a scalar or vector
      caller should cast it to (double **) for an array
-   for per-atom or local data, returns a pointer to the
+   for per-atom or local vector/array data, returns a pointer to the
      compute's internal data structure for the entity
      caller should cast it to (double *) for a vector
      caller should cast it to (double **) for an array
+   for local data, accessing scalar data for the compute (type = 0),
+   returns a pointer that should be cast to (int *) which points to
+   an int with the number of local rows, i.e. the length of the local array.
    returns a void pointer to the compute's internal data structure
      for the entity which the caller can cast to the proper data type
    returns a NULL if id is not recognized or style/type not supported
@@ -541,6 +544,11 @@ void *lammps_extract_compute(void *ptr, char *id, int style, int type)
 
     if (style == 2) {
       if (!compute->local_flag) return NULL;
+      if (type == 0) {
+        if (compute->invoked_local != lmp->update->ntimestep)
+          compute->compute_local();
+        return (void *) &compute->size_local_rows;
+      }
       if (type == 1) {
         if (compute->invoked_local != lmp->update->ntimestep)
           compute->compute_local();