(radiation) signal speed used in tcour should include radiation pressure

danieljprice · Oct 12, 2023 · 09532f9 · 09532f9
1 parent 7b19d4c
commit 09532f9
Showing 1 changed file with 4 additions and 0 deletions.
diff --git 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