diff --git a/pySD/thermobinary_src/thermogen.py b/pySD/thermobinary_src/thermogen.py
index 8df3f3005..aa02d2ae4 100644
--- a/pySD/thermobinary_src/thermogen.py
+++ b/pySD/thermobinary_src/thermogen.py
@@ -399,7 +399,6 @@ def generate_qvap(self, zfulls, xfulls, PRESS, TEMP):
 
   def generate_thermo(self, gbxbounds, ndims, ntime):
 
-    ngridboxes = int(np.prod(ndims))
     zfulls, xfulls, yfulls = rgrid.fullcoords_forallgridboxes(gbxbounds,
                                                               ndims)
     PRESS, TEMP = self.hydrostatic_adiabatic_thermo(zfulls)
@@ -425,7 +424,6 @@ class ConstHydrostaticLapseRates:
 
   def generate_thermo(self, gbxbounds, ndims, ntime):
 
-    ngridboxes = int(np.prod(ndims))
     zfulls, xfulls, yfulls = rgrid.fullcoords_forallgridboxes(gbxbounds,
                                                               ndims)
     PRESS, TEMP = self.hydrostatic_adiabatic_thermo(zfulls)