test_radialI: save diag fig

test/onedim_test.py

@@ -4,14 +4,14 @@
 from mpol.onedim import radialI, radialV
 from mpol.plot import plot_image
 
-def test_radialI(mock_1d_image_model):
+def test_radialI(mock_1d_image_model, tmp_path):
     # obtain a 1d radial brightness profile I(r) from an image    
 
     r, i, i2d, _, _, geom, coords = mock_1d_image_model
 
     rtest, itest = radialI(i2d, coords, geom, bins=None)
 
-    _, ax = plt.subplots(ncols=2)
+    fig, ax = plt.subplots(ncols=2)
 
     plot_image(i2d, extent=coords.img_ext, ax=ax[0], clab='Jy sr$^{-2}$')
 
@@ -20,6 +20,9 @@ def test_radialI(mock_1d_image_model):
     ax[1].plot(r, i, 'k', label='truth')
     ax[1].plot(rtest, itest, 'r.-', label='result')
 
+    fig.savefig(tmp_path / "test_radialI.png", dpi=300)
+    plt.close("all")
+
     expected = [
         5.79780326e+10, 2.47990375e+10, 4.19794053e+09, 1.63165616e+10, 
         2.56197452e+10, 1.86014523e+10, 1.39800643e+10, 1.14935415e+10,