From e1df93e9a6f6588bc0a6850b282c9e6b857f59b6 Mon Sep 17 00:00:00 2001
From: Andrew Lock <69500916+andrewjlock@users.noreply.github.com>
Date: Wed, 30 Oct 2024 14:02:47 +1000
Subject: [PATCH] fix(OPM): oblique_rho2_rho1 function

Add missing brackets to oblique shock state relations
---
 src/pysagas/cfd/oblique_prandtl_meyer.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/pysagas/cfd/oblique_prandtl_meyer.py b/src/pysagas/cfd/oblique_prandtl_meyer.py
index fc3f5c5..51f2535 100644
--- a/src/pysagas/cfd/oblique_prandtl_meyer.py
+++ b/src/pysagas/cfd/oblique_prandtl_meyer.py
@@ -507,7 +507,7 @@ def oblique_rho2_rho1(M1: float, beta: float, gamma: float = 1.4):
         Peter Jacobs
         """
         M1n = M1 * abs(np.sin(beta))
-        rho2_rho1 = (gamma + 1) * M1n**2 / 2 + (gamma - 1) * M1n**2
+        rho2_rho1 = (gamma + 1) * M1n**2 / (2 + (gamma - 1) * M1n**2)
         return rho2_rho1
 
     @staticmethod