diff --git a/examples/dpa/plot_dpa_resid.py b/examples/dpa/plot_dpa_resid.py
index 70aba8a2..fa9adca4 100644
--- a/examples/dpa/plot_dpa_resid.py
+++ b/examples/dpa/plot_dpa_resid.py
@@ -1,17 +1,19 @@
 # Licensed under a 3-clause BSD style license - see LICENSE.rst
-import xija
-import numpy as np
 import matplotlib.pyplot as plt
-from Ska.Matplotlib import pointpair
+import numpy as np
+from ska_matplotlib import pointpair
 
-start = '2010:001'
-stop = '2011:345'
+import xija
+from xija.get_model_spec import get_xija_model_spec
+
+start = "2025:001"
+stop = "2025:030"
 
-msid = '1dpamzt'
-model_spec = 'dpa.json'
+msid = "1dpamzt"
+model_spec, _ = get_xija_model_spec("dpa")
 
-model = xija.ThermalModel('dpa', start=start, stop=stop,
-                          model_spec=model_spec)
+model = xija.ThermalModel("dpa", start=start, stop=stop, model_spec=model_spec)
+model.comp["dpa0"].set_data(0.0)
 model.make()
 model.calc()
 
@@ -21,17 +23,17 @@
 xscatter = np.random.uniform(-0.2, 0.2, size=len(dpa.dvals))
 yscatter = np.random.uniform(-0.2, 0.2, size=len(dpa.dvals))
 plt.clf()
-plt.plot(dpa.dvals + xscatter, resid + yscatter, '.', ms=1.0, alpha=1)
-plt.xlabel('{} telemetry (degC)'.format(msid.upper()))
-plt.ylabel('Data - Model (degC)')
-plt.title('Residual vs. Data ({} - {})'.format(start, stop))
+plt.plot(dpa.dvals + xscatter, resid + yscatter, ".", ms=1.0, alpha=1)
+plt.xlabel("{} telemetry (degC)".format(msid.upper()))
+plt.ylabel("Data - Model (degC)")
+plt.title("Residual vs. Data ({} - {})".format(start, stop))
 
-bins = np.arange(6, 26.1, 2.0)
+bins = np.arange(10, 35, 2.0)
 r1 = []
 r99 = []
 ns = []
 xs = []
-for x0, x1 in zip(bins[:-1], bins[1:]):
+for x0, x1 in zip(bins[:-1], bins[1:], strict=True):
     ok = (dpa.dvals >= x0) & (dpa.dvals < x1)
     val1, val99 = np.percentile(resid[ok], [1, 99])
     xs.append((x0 + x1) / 2)
@@ -43,18 +45,22 @@
 r1pp = pointpair(r1)
 r99pp = pointpair(r99)
 
-plt.plot(xspp, r1pp, '-r')
-plt.plot(xspp, r99pp, '-r', label='1% and 99% limits')
+plt.plot(xspp, r1pp, "-r")
+plt.plot(xspp, r99pp, "-r", label="1% and 99% limits")
 plt.grid()
 plt.ylim(-8, 14)
 plt.xlim(5, 31)
 
-plt.plot([5, 31], [3.5, 3.5], 'g--', alpha=1, label='+/- 3.5 degC')
-plt.plot([5, 31], [-3.5, -3.5], 'g--', alpha=1)
-for x, n, y in zip(xs, ns, r99):
-    plt.text(x, max(y + 1, 5), 'N={}'.format(n),
-         rotation='vertical', va='bottom', ha='center')
-
-plt.legend(loc='upper right')
+plt.plot([5, 31], [3.5, 3.5], "g--", alpha=1, label="+/- 3.5 degC")
+plt.plot([5, 31], [-3.5, -3.5], "g--", alpha=1)
+for x, n, y in zip(xs, ns, r99, strict=True):
+    plt.text(
+        x,
+        max(y + 1, 5),
+        "N={}".format(n),
+        rotation="vertical",
+        va="bottom",
+        ha="center",
+    )
 
-plt.savefig('dpa_resid_{}_{}.png'.format(start, stop))
+plt.legend(loc="upper right")