From 09532f95a96161ef1ecdd0edef8b2c857766d907 Mon Sep 17 00:00:00 2001
From: Mike Lau <mike.lau@monash.edu>
Date: Thu, 12 Oct 2023 21:18:09 +1100
Subject: [PATCH] (radiation) signal speed used in tcour should include
 radiation pressure

---
 src/main/force.F90 | 4 ++++
 1 file changed, 4 insertions(+)

diff --git a/src/main/force.F90 b/src/main/force.F90
index 2e7eb3f58..bff01b4c0 100644
--- a/src/main/force.F90
+++ b/src/main/force.F90
@@ -2043,6 +2043,10 @@ subroutine get_stress(pri,spsoundi,rhoi,rho1i,xi,yi,zi, &
 !
     pro2i  = (pri + radPi)*rho1i*rho1i + stressiso
     vwavei = spsoundi
+
+    if (do_radiation) then
+       vwavei = sqrt(vwavei*vwavei + 4.*radPi/(3.*rhoi))  ! Commercon et al. (2011)
+    endif
  endif
 
 end subroutine get_stress