From 90fd2874de4f64961009b628d3c1bd3bf9332e9e Mon Sep 17 00:00:00 2001
From: Roy Smart <roytsmart@gmail.com>
Date: Mon, 4 Nov 2024 16:50:21 -0700
Subject: [PATCH] background adjustments

---
 iris/sg/background/_background.py | 26 +++++++++++++-------------
 1 file changed, 13 insertions(+), 13 deletions(-)

diff --git a/iris/sg/background/_background.py b/iris/sg/background/_background.py
index 5d16401..89312ec 100644
--- a/iris/sg/background/_background.py
+++ b/iris/sg/background/_background.py
@@ -409,21 +409,21 @@ def function(x: na.CartesianNdVectorArray):
         guess=guess,
         step_size=na.CartesianNdVectorArray.from_components(
             components=dict(
-                amplitude=0.01 * u.DN / u.s,
-                shift=0.01 * (u.km / u.s) ** 2 / (u.DN / u.s),
-                width=0.01 * (u.km / u.s) ** 2 / (u.DN / u.s),
-                bias=0.001 * u.DN / u.s,
-                kappa=0.0001 / (u.DN / u.s),
-                slope=1e-6 * u.DN / u.s / (u.km / u.s) ** 2,
+                amplitude=0.001 * u.DN / u.s,
+                shift=0.1 * (u.km / u.s) ** 2 / (u.DN / u.s),
+                width=0.1 * (u.km / u.s) ** 2 / (u.DN / u.s),
+                bias=0.0001 * u.DN / u.s,
+                kappa=0.001 / (u.DN / u.s),
+                slope=1e-7 * u.DN / u.s / (u.km / u.s) ** 2,
             ),
         ),
         momentum=0.9,
         min_gradient=na.CartesianNdVectorArray.from_components(
             components=dict(
                 amplitude=0.5,
-                shift=0.5 * u.DN / u.s / (u.km / u.s),
-                width=0.5 * u.DN / u.s / (u.km / u.s),
-                kappa=0.5 * u.DN / u.s,
+                shift=0.05 * u.DN / u.s / (u.km / u.s),
+                width=0.05 * u.DN / u.s / (u.km / u.s),
+                kappa=0.05 * u.DN / u.s,
                 bias=1,
                 slope=110 * (u.km / u.s),
             )
@@ -500,8 +500,8 @@ def subtract_spectral_line(
                 bg.inputs.wavelength,
                 bg.inputs.position.y,
                 C=bg.outputs.value,
-                vmin=-5,
-                vmax=+5,
+                vmin=-1,
+                vmax=+1,
             )
             ax.set_xlabel(f"wavelength ({ax.get_xlabel()})")
             ax.set_ylabel(f"helioprojective $y$ ({ax.get_ylabel()})")
@@ -748,7 +748,7 @@ def estimate(
                 constrained_layout=True,
             )
             mappable = plt.cm.ScalarMappable(
-                norm=plt.Normalize(vmin=0, vmax=5),
+                norm=plt.Normalize(vmin=0, vmax=1),
             )
             ax[0].set_title("original")
             na.plt.pcolormesh(
@@ -804,7 +804,7 @@ def estimate(
             )
             ax.set_xlabel(f"Doppler velocity ({ax.get_xlabel()})")
             ax.set_ylabel(f"median spectral radiance ({ax.get_ylabel()})")
-            ax.set_ylim(top=10)
+            ax.set_ylim(top=1)
             ax.legend()
     """
     avg = average(