diff --git a/pyro/incompressible/simulation.py b/pyro/incompressible/simulation.py
index 7c2e4930e..d79da2578 100644
--- a/pyro/incompressible/simulation.py
+++ b/pyro/incompressible/simulation.py
@@ -427,7 +427,7 @@ def dovis(self):
             0.5*(u.ip(1) - u.ip(-1))/myg.dx + \
             0.5*(v.jp(1) - v.jp(-1))/myg.dy
 
-        _, axes = plt.subplots(nrows=2, ncols=2, num=1)
+        _, axes = plt.subplots(nrows=2, ncols=2, num=1, clear=True)
         plt.subplots_adjust(hspace=0.25)
 
         fields = [u, v, vort, divU]
diff --git a/pyro/lm_atm/simulation.py b/pyro/lm_atm/simulation.py
index 04ac8e7d1..3f3b4b749 100644
--- a/pyro/lm_atm/simulation.py
+++ b/pyro/lm_atm/simulation.py
@@ -643,7 +643,7 @@ def dovis(self):
 
         vort.v()[:, :] = dv - du
 
-        _, axes = plt.subplots(nrows=2, ncols=2, num=1)
+        _, axes = plt.subplots(nrows=2, ncols=2, num=1, clear=True)
         plt.subplots_adjust(hspace=0.25)
 
         fields = [rho, magvel, vort, rhoprime]