# Reference example: minimal steady-state solve.
#
# This input is intentionally tiny and is useful when you want the smallest
# possible example of `steady_state_solve`.

# Species and their initial guess for the steady-state solver.
components {
  A* {phase=surface, init=0.2},
  B* {phase=surface, init=0.8}
}

# A reversible interconversion with zero activation barriers.
reactions {
  {A*}=>{B*} @ ArrheniusDefault(Vf=1.0, Vb=1.0, Eaf=0.0, Eab=0.0)
}

# Direct steady-state solve using KINSOL.
simulation {
  eq @ steady_state_solve(
    model=thermal,
    solver={method=kinsol_newton, fnormtol=1e-12, scsteptol=1e-12, max_iter=50, strategy=line_search},
    conditions=[{T=600}, {T=650}, {T=700}]
  )
}

# Output files for the steady-state result.
output {
  formats=[tabdelim, hdf5],
  folder="results_steady_state"
}
