diff --git a/src/USER-CGDNA/fix_nve_dot.cpp b/src/USER-CGDNA/fix_nve_dot.cpp
index 9a90af365d..68ea1999ab 100644
--- a/src/USER-CGDNA/fix_nve_dot.cpp
+++ b/src/USER-CGDNA/fix_nve_dot.cpp
@@ -146,7 +146,7 @@ void FixNVEDot::initial_integrate(int vflag)
 void FixNVEDot::final_integrate()
 {
 
-  double *shape,*quat;
+  double *quat;
   double fquat[4],conjqm[4];
   double conjqm_dot_quat;
 
@@ -171,7 +171,6 @@ void FixNVEDot::final_integrate()
 
       dthlfm = dthlf / rmass[i];
       quat = bonus[ellipsoid[i]].quat;
-      shape = bonus[ellipsoid[i]].shape;
 
       // update momentum
       v[i][0] += dthlfm * f[i][0];
diff --git a/src/USER-CGDNA/fix_nve_dotc_langevin.cpp b/src/USER-CGDNA/fix_nve_dotc_langevin.cpp
index bc3cfe09f8..53cf4b08b9 100644
--- a/src/USER-CGDNA/fix_nve_dotc_langevin.cpp
+++ b/src/USER-CGDNA/fix_nve_dotc_langevin.cpp
@@ -263,7 +263,7 @@ void FixNVEDotcLangevin::initial_integrate(int vflag)
 void FixNVEDotcLangevin::final_integrate()
 {
 
-  double *shape,*quat;
+  double *quat;
   double fquat[4],conjqm[4];
   double conjqm_dot_quat;
 
@@ -288,7 +288,6 @@ void FixNVEDotcLangevin::final_integrate()
 
       dthlfm = dthlf / rmass[i];
       quat = bonus[ellipsoid[i]].quat;
-      shape = bonus[ellipsoid[i]].shape;
 
       // update momentum by 1/2 step
       v[i][0] += dthlfm * f[i][0];