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:5438: 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
    },
)
  • plot5 legacyFit Trot Tvib
  • plot5 legacyFit Trot Tvib
/home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/latest/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/latest/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/latest/radis/test/files/cdsd_hitemp_09_fragment.txt'], 'format': 'cdsd-hitemp', 'levelsfmt': 'radis'}



/home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/latest/radis/misc/warning.py:434: MissingReferenceWarning: Missing doi for CDSD-HITEMP. Use HITEMP-2010?
  warnings.warn(WarningType(message))
0.10s - Loaded database
Calculating Non-Equilibrium Spectrum
Physical Conditions
----------------------------------------
   Tgas                 491 K
   Trot                 491 K
   Tvib                 517,517,2641 K
   isotope              1,2
   medium               vacuum
   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_poly
   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
   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%
(0s)    5.3%
(1s)    7.9%
(1s)    10.5%
(1s)    13.2%
(2s)    15.8%
(2s)    18.4%
(2s)    21.1%
(2s)    23.7%
(3s)    26.3%
(3s)    28.9%
(3s)    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%
(4s)    60.5%
(4s)    63.2%
(4s)    65.8%
(4s)    68.4%
(4s)    71.1%
(4s)    73.7%
(4s)    76.3%
(4s)    78.9%
(4s)    81.6%
(4s)    84.2%
(4s)    86.8%
(4s)    89.5%
(4s)    92.1%
(4s)    94.7%
(4s)    97.4%
(4s)    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%
(2s)    15.8%
(2s)    18.4%
(2s)    21.1%
(2s)    23.7%
(2s)    26.3%
(3s)    28.9%
(3s)    31.6%
(3s)    34.2%
(3s)    36.8%
(3s)    39.5%
(3s)    42.1%
(3s)    44.7%
(3s)    47.4%
(3s)    50.0%
(4s)    52.6%
(4s)    55.3%
(4s)    57.9%
(4s)    60.5%
(4s)    63.2%
(4s)    65.8%
(4s)    68.4%
(4s)    71.1%
(4s)    73.7%
(4s)    76.3%
(4s)    78.9%
(4s)    81.6%
(4s)    84.2%
(4s)    86.8%
(4s)    89.5%
(4s)    92.1%
(4s)    94.7%
(4s)    97.4%
(4s)    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:434: NegativeEnergiesWarning: There are negative rotational energies in the database
  warnings.warn(WarningType(message))
13.30s - Spectrum calculated
/home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/latest/radis/misc/warning.py:434: NegativeEnergiesWarning: There are negative rotational energies in the database
  warnings.warn(WarningType(message))
------------------------------
TYPICAL FIT CALCULATION TIME:
Fit (in progress) profiler :
    spectrum_calculation      0.083s β–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆ
        check_line_databank                  0.001s
        check_non_eq_param                   0.000s
        reinitialize                         0.007s β–ˆ
            copy_database                        0.005s
            memory_usage_warning                 0.002s
            reset_population                     0.000s
        calc_noneq_population_multiTvib      0.053s β–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆ
            part_function                        0.049s β–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆβ–ˆ
            others                               0.004s
        scaled_non_eq_linestrength           0.002s
        scaled_S0                            0.001s
        calc_emission_integral               0.002s
        applied_linestrength_cutoff          0.002s
        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.003s
------------------------------
/home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/latest/radis/misc/warning.py:434: NegativeEnergiesWarning: There are negative rotational energies in the database
  warnings.warn(WarningType(message))
/home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/latest/radis/misc/warning.py:434: 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/latest/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

Total running time of the script: (0 minutes 44.472 seconds)