Merge pull request #474 from danieljprice/tcourant

Signal speed used in tcour should include radiation pressure
danieljprice · Oct 23, 2023 · 420f5a6 · 420f5a6
2 parents cef5581 + 0d3f89f
commit 420f5a6
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