diff --git a/docs/source/shape_editor_install.rst b/docs/source/shape_editor_install.rst
index 232278cb..b568ddf5 100644
--- a/docs/source/shape_editor_install.rst
+++ b/docs/source/shape_editor_install.rst
@@ -87,6 +87,7 @@ https://blfauger.gitlabpages.inria.fr/nice/install.html.
     cp Makefile.TEMPLATE Makefile
     make -j
     make nice_imas_inv_muscle3
+    make nice_imas_dir_muscle3
 
-If this was successful, you should find the program ``nice_imas_inv_muscle3`` in the
-``nice/run`` folder.
+If this was successful, you should find the programs ``nice_imas_inv_muscle3`` and
+``nice_imas_dir_muscle3`` in the ``nice/run`` folder.
diff --git a/waveform_editor/gui/shape_editor.py b/waveform_editor/gui/shape_editor.py
index e7b6228f..20001580 100644
--- a/waveform_editor/gui/shape_editor.py
+++ b/waveform_editor/gui/shape_editor.py
@@ -1,4 +1,5 @@
 import importlib.resources
+import logging
 import xml.etree.ElementTree as ET
 
 import imas
@@ -14,6 +15,8 @@
 from waveform_editor.shape_editor.plasma_properties import PlasmaProperties
 from waveform_editor.shape_editor.plasma_shape import PlasmaShape
 
+logger = logging.getLogger(__name__)
+
 
 def _reactive_title(title, is_valid):
     return title if is_valid else f"{title} ⚠️"
@@ -37,40 +40,52 @@ def __init__(self):
         self.plasma_properties = PlasmaProperties()
         self.coil_currents = CoilCurrents()
         self.nice_plotter = NicePlotter(
-            self.communicator, self.plasma_shape, self.plasma_properties
+            communicator=self.communicator,
+            plasma_shape=self.plasma_shape,
+            plasma_properties=self.plasma_properties,
         )
         self.nice_settings = settings.nice
 
-        self.xml_params = ET.fromstring(
+        self.xml_params_inv = ET.fromstring(
+            importlib.resources.files("waveform_editor.shape_editor.xml_param")
+            .joinpath("inverse_param.xml")
+            .read_text()
+        )
+        self.xml_params_dir = ET.fromstring(
             importlib.resources.files("waveform_editor.shape_editor.xml_param")
-            .joinpath("param.xml")
+            .joinpath("direct_param.xml")
             .read_text()
         )
 
         # UI Configuration
         button_start = pn.widgets.Button(name="Run", on_click=self.submit)
         button_start.disabled = (
-            self.plasma_shape.param.has_shape.rx.not_()
+            (
+                self.plasma_shape.param.has_shape.rx.not_()
+                & self.nice_settings.param.is_inverse_mode.rx()
+            )
             | self.plasma_properties.param.has_properties.rx.not_()
-            | param.rx(self.nice_settings.required_params_filled).rx.not_()
+            | self.nice_settings.param.are_required_filled.rx.not_()
         )
         button_stop = pn.widgets.Button(name="Stop", on_click=self.stop_nice)
-        buttons = pn.Row(button_start, button_stop)
+        nice_mode_radio = pn.widgets.RadioBoxGroup.from_param(
+            self.nice_settings.param.mode, inline=True, margin=(15, 20, 0, 20)
+        )
+        buttons = pn.Row(button_start, button_stop, nice_mode_radio)
 
         # Accordion does not allow dynamic titles, so use separate card for each option
         options = pn.Column(
             self._create_card(
                 self.nice_settings.panel,
                 "NICE Configuration",
-                is_valid=param.rx(self.nice_settings.required_params_filled),
-            ),
-            self._create_card(
-                pn.Param(self.nice_plotter, show_name=False), "Plotting Parameters"
+                is_valid=self.nice_settings.param.are_required_filled.rx(),
             ),
+            self._create_card(self.nice_plotter, "Plotting Parameters"),
             self._create_card(
                 self.plasma_shape,
                 "Plasma Shape",
                 is_valid=self.plasma_shape.param.has_shape,
+                visible=self.nice_settings.param.is_inverse_mode.rx(),
             ),
             self._create_card(
                 pn.Column(self.plasma_properties, self.nice_plotter.profiles_pane),
@@ -91,7 +106,7 @@ def __init__(self):
             ),
         )
 
-    def _create_card(self, panel_object, title, is_valid=None):
+    def _create_card(self, panel_object, title, is_valid=None, visible=True):
         """Create a collapsed card containing a panel object and a title.
 
         Args:
@@ -99,6 +114,7 @@ def _create_card(self, panel_object, title, is_valid=None):
             title: The title to give the card.
             is_valid: If supplied, binds the card title to update reactively using
                 `_reactive_title`.
+            visible: Whether the card is visible.
         """
         if is_valid:
             title = param.bind(_reactive_title, title=title, is_valid=is_valid)
@@ -107,6 +123,7 @@ def _create_card(self, panel_object, title, is_valid=None):
             title=title,
             sizing_mode="stretch_width",
             collapsed=True,
+            visible=visible,
         )
         return card
 
@@ -169,9 +186,10 @@ def _create_equilibrium(self):
         equilibrium.time_slice.resize(1)
         equilibrium.vacuum_toroidal_field.b0.resize(1)
 
-        # Fill plasma shape
-        equilibrium.time_slice[0].boundary.outline.r = self.plasma_shape.outline_r
-        equilibrium.time_slice[0].boundary.outline.z = self.plasma_shape.outline_z
+        # Only fill plasma shape for NICE inverse mode
+        if self.nice_settings.is_inverse_mode:
+            equilibrium.time_slice[0].boundary.outline.r = self.plasma_shape.outline_r
+            equilibrium.time_slice[0].boundary.outline.z = self.plasma_shape.outline_z
 
         # Fill plasma properties
         equilibrium.vacuum_toroidal_field.r0 = self.plasma_properties.r0
