diff --git a/final_plots.py b/final_plots.py
index 95ce3b32fb5e7861d5bf2beb59a3e6c727da89ef..94533c77330445ae65a936c482f8b3526156e656 100644
--- a/final_plots.py
+++ b/final_plots.py
@@ -401,6 +401,22 @@ def export_interp(
             results[f"gdc_{method}_{order}{suffix}_ispline_old"] = bisplev(y_arr, x_arr, idc_tck_old, dy=1)
         except:
             settings.logger.exception(f"in spline interpolation of Idc for {order}{suffix} {method} (old)")
+        try:
+            results[f"idc_{method}_{order}{suffix}_old"] = griddata(
+                    xy_data,
+                    g_data.dc_current,
+                    mesh,
+                    method="cubic")
+        except:
+            settings.logger.exception(f"in griddata interpolation of Idc for {order}{suffix} {method} (old)")
+        try:
+            results[f"iac_{method}_{order}{suffix}_old"] = griddata(
+                    xy_data,
+                    2*g_data.ac_current_abs,
+                    mesh,
+                    method="cubic")
+        except:
+            settings.logger.exception(f"in griddata interpolation of Iac for {order}{suffix} {method} (old)")
         if has_current:
             try:
                 results[f"idc_{method}_{order}{suffix}"] = griddata(
diff --git a/gen_data.py b/gen_data.py
index be69e3e9358b61d16ea5b4d7e8e4e9ab87626372..da9b9ad14ed20a1b4a793497b172c1e5bcba215b 100644
--- a/gen_data.py
+++ b/gen_data.py
@@ -244,7 +244,7 @@ python -m frtrg_kondo.gen_data \\
             type=int, default=0,
             help = "Set initial condition for gammaL to 0")
     method_group.add_argument("--improved_initial_conditions", metavar="bool",
-            type=int, default=0,
+            type=int, default=1,
             help = "Set initial condition for dgammaL/dE to linear response estimate")
     method_group.add_argument("--include_Ga", metavar="int", type=int,
             default=1, choices=(0,1),
diff --git a/gen_pulse_data.py b/gen_pulse_data.py
index 3fca9e9959c024fc46d3306aea437be3082ae4fc..54413d2ff7967dabe28122a0bf98b101aaa2e72f 100644
--- a/gen_pulse_data.py
+++ b/gen_pulse_data.py
@@ -67,8 +67,12 @@ def main():
     method_group = parser.add_argument_group(title="Method")
     method_group.add_argument("--method", type=str, required=True, choices=('J', 'mu'),
             help = "J: include all time dependence in coupling by unitary transformation.\nmu: describe time dependence by Floquet matrix for chemical potentials.")
-    method_group.add_argument("--simplified_initial_conditions", metavar="bool", type=bool, default=False,
+    method_group.add_argument("--simplified_initial_conditions", metavar="bool",
+            type=int, default=0,
             help = "Set initial condition for gammaL to 0")
+    method_group.add_argument("--improved_initial_conditions", metavar="bool",
+            type=int, default=1,
+            help = "Set initial condition for dgammaL/dE to linear response estimate")
     method_group.add_argument("--include_Ga", metavar="bool", type=bool, default=False,
             help = "include vertex parameter Ga in RG equations")
     method_group.add_argument("--solve_integral_exactly", metavar="bool", type=bool, default=False,