From da80977954ce38820ed7135a868b5c2bee130488 Mon Sep 17 00:00:00 2001
From: doreenfan <35239173+doreenfan@users.noreply.github.com>
Date: Sun, 31 May 2020 09:05:50 -0700
Subject: [PATCH] bug fix for SDC_EVOLVE_ENTHALPY (#317)

---
 integration/VODE/vode_type_simplified_sdc.F90 | 4 ++--
 interfaces/burn_type.F90                      | 5 +++++
 2 files changed, 7 insertions(+), 2 deletions(-)

diff --git a/integration/VODE/vode_type_simplified_sdc.F90 b/integration/VODE/vode_type_simplified_sdc.F90
index fa12392bc7..b2881a90c0 100644
--- a/integration/VODE/vode_type_simplified_sdc.F90
+++ b/integration/VODE/vode_type_simplified_sdc.F90
@@ -544,8 +544,8 @@ subroutine jac_to_vode(time, jac_react, vode_state, jac)
 
     eos_state % rho = vode_state % rpar(irp_SRHO)
     eos_state % T = 1.e4   ! initial guess
-    eos_state % xn(:) = y(SFS:SFS-1+nspec)/vode_state % rpar(irp_SRHO)
-    eos_state % h = y(SENTH)/vode_state % rpar(irp_SRHO)
+    eos_state % xn(:) = vode_state % y(SFS:SFS-1+nspec)/vode_state % rpar(irp_SRHO)
+    eos_state % h = vode_state % y(SENTH)/vode_state % rpar(irp_SRHO)
 
     call eos(eos_input_rh, eos_state)
 
diff --git a/interfaces/burn_type.F90 b/interfaces/burn_type.F90
index 72d5254a93..427251bee2 100644
--- a/interfaces/burn_type.F90
+++ b/interfaces/burn_type.F90
@@ -44,6 +44,11 @@ module burn_type_module
     real(rt) :: aux(naux)
 #endif
 
+#if SDC_EVOLVE_ENTHALPY
+    ! make pressure available to RHS
+    real(rt) :: p0
+#endif
+    
     real(rt) :: cv
     real(rt) :: cp
     real(rt) :: y_e