Legacy #2: non-equilibrium CO2 (Tvib_12, Tvib_3, Trot)

A 3 temperature fitting example . reproducing the validation case of Klarenaar 2017 [1], who calculated a transmittance spectrum from the initial data of Dang 1973 [2], with a 1 rotational temperature + 3 vibrational temperature (Treanor distributions) model.

|CO2| Energies are calculated from Dunham developments in an uncoupled harmonic oscillator - rigid rotor model. The example is based on one of RADIS validation cases. It makes use of the RADIS Spectrum class and the associated compare and load functions

Typical output is similar to the radis-examples Multi-temperature fit :

The method requires a fitting model. An example model is provided in radis.tools.fitting : Tvib12Tvib3Trot_NonLTEModel(). Other models can be used, such as in the one-temperature fit example

More advanced tools for interactive fitting of multi-dimensional, multi-slabs spectra can be found in fitroom. Finally, the GPU-accelerated example shows how to obtain real-time interactive spectra.

[1]

Klarenaar et al 2017, “Time evolution of vibrational temperatures in a CO2 glow discharge measured with infrared absorption spectroscopy” doi/10.1088/1361-6595/aa902e

[2]

Dang et al 1982, “Detailed vibrational population distributions in a CO2 laser discharge as measured with a tunable diode laser” doi/10.1007/BF00694640

from os.path import join

from radis import Spectrum, SpectrumFactory

Get Fitted Data

from radis.test.utils import getValidationCase, setup_test_line_databases
from radis.tools.fitting import Tvib12Tvib3Trot_NonLTEModel

setup_test_line_databases()
# Data from Dang, adapted by Klarenaar, digitized by us
s_exp = Spectrum.from_txt(
    getValidationCase(
        join("test_CO2_3Tvib_vs_klarenaar_data", "klarenaar_2017_digitized_data.csv")
    ),
    "transmittance_noslit",
    wunit="cm-1",
    unit="",
    delimiter=",",
    name="Klarenaar 2017",
)

