add magneto disk

YihanWangAstro · Oct 2, 2023 · 4c52c10 · 4c52c10
1 parent fb28793
commit 4c52c10
Show file tree

Hide file tree

Showing 4 changed files with 626 additions and 4 deletions.
diff --git a/src/interaction/disk-capture.hpp b/src/interaction/disk-capture.hpp
@@ -87,7 +87,7 @@ namespace hub::force {
             auto vabs = sqrt(v2);
             auto cs = Omega * H;
             auto Mach = vabs / cs;
-            auto logR = 3.0;
+            auto logR = 8.0;
             // std::max(log(m[0] / m[i]), 3.0);
 
             double I = std::max(logR + offset, 1.0);