@@ -194,14 +212,21 @@ async def submit(self, event=None):
         description IDSs and an input equilibrium IDS."""
 
         self.coil_currents.fill_pf_active(self.pf_active)
+        if self.nice_settings.is_direct_mode:
+            xml_params = self.xml_params_dir
+        else:
+            xml_params = self.xml_params_inv
+            self.coil_currents.update_fixed_coils_in_xml(xml_params)
+
         # Update XML parameters:
-        self.coil_currents.update_fixed_coils_in_xml(self.xml_params)
-        self.xml_params.find("verbose").text = str(self.nice_settings.verbose)
+        xml_params.find("verbose").text = str(self.nice_settings.verbose)
         equilibrium = self._create_equilibrium()
         if not self.communicator.running:
-            await self.communicator.run()
+            await self.communicator.run(
+                is_direct_mode=self.nice_settings.is_direct_mode
+            )
         await self.communicator.submit(
-            ET.tostring(self.xml_params, encoding="unicode"),
+            ET.tostring(xml_params, encoding="unicode"),
             equilibrium.serialize(),
             self.pf_active.serialize(),
             self.pf_passive.serialize(),
@@ -210,7 +235,9 @@ async def submit(self, event=None):
         )
         self.coil_currents.sync_ui_with_pf_active(self.communicator.pf_active)
 
-    async def stop_nice(self, event):
+    @param.depends("nice_settings.mode", watch=True)
+    async def stop_nice(self, event=None):
+        logger.info("Stopping NICE...")
         await self.communicator.close()
 
     def __panel__(self):
diff --git a/waveform_editor/settings.py b/waveform_editor/settings.py
index ab09f606..ce030efa 100644
--- a/waveform_editor/settings.py
+++ b/waveform_editor/settings.py
@@ -16,17 +16,23 @@
 
 
 class NiceSettings(param.Parameterized):
-    REQUIRED = (
-        "executable",
+    INVERSE_MODE = "NICE Inverse"
+    DIRECT_MODE = "NICE Direct"
+
+    BASE_REQUIRED = (
         "md_pf_active",
         "md_pf_passive",
         "md_wall",
         "md_iron_core",
     )
-    executable = param.String(
-        label="NICE executable path",
+    inv_executable = param.String(
+        label="NICE inverse executable path",
         doc="Path to NICE inverse IMAS MUSCLE3 executable",
     )
+    dir_executable = param.String(
+        label="NICE direct executable path",
+        doc="Path to NICE direct IMAS MUSCLE3 executable",
+    )
     environment = param.Dict(
         default={},
         label="NICE environment variables",
@@ -37,10 +43,32 @@ class NiceSettings(param.Parameterized):
     md_wall = param.String(label="'wall' machine description URI")
     md_iron_core = param.String(label="'iron_core' machine description URI")
     verbose = param.Integer(label="NICE verbosity (set to 1 for more verbose output)")
+    mode = param.Selector(
+        objects=[INVERSE_MODE, DIRECT_MODE], default=INVERSE_MODE, precedence=-1
+    )
+    are_required_filled = param.Boolean(precedence=-1)
+    is_direct_mode = param.Boolean(precedence=-1)
+    is_inverse_mode = param.Boolean(precedence=-1)
+
+    @param.depends("mode", watch=True, on_init=True)
+    def set_mode_flags(self):
+        self.is_direct_mode = self.mode == self.DIRECT_MODE
+        self.is_inverse_mode = self.mode == self.INVERSE_MODE
+
+    @param.depends(
+        *BASE_REQUIRED, "inv_executable", "dir_executable", "mode", watch=True
+    )
+    def check_required_params_filled(self):
+        base_ready = all(getattr(self, p) for p in self.BASE_REQUIRED)
 
-    @param.depends(*REQUIRED)
-    def required_params_filled(self):
-        return all(getattr(self, required) for required in self.REQUIRED)
+        if not base_ready:
+            self.are_required_filled = False
+            return
+
+        if self.mode == self.INVERSE_MODE:
+            self.are_required_filled = bool(self.inv_executable)
+        else:
+            self.are_required_filled = bool(self.dir_executable)
 
     def apply_settings(self, params):
         """Update parameters from a dictionary, skipping unknown keys."""
@@ -51,24 +79,33 @@ def apply_settings(self, params):
         self.param.update(**params)
 
     def to_dict(self):
-        """Returns a dictionary representation of current parameter values."""
-        return {p: getattr(self, p) for p in self.param if p != "name"}
+        """Returns a dictionary representation of current parameter values, excluding
+        params with a precendence of -1."""
+        result = {}
+        for p in self.param:
+            param_obj = self.param[p]
+            if p != "name" and param_obj.precedence != -1:
+                result[p] = getattr(self, p)
+        return result
 
     def panel(self):
         items = []
 
         for p in self.param:
             if p == "name":
-                pass
-            elif p in self.REQUIRED:
-                items.append(
-                    pn.Row(
-                        pn.Param(self.param[p], show_name=False),
-                        WarningIndicator(visible=self.param[p].rx.not_()),
-                    )
-                )
-            else:
-                items.append(pn.Param(self.param[p], show_name=False))
+                continue
+
+            # Add warning indicator if required parameter is not filled
+            is_inv_required = p == "inv_executable" and self.is_inverse_mode
+            is_dir_required = p == "dir_executable" and self.is_direct_mode
+            is_base_required = p in self.BASE_REQUIRED
+
+            row_content = [pn.Param(self.param[p], show_name=False)]
+            if is_inv_required or is_dir_required or is_base_required:
+                warning = WarningIndicator(visible=self.param[p].rx.not_())
+                row_content.append(warning)
+
+            items.append(pn.Row(*row_content))
 
         return pn.Column(*items)
 
diff --git a/waveform_editor/shape_editor/coil_currents.py b/waveform_editor/shape_editor/coil_currents.py
index d561b3e5..0f009d14 100644
--- a/waveform_editor/shape_editor/coil_currents.py
+++ b/waveform_editor/shape_editor/coil_currents.py
@@ -5,14 +5,17 @@
 import param
 from panel.viewable import Viewer
 
+from waveform_editor.settings import settings
+
 
 class CoilCurrents(Viewer):
     coil_ui = param.List(
         doc="List of tuples containing the checkboxes and sliders for the coil currents"
     )
 
-    def __init__(self):
-        super().__init__()
+    def __init__(self, **params):
+        super().__init__(**params)
+        self.nice_settings = settings.nice
         self.sliders_ui = pn.Column(visible=self.param.coil_ui.rx.bool())
         guide_message = pn.pane.Markdown(
             "_To fix a coil to a specific current, enable the checkbox and provide "
@@ -44,13 +47,15 @@ def create_ui(self, pf_active):
         new_coil_ui = []
         for coil in pf_active.coil:
             coil_current = coil.current
-            checkbox = pn.widgets.Checkbox(value=False, margin=(30, 10, 10, 10))
+            checkbox = pn.widgets.Checkbox(
+                margin=(30, 10, 10, 10),
+                disabled=self.nice_settings.param.is_direct_mode,
+            )
             slider = pn.widgets.EditableFloatSlider(
                 name=f"{coil.name} Current [{coil_current.metadata.units}]",
                 value=coil_current.data[0] if coil_current.data.has_value else 0.0,
                 start=-5e4,
                 end=5e4,
-                disabled=checkbox.param.value.rx.not_(),
                 format="0",
                 width=450,
             )
@@ -67,9 +72,8 @@ def fill_pf_active(self, pf_active):
             pf_active: pf_active IDS to update the coil currents for.
         """
         for i, coil_ui in enumerate(self.coil_ui):
-            checkbox, slider = coil_ui.objects
-            if checkbox.value:
-                pf_active.coil[i].current.data = np.array([slider.value])
+            _, slider = coil_ui.objects
+            pf_active.coil[i].current.data = np.array([slider.value])
 
     def sync_ui_with_pf_active(self, pf_active):
         """Synchronize UI sliders with the current values from the pf_active IDS.
diff --git a/waveform_editor/shape_editor/nice_integration.py b/waveform_editor/shape_editor/nice_integration.py
index f117171c..d4345ce8 100644
--- a/waveform_editor/shape_editor/nice_integration.py
+++ b/waveform_editor/shape_editor/nice_integration.py
@@ -1,4 +1,4 @@
-"""Module that enables running the NICE inverse equilibrium solver through MUSCLE3."""
+"""Module that enables running the NICE equilibrium solver through MUSCLE3."""
 
 import asyncio
 import multiprocessing
@@ -18,7 +18,8 @@
 
 from waveform_editor.settings import settings
 
-_muscle3_configuration = """
+# YMMSL configuration for NICE inverse mode
+_muscle3_inv_configuration = """
 ymmsl_version: v0.1
 model:
     name: shape_editor
@@ -42,9 +43,33 @@
     nice_inv.xml_path: {xml_path}
 """
 
+# YMMSL configuration for NICE direct mode
+_muscle3_dir_configuration = """
+ymmsl_version: v0.1
+model:
+    name: shape_editor
+    components:
+        shape_editor:
+            implementation: shape_editor
+        nice_dir:
+            implementation: nice_dir
+
+    conduits:
+        shape_editor.equilibrium_out: nice_dir.equilibrium_in
+        shape_editor.pf_passive_out: nice_dir.pf_passive_in
+        shape_editor.pf_active_out: nice_dir.pf_active_in
+        shape_editor.iron_core_out: nice_dir.iron_core_in
+        shape_editor.wall_out: nice_dir.wall_in
+        nice_dir.equilibrium_out: shape_editor.equilibrium_in
+
+settings:
+    muscle_profile_level: none  # Disable profiling
+    nice_dir.xml_path: {xml_path}
+"""
+
 
 class NiceIntegration(param.Parameterized):