/home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/master/radis/spectrum/spectrum.py:5425: UserWarning: Wavespace is not evenly spaced (0.000%) for transmittance_noslit. This may create problems if later convolving with slit function (`s.apply_slit()`). You can use `s.resample_even()`
  warn(

Calculate

sf = SpectrumFactory(
    2284.2,
    2284.6,
    wstep=0.001,  # cm-1
    pressure=20 * 1e-3,  # bar
    db_use_cached=True,
    lvl_use_cached=True,
    cutoff=1e-25,
    isotope="1,2",
    path_length=10,  # cm-1
    mole_fraction=0.1 * 28.97 / 44.07,
    truncation=1,  # cm-1
    medium="vacuum",
    export_populations=None,  # 'vib',
    # parsum_mode="tabulation"
)
sf.warnings["MissingSelfBroadeningWarning"] = "ignore"
sf.warnings["PerformanceWarning"] = "ignore"
sf.load_databank("HITEMP-CO2-TEST")

s_best, best = sf.fit_legacy(
    s_exp.take("transmittance_noslit"),
    model=Tvib12Tvib3Trot_NonLTEModel,
    fit_parameters={
        "T12": 517,
        "T3": 2641,
        "Trot": 491,
    },
    bounds={"T12": [300, 2000], "T3": [300, 5000], "Trot": [300, 2000]},
    fixed_parameters={"vib_distribution": "treanor"},
    plot=True,
    solver_options={
        "method": "TNC",
        "maxiter": 80,  # 👈 increase to let the fit converge
    },
)

• 
• 

/home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/master/examples/3_Fitting/plot5_legacyFit_Trot-Tvib.py:66: DeprecationWarning: 'db_use_cached' removed from SpectrumFactory init and moved in load/fetch_databank()
  sf = SpectrumFactory(
/home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/master/examples/3_Fitting/plot5_legacyFit_Trot-Tvib.py:66: DeprecationWarning: 'lvl_use_cached' removed from SpectrumFactory init and moved in load/fetch_databank()
  sf = SpectrumFactory(
Using database: HITEMP-CO2-TEST
'HITEMP-CO2-TEST':
{'info': 'HITEMP-2010, CO2, 3 main isotope (CO2-626, 636, 628), 2283.7-2285.1 cm-1', 'path': ['/home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/master/radis/test/files/cdsd_hitemp_09_fragment.txt'], 'format': 'cdsd-hitemp', 'levelsfmt': 'radis'}



/home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/master/radis/misc/warning.py:427: MissingReferenceWarning: Missing doi for CDSD-HITEMP. Use HITEMP-2010?
  warnings.warn(WarningType(message))
Calculating Non-Equilibrium Spectrum
Physical Conditions
----------------------------------------
   Tgas                 491 K
   Trot                 491 K
   Tvib                 517,517,2641 K
   isotope              1,2
   mole_fraction        0.06573632856818698
   path_length          10 cm
   pressure             0.02 bar
   rot_distribution     boltzmann
   self_absorption      True
   species              CO2
   state                X
   vib_distribution     treanor
   wavenum_max          2284.6000 cm-1
   wavenum_min          2284.2000 cm-1
Computation Parameters
----------------------------------------
   Tref                 296 K
   add_at_used
   broadening_method    voigt
   cutoff               1e-25 cm-1/(#.cm-2)
   dbformat             cdsd-hitemp
   dbpath               /home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/master/radis/test/files/cdsd_hitemp...
   diluent              air
   folding_thresh       1e-06
   include_neighbouring_lines  True
   isatom               False
   isneutral            None
   lbfunc               None
   memory_mapping_engine  auto
   neighbour_lines      0 cm-1
   optimization         simple
   parsum_mode          full summation
   pfsource             default
   potential_lowering   None
   pseudo_continuum_threshold  0
   sparse_ldm           auto
   truncation           1 cm-1
   waveunit             cm-1
   wstep                0.001 cm-1
   zero_padding         -1
----------------------------------------
Calculating energy levels with Dunham expansion for CO2(X1Σu+)(iso1)

(0s)    0.0%
(0s)    2.6%
(1s)    5.3%
(1s)    7.9%
(1s)    10.5%
(1s)    13.2%
(2s)    15.8%
(2s)    18.4%
(2s)    21.1%
(2s)    23.7%
(2s)    26.3%
(2s)    28.9%
(2s)    31.6%
(2s)    34.2%
(2s)    36.8%
(3s)    39.5%
(3s)    42.1%
(3s)    44.7%
(3s)    47.4%
(3s)    50.0%
(3s)    52.6%
(3s)    55.3%
(3s)    57.9%
(3s)    60.5%
(3s)    63.2%
(3s)    65.8%
(3s)    68.4%
(3s)    71.1%
(3s)    73.7%
(3s)    76.3%
(3s)    78.9%
(3s)    81.6%
(3s)    84.2%
(3s)    86.8%
(3s)    89.5%
(3s)    92.1%
(3s)    94.7%
(3s)    97.4%
(3s)    100.0%
Database generated up to v1=37, v2=65, v3=22, J=344
Calculating energy levels with Dunham expansion for CO2(X1Σu+)(iso2)

(0s)    0.0%
(0s)    2.6%
(1s)    5.3%
(1s)    7.9%
(1s)    10.5%
(1s)    13.2%
(1s)    15.8%
(2s)    18.4%
(2s)    21.1%
(2s)    23.7%
(2s)    26.3%
(2s)    28.9%
(2s)    31.6%
(3s)    34.2%
(3s)    36.8%
(3s)    39.5%
(3s)    42.1%
(3s)    44.7%
(3s)    47.4%
(3s)    50.0%
(3s)    52.6%
(3s)    55.3%
(3s)    57.9%
(3s)    60.5%
(3s)    63.2%
(3s)    65.8%
(3s)    68.4%
(3s)    71.1%
(3s)    73.7%
(3s)    76.3%
(3s)    78.9%
(3s)    81.6%
(3s)    84.2%
(3s)    86.8%
(3s)    89.5%
(3s)    92.1%
(3s)    94.7%
(3s)    97.4%
(3s)    100.0%
Database generated up to v1=37, v2=67, v3=23, J=344
Fetching Evib & Erot. If using this code several times you should consider updating the database directly. See functions in factory.py
/home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/master/radis/misc/warning.py:427: NegativeEnergiesWarning: There are negative rotational energies in the database
  warnings.warn(WarningType(message))
12.29s - Spectrum calculated
/home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/master/radis/misc/warning.py:427: NegativeEnergiesWarning: There are negative rotational energies in the database
  warnings.warn(WarningType(message))
------------------------------
TYPICAL FIT CALCULATION TIME:
Fit (in progress) profiler :
    spectrum_calculation      0.084s ████████████████
        check_line_databank                  0.001s
        check_non_eq_param                   0.000s
        reinitialize                         0.004s
            copy_database                        0.001s
            memory_usage_warning                 0.002s
            reset_population                     0.000s
        calc_noneq_population_multiTvib      0.058s ██████████
            part_function                        0.053s ██████████
            others                               0.005s
        scaled_non_eq_linestrength           0.002s
        scaled_S0                            0.001s
        calc_emission_integral               0.002s
        applied_linestrength_cutoff          0.001s
        calc_hwhm                            0.004s
        calc_lineshift                       0.001s
        generate_wavenumber_arrays           0.000s
        calc_line_broadening                 0.012s ██
            precompute_LDM_lineshapes            0.001s
            LDM_Initialized_vectors              0.000s
            LDM_closest_matching_line            0.000s
            LDM_Distribute_lines                 0.008s █
            LDM_convolve                         0.002s
        calc_other_spectral_quan             0.000s
        generate_spectrum_obj                0.000s
        others                               0.004s
------------------------------
/home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/master/radis/misc/warning.py:427: NegativeEnergiesWarning: There are negative rotational energies in the database
  warnings.warn(WarningType(message))
/home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/master/radis/misc/warning.py:427: NegativeEnergiesWarning: There are negative rotational energies in the database
  warnings.warn(WarningType(message))

Now starting the fitting process:
---------------------------------

/home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/master/radis/tools/fitting.py:480: OptimizeWarning: Unknown solver options: maxiter
  best = minimize(
T12=1150.0,T3=2650.0,Trot=1150.0, Residual: 0.0143 🏆
T12=1170.0,T3=2650.0,Trot=1150.0, Residual: 0.0144
T12=1150.0,T3=2670.0,Trot=1150.0, Residual: 0.0143 🏆
T12=1150.0,T3=2650.0,Trot=1170.0, Residual: 0.0144
T12=1150.0,T3=2650.0,Trot=1150.0, Residual: 0.0143
T12=1170.0,T3=2650.0,Trot=1150.0, Residual: 0.0144
T12=1150.0,T3=2670.0,Trot=1150.0, Residual: 0.0143 🏆
T12=1150.0,T3=2650.0,Trot=1170.0, Residual: 0.0144
T12=791.0,T3=5000.0,Trot=464.0, Residual: 0.0127 🏆
T12=811.0,T3=5000.0,Trot=464.0, Residual: 0.0128
T12=791.0,T3=4980.0,Trot=464.0, Residual: 0.0125 🏆
T12=791.0,T3=5000.0,Trot=484.0, Residual: 0.0105 🏆
T12=970.0,T3=3825.0,Trot=807.0, Residual: 0.0078 🏆
T12=990.0,T3=3825.0,Trot=807.0, Residual: 0.0080
T12=970.0,T3=3845.0,Trot=807.0, Residual: 0.0078 🏆
T12=970.0,T3=3825.0,Trot=827.0, Residual: 0.0080
T12=880.0,T3=4412.0,Trot=635.0, Residual: 0.0059 🏆
T12=900.0,T3=4412.0,Trot=635.0, Residual: 0.0061
T12=880.0,T3=4432.0,Trot=635.0, Residual: 0.0059
T12=880.0,T3=4412.0,Trot=655.0, Residual: 0.0058 🏆
T12=880.0,T3=4413.0,Trot=635.0, Residual: 0.0059
T12=900.0,T3=4413.0,Trot=635.0, Residual: 0.0061
T12=880.0,T3=4433.0,Trot=635.0, Residual: 0.0059
T12=880.0,T3=4413.0,Trot=655.0, Residual: 0.0058 🏆
T12=604.0,T3=4938.0,Trot=546.0, Residual: 0.0054 🏆
T12=624.0,T3=4938.0,Trot=546.0, Residual: 0.0055
T12=604.0,T3=4958.0,Trot=546.0, Residual: 0.0055
T12=604.0,T3=4938.0,Trot=566.0, Residual: 0.0050 🏆
T12=742.0,T3=4675.0,Trot=590.0, Residual: 0.0050
T12=762.0,T3=4675.0,Trot=590.0, Residual: 0.0052
T12=742.0,T3=4695.0,Trot=590.0, Residual: 0.0050
T12=742.0,T3=4675.0,Trot=610.0, Residual: 0.0048 🏆
T12=673.0,T3=4807.0,Trot=568.0, Residual: 0.0050
T12=693.0,T3=4807.0,Trot=568.0, Residual: 0.0051
T12=673.0,T3=4827.0,Trot=568.0, Residual: 0.0050
T12=673.0,T3=4807.0,Trot=588.0, Residual: 0.0047 🏆
T12=673.0,T3=4807.0,Trot=568.0, Residual: 0.0050
T12=693.0,T3=4807.0,Trot=568.0, Residual: 0.0051
T12=673.0,T3=4827.0,Trot=568.0, Residual: 0.0050
T12=673.0,T3=4807.0,Trot=588.0, Residual: 0.0047 🏆
T12=389.0,T3=4894.0,Trot=793.0, Residual: 0.0050
T12=409.0,T3=4894.0,Trot=793.0, Residual: 0.0048
T12=389.0,T3=4914.0,Trot=793.0, Residual: 0.0050
T12=389.0,T3=4894.0,Trot=813.0, Residual: 0.0051
T12=531.0,T3=4850.0,Trot=680.0, Residual: 0.0038 🏆
T12=551.0,T3=4850.0,Trot=680.0, Residual: 0.0037 🏆
T12=531.0,T3=4870.0,Trot=680.0, Residual: 0.0038
T12=531.0,T3=4850.0,Trot=700.0, Residual: 0.0037
T12=531.0,T3=4850.0,Trot=680.0, Residual: 0.0038
T12=551.0,T3=4850.0,Trot=680.0, Residual: 0.0037 🏆
T12=531.0,T3=4870.0,Trot=680.0, Residual: 0.0038
T12=531.0,T3=4850.0,Trot=700.0, Residual: 0.0037
T12=571.0,T3=4833.0,Trot=715.0, Residual: 0.0036 🏆
T12=591.0,T3=4833.0,Trot=715.0, Residual: 0.0036
T12=571.0,T3=4853.0,Trot=715.0, Residual: 0.0036
T12=571.0,T3=4833.0,Trot=735.0, Residual: 0.0036
T12=571.0,T3=4833.0,Trot=715.0, Residual: 0.0036 🏆
T12=591.0,T3=4833.0,Trot=715.0, Residual: 0.0036
T12=571.0,T3=4853.0,Trot=715.0, Residual: 0.0036
T12=571.0,T3=4833.0,Trot=735.0, Residual: 0.0036
T12=654.0,T3=3333.0,Trot=631.0, Residual: 0.0033 🏆
T12=674.0,T3=3333.0,Trot=631.0, Residual: 0.0035
T12=654.0,T3=3353.0,Trot=631.0, Residual: 0.0033 🏆
T12=654.0,T3=3333.0,Trot=651.0, Residual: 0.0034
T12=612.0,T3=4083.0,Trot=673.0, Residual: 0.0031 🏆
T12=632.0,T3=4083.0,Trot=673.0, Residual: 0.0031
T12=612.0,T3=4103.0,Trot=673.0, Residual: 0.0031
T12=612.0,T3=4083.0,Trot=693.0, Residual: 0.0031
T12=633.0,T3=3708.0,Trot=652.0, Residual: 0.0030 🏆
T12=653.0,T3=3708.0,Trot=652.0, Residual: 0.0032
T12=633.0,T3=3728.0,Trot=652.0, Residual: 0.0030
T12=633.0,T3=3708.0,Trot=672.0, Residual: 0.0031
T12=633.0,T3=3708.0,Trot=652.0, Residual: 0.0030 🏆
T12=653.0,T3=3708.0,Trot=652.0, Residual: 0.0032
T12=633.0,T3=3728.0,Trot=652.0, Residual: 0.0030
T12=633.0,T3=3708.0,Trot=672.0, Residual: 0.0031
T12=475.0,T3=2044.0,Trot=425.0, Residual: 0.0042
T12=495.0,T3=2044.0,Trot=425.0, Residual: 0.0043
T12=475.0,T3=2064.0,Trot=425.0, Residual: 0.0040
T12=475.0,T3=2044.0,Trot=445.0, Residual: 0.0046
T12=554.0,T3=2876.0,Trot=538.0, Residual: 0.0025 🏆
T12=574.0,T3=2876.0,Trot=538.0, Residual: 0.0027
T12=554.0,T3=2896.0,Trot=538.0, Residual: 0.0025
T12=554.0,T3=2876.0,Trot=558.0, Residual: 0.0026
T12=515.0,T3=2460.0,Trot=482.0, Residual: 0.0026
T12=535.0,T3=2460.0,Trot=482.0, Residual: 0.0027
T12=515.0,T3=2480.0,Trot=482.0, Residual: 0.0025
T12=515.0,T3=2460.0,Trot=502.0, Residual: 0.0029
T12=535.0,T3=2668.0,Trot=510.0, Residual: 0.0025 🏆
T12=555.0,T3=2668.0,Trot=510.0, Residual: 0.0026
T12=535.0,T3=2688.0,Trot=510.0, Residual: 0.0024 🏆
T12=535.0,T3=2668.0,Trot=530.0, Residual: 0.0026
T12=525.0,T3=2564.0,Trot=496.0, Residual: 0.0025
T12=545.0,T3=2564.0,Trot=496.0, Residual: 0.0026
T12=525.0,T3=2584.0,Trot=496.0, Residual: 0.0025
T12=525.0,T3=2564.0,Trot=516.0, Residual: 0.0027
T12=530.0,T3=2616.0,Trot=503.0, Residual: 0.0025
T12=550.0,T3=2616.0,Trot=503.0, Residual: 0.0026
T12=530.0,T3=2636.0,Trot=503.0, Residual: 0.0024 🏆
T12=530.0,T3=2616.0,Trot=523.0, Residual: 0.0027
T12=532.0,T3=2642.0,Trot=507.0, Residual: 0.0025
T12=552.0,T3=2642.0,Trot=507.0, Residual: 0.0026
T12=532.0,T3=2662.0,Trot=507.0, Residual: 0.0024 🏆
T12=532.0,T3=2642.0,Trot=527.0, Residual: 0.0027
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2648.0,Trot=507.0, Residual: 0.0025
T12=553.0,T3=2648.0,Trot=507.0, Residual: 0.0026
T12=533.0,T3=2668.0,Trot=507.0, Residual: 0.0024
T12=533.0,T3=2648.0,Trot=527.0, Residual: 0.0026
T12=533.0,T3=2652.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2652.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2672.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2652.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2653.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2653.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2673.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2653.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2654.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2654.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2674.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2654.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2654.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2654.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2674.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2654.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=515.0,T3=2463.0,Trot=488.0, Residual: 0.0026
T12=535.0,T3=2463.0,Trot=488.0, Residual: 0.0027
T12=515.0,T3=2483.0,Trot=488.0, Residual: 0.0026
T12=515.0,T3=2463.0,Trot=508.0, Residual: 0.0030
T12=524.0,T3=2559.0,Trot=498.0, Residual: 0.0025
T12=544.0,T3=2559.0,Trot=498.0, Residual: 0.0026
T12=524.0,T3=2579.0,Trot=498.0, Residual: 0.0025
T12=524.0,T3=2559.0,Trot=518.0, Residual: 0.0028
T12=529.0,T3=2607.0,Trot=503.0, Residual: 0.0025
T12=549.0,T3=2607.0,Trot=503.0, Residual: 0.0026
T12=529.0,T3=2627.0,Trot=503.0, Residual: 0.0024
T12=529.0,T3=2607.0,Trot=523.0, Residual: 0.0027
T12=531.0,T3=2631.0,Trot=506.0, Residual: 0.0025
T12=551.0,T3=2631.0,Trot=506.0, Residual: 0.0026
T12=531.0,T3=2651.0,Trot=506.0, Residual: 0.0024
T12=531.0,T3=2631.0,Trot=526.0, Residual: 0.0027
T12=532.0,T3=2643.0,Trot=507.0, Residual: 0.0025
T12=552.0,T3=2643.0,Trot=507.0, Residual: 0.0026
T12=532.0,T3=2663.0,Trot=507.0, Residual: 0.0024
T12=532.0,T3=2643.0,Trot=527.0, Residual: 0.0027
T12=533.0,T3=2649.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2649.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2669.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2649.0,Trot=528.0, Residual: 0.0027
T12=533.0,T3=2652.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2652.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2672.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2652.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2653.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2653.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2673.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2653.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2654.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2654.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2674.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2654.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2654.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2654.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2674.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2654.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
T12=533.0,T3=2655.0,Trot=508.0, Residual: 0.0025
T12=553.0,T3=2655.0,Trot=508.0, Residual: 0.0026
T12=533.0,T3=2675.0,Trot=508.0, Residual: 0.0024
T12=533.0,T3=2655.0,Trot=528.0, Residual: 0.0026
Init ['T12', 'T3', 'Trot'] = [1150. 2650. 1150.]['', '', '']
Final ['T12', 'T3', 'Trot'] = [ 533. 2655.  508.]['', '', '']
Best ['T12', 'T3', 'Trot'] = [ 532.11786971 2661.78385922  506.54076838]['', '', ''] reached at iteration 104/368