From d00a530d2d299487f2eab6f2daae9258ef5d394c Mon Sep 17 00:00:00 2001
From: Martin Diehl <mail@martin-diehl.net>
Date: Sun, 4 Jun 2023 10:35:20 +0200
Subject: [PATCH] f_phi factor corrected

---
 PRIVATE              | 2 +-
 src/phase_damage.f90 | 2 +-
 2 files changed, 2 insertions(+), 2 deletions(-)

diff --git a/PRIVATE b/PRIVATE
index 22a23a9d5..4cd6c7350 160000
--- a/PRIVATE
+++ b/PRIVATE
@@ -1 +1 @@
-Subproject commit 22a23a9d5939d49d9d277c7066d9b68003a33324
+Subproject commit 4cd6c7350b0a9d4ad3efcb5fe6c6cfffa99c426f
diff --git a/src/phase_damage.f90 b/src/phase_damage.f90
index df5e00575..af0a8fb7e 100644
--- a/src/phase_damage.f90
+++ b/src/phase_damage.f90
@@ -210,7 +210,7 @@ module function phase_f_phi(phi,co,ce) result(f)
   select case(phase_damage(ph))
     case(DAMAGE_ISOBRITTLE_ID,DAMAGE_ANISOBRITTLE_ID)
       f = 1.0_pReal &
-        - phi*damageState(ph)%state(1,en)
+        - 2.0_pREAL * phi*damageState(ph)%state(1,en)
     case default
       f = 0.0_pReal
   end select