From 140c60bf3218d6e39264642891888525d6ddf33e Mon Sep 17 00:00:00 2001
From: Javier Emilio Alfonso Ramos
 <58225254+javialra97@users.noreply.github.com>
Date: Tue, 5 Nov 2024 15:09:49 +0100
Subject: [PATCH] Minor improvements for generation of dataset of reaction
 profiles. (#359)

* tight criteria

* csv output

* single point refinement bool

* bug #358

* Revert "csv output"

This reverts commit 1f4c6e027dc19a159db344a7f4cb7aa6883c7096.

* Revert "tight criteria"

This reverts commit 675b82ec32274719a59571b92c2b3c67ae90dd0f.

* template

* update changelog and readme

* parentheses

* provisional patch

---------

Co-authored-by: t-young31 <tom.young.316@gmail.com>
---
 README.md                    | 1 +
 autode/reactions/reaction.py | 7 ++++++-
 autode/wrappers/G09.py       | 4 +++-
 doc/changelog.rst            | 4 ++--
 4 files changed, 12 insertions(+), 4 deletions(-)

diff --git a/README.md b/README.md
index 558dbe879..7b2b0e17b 100644
--- a/README.md
+++ b/README.md
@@ -97,3 +97,4 @@ If **autodE** is used in a publication please consider citing the [paper](https:
 - Daniel Hollas ([@danielhollas](https://github.com/danielhollas))
 - Nils Heunemann ([@nilsheunemann](https://github.com/NilsHeunemann))
 - Sijie Fu ([@sijiefu](https://github.com/SijieFu))
+- Javier Alfonso ([@javialra97](https://github.com/javialra97))
diff --git a/autode/reactions/reaction.py b/autode/reactions/reaction.py
index e24899bce..d0ae649df 100644
--- a/autode/reactions/reaction.py
+++ b/autode/reactions/reaction.py
@@ -119,6 +119,7 @@ def __str__(self):
     def calculate_reaction_profile(
         self,
         units: Union["Unit", str] = "kcal mol-1",
+        single_point_refinement: bool = True,
         with_complexes: bool = False,
         free_energy: bool = False,
         enthalpy: bool = False,
@@ -132,6 +133,9 @@ def calculate_reaction_profile(
         Keyword Arguments:
             units (autode.units.Unit | str):
 
+            single_point_refinement (bool): Calculate single point at a
+                                            higher level of theory
+
             with_complexes (bool): Calculate the lowest energy conformers
                                    of the reactant and product complexes
 
@@ -161,7 +165,8 @@ def calculate(reaction):
                 reaction.calculate_complexes()
             if free_energy or enthalpy:
                 reaction.calculate_thermochemical_cont()
-            reaction.calculate_single_points()
+            if single_point_refinement:
+                reaction.calculate_single_points()
             reaction.print_output()
             return None
 
diff --git a/autode/wrappers/G09.py b/autode/wrappers/G09.py
index 684521fa6..b2cfcbf05 100644
--- a/autode/wrappers/G09.py
+++ b/autode/wrappers/G09.py
@@ -546,7 +546,9 @@ def coordinates_from(self, calc: "CalculationExecutor") -> Coordinates:
         coords: List[List[float]] = []
 
         for i, line in enumerate(calc.output.file_lines):
-            if "Input orientation" in line:
+            if "Input orientation" in line or (
+                "Standard orientation" in line and len(coords) == 0
+            ):
                 coords.clear()
                 xyz_lines = calc.output.file_lines[
                     i + 5 : i + 5 + calc.molecule.n_atoms
diff --git a/doc/changelog.rst b/doc/changelog.rst
index 031d97403..38d9b3e3c 100644
--- a/doc/changelog.rst
+++ b/doc/changelog.rst
@@ -7,11 +7,11 @@ Changelog
 
 Functionality improvements
 **************************
-- ...
+- Adds a boolean option for calculating single points refinements
 
 Bug Fixes
 *********
-- ...
+- Fixes coordinate extraction in some G16 output files
 
 
 1.4.4