Merge pull request #28 from ipqa-research/dev

fixed PR78
ipqa-research · Feb 20, 2024 · ed9c128 · ed9c128
2 parents 5516984 + 74310a0
commit ed9c128
Showing 1 changed file with 2 additions and 5 deletions.
diff --git a/src/models/residual_helmholtz/cubic/implementations/implementations.f90 b/src/models/residual_helmholtz/cubic/implementations/implementations.f90
@@ -125,12 +125,9 @@ type(CubicEoS) function PengRobinson78(tc, pc, w, kij, lij) result(model)
         composition%pc = pc
         composition%w = w
 
-        alpha%k = 0.37464_pr &
-                 + 1.54226_pr * composition%w &
-                 - 0.26993_pr * composition%w**2
-
+        allocate(alpha%k(nc))
         where (composition%w <=0.491)
-              alpha%k = 0.37464 + 1.54226 * composition%w(i) - 0.26992 * composition%w**2
+              alpha%k = 0.37464 + 1.54226 * composition%w - 0.26992 * composition%w**2
         elsewhere
               alpha%k = 0.379642 + 1.48503 * composition%w - 0.164423 * composition%w**2 + 0.016666 * composition%w**3
         end where