Note

You can download below the full example code and run it with πŸ”¬ Radis-Lab,

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 :

https://raw.githubusercontent.com/radis/radis-examples/master/docs/multi-temperature-fit.gif

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.

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/latest/radis/spectrum/spectrum.py:4921: 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()`

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
    },
)
  • plot5 legacyFit Trot Tvib
  • plot5 legacyFit Trot Tvib
/home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/latest/radis/db/molparam.py:252: FutureWarning:

The 'delim_whitespace' keyword in pd.read_csv is deprecated and will be removed in a future version. Use ``sep='\s+'`` instead

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/latest/radis/test/files/cdsd_hitemp_09_fragment.txt'], 'format': 'cdsd-hitemp', 'parfuncfmt': 'hapi', 'levelsfmt': 'radis'}



/home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/latest/radis/misc/warning.py:427: MissingReferenceWarning:

Missing doi for CDSD-HITEMP. Use HITEMP-2010?

Calculating Non-Equilibrium Spectrum
Physical Conditions
----------------------------------------
   Tgas                 491 K
   Trot                 491 K
   Tvib                 517,517,2641 K
   isotope              1,2
   mole_fraction        0.06573632856818698
   molecule             CO2
   path_length          10 cm
   pressure             0.02 bar
   rot_distribution     boltzmann
   self_absorption      True
   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/latest/radis/test/files/cdsd_hitemp...
   diluent              air
   folding_thresh       1e-06
   include_neighbouring_lines  True
   memory_mapping_engine  auto
   neighbour_lines      0 cm-1
   optimization         simple
   parfuncfmt           hapi
   parsum_mode          full summation
   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%
(2s)    7.9%
(2s)    10.5%
(3s)    13.2%
(3s)    15.8%
(3s)    18.4%
(4s)    21.1%
(4s)    23.7%
(4s)    26.3%
(5s)    28.9%
(5s)    31.6%
(5s)    34.2%
(5s)    36.8%
(6s)    39.5%
(6s)    42.1%
(6s)    44.7%
(6s)    47.4%
(6s)    50.0%
(6s)    52.6%
(6s)    55.3%
(6s)    57.9%
(6s)    60.5%
(6s)    63.2%
(6s)    65.8%
(6s)    68.4%
(6s)    71.1%
(6s)    73.7%
(7s)    76.3%
(7s)    78.9%
(7s)    81.6%
(7s)    84.2%
(7s)    86.8%
(7s)    89.5%
(7s)    92.1%
(7s)    94.7%
(7s)    97.4%
(7s)    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%
(1s)    2.6%
(1s)    5.3%
(2s)    7.9%
(2s)    10.5%
(3s)    13.2%
(3s)    15.8%
(4s)    18.4%
(4s)    21.1%
(4s)    23.7%
(5s)    26.3%
(5s)    28.9%
(5s)    31.6%
(5s)    34.2%
(5s)    36.8%
(5s)    39.5%
(6s)    42.1%
(6s)    44.7%
(6s)    47.4%
(6s)    50.0%
(6s)    52.6%
(6s)    55.3%
(6s)    57.9%
(7s)    60.5%
(7s)    63.2%
(7s)    65.8%
(7s)    68.4%
(7s)    71.1%
(7s)    73.7%
(7s)    76.3%
(7s)    78.9%
(7s)    81.6%
(7s)    84.2%
(7s)    86.8%
(7s)    89.5%
(7s)    92.1%
(7s)    94.7%
(7s)    97.4%
(7s)    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/latest/radis/misc/warning.py:427: NegativeEnergiesWarning:

There are negative rotational energies in the database

22.78s - Spectrum calculated
/home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/latest/radis/misc/warning.py:427: NegativeEnergiesWarning:

There are negative rotational energies in the database

------------------------------
TYPICAL FIT CALCULATION TIME:
Fit (in progress) profiler :
    spectrum_calculation      0.106s β–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆ
        check_line_databank                  0.002s
        check_non_eq_param                   0.000s
        reinitialize                         0.005s
            copy_database                        0.001s
            memory_usage_warning                 0.003s
            reset_population                     0.000s
        calc_noneq_population_multiTvib      0.071s β–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆ
            part_function                        0.066s β–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆ
            others                               0.006s
        scaled_non_eq_linestrength           0.003s
        calc_emission_integral               0.003s
        applied_linestrength_cutoff          0.001s
        calc_lineshift                       0.001s
        calc_hwhm                            0.005s
        generate_wavenumber_arrays           0.000s
        calc_line_broadening                 0.014s β–ˆβ–ˆ
            precompute_LDM_lineshapes            0.001s
            LDM_Initialized_vectors              0.000s
            LDM_closest_matching_line            0.000s
            LDM_Distribute_lines                 0.009s β–ˆ
            LDM_convolve                         0.003s
        calc_other_spectral_quan             0.000s
        generate_spectrum_obj                0.000s
        others                               0.007s
------------------------------
/home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/latest/radis/misc/warning.py:427: NegativeEnergiesWarning:

There are negative rotational energies in the database

/home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/latest/radis/misc/warning.py:427: NegativeEnergiesWarning:

There are negative rotational energies in the database


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

/home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/latest/radis/tools/fitting.py:486: OptimizeWarning:

Unknown solver options: maxiter

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=782.0,T3=5000.0,Trot=459.0, Residual: 0.0133  πŸ†
T12=802.0,T3=5000.0,Trot=459.0, Residual: 0.0134
T12=782.0,T3=4980.0,Trot=459.0, Residual: 0.0131  πŸ†
T12=782.0,T3=5000.0,Trot=479.0, Residual: 0.0110  πŸ†
T12=966.0,T3=3825.0,Trot=804.0, Residual: 0.0077  πŸ†
T12=986.0,T3=3825.0,Trot=804.0, Residual: 0.0080
T12=966.0,T3=3845.0,Trot=804.0, Residual: 0.0077  πŸ†
T12=966.0,T3=3825.0,Trot=824.0, Residual: 0.0079
T12=874.0,T3=4412.0,Trot=632.0, Residual: 0.0059  πŸ†
T12=894.0,T3=4412.0,Trot=632.0, Residual: 0.0061
T12=874.0,T3=4432.0,Trot=632.0, Residual: 0.0059
T12=874.0,T3=4412.0,Trot=652.0, Residual: 0.0058  πŸ†
T12=874.0,T3=4413.0,Trot=632.0, Residual: 0.0059
T12=894.0,T3=4413.0,Trot=632.0, Residual: 0.0061
T12=874.0,T3=4433.0,Trot=632.0, Residual: 0.0059
T12=874.0,T3=4413.0,Trot=652.0, Residual: 0.0058  πŸ†
T12=576.0,T3=4971.0,Trot=546.0, Residual: 0.0056  πŸ†
T12=596.0,T3=4971.0,Trot=546.0, Residual: 0.0056
T12=576.0,T3=4991.0,Trot=546.0, Residual: 0.0056
T12=576.0,T3=4971.0,Trot=566.0, Residual: 0.0051  πŸ†
T12=725.0,T3=4692.0,Trot=589.0, Residual: 0.0050  πŸ†
T12=745.0,T3=4692.0,Trot=589.0, Residual: 0.0051
T12=725.0,T3=4712.0,Trot=589.0, Residual: 0.0050
T12=725.0,T3=4692.0,Trot=609.0, Residual: 0.0048  πŸ†
T12=650.0,T3=4832.0,Trot=567.0, Residual: 0.0050
T12=670.0,T3=4832.0,Trot=567.0, Residual: 0.0051
T12=650.0,T3=4852.0,Trot=567.0, Residual: 0.0050
T12=650.0,T3=4832.0,Trot=587.0, Residual: 0.0047  πŸ†
T12=688.0,T3=4762.0,Trot=578.0, Residual: 0.0049
T12=708.0,T3=4762.0,Trot=578.0, Residual: 0.0051
T12=688.0,T3=4782.0,Trot=578.0, Residual: 0.0050
T12=688.0,T3=4762.0,Trot=598.0, Residual: 0.0047  πŸ†
T12=688.0,T3=4762.0,Trot=578.0, Residual: 0.0049
T12=708.0,T3=4762.0,Trot=578.0, Residual: 0.0051
T12=688.0,T3=4782.0,Trot=578.0, Residual: 0.0050
T12=688.0,T3=4762.0,Trot=598.0, Residual: 0.0047  πŸ†
T12=315.0,T3=5000.0,Trot=850.0, Residual: 0.0059
T12=335.0,T3=5000.0,Trot=850.0, Residual: 0.0057
T12=315.0,T3=4980.0,Trot=850.0, Residual: 0.0059
T12=315.0,T3=5000.0,Trot=870.0, Residual: 0.0060
T12=501.0,T3=4881.0,Trot=714.0, Residual: 0.0039  πŸ†
T12=521.0,T3=4881.0,Trot=714.0, Residual: 0.0038  πŸ†
T12=501.0,T3=4901.0,Trot=714.0, Residual: 0.0039
T12=501.0,T3=4881.0,Trot=734.0, Residual: 0.0039
T12=594.0,T3=4821.0,Trot=646.0, Residual: 0.0039
T12=614.0,T3=4821.0,Trot=646.0, Residual: 0.0040
T12=594.0,T3=4841.0,Trot=646.0, Residual: 0.0040
T12=594.0,T3=4821.0,Trot=666.0, Residual: 0.0038  πŸ†
T12=548.0,T3=4851.0,Trot=680.0, Residual: 0.0038  πŸ†
T12=568.0,T3=4851.0,Trot=680.0, Residual: 0.0038  πŸ†
T12=548.0,T3=4871.0,Trot=680.0, Residual: 0.0038
T12=548.0,T3=4851.0,Trot=700.0, Residual: 0.0037  πŸ†
T12=548.0,T3=4851.0,Trot=680.0, Residual: 0.0038
T12=568.0,T3=4851.0,Trot=680.0, Residual: 0.0038
T12=548.0,T3=4871.0,Trot=680.0, Residual: 0.0038
T12=548.0,T3=4851.0,Trot=700.0, Residual: 0.0037  πŸ†
T12=584.0,T3=4793.0,Trot=714.0, Residual: 0.0036  πŸ†
T12=604.0,T3=4793.0,Trot=714.0, Residual: 0.0037
T12=584.0,T3=4813.0,Trot=714.0, Residual: 0.0036
T12=584.0,T3=4793.0,Trot=734.0, Residual: 0.0036
T12=584.0,T3=4793.0,Trot=714.0, Residual: 0.0036  πŸ†
T12=604.0,T3=4793.0,Trot=714.0, Residual: 0.0037
T12=584.0,T3=4813.0,Trot=714.0, Residual: 0.0036
T12=584.0,T3=4793.0,Trot=734.0, Residual: 0.0036
T12=584.0,T3=3812.0,Trot=702.0, Residual: 0.0032  πŸ†
T12=604.0,T3=3812.0,Trot=702.0, Residual: 0.0032  πŸ†
T12=584.0,T3=3832.0,Trot=702.0, Residual: 0.0032
T12=584.0,T3=3812.0,Trot=722.0, Residual: 0.0034
T12=584.0,T3=4302.0,Trot=708.0, Residual: 0.0032
T12=604.0,T3=4302.0,Trot=708.0, Residual: 0.0032
T12=584.0,T3=4322.0,Trot=708.0, Residual: 0.0033
T12=584.0,T3=4302.0,Trot=728.0, Residual: 0.0033
T12=584.0,T3=4057.0,Trot=705.0, Residual: 0.0032  πŸ†
T12=604.0,T3=4057.0,Trot=705.0, Residual: 0.0032  πŸ†
T12=584.0,T3=4077.0,Trot=705.0, Residual: 0.0032
T12=584.0,T3=4057.0,Trot=725.0, Residual: 0.0033
T12=584.0,T3=4057.0,Trot=705.0, Residual: 0.0032
T12=604.0,T3=4057.0,Trot=705.0, Residual: 0.0032  πŸ†
T12=584.0,T3=4077.0,Trot=705.0, Residual: 0.0032
T12=584.0,T3=4057.0,Trot=725.0, Residual: 0.0033
T12=599.0,T3=2565.0,Trot=470.0, Residual: 0.0038
T12=619.0,T3=2565.0,Trot=470.0, Residual: 0.0042
T12=599.0,T3=2585.0,Trot=470.0, Residual: 0.0038
T12=599.0,T3=2565.0,Trot=490.0, Residual: 0.0035
T12=591.0,T3=3311.0,Trot=587.0, Residual: 0.0028  πŸ†
T12=611.0,T3=3311.0,Trot=587.0, Residual: 0.0029
T12=591.0,T3=3331.0,Trot=587.0, Residual: 0.0028
T12=591.0,T3=3311.0,Trot=607.0, Residual: 0.0028  πŸ†
T12=595.0,T3=2938.0,Trot=529.0, Residual: 0.0030
T12=615.0,T3=2938.0,Trot=529.0, Residual: 0.0033
T12=595.0,T3=2958.0,Trot=529.0, Residual: 0.0030
T12=595.0,T3=2938.0,Trot=549.0, Residual: 0.0029
T12=593.0,T3=3124.0,Trot=558.0, Residual: 0.0029
T12=613.0,T3=3124.0,Trot=558.0, Residual: 0.0031
T12=593.0,T3=3144.0,Trot=558.0, Residual: 0.0029
T12=593.0,T3=3124.0,Trot=578.0, Residual: 0.0028  πŸ†
T12=592.0,T3=3218.0,Trot=573.0, Residual: 0.0028
T12=612.0,T3=3218.0,Trot=573.0, Residual: 0.0030
T12=592.0,T3=3238.0,Trot=573.0, Residual: 0.0028
T12=592.0,T3=3218.0,Trot=593.0, Residual: 0.0028  πŸ†
T12=592.0,T3=3264.0,Trot=580.0, Residual: 0.0028
T12=612.0,T3=3264.0,Trot=580.0, Residual: 0.0030
T12=592.0,T3=3284.0,Trot=580.0, Residual: 0.0028
T12=592.0,T3=3264.0,Trot=600.0, Residual: 0.0028
T12=592.0,T3=3288.0,Trot=584.0, Residual: 0.0028
T12=612.0,T3=3288.0,Trot=584.0, Residual: 0.0029
T12=592.0,T3=3308.0,Trot=584.0, Residual: 0.0028
T12=592.0,T3=3288.0,Trot=604.0, Residual: 0.0028
T12=592.0,T3=3299.0,Trot=586.0, Residual: 0.0028
T12=612.0,T3=3299.0,Trot=586.0, Residual: 0.0029
T12=592.0,T3=3319.0,Trot=586.0, Residual: 0.0028
T12=592.0,T3=3299.0,Trot=606.0, Residual: 0.0028
T12=592.0,T3=3305.0,Trot=587.0, Residual: 0.0028
T12=612.0,T3=3305.0,Trot=587.0, Residual: 0.0029
T12=592.0,T3=3325.0,Trot=587.0, Residual: 0.0028
T12=592.0,T3=3305.0,Trot=607.0, Residual: 0.0028
T12=591.0,T3=3308.0,Trot=587.0, Residual: 0.0028
T12=611.0,T3=3308.0,Trot=587.0, Residual: 0.0029
T12=591.0,T3=3328.0,Trot=587.0, Residual: 0.0028
T12=591.0,T3=3308.0,Trot=607.0, Residual: 0.0028
T12=591.0,T3=3309.0,Trot=587.0, Residual: 0.0028
T12=611.0,T3=3309.0,Trot=587.0, Residual: 0.0029
T12=591.0,T3=3329.0,Trot=587.0, Residual: 0.0028
T12=591.0,T3=3309.0,Trot=607.0, Residual: 0.0028
T12=591.0,T3=3310.0,Trot=587.0, Residual: 0.0028
T12=611.0,T3=3310.0,Trot=587.0, Residual: 0.0029
T12=591.0,T3=3330.0,Trot=587.0, Residual: 0.0028
T12=591.0,T3=3310.0,Trot=607.0, Residual: 0.0028
T12=591.0,T3=3311.0,Trot=587.0, Residual: 0.0028
T12=611.0,T3=3311.0,Trot=587.0, Residual: 0.0029
T12=591.0,T3=3331.0,Trot=587.0, Residual: 0.0028
T12=591.0,T3=3311.0,Trot=607.0, Residual: 0.0028
T12=591.0,T3=3311.0,Trot=587.0, Residual: 0.0028
T12=611.0,T3=3311.0,Trot=587.0, Residual: 0.0029
T12=591.0,T3=3331.0,Trot=587.0, Residual: 0.0028
T12=591.0,T3=3311.0,Trot=607.0, Residual: 0.0028
T12=591.0,T3=3311.0,Trot=587.0, Residual: 0.0028
T12=611.0,T3=3311.0,Trot=587.0, Residual: 0.0029
T12=591.0,T3=3331.0,Trot=587.0, Residual: 0.0028
T12=591.0,T3=3311.0,Trot=607.0, Residual: 0.0028
T12=591.0,T3=3311.0,Trot=587.0, Residual: 0.0028
T12=611.0,T3=3311.0,Trot=587.0, Residual: 0.0029
T12=591.0,T3=3331.0,Trot=587.0, Residual: 0.0028
T12=591.0,T3=3311.0,Trot=607.0, Residual: 0.0028
T12=591.0,T3=3311.0,Trot=587.0, Residual: 0.0028
T12=611.0,T3=3311.0,Trot=587.0, Residual: 0.0029
T12=591.0,T3=3331.0,Trot=587.0, Residual: 0.0028
T12=591.0,T3=3311.0,Trot=607.0, Residual: 0.0028
T12=591.0,T3=3311.0,Trot=587.0, Residual: 0.0028
T12=611.0,T3=3311.0,Trot=587.0, Residual: 0.0029
T12=591.0,T3=3331.0,Trot=587.0, Residual: 0.0028
T12=591.0,T3=3311.0,Trot=607.0, Residual: 0.0028
T12=591.0,T3=3311.0,Trot=587.0, Residual: 0.0028
T12=611.0,T3=3311.0,Trot=587.0, Residual: 0.0029
T12=591.0,T3=3331.0,Trot=587.0, Residual: 0.0028
T12=591.0,T3=3311.0,Trot=607.0, Residual: 0.0028
T12=591.0,T3=3311.0,Trot=587.0, Residual: 0.0028
T12=611.0,T3=3311.0,Trot=587.0, Residual: 0.0029
T12=591.0,T3=3331.0,Trot=587.0, Residual: 0.0028
T12=591.0,T3=3311.0,Trot=607.0, Residual: 0.0028
T12=591.0,T3=3311.0,Trot=587.0, Residual: 0.0028
T12=611.0,T3=3311.0,Trot=587.0, Residual: 0.0029
T12=591.0,T3=3331.0,Trot=587.0, Residual: 0.0028
T12=591.0,T3=3311.0,Trot=607.0, Residual: 0.0028
T12=591.0,T3=3311.0,Trot=587.0, Residual: 0.0028
T12=611.0,T3=3311.0,Trot=587.0, Residual: 0.0029
T12=591.0,T3=3331.0,Trot=587.0, Residual: 0.0028
T12=591.0,T3=3311.0,Trot=607.0, Residual: 0.0028
T12=591.0,T3=3311.0,Trot=587.0, Residual: 0.0028
T12=611.0,T3=3311.0,Trot=587.0, Residual: 0.0029
T12=591.0,T3=3331.0,Trot=587.0, Residual: 0.0028
T12=591.0,T3=3311.0,Trot=607.0, Residual: 0.0028
T12=591.0,T3=3311.0,Trot=587.0, Residual: 0.0028
T12=611.0,T3=3311.0,Trot=587.0, Residual: 0.0029
T12=591.0,T3=3331.0,Trot=587.0, Residual: 0.0028
T12=591.0,T3=3311.0,Trot=607.0, Residual: 0.0028
T12=591.0,T3=3311.0,Trot=587.0, Residual: 0.0028
T12=611.0,T3=3311.0,Trot=587.0, Residual: 0.0029
T12=591.0,T3=3331.0,Trot=587.0, Residual: 0.0028
T12=591.0,T3=3311.0,Trot=607.0, Residual: 0.0028
T12=591.0,T3=3311.0,Trot=587.0, Residual: 0.0028
T12=611.0,T3=3311.0,Trot=587.0, Residual: 0.0029
T12=591.0,T3=3331.0,Trot=587.0, Residual: 0.0028
T12=591.0,T3=3311.0,Trot=607.0, Residual: 0.0028
T12=434.0,T3=1587.0,Trot=300.0, Residual: 0.0083
T12=454.0,T3=1587.0,Trot=300.0, Residual: 0.0083
T12=434.0,T3=1607.0,Trot=300.0, Residual: 0.0080
T12=434.0,T3=1587.0,Trot=320.0, Residual: 0.0084
T12=513.0,T3=2449.0,Trot=444.0, Residual: 0.0029
T12=533.0,T3=2449.0,Trot=444.0, Residual: 0.0031
T12=513.0,T3=2469.0,Trot=444.0, Residual: 0.0030
T12=513.0,T3=2449.0,Trot=464.0, Residual: 0.0026  πŸ†
T12=552.0,T3=2880.0,Trot=516.0, Residual: 0.0027
T12=572.0,T3=2880.0,Trot=516.0, Residual: 0.0029
T12=552.0,T3=2900.0,Trot=516.0, Residual: 0.0027
T12=552.0,T3=2880.0,Trot=536.0, Residual: 0.0025  πŸ†
T12=532.0,T3=2664.0,Trot=480.0, Residual: 0.0027
T12=552.0,T3=2664.0,Trot=480.0, Residual: 0.0029
T12=532.0,T3=2684.0,Trot=480.0, Residual: 0.0028
T12=532.0,T3=2664.0,Trot=500.0, Residual: 0.0025  πŸ†
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=537.0,T3=2718.0,Trot=489.0, Residual: 0.0027
T12=557.0,T3=2718.0,Trot=489.0, Residual: 0.0029
T12=537.0,T3=2738.0,Trot=489.0, Residual: 0.0027
T12=537.0,T3=2718.0,Trot=509.0, Residual: 0.0025
T12=540.0,T3=2745.0,Trot=493.0, Residual: 0.0027
T12=560.0,T3=2745.0,Trot=493.0, Residual: 0.0029
T12=540.0,T3=2765.0,Trot=493.0, Residual: 0.0027
T12=540.0,T3=2745.0,Trot=513.0, Residual: 0.0025
T12=541.0,T3=2759.0,Trot=495.0, Residual: 0.0027
T12=561.0,T3=2759.0,Trot=495.0, Residual: 0.0029
T12=541.0,T3=2779.0,Trot=495.0, Residual: 0.0027
T12=541.0,T3=2759.0,Trot=515.0, Residual: 0.0025
T12=542.0,T3=2765.0,Trot=496.0, Residual: 0.0027
T12=562.0,T3=2765.0,Trot=496.0, Residual: 0.0029
T12=542.0,T3=2785.0,Trot=496.0, Residual: 0.0027
T12=542.0,T3=2765.0,Trot=516.0, Residual: 0.0025
T12=542.0,T3=2769.0,Trot=497.0, Residual: 0.0027
T12=562.0,T3=2769.0,Trot=497.0, Residual: 0.0029
T12=542.0,T3=2789.0,Trot=497.0, Residual: 0.0027
T12=542.0,T3=2769.0,Trot=517.0, Residual: 0.0025
T12=542.0,T3=2771.0,Trot=497.0, Residual: 0.0027
T12=562.0,T3=2771.0,Trot=497.0, Residual: 0.0029
T12=542.0,T3=2791.0,Trot=497.0, Residual: 0.0027
T12=542.0,T3=2771.0,Trot=517.0, Residual: 0.0025
T12=542.0,T3=2771.0,Trot=497.0, Residual: 0.0027
T12=562.0,T3=2771.0,Trot=497.0, Residual: 0.0029
T12=542.0,T3=2791.0,Trot=497.0, Residual: 0.0027
T12=542.0,T3=2771.0,Trot=517.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=490.0,T3=2299.0,Trot=407.0, Residual: 0.0035
T12=510.0,T3=2299.0,Trot=407.0, Residual: 0.0037
T12=490.0,T3=2319.0,Trot=407.0, Residual: 0.0036
T12=490.0,T3=2299.0,Trot=427.0, Residual: 0.0031
T12=516.0,T3=2536.0,Trot=452.0, Residual: 0.0029
T12=536.0,T3=2536.0,Trot=452.0, Residual: 0.0031
T12=516.0,T3=2556.0,Trot=452.0, Residual: 0.0030
T12=516.0,T3=2536.0,Trot=472.0, Residual: 0.0025
T12=529.0,T3=2654.0,Trot=475.0, Residual: 0.0028
T12=549.0,T3=2654.0,Trot=475.0, Residual: 0.0029
T12=529.0,T3=2674.0,Trot=475.0, Residual: 0.0028
T12=529.0,T3=2654.0,Trot=495.0, Residual: 0.0025
T12=536.0,T3=2713.0,Trot=486.0, Residual: 0.0027
T12=556.0,T3=2713.0,Trot=486.0, Residual: 0.0029
T12=536.0,T3=2733.0,Trot=486.0, Residual: 0.0028
T12=536.0,T3=2713.0,Trot=506.0, Residual: 0.0025
T12=539.0,T3=2743.0,Trot=492.0, Residual: 0.0027
T12=559.0,T3=2743.0,Trot=492.0, Residual: 0.0029
T12=539.0,T3=2763.0,Trot=492.0, Residual: 0.0027
T12=539.0,T3=2743.0,Trot=512.0, Residual: 0.0025
T12=541.0,T3=2757.0,Trot=495.0, Residual: 0.0027
T12=561.0,T3=2757.0,Trot=495.0, Residual: 0.0029
T12=541.0,T3=2777.0,Trot=495.0, Residual: 0.0027
T12=541.0,T3=2757.0,Trot=515.0, Residual: 0.0025
T12=541.0,T3=2765.0,Trot=496.0, Residual: 0.0027
T12=561.0,T3=2765.0,Trot=496.0, Residual: 0.0029
T12=541.0,T3=2785.0,Trot=496.0, Residual: 0.0027
T12=541.0,T3=2765.0,Trot=516.0, Residual: 0.0025
T12=542.0,T3=2768.0,Trot=497.0, Residual: 0.0027
T12=562.0,T3=2768.0,Trot=497.0, Residual: 0.0029
T12=542.0,T3=2788.0,Trot=497.0, Residual: 0.0027
T12=542.0,T3=2768.0,Trot=517.0, Residual: 0.0025
T12=542.0,T3=2770.0,Trot=497.0, Residual: 0.0027
T12=562.0,T3=2770.0,Trot=497.0, Residual: 0.0029
T12=542.0,T3=2790.0,Trot=497.0, Residual: 0.0027
T12=542.0,T3=2770.0,Trot=517.0, Residual: 0.0025
T12=542.0,T3=2771.0,Trot=497.0, Residual: 0.0027
T12=562.0,T3=2771.0,Trot=497.0, Residual: 0.0029
T12=542.0,T3=2791.0,Trot=497.0, Residual: 0.0027
T12=542.0,T3=2771.0,Trot=517.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
T12=542.0,T3=2772.0,Trot=498.0, Residual: 0.0027
T12=562.0,T3=2772.0,Trot=498.0, Residual: 0.0029
T12=542.0,T3=2792.0,Trot=498.0, Residual: 0.0027
T12=542.0,T3=2772.0,Trot=518.0, Residual: 0.0025
Best ['T12', 'T3', 'Trot'] = [ 532.40831841 2664.44122594  499.63988523]['', '', ''] reached at iteration 209/404

Total running time of the script: (1 minutes 7.698 seconds)