-    """Main API for running NICE, submitting inverse problems and getting the resulting
+    """Main API for running NICE, submitting problems and getting the resulting
     equilibrium back.
     """
 
@@ -118,7 +143,7 @@ async def close(self):
         self.closing = self.running = False
         self._update_state()
 
-    async def run(self):
+    async def run(self, is_direct_mode=False):
         """Start NICE and the controlling processes."""
         if self.running:
             raise RuntimeError("Already running!")
@@ -130,7 +155,7 @@ async def run(self):
         self.manager_pipe, pipe = multiprocessing.Pipe()
         self.manager = multiprocessing.Process(
             target=run_muscle_manager,
-            args=[pipe, self.xml_config_file.name],
+            args=[pipe, self.xml_config_file.name, is_direct_mode],
             name="MUSCLE Manager",
         )
         self.manager.start()
@@ -140,7 +165,7 @@ async def run(self):
         self.communicator_pipe, pipe = multiprocessing.Pipe()
         self.communicator = multiprocessing.Process(
             target=run_muscle3_communicator,
-            args=[manager_location, pipe],
+            args=[manager_location, pipe, is_direct_mode],
             name="NICE Communicator",
         )
         self.communicator.start()
@@ -148,15 +173,22 @@ async def run(self):
         # NICE
         nice_env = os.environ.copy()
         nice_env.update(settings.nice.environment)
-        nice_env["MUSCLE_INSTANCE"] = "nice_inv"
         nice_env["MUSCLE_MANAGER"] = manager_location
 
-        self.terminal.write(f"{os.getcwd()}$ {settings.nice.executable}\n")
+        if is_direct_mode:
+            executable = settings.nice.dir_executable
+            nice_env["MUSCLE_INSTANCE"] = "nice_dir"
+        else:
+            executable = settings.nice.inv_executable
+            nice_env["MUSCLE_INSTANCE"] = "nice_inv"
+
+        self.terminal.write(f"{os.getcwd()}$ {executable}\n")
+
         loop = asyncio.get_running_loop()
         try:
             self.nice_transport, self.nice_protocol = await loop.subprocess_exec(
                 self.create_communicator_protocol,
-                settings.nice.executable,
+                executable,
                 env=nice_env,
                 stdin=subprocess.DEVNULL,
             )
@@ -267,6 +299,7 @@ def pipe_data_received(self, fd, data):
 def run_muscle3_communicator(
     server_location: str,
     pipe: multiprocessing.connection.Connection,
+    is_direct_mode: bool,
 ):
     """Run MUSCLE3 actor for communicating with NICE."""
     os.environ["MUSCLE_INSTANCE"] = "shape_editor"
@@ -280,8 +313,11 @@ def run_muscle3_communicator(
             "wall_out",
             "iron_core_out",
         ],
-        Operator.S: ["equilibrium_in", "pf_active_in"],
+        Operator.S: ["equilibrium_in"],
     }
+    if not is_direct_mode:
+        ports[Operator.S].append("pf_active_in")
+
     instance = Instance(ports)
 
     while instance.reuse_instance():
@@ -299,13 +335,19 @@ def run_muscle3_communicator(
 
             # Wait for nice to process
             eq = instance.receive("equilibrium_in").data
-            pfa = instance.receive("pf_active_in").data
+            if not is_direct_mode:
+                pfa = instance.receive("pf_active_in").data
             pipe.send((eq, pfa))
 
 
-def run_muscle_manager(pipe: multiprocessing.connection.Connection, xml_path: str):
-    """Run the muscle_manager"""
-    config = ymmsl.load(_muscle3_configuration.format(xml_path=xml_path))
+def run_muscle_manager(
+    pipe: multiprocessing.connection.Connection, xml_path: str, is_direct_mode: bool
+):
+    """Run the muscle_manager with a given configuration."""
+    config_str = (
+        _muscle3_dir_configuration if is_direct_mode else _muscle3_inv_configuration
+    )
+    config = ymmsl.load(config_str.format(xml_path=xml_path))
     manager = Manager(config)
     server_location = manager.get_server_location()
     pipe.send(server_location)
diff --git a/waveform_editor/shape_editor/nice_plotter.py b/waveform_editor/shape_editor/nice_plotter.py
index d96caaec..f9f2b021 100644
--- a/waveform_editor/shape_editor/nice_plotter.py
+++ b/waveform_editor/shape_editor/nice_plotter.py
@@ -8,7 +8,9 @@
 import param
 import scipy
 from imas.ids_toplevel import IDSToplevel
+from panel.viewable import Viewer
 
+from waveform_editor.settings import NiceSettings, settings
 from waveform_editor.shape_editor.nice_integration import NiceIntegration
 from waveform_editor.shape_editor.plasma_properties import PlasmaProperties
 from waveform_editor.shape_editor.plasma_shape import PlasmaShape
@@ -17,13 +19,14 @@
 logger = logging.getLogger(__name__)
 
 
-class NicePlotter(param.Parameterized):
+class NicePlotter(Viewer):
     # Input data, use negative precedence to hide from the UI
     communicator = param.ClassSelector(class_=NiceIntegration, precedence=-1)
     wall = param.ClassSelector(class_=IDSToplevel, precedence=-1)
     pf_active = param.ClassSelector(class_=IDSToplevel, precedence=-1)
     plasma_shape = param.ClassSelector(class_=PlasmaShape, precedence=-1)
     plasma_properties = param.ClassSelector(class_=PlasmaProperties, precedence=-1)
+    nice_settings = param.ClassSelector(class_=NiceSettings, precedence=-1)
 
     # Plot parameters
     show_contour = param.Boolean(default=True, label="Show contour lines")
@@ -45,13 +48,8 @@ class NicePlotter(param.Parameterized):
     PROFILE_WIDTH = 350
     PROFILE_HEIGHT = 350
 
-    def __init__(self, communicator, plasma_shape, plasma_properties, **params):
-        super().__init__(
-            **params,
-            communicator=communicator,
-            plasma_shape=plasma_shape,
-            plasma_properties=plasma_properties,
-        )
+    def __init__(self, **params):
+        super().__init__(**params)
         self.DEFAULT_OPTS = hv.opts.Overlay(
             xlim=(0, 13),
             ylim=(-10, 10),
@@ -59,6 +57,7 @@ def __init__(self, communicator, plasma_shape, plasma_properties, **params):
             xlabel="r [m]",
             ylabel="z [m]",
         )
+        self.nice_settings = settings.nice
         self.CONTOUR_OPTS = hv.opts.Contours(
             cmap="viridis",
             colorbar=True,
@@ -90,6 +89,15 @@ def __init__(self, communicator, plasma_shape, plasma_properties, **params):
         self.profiles_pane = pn.pane.HoloViews(
             profiles_plot, width=self.PROFILE_WIDTH, height=self.PROFILE_HEIGHT
         )
+        self.panel_layout = pn.Param(
+            self.param,
+            show_name=False,
+            widgets={
+                "show_desired_shape": {
+                    "visible": self.nice_settings.param.is_inverse_mode
+                }
+            },
+        )
 
     @pn.depends("plasma_properties.profile_updated")
     def _plot_profiles(self):
@@ -124,9 +132,15 @@ def _plot_profiles(self):
             hv.opts.Overlay(title="Plasma Profiles"), hv.opts.Curve(framewise=True)
         )
 
-    @pn.depends("plasma_shape.shape_updated", "show_desired_shape")
+    @pn.depends(
+        "plasma_shape.shape_updated", "show_desired_shape", "nice_settings.mode"
+    )
     def _plot_plasma_shape(self):
-        if not self.show_desired_shape or not self.plasma_shape.has_shape:
+        if (
+            self.nice_settings.is_direct_mode
+            or not self.show_desired_shape
+            or not self.plasma_shape.has_shape
+        ):
             return hv.Overlay([hv.Curve([]).opts(self.DESIRED_SHAPE_OPTS)])
 
         r = self.plasma_shape.outline_r
@@ -388,3 +402,6 @@ def _plot_xo_points(self):
             hover_tooltips=[("", "X-point")],
         )
         return o_scatter * x_scatter
+
+    def __panel__(self):
+        return self.panel_layout
diff --git a/waveform_editor/shape_editor/xml_param/direct_param.xml b/waveform_editor/shape_editor/xml_param/direct_param.xml
new file mode 100644
index 00000000..2bb5615d
--- /dev/null
+++ b/waveform_editor/shape_editor/xml_param/direct_param.xml
@@ -0,0 +1,262 @@
+<?xml version="1.0"?>
+<!-- Taken from: https://gitlab.inria.fr/blfauger/nice
+     Path: nice/run/iwrap/param/dir/param.xml
+     Author: Blaise Faugeras
+     License: LGPL -->
+<parameters>
+  <verbose>1</verbose>
+  <!-- overwritten by shape editor -->
+  <outputTxtFiles>0</outputTxtFiles>
+  <outputRestartFiles>0</outputRestartFiles>
+  <algoMode>11</algoMode>
+  <algoMeshDomain>2</algoMeshDomain>
+  <algoVacTHLRSolver>0</algoVacTHLRSolver>
+  <algoVacTHonly>0</algoVacTHonly>
+  <algoDirect>1</algoDirect>
+  <algoIpFreeScaleInitLambda>0</algoIpFreeScaleInitLambda>
+  <algoRecon>2</algoRecon>
+  <algoReconUseInputIp>1</algoReconUseInputIp>
+  <algoReconInterf>0</algoReconInterf>
+  <algoReconPolar>0</algoReconPolar>
+  <algoReconPressure>0</algoReconPressure>
+  <algoReconMSE>0</algoReconMSE>
+  <algoStartFromScratch>1</algoStartFromScratch>
+  <algoStartFromScratchReconAB>1</algoStartFromScratchReconAB>
+  <algoStartPsiFromInData>0</algoStartPsiFromInData>
+  <algoWithIp>1</algoWithIp>
+  <iterMaxDir>40</iterMaxDir>
+  <epsStopDir>1.0e-10</epsStopDir>
+  <iterMaxDirInitRecon>2</iterMaxDirInitRecon>
+  <iterMaxDirInitInv>0</iterMaxDirInitInv>
+  <iterMaxRecon>20</iterMaxRecon>
+  <epsStopRecon>1.0e-10</epsStopRecon>
+  <iterMaxInv>40</iterMaxInv>
+  <epsStopInv>1.0e-10</epsStopInv>
+  <iterMaxEvo>50</iterMaxEvo>
+  <epsStopEvo>1.0e-10</epsStopEvo>
+  <useDefaultDirectInit>1</useDefaultDirectInit>
+  <radiusInit>2.25</radiusInit>
+  <rcenterInit>6.2</rcenterInit>
+  <zcenterInit>0.0</zcenterInit>
+  <useMobileWall>0</useMobileWall>
+  <refinedZone1>1</refinedZone1>
+  <nPointsZone1>4</nPointsZone1>
+  <rPointsZone1>4.45 5.39 5.78 4.70</rPointsZone1>
+  <zPointsZone1>-3.5 -3.86 -3.20 -2.81</zPointsZone1>
+  <rInPointZone1>5.1</rInPointZone1>
+  <zInPointZone1>-3.3</zInPointZone1>
+  <zone1TriangleArea>0.0005</zone1TriangleArea>
+  <refinedZone2>0</refinedZone2>
+  <nPointsZone2>4</nPointsZone2>
+  <rPointsZone2>2.44 1.91 1.94 2.44</rPointsZone2>
+  <zPointsZone2>0.79 0.573 0.49 0.66</zPointsZone2>
+  <rInPointZone2>2.18</rInPointZone2>
+  <zInPointZone2>0.6</zInPointZone2>
+  <zone2TriangleArea>0.00003</zone2TriangleArea>
+  <refinedZone3>1</refinedZone3>
+  <nPointsZone3>10</nPointsZone3>
+  <rPointsZone3>6.0 6.2 6.2 6.4 6.6 6.6 6.4 6.4 6.2 6.0</rPointsZone3>
+  <zPointsZone3>0.2 0.2 0.4 0.2 0.2 1.0 1.0 0.8 1.0 1.0</zPointsZone3>
+  <rInPointZone3>6.1</rInPointZone3>
+  <zInPointZone3>0.3</zInPointZone3>
+  <zone3TriangleArea>0.0005</zone3TriangleArea>
+  <vacuumTriangleArea>0.002</vacuumTriangleArea>
+  <airTriangleArea>0.04</airTriangleArea>
+  <airInTriangleArea>0.01</airInTriangleArea>
+  <ironTriangleArea>0.005</ironTriangleArea>
+  <coilTriangleArea>0.005</coilTriangleArea>
+  <vesselTriangleArea>0.005</vesselTriangleArea>
+  <passiveStructTriangleArea>0.005</passiveStructTriangleArea>
+  <useNewCOCOSManager>1</useNewCOCOSManager>
+  <inoutCOCOS>17</inoutCOCOS>
+  <change_coil_sign>1</change_coil_sign>
+  <change_psimeas_sign>0</change_psimeas_sign>
+  <change_bmeas_sign>0</change_bmeas_sign>
+  <change_BpSensorAngle_sign>0</change_BpSensorAngle_sign>
+  <change_polarmeas_sign>0</change_polarmeas_sign>
+  <change_dirichlet_sign>0</change_dirichlet_sign>
+  <change_dpressure_dpsi_sign>0</change_dpressure_dpsi_sign>
+  <change_f_df_dpsi_sign>0</change_f_df_dpsi_sign>
+  <change_TPI>0</change_TPI>
+  <change_psi_sign>0</change_psi_sign>
+  <inCOCOS>17</inCOCOS>
+  <outCOCOS>17</outCOCOS>
+  <n_profile_1d>50</n_profile_1d>
+  <allErrorBar>0</allErrorBar>
+  <r0>6.2</r0>
+  <!-- plasma current density  -->
+  <Afunc>0</Afunc>
+  <Aknots_n>11</Aknots_n>
+  <Aknots>0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1</Aknots>
+  <Aabg>2 0.6 1.4</Aabg>
+  <Bfunc>0</Bfunc>
+  <Bknots_n>11</Bknots_n>
+  <Bknots>0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1</Bknots>
+  <Babg>2 0.4 1.4</Babg>
+  <useInputAB>1</useInputAB>
+  <!-- electron density  -->
+  <Cfunc>0</Cfunc>
+  <Cknots_n>11</Cknots_n>
+  <Cknots>0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1</Cknots>
+  <lambdaNe>1e20</lambdaNe>
+  <Cparam_n>4</Cparam_n>
+  <Cparam>3 1 0 1</Cparam>
+  <!-- electron temperature  -->
+  <Dfunc>0</Dfunc>
+  <Dknots_n>11</Dknots_n>
+  <Dknots>0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1</Dknots>
+  <lambdaTe>25</lambdaTe>
+  <agmM_Te>1 1.5 0 1</agmM_Te>
+  <!--   -->
+  <epsA>1.0e-1</epsA>
+  <sA>1 1 1 1 1 1 1 1 1 1 1</sA>
+  <penalA>1.0e12</penalA>
+  <penalNegativeA>0.0e10</penalNegativeA>
+  <epsB>1.0e-1</epsB>
+  <sB>1 1 1 1 1 1 1 1 1 1 1</sB>
+  <penalB>1.0e12</penalB>
+  <epsC>1.0e-1</epsC>
+  <sC>1 1 1 1 1 1 1 1 1 1 1</sC>
+  <penalC>1.0e12</penalC>
+  <epsCparam>1.0e3 1.0e3 1.0e0 1.0e12 1.0e0</epsCparam>
+  <epsD>1.0e-1</epsD>
+  <sD>1 1 1 1 1 1 1 1 1 1 1 1</sD>
+  <penalD>0.0e20</penalD>
+  <penalD0>0.0e0</penalD0>
+  <!-- magnetics -->
+  <neumannPercent>0.1</neumannPercent>
+  <abserrB>1e-3</abserrB>
+  <relerrB>0.01</relerrB>
+  <force_default_errB>0</force_default_errB>
+  <n_B_special_weight>0</n_B_special_weight>
+  <B_special_weight_index>0</B_special_weight_index>
+  <B_special_weight>0</B_special_weight>
+  <abserrF>0.01</abserrF>
+  <relerrF>0.01</relerrF>
+  <force_default_errF>0</force_default_errF>
+  <n_F_special_weight>0</n_F_special_weight>
+  <F_special_weight_index>0</F_special_weight_index>
+  <F_special_weight>0</F_special_weight>
+  <add_fictitious_full_loop_if_not_exist>0</add_fictitious_full_loop_if_not_exist>
+  <!-- interferometry -->
+  <chord_dl>0.01</chord_dl>
+  <constructLOS>0</constructLOS>
+  <abserrNe>1e18</abserrNe>
+  <relerrNe>0.1</relerrNe>
+  <force_default_errNe>0</force_default_errNe>
+  <n_Ne_special_weight>0</n_Ne_special_weight>
+  <Ne_special_weight_index>0 1 2 3 4 5 6 7 8 9 10 11 12 13 14</Ne_special_weight_index>
+  <!--Ne_special_weight>1 1 1 1 1 1 1 1 1 1 1 1 1 1 1</Ne_special_weight-->
+  <Ne_special_weight>0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1</Ne_special_weight>
+  <!--Ne_special_weight>0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01</Ne_special_weight-->
+  <!-- polarimetry -->
+  <wavelength>0.195</wavelength>
+  <C1>2.42e-20</C1>
+  <C3>5.23e-19</C3>
+  <abserrPolar>0.001</abserrPolar>
+  <relerrPolar>0.1</relerrPolar>
+  <force_default_errPolar>0</force_default_errPolar>
+  <n_Polar_special_weight>0</n_Polar_special_weight>
+  <Polar_special_weight_index>0 1 2 3 4 5 6 7 8 9 10 11 12 13 14</Polar_special_weight_index>
+  <!--Polar_special_weight>1 1 1 1 1 1 1 1 1 1 1 1 1 1 1</Polar_special_weight-->
+  <Polar_special_weight>0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1</Polar_special_weight>
+  <abserrS>0.001</abserrS>
+  <relerrS>0.1</relerrS>
+  <force_default_errS>0</force_default_errS>
+  <n_S_special_weight>0</n_S_special_weight>
+  <S_special_weight_index>0 1 2 3 4 5 6 7 8 9 10 11 12 13 14</S_special_weight_index>
+  <!--S_special_weight>0</S_special_weight-->
+  <!--S_special_weight>0.1 0 0 0 0 0 0 0 0 0 0 0 0 0 0</S_special_weight-->
+  <S_special_weight>0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1</S_special_weight>
+  <!--S_special_weight>0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01</S_special_weight-->
+  <!--S_special_weight>0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001</S_special_weight-->
+  <S012_weight>1 1 1</S012_weight>
+  <!-- pressure -->
+  <abserrPressure>1.0e3</abserrPressure>
+  <relerrPressure>0.1</relerrPressure>
+  <force_default_errPressure>0</force_default_errPressure>
+  <n_Pressure_special_weight>0</n_Pressure_special_weight>
+  <Pressure_special_weight_index>0</Pressure_special_weight_index>
+  <Pressure_special_weight>0</Pressure_special_weight>
+  <!-- mse -->
+  <abserrMSE>0.001</abserrMSE>
+  <relerrMSE>0.1</relerrMSE>
+  <force_default_errMSE>0</force_default_errMSE>
+  <n_MSE_special_weight>25</n_MSE_special_weight>
+  <MSE_special_weight_index>0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24</MSE_special_weight_index>
+  <MSE_special_weight>0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 0 0 1 1 1 0 0</MSE_special_weight>
+  <!-- vacth -->
+  <n_points_eqx_contour>13</n_points_eqx_contour>
+  <!--r_eqx_contour>3.88 4.00 5.55 7.36 8.04 8.45 8.66 8.44 7.62 5.80 4.87 4.18 3.88</r_eqx_contour-->
+  <!--z_eqx_contour>-2.53 -3.92 -4.65 -2.45 -1.42 -0.47 0.66 1.65 3.21 4.65 4.85 4.39 3.58</z_eqx_contour-->
+  <r_eqx_contour>3.88 4.00 5.55 7.36 8.0 8.45 8.66 8.44 7.62 5.80 4.87 4.18 3.88</r_eqx_contour>
+  <z_eqx_contour>-2.53 -3.92 -4.65 -2.45 -1.4 -0.47 0.66 1.65 3.21 4.65 4.85 4.39 3.58</z_eqx_contour>
+  <max_segment_length_eqx_contour>0.05</max_segment_length_eqx_contour>
+  <vacth_contour_choice>0</vacth_contour_choice>
+  <deltar>0.15</deltar>
+  <deltaz>0.15</deltaz>
+  <d_filament>0.05</d_filament>
+  <use_pfcoils>1</use_pfcoils>
+  <use_same_order>0</use_same_order>
+  <nqc>8</nqc>
+  <qc>0 1 2 3 4 5 6 7</qc>
+  <nqs>7</nqs>
+  <qs>1 2 3 4 5 6 7</qs>
+  <npc>8</npc>
+  <pc>0 1 2 3 4 5 6 7</pc>
+  <nps>7</nps>
+  <ps>1 2 3 4 5 6 7</ps>
+  <manip>1</manip>
+  <maxdist>0.05</maxdist>
+  <vacth_regularization>3</vacth_regularization>
+  <vacth_fracofdeltar_regularization>0.1</vacth_fracofdeltar_regularization>
+  <vacth_eps_regularization>1.0e-6</vacth_eps_regularization>
+  <fep_eps_regularization>1.0e-3</fep_eps_regularization>
+  <compute_plasma_boundary>1</compute_plasma_boundary>
+  <n_plasma_boundary_points>60</n_plasma_boundary_points>
+  <output_vacth_files>1</output_vacth_files>
+  <!-- static inverse -->
+  <parametric_bnd>0</parametric_bnd>
+  <a>1.9</a>
+  <center_r>6.2</center_r>
+  <center_z>0.545</center_z>
+  <kappa>1.8</kappa>
+  <delta>0.43</delta>
+  <rx>5.089</rx>
+  <zx>-3.346</zx>
+  <n_desired_bnd_points>96</n_desired_bnd_points>
+  <use_desired_psib>0</use_desired_psib>
+  <!-- InverseSQP2  -->
+  <abserrIg>3e7</abserrIg>
+  <relerrIg>0</relerrIg>
+  <n_coil_group_index>14</n_coil_group_index>
+  <coil_group_index>0 1 2 2 3 4 5 6 7 8 9 10 11 12</coil_group_index>
+  <n_group_current_index>13</n_group_current_index>
+  <group_current_index>0 1 2 4 5 6 7 8 9 10 11 12 13</group_current_index>
+  <n_group_fixed_index>2</n_group_fixed_index>
+  <group_fixed_index>11 12</group_fixed_index>
+  <n_group_penalized_to_zero_index>0</n_group_penalized_to_zero_index>
+  <group_penalized_to_zero_index>0 1 2 3 4 5 6 7 8 9 10</group_penalized_to_zero_index>
+  <!-- evolution inverse -->
+  <!-- same basis representation for all voltages  -->
+  <Vfunc>2</Vfunc>
+  <Vknots_n>6</Vknots_n>
+  <Vknots>0 0.2 0.4 0.6 0.8 1</Vknots>
+  <ws>1.0</ws>
+  <wR1>1.0e-3</wR1>
+  <wPVmax>0.0e12</wPVmax>
+  <wPVmin>0.0e12</wPVmin>
+  <wPinduced>1.0e-3</wPinduced>
+  <!-- for Clough Tocher  -->
+  <ct_profile_method>0</ct_profile_method>
+  <ct_approx_n_by_triangle>5</ct_approx_n_by_triangle>
+  <!-- control outputs, interpolated magaxis, xpoint, gaps  -->
+  <outputForControl>0</outputForControl>
+  <rd_basis_choice>2</rd_basis_choice>
+  <rd_n_dof>10</rd_n_dof>
+  <rd_basis_choice_petajnib>1</rd_basis_choice_petajnib>
+  <rd_n_dof_petajnib>20</rd_n_dof_petajnib>
+  <dt>1.0e-3</dt>
+  <!-- inverse (psi, theta) grid  -->
+  <n_theta_grid>101</n_theta_grid>
+</parameters>
diff --git a/waveform_editor/shape_editor/xml_param/inverse_param.xml b/waveform_editor/shape_editor/xml_param/inverse_param.xml
new file mode 100644
index 00000000..864f91a5
--- /dev/null
+++ b/waveform_editor/shape_editor/xml_param/inverse_param.xml
@@ -0,0 +1,267 @@
+<?xml version="1.0"?>
+<!-- Taken from: https://gitlab.inria.fr/blfauger/nice
+     Path: nice/run/iwrap/param/inv/iter/param.xml
+     Author: Blaise Faugeras
+     License: LGPL -->
+<parameters>
+  <verbose>1</verbose>
+  <!-- overwritten by shape editor -->
+  <outputTxtFiles>0</outputTxtFiles>
+  <outputRestartFiles>0</outputRestartFiles>
+  <algoMode>31</algoMode>
+  <algoMeshDomain>2</algoMeshDomain>
+  <algoVacTHLRSolver>0</algoVacTHLRSolver>
+  <algoVacTHonly>0</algoVacTHonly>
+  <algoDirect>1</algoDirect>
+  <algoIpFreeScaleInitLambda>0</algoIpFreeScaleInitLambda>
+  <algoRecon>2</algoRecon>
+  <algoReconUseInputIp>1</algoReconUseInputIp>
+  <algoReconInterf>0</algoReconInterf>
+  <algoReconPolar>0</algoReconPolar>
+  <algoReconPressure>0</algoReconPressure>
+  <algoReconMSE>0</algoReconMSE>
+  <algoStartFromScratch>1</algoStartFromScratch>
+  <algoStartFromScratchReconAB>1</algoStartFromScratchReconAB>
+  <algoStartPsiFromInData>0</algoStartPsiFromInData>
+  <algoWithIp>1</algoWithIp>
+  <iterMaxDir>40</iterMaxDir>
+  <epsStopDir>1.0e-10</epsStopDir>
+  <iterMaxDirInitRecon>2</iterMaxDirInitRecon>
+  <iterMaxDirInitInv>0</iterMaxDirInitInv>
+  <iterMaxRecon>20</iterMaxRecon>
+  <epsStopRecon>1.0e-10</epsStopRecon>
+  <iterMaxInv>40</iterMaxInv>
+  <epsStopInv>1.0e-10</epsStopInv>
+  <iterMaxEvo>50</iterMaxEvo>
+  <epsStopEvo>1.0e-10</epsStopEvo>
+  <useDefaultDirectInit>1</useDefaultDirectInit>
+  <radiusInit>2.25</radiusInit>
+  <rcenterInit>6.2</rcenterInit>
+  <zcenterInit>0.0</zcenterInit>
+  <useMobileWall>0</useMobileWall>
+  <refinedZone1>1</refinedZone1>
+  <nPointsZone1>4</nPointsZone1>
+  <rPointsZone1>4.45 5.39 5.78 4.70</rPointsZone1>
+  <zPointsZone1>-3.5 -3.86 -3.20 -2.81</zPointsZone1>
+  <rInPointZone1>5.1</rInPointZone1>
+  <zInPointZone1>-3.3</zInPointZone1>
+  <zone1TriangleArea>0.0005</zone1TriangleArea>
+  <refinedZone2>0</refinedZone2>
+  <nPointsZone2>4</nPointsZone2>
+  <rPointsZone2>2.44 1.91 1.94 2.44</rPointsZone2>
+  <zPointsZone2>0.79 0.573 0.49 0.66</zPointsZone2>
+  <rInPointZone2>2.18</rInPointZone2>
+  <zInPointZone2>0.6</zInPointZone2>
+  <zone2TriangleArea>0.00003</zone2TriangleArea>
+  <refinedZone3>1</refinedZone3>
+  <nPointsZone3>10</nPointsZone3>
+  <rPointsZone3>6.0 6.2 6.2 6.4 6.6 6.6 6.4 6.4 6.2 6.0</rPointsZone3>
+  <zPointsZone3>0.2 0.2 0.4 0.2 0.2 1.0 1.0 0.8 1.0 1.0</zPointsZone3>
+  <rInPointZone3>6.1</rInPointZone3>
+  <zInPointZone3>0.3</zInPointZone3>
+  <zone3TriangleArea>0.0005</zone3TriangleArea>
+  <vacuumTriangleArea>0.002</vacuumTriangleArea>
+  <airTriangleArea>0.04</airTriangleArea>
+  <airInTriangleArea>0.01</airInTriangleArea>
+  <ironTriangleArea>0.005</ironTriangleArea>
+  <coilTriangleArea>0.005</coilTriangleArea>
+  <vesselTriangleArea>0.005</vesselTriangleArea>
+  <passiveStructTriangleArea>0.005</passiveStructTriangleArea>
+  <useNewCOCOSManager>1</useNewCOCOSManager>
+  <inoutCOCOS>17</inoutCOCOS>
+  <change_coil_sign>1</change_coil_sign>
+  <change_psimeas_sign>0</change_psimeas_sign>
+  <change_bmeas_sign>0</change_bmeas_sign>
+  <change_BpSensorAngle_sign>0</change_BpSensorAngle_sign>
+  <change_polarmeas_sign>0</change_polarmeas_sign>
+  <change_dirichlet_sign>0</change_dirichlet_sign>
+  <change_dpressure_dpsi_sign>0</change_dpressure_dpsi_sign>
+  <change_f_df_dpsi_sign>0</change_f_df_dpsi_sign>
+  <change_TPI>0</change_TPI>
+  <change_psi_sign>0</change_psi_sign>
+  <inCOCOS>17</inCOCOS>
+  <outCOCOS>17</outCOCOS>
+  <profile_method>1</profile_method>
+  <n_profile_1d>50</n_profile_1d>
+  <allErrorBar>0</allErrorBar>
+  <r0>6.2</r0>
+  <!-- plasma current density  -->
+  <Afunc>0</Afunc>
+  <Aknots_n>11</Aknots_n>
+  <Aknots>0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1</Aknots>
+  <Aabg>2 0.6 1.4</Aabg>
+  <Bfunc>0</Bfunc>
+  <Bknots_n>11</Bknots_n>
+  <Bknots>0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1</Bknots>
+  <Babg>2 0.4 1.4</Babg>
+  <useInputAB>1</useInputAB>
+  <!-- electron density  -->
+  <Cfunc>0</Cfunc>
+  <Cknots_n>11</Cknots_n>
+  <Cknots>0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1</Cknots>
+  <lambdaNe>1e20</lambdaNe>
+  <Cparam_n>4</Cparam_n>
+  <Cparam>3 1 0 1</Cparam>
+  <!-- electron temperature  -->
+  <Dfunc>0</Dfunc>
+  <Dknots_n>11</Dknots_n>
+  <Dknots>0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1</Dknots>
+  <lambdaTe>25</lambdaTe>
+  <agmM_Te>1 1.5 0 1</agmM_Te>
+  <!--   -->
+  <epsA>1.0e-1</epsA>
+  <sA>1 1 1 1 1 1 1 1 1 1 1</sA>
+  <penalA>1.0e12</penalA>
+  <penalNegativeA>0.0e10</penalNegativeA>
+  <epsB>1.0e-1</epsB>
+  <sB>1 1 1 1 1 1 1 1 1 1 1</sB>
+  <penalB>1.0e12</penalB>
+  <epsC>1.0e-1</epsC>
+  <sC>1 1 1 1 1 1 1 1 1 1 1</sC>
+  <penalC>1.0e12</penalC>
+  <epsCparam>1.0e3 1.0e3 1.0e0 1.0e12 1.0e0</epsCparam>
+  <epsD>1.0e-1</epsD>
+  <sD>1 1 1 1 1 1 1 1 1 1 1 1</sD>
+  <penalD>0.0e20</penalD>
+  <penalD0>0.0e0</penalD0>
+  <!-- magnetics -->
+  <neumannPercent>0.1</neumannPercent>
+  <abserrB>1e-3</abserrB>
+  <relerrB>0.01</relerrB>
+  <force_default_errB>0</force_default_errB>
+  <n_B_special_weight>0</n_B_special_weight>
+  <B_special_weight_index>0</B_special_weight_index>
+  <B_special_weight>0</B_special_weight>
+  <abserrF>0.01</abserrF>
+  <relerrF>0.01</relerrF>
+  <force_default_errF>0</force_default_errF>
+  <n_F_special_weight>0</n_F_special_weight>
+  <F_special_weight_index>0</F_special_weight_index>
+  <F_special_weight>0</F_special_weight>
+  <add_fictitious_full_loop_if_not_exist>0</add_fictitious_full_loop_if_not_exist>
+  <!-- interferometry -->
+  <chord_dl>0.01</chord_dl>
+  <constructLOS>0</constructLOS>
+  <abserrNe>1e18</abserrNe>
+  <relerrNe>0.1</relerrNe>
+  <force_default_errNe>0</force_default_errNe>
+  <n_Ne_special_weight>0</n_Ne_special_weight>
+  <Ne_special_weight_index>0 1 2 3 4 5 6 7 8 9 10 11 12 13 14</Ne_special_weight_index>
+  <!--Ne_special_weight>1 1 1 1 1 1 1 1 1 1 1 1 1 1 1</Ne_special_weight-->
+  <Ne_special_weight>0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1</Ne_special_weight>
+  <!--Ne_special_weight>0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01</Ne_special_weight-->
+  <!-- polarimetry -->
+  <wavelength>0.195</wavelength>
+  <C1>2.42e-20</C1>
+  <C3>5.23e-19</C3>
+  <abserrPolar>0.001</abserrPolar>
+  <relerrPolar>0.1</relerrPolar>
+  <force_default_errPolar>0</force_default_errPolar>
+  <n_Polar_special_weight>0</n_Polar_special_weight>
+  <Polar_special_weight_index>0 1 2 3 4 5 6 7 8 9 10 11 12 13 14</Polar_special_weight_index>
+  <!--Polar_special_weight>1 1 1 1 1 1 1 1 1 1 1 1 1 1 1</Polar_special_weight-->
+  <Polar_special_weight>0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1</Polar_special_weight>
+  <abserrS>0.001</abserrS>
+  <relerrS>0.1</relerrS>
+  <force_default_errS>0</force_default_errS>
+  <n_S_special_weight>0</n_S_special_weight>
+  <S_special_weight_index>0 1 2 3 4 5 6 7 8 9 10 11 12 13 14</S_special_weight_index>
+  <!--S_special_weight>0</S_special_weight-->
+  <!--S_special_weight>0.1 0 0 0 0 0 0 0 0 0 0 0 0 0 0</S_special_weight-->
+  <S_special_weight>0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1</S_special_weight>
+  <!--S_special_weight>0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01</S_special_weight-->
+  <!--S_special_weight>0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001</S_special_weight-->
+  <S012_weight>1 1 1</S012_weight>
+  <!-- pressure -->
+  <abserrPressure>1.0e3</abserrPressure>
+  <relerrPressure>0.1</relerrPressure>
+  <force_default_errPressure>0</force_default_errPressure>
+  <n_Pressure_special_weight>0</n_Pressure_special_weight>
+  <Pressure_special_weight_index>0</Pressure_special_weight_index>
+  <Pressure_special_weight>0</Pressure_special_weight>
+  <!-- mse -->
+  <abserrMSE>0.001</abserrMSE>
+  <relerrMSE>0.1</relerrMSE>
+  <force_default_errMSE>0</force_default_errMSE>
+  <n_MSE_special_weight>25</n_MSE_special_weight>
+  <MSE_special_weight_index>0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24</MSE_special_weight_index>
+  <MSE_special_weight>0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 0 0 1 1 1 0 0</MSE_special_weight>
+  <!-- vacth -->
+  <n_points_eqx_contour>13</n_points_eqx_contour>
+  <!--r_eqx_contour>3.88 4.00 5.55 7.36 8.04 8.45 8.66 8.44 7.62 5.80 4.87 4.18 3.88</r_eqx_contour-->
+  <!--z_eqx_contour>-2.53 -3.92 -4.65 -2.45 -1.42 -0.47 0.66 1.65 3.21 4.65 4.85 4.39 3.58</z_eqx_contour-->
+  <r_eqx_contour>3.88 4.00 5.55 7.36 8.0 8.45 8.66 8.44 7.62 5.80 4.87 4.18 3.88</r_eqx_contour>
+  <z_eqx_contour>-2.53 -3.92 -4.65 -2.45 -1.4 -0.47 0.66 1.65 3.21 4.65 4.85 4.39 3.58</z_eqx_contour>
+  <max_segment_length_eqx_contour>0.05</max_segment_length_eqx_contour>
+  <vacth_contour_choice>0</vacth_contour_choice>
+  <deltar>0.15</deltar>
+  <deltaz>0.15</deltaz>
+  <d_filament>0.05</d_filament>
+  <use_pfcoils>1</use_pfcoils>
+  <use_same_order>0</use_same_order>
+  <nqc>8</nqc>
+  <qc>0 1 2 3 4 5 6 7</qc>
+  <nqs>7</nqs>
+  <qs>1 2 3 4 5 6 7</qs>
+  <npc>8</npc>
+  <pc>0 1 2 3 4 5 6 7</pc>
+  <nps>7</nps>
+  <ps>1 2 3 4 5 6 7</ps>
+  <manip>1</manip>
+  <maxdist>0.05</maxdist>
+  <vacth_regularization>3</vacth_regularization>
+  <vacth_fracofdeltar_regularization>0.1</vacth_fracofdeltar_regularization>
+  <vacth_eps_regularization>1.0e-6</vacth_eps_regularization>
+  <fep_eps_regularization>1.0e-3</fep_eps_regularization>
+  <compute_plasma_boundary>1</compute_plasma_boundary>
+  <n_plasma_boundary_points>60</n_plasma_boundary_points>
+  <output_vacth_files>1</output_vacth_files>
+  <!-- static inverse -->
+  <parametric_bnd>0</parametric_bnd>
+  <a>1.9</a>
+  <center_r>6.2</center_r>
+  <center_z>0.545</center_z>
+  <kappa>1.8</kappa>
+  <delta>0.43</delta>
+  <rx>5.089</rx>
+  <zx>-3.346</zx>
+  <n_desired_bnd_points>96</n_desired_bnd_points>
+  <use_desired_psib>0</use_desired_psib>
+  <!-- InverseSQP2  -->
+  <!-- Inverse old with turns  -->
+  <!--abserrIg>3e7</abserrIg-->
+  <!-- Inverse new  -->
+  <abserrIg>3e5</abserrIg>
+  <relerrIg>0</relerrIg>
+  <n_coil_group_index>14</n_coil_group_index>
+  <coil_group_index>0 1 2 2 3 4 5 6 7 8 9 10 11 12</coil_group_index>
+  <n_group_current_index>13</n_group_current_index>
+  <group_current_index>0 1 2 4 5 6 7 8 9 10 11 12 13</group_current_index>
+  <!-- Fixed coil currents are updated by the shape editor -->
+  <n_group_fixed_index>2</n_group_fixed_index>
+  <group_fixed_index>11 12</group_fixed_index>
+  <n_group_penalized_to_zero_index>0</n_group_penalized_to_zero_index>
+  <group_penalized_to_zero_index>0 1 2 3 4 5 6 7 8 9 10</group_penalized_to_zero_index>
+  <!-- evolution inverse -->
+  <!-- same basis representation for all voltages  -->
+  <Vfunc>2</Vfunc>
+  <Vknots_n>6</Vknots_n>
+  <Vknots>0 0.2 0.4 0.6 0.8 1</Vknots>
+  <ws>1.0</ws>
+  <wR1>1.0e-3</wR1>
+  <wPVmax>0.0e12</wPVmax>
+  <wPVmin>0.0e12</wPVmin>
+  <wPinduced>1.0e-3</wPinduced>
+  <!-- for Clough Tocher  -->
+  <ct_profile_method>0</ct_profile_method>
+  <ct_approx_n_by_triangle>5</ct_approx_n_by_triangle>
+  <!-- control outputs, interpolated magaxis, xpoint, gaps  -->
+  <outputForControl>0</outputForControl>
+  <rd_basis_choice>2</rd_basis_choice>
+  <rd_n_dof>10</rd_n_dof>
+  <rd_basis_choice_petajnib>1</rd_basis_choice_petajnib>
+  <rd_n_dof_petajnib>20</rd_n_dof_petajnib>
+  <dt>1.0e-3</dt>
+  <!-- inverse (psi, theta) grid  -->
+  <n_theta_grid>101</n_theta_grid>
+</parameters>
diff --git a/waveform_editor/shape_editor/xml_param/param.xml b/waveform_editor/shape_editor/xml_param/param.xml
deleted file mode 100644
index d8a25eb8..00000000
--- a/waveform_editor/shape_editor/xml_param/param.xml
+++ /dev/null
@@ -1,323 +0,0 @@
-<?xml version="1.0"?>
-<!-- Taken from: https://gitlab.inria.fr/blfauger/nice
-     Path: nice/run/iwrap/param/inv/iter/param.xml
-     Author: Blaise Faugeras
-     License: LGPL -->
-<parameters>
-    <verbose>1</verbose>  <!-- overwritten by shape editor -->
-    <outputTxtFiles>0</outputTxtFiles>
-    <outputRestartFiles>0</outputRestartFiles>
-    <algoMode>31</algoMode>
-    <algoMeshDomain>2</algoMeshDomain>
-    <algoVacTHLRSolver>0</algoVacTHLRSolver>
-    <algoVacTHonly>0</algoVacTHonly>
-    <algoDirect>1</algoDirect>
-    <algoIpFreeScaleInitLambda>0</algoIpFreeScaleInitLambda>
-    <algoRecon>2</algoRecon>
-    <algoReconUseInputIp>1</algoReconUseInputIp>
-    <algoReconInterf>0</algoReconInterf>
-    <algoReconPolar>0</algoReconPolar>
-    <algoReconPressure>0</algoReconPressure>
-    <algoReconMSE>0</algoReconMSE>
-    <algoStartFromScratch>1</algoStartFromScratch>
-    <algoStartFromScratchReconAB>1</algoStartFromScratchReconAB>
-    <algoStartPsiFromInData>0</algoStartPsiFromInData>
-    <algoWithIp>1</algoWithIp>
-    <iterMaxDir>40</iterMaxDir>
-    <epsStopDir>1.0e-10</epsStopDir>
-    <iterMaxDirInitRecon>2</iterMaxDirInitRecon>
-    <iterMaxDirInitInv>0</iterMaxDirInitInv>
-    <iterMaxRecon>20</iterMaxRecon>
-    <epsStopRecon>1.0e-10</epsStopRecon>
-    <iterMaxInv>40</iterMaxInv>
-    <epsStopInv>1.0e-10</epsStopInv>
-    <iterMaxEvo>50</iterMaxEvo>
-    <epsStopEvo>1.0e-10</epsStopEvo>
-    <useDefaultDirectInit>1</useDefaultDirectInit>
-    <radiusInit>2.25</radiusInit>
-    <rcenterInit>6.2</rcenterInit>
-    <zcenterInit>0.0</zcenterInit>
-
-    <useMobileWall>0</useMobileWall>
-
-    <refinedZone1>1</refinedZone1>
-    <nPointsZone1>4</nPointsZone1>
-    <rPointsZone1>4.45 5.39 5.78 4.70</rPointsZone1>
-    <zPointsZone1>-3.5 -3.86 -3.20 -2.81</zPointsZone1>
-    <rInPointZone1>5.1</rInPointZone1>
-    <zInPointZone1>-3.3</zInPointZone1>
-    <zone1TriangleArea>0.0005</zone1TriangleArea>
-    
-    <refinedZone2>0</refinedZone2>
-    <nPointsZone2>4</nPointsZone2>
-    <rPointsZone2>2.44 1.91 1.94 2.44</rPointsZone2>
-    <zPointsZone2>0.79 0.573 0.49 0.66</zPointsZone2>
-    <rInPointZone2>2.18</rInPointZone2>
-    <zInPointZone2>0.6</zInPointZone2>
-    <zone2TriangleArea>0.00003</zone2TriangleArea>
-    
-    <refinedZone3>1</refinedZone3>
-    <nPointsZone3>10</nPointsZone3>
-    <rPointsZone3>6.0 6.2 6.2 6.4 6.6 6.6 6.4 6.4 6.2 6.0</rPointsZone3>
-    <zPointsZone3>0.2 0.2 0.4 0.2 0.2 1.0 1.0 0.8 1.0 1.0</zPointsZone3>
-    <rInPointZone3>6.1</rInPointZone3>
-    <zInPointZone3>0.3</zInPointZone3>
-    <zone3TriangleArea>0.0005</zone3TriangleArea>
-
-    <vacuumTriangleArea>0.002</vacuumTriangleArea>
-    <airTriangleArea>0.04</airTriangleArea>
-    <airInTriangleArea>0.01</airInTriangleArea>
-    <ironTriangleArea>0.005</ironTriangleArea>
-    <coilTriangleArea>0.005</coilTriangleArea>
-    <vesselTriangleArea>0.005</vesselTriangleArea>
-    <passiveStructTriangleArea>0.005</passiveStructTriangleArea>
-
-    <useNewCOCOSManager>1</useNewCOCOSManager>
-    <inoutCOCOS>17</inoutCOCOS>
-
-    <change_coil_sign>1</change_coil_sign>
-    <change_psimeas_sign>0</change_psimeas_sign>
-    <change_bmeas_sign>0</change_bmeas_sign>
-    <change_BpSensorAngle_sign>0</change_BpSensorAngle_sign>
-    <change_polarmeas_sign>0</change_polarmeas_sign>
-    <change_dirichlet_sign>0</change_dirichlet_sign>
-    <change_dpressure_dpsi_sign>0</change_dpressure_dpsi_sign>
-    <change_f_df_dpsi_sign>0</change_f_df_dpsi_sign>
-    <change_TPI>0</change_TPI>
-    <change_psi_sign>0</change_psi_sign>
-
-    <inCOCOS>17</inCOCOS>
-    <outCOCOS>17</outCOCOS>
-    
-    <profile_method>1</profile_method>
-
-    <n_profile_1d>50</n_profile_1d>
-    <allErrorBar>0</allErrorBar>
-
-    <r0>6.2</r0>
-
-<!-- plasma current density  -->
-    <Afunc>0</Afunc>
-    <Aknots_n>11</Aknots_n>
-    <Aknots>0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1</Aknots>
-    <Aabg>2 0.6 1.4</Aabg>
-
-    <Bfunc>0</Bfunc>
-    <Bknots_n>11</Bknots_n>
-    <Bknots>0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1</Bknots>
-    <Babg>2 0.4 1.4</Babg>
-
-    <useInputAB>1</useInputAB>
-
-<!-- electron density  -->
-    <Cfunc>0</Cfunc>
-    <Cknots_n>11</Cknots_n>
-    <Cknots>0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1</Cknots>
-    <lambdaNe>1e20</lambdaNe>
-    <Cparam_n>4</Cparam_n>
-    <Cparam>3 1 0 1</Cparam>
-
-<!-- electron temperature  -->
-    <Dfunc>0</Dfunc>
-    <Dknots_n>11</Dknots_n>
-    <Dknots>0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1</Dknots>
-    <lambdaTe>25</lambdaTe>
-    <agmM_Te>1 1.5 0 1</agmM_Te>
-
-
-<!--   -->
-    <epsA>1.0e-1</epsA>
-    <sA>1 1 1 1 1 1 1 1 1 1 1</sA>
-    <penalA>1.0e12</penalA>
-    <penalNegativeA>0.0e10</penalNegativeA>
-
-    <epsB>1.0e-1</epsB>
-    <sB>1 1 1 1 1 1 1 1 1 1 1</sB>
-    <penalB>1.0e12</penalB>
-
-    <epsC>1.0e-1</epsC>
-    <sC>1 1 1 1 1 1 1 1 1 1 1</sC>
-    <penalC>1.0e12</penalC>
-    <epsCparam>1.0e3 1.0e3 1.0e0 1.0e12 1.0e0</epsCparam>
-
-    <epsD>1.0e-1</epsD>
-    <sD>1 1 1 1 1 1 1 1 1 1 1 1</sD>
-    <penalD>0.0e20</penalD>
-    <penalD0>0.0e0</penalD0>
-
-<!-- magnetics -->
-    <neumannPercent>0.1</neumannPercent>
-
-    <abserrB>1e-3</abserrB>
-    <relerrB>0.01</relerrB>
-    <force_default_errB>0</force_default_errB>
-    <n_B_special_weight>0</n_B_special_weight>
-    <B_special_weight_index>0</B_special_weight_index>
-    <B_special_weight>0</B_special_weight>
-
-    <abserrF>0.01</abserrF>
-    <relerrF>0.01</relerrF>
-    <force_default_errF>0</force_default_errF>
-    <n_F_special_weight>0</n_F_special_weight>
-    <F_special_weight_index>0</F_special_weight_index>
-    <F_special_weight>0</F_special_weight>
-    <add_fictitious_full_loop_if_not_exist>0</add_fictitious_full_loop_if_not_exist>
-
-
-<!-- interferometry -->
-    <chord_dl>0.01</chord_dl>
-    <constructLOS>0</constructLOS>
-
-    <abserrNe>1e18</abserrNe>
-    <relerrNe>0.1</relerrNe>
-    <force_default_errNe>0</force_default_errNe>
-    <n_Ne_special_weight>0</n_Ne_special_weight>
-    <Ne_special_weight_index>0 1 2 3 4 5 6 7 8 9 10 11 12 13 14</Ne_special_weight_index>
-    <!--Ne_special_weight>1 1 1 1 1 1 1 1 1 1 1 1 1 1 1</Ne_special_weight-->
-    <Ne_special_weight>0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1</Ne_special_weight>
-    <!--Ne_special_weight>0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01</Ne_special_weight-->
-
-<!-- polarimetry -->
-    <wavelength>0.195</wavelength>
-    <C1>2.42e-20</C1>
-    <C3>5.23e-19</C3>
-
-    <abserrPolar>0.001</abserrPolar>
-    <relerrPolar>0.1</relerrPolar>
-    <force_default_errPolar>0</force_default_errPolar>
-    <n_Polar_special_weight>0</n_Polar_special_weight>
-    <Polar_special_weight_index>0 1 2 3 4 5 6 7 8 9 10 11 12 13 14</Polar_special_weight_index>
-    <!--Polar_special_weight>1 1 1 1 1 1 1 1 1 1 1 1 1 1 1</Polar_special_weight-->
-    <Polar_special_weight>0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1</Polar_special_weight>
-
-    <abserrS>0.001</abserrS>
-    <relerrS>0.1</relerrS>
-    <force_default_errS>0</force_default_errS>
-    <n_S_special_weight>0</n_S_special_weight>
-    <S_special_weight_index>0 1 2 3 4 5 6 7 8 9 10 11 12 13 14</S_special_weight_index>
-    <!--S_special_weight>0</S_special_weight-->
-    <!--S_special_weight>0.1 0 0 0 0 0 0 0 0 0 0 0 0 0 0</S_special_weight-->
-    <S_special_weight>0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1</S_special_weight>
-    <!--S_special_weight>0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01</S_special_weight-->
-    <!--S_special_weight>0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001</S_special_weight-->
-    <S012_weight>1 1 1</S012_weight>
-
-<!-- pressure -->
-    <abserrPressure>1.0e3</abserrPressure>
-    <relerrPressure>0.1</relerrPressure>
-    <force_default_errPressure>0</force_default_errPressure>
-    <n_Pressure_special_weight>0</n_Pressure_special_weight>
-    <Pressure_special_weight_index>0</Pressure_special_weight_index>
-    <Pressure_special_weight>0</Pressure_special_weight>
-
-<!-- mse -->
-    <abserrMSE>0.001</abserrMSE>
-    <relerrMSE>0.1</relerrMSE>
-    <force_default_errMSE>0</force_default_errMSE>
-    <n_MSE_special_weight>25</n_MSE_special_weight>
-    <MSE_special_weight_index>0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24</MSE_special_weight_index>
-    <MSE_special_weight>0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 0 0 1 1 1 0 0</MSE_special_weight>
-
-
-
-<!-- vacth --> 
-    <n_points_eqx_contour>13</n_points_eqx_contour>
-    <!--r_eqx_contour>3.88 4.00 5.55 7.36 8.04 8.45 8.66 8.44 7.62 5.80 4.87 4.18 3.88</r_eqx_contour-->
-    <!--z_eqx_contour>-2.53 -3.92 -4.65 -2.45 -1.42 -0.47 0.66 1.65 3.21 4.65 4.85 4.39 3.58</z_eqx_contour-->
-    <r_eqx_contour>3.88 4.00 5.55 7.36 8.0 8.45 8.66 8.44 7.62 5.80 4.87 4.18 3.88</r_eqx_contour>
-    <z_eqx_contour>-2.53 -3.92 -4.65 -2.45 -1.4 -0.47 0.66 1.65 3.21 4.65 4.85 4.39 3.58</z_eqx_contour>
-
-    <max_segment_length_eqx_contour>0.05</max_segment_length_eqx_contour>
-    <vacth_contour_choice>0</vacth_contour_choice>
-
-    <deltar>0.15</deltar>
-    <deltaz>0.15</deltaz>
-    <d_filament>0.05</d_filament>
-
-    <use_pfcoils>1</use_pfcoils>
-
-    <use_same_order>0</use_same_order>
-    <nqc>8</nqc>      
-    <qc>0 1 2 3 4 5 6 7</qc> 
-    <nqs>7</nqs>      
-    <qs>1 2 3 4 5 6 7</qs>     
-    <npc>8</npc>
-    <pc>0 1 2 3 4 5 6 7</pc>
-    <nps>7</nps>
-    <ps>1 2 3 4 5 6 7</ps> 
-    <manip>1</manip>      
-    <maxdist>0.05</maxdist>    
-         
-    <vacth_regularization>3</vacth_regularization>     
-    <vacth_fracofdeltar_regularization>0.1</vacth_fracofdeltar_regularization>   
-    <vacth_eps_regularization>1.0e-6</vacth_eps_regularization>   
-    <fep_eps_regularization>1.0e-3</fep_eps_regularization>   
-
-    <compute_plasma_boundary>1</compute_plasma_boundary>      
-    <n_plasma_boundary_points>60</n_plasma_boundary_points>      
-    <output_vacth_files>1</output_vacth_files>      
-
-<!-- static inverse -->
-    <parametric_bnd>0</parametric_bnd>
-    <a>1.9</a>
-    <center_r>6.2</center_r>
-    <center_z>0.545</center_z>
-    <kappa>1.8</kappa>
-    <delta>0.43</delta>
-    <rx>5.089</rx>
-    <zx>-3.346</zx>
-    <n_desired_bnd_points>96</n_desired_bnd_points>
-
-    <use_desired_psib>0</use_desired_psib>
-    
-    <!-- InverseSQP2  -->
-    <!-- Inverse old with turns  -->
-    <!--abserrIg>3e7</abserrIg-->
-    <!-- Inverse new  -->
-    <abserrIg>3e5</abserrIg>
-    <relerrIg>0</relerrIg>
-    <n_coil_group_index>14</n_coil_group_index>
-    <coil_group_index>0 1 2 2 3 4 5 6 7 8 9 10 11 12</coil_group_index>
-    <n_group_current_index>13</n_group_current_index>
-    <group_current_index>0 1 2 4 5 6 7 8 9 10 11 12 13</group_current_index>
-
-    <!-- Fixed coil currents are updated by the shape editor -->
-    <n_group_fixed_index>2</n_group_fixed_index>
-    <group_fixed_index>11 12</group_fixed_index>
-
-    <n_group_penalized_to_zero_index>0</n_group_penalized_to_zero_index>
-    <group_penalized_to_zero_index>0 1 2 3 4 5 6 7 8 9 10</group_penalized_to_zero_index>
-
-
-<!-- evolution inverse -->
-    
-<!-- same basis representation for all voltages  -->
-    <Vfunc>2</Vfunc>
-    <Vknots_n>6</Vknots_n>
-    <Vknots>0 0.2 0.4 0.6 0.8 1</Vknots>
-    
-    <ws>1.0</ws>
-    <wR1>1.0e-3</wR1>
-    <wPVmax>0.0e12</wPVmax>
-    <wPVmin>0.0e12</wPVmin>
-    <wPinduced>1.0e-3</wPinduced>
-
-<!-- for Clough Tocher  -->
-    <ct_profile_method>0</ct_profile_method>
-    <ct_approx_n_by_triangle>5</ct_approx_n_by_triangle>
-
-<!-- control outputs, interpolated magaxis, xpoint, gaps  -->
-    <outputForControl>0</outputForControl>
-
-   <rd_basis_choice>2</rd_basis_choice>
-   <rd_n_dof>10</rd_n_dof> 
-    <rd_basis_choice_petajnib>1</rd_basis_choice_petajnib>
-    <rd_n_dof_petajnib>20</rd_n_dof_petajnib>
-    
-    <dt>1.0e-3</dt>
-
-<!-- inverse (psi, theta) grid  -->
-    <n_theta_grid>101</n_theta_grid>
-
-</parameters>
-