Ensuring that an independent doublet processes grids the same as the Java-code run from utc

import numpy as np

import xarray as xr

def calculate_performance(doublet, input_params: dict, hydrocarbons: float, depth: float, thickness: float, porosity: float, ntg: float, temperature: float, transmissivity: float, transmissivity_with_ntg: float) -> dict:

    if hydrocarbons == 0.0:

def calculate_performance(hydrocarbons: float, depth: float, thickness: float, porosity: float, ntg: float, temperature: float, transmissivity: float, transmissivity_with_ntg: float, doublet=None,

                          input_params: dict = None, return_dict: bool = False):

    output_values ={"power": 0.0,

                "heat_pump_power": 0.0,

                "capex": 0.0,

                "flow_rate": 0.0,

                "welld": 0.0,

                }

        return output_values

    if hydrocarbons != 0.0:

        set_doublet_parameters(doublet, transmissivity_with_ntg, depth, porosity, ntg, temperature, input_params["use_stimulation"], input_params["stimKhMax"], input_params["surface_temperature"],

                               input_params["return_temperature"], input_params["use_heat_pump"], input_params["max_cooling_temperature_range"], input_params["hp_minimum_injection_temperature"])

        # Reset doublet variables

        doublet.setProjectVariables(False, 0.0)

    if return_dict:

        return output_values

    else:

        return output_values["power"], output_values["heat_pump_power"], output_values["capex"], output_values["opex"], output_values["utc"], output_values["npv"], output_values["hprod"], output_values["cop"], output_values["cophp"], output_values["pres"], output_values["flow_rate"], output_values["welld"]

def set_doublet_parameters(doublet, transmissivity_with_ntg, depth, porosity, ntg, temperature, useStimulation, stimKhMax, surface_temperature, return_temperature, use_heat_pump, max_cooling_temperature_range,