# Data sources
database(
    thermoLibraries = ['primaryThermoLibrary'],
    reactionLibraries = [],
    seedMechanisms = [],
    kineticsDepositories = ['training'],
    kineticsFamilies = 'default',
    kineticsEstimator = 'rate rules',
)

# Constraints on generated species
generatedSpeciesConstraints(
    maximumRadicalElectrons = 3,
)

# List of species
species(
    label='reactant1',
    reactive=True,
    structure=SMILES("O=C(CCl)NCC1CCC(O)CC1"),
)

species(
    label='reactant2',
    reactive=True,
    structure=SMILES("Cc1n[nH]cc1B1OC(C)(C)C(C)(C)O1"),
)
species(
    label='water',
    reactive=True,
    structure=SMILES("O"),
)
# Reaction systems
liquidReactor(
    temperature=(300,'K'),
    initialConcentrations={
        "reactant1": (6.154e-3,'mol/cm^3'),
        "reactant2": (6.154e-3,'mol/cm^3'),
        "water": (12.154e-3,'mol/cm^3')
    },
    terminationTime=(5,'s'),
)

solvation(
	solvent='water'
)

simulator(
    atol=1e-16,
    rtol=1e-8,
)

model(
    toleranceKeepInEdge=1E-9,
    toleranceMoveToCore=0.001,
    toleranceInterruptSimulation=0.1,
    maximumEdgeSpecies=100000
)

options(
    units='si',
    generateOutputHTML=True,
    generatePlots=True,
    saveEdgeSpecies=True,
    saveSimulationProfiles=True,
)