From ee15064f41254deae9cd942ac859b99fe97cfd6b Mon Sep 17 00:00:00 2001
From: Thomas Dauser <thomas.dauser@sternwarte.uni-erlangen.de>
Date: Wed, 15 Feb 2023 07:35:53 +0100
Subject: [PATCH] adds new function to calculate the grav radius

---
 src/slang/physics/gravitational_radius.sl | 19 +++++++++++++++++++
 1 file changed, 19 insertions(+)
 create mode 100644 src/slang/physics/gravitational_radius.sl

diff --git a/src/slang/physics/gravitational_radius.sl b/src/slang/physics/gravitational_radius.sl
new file mode 100644
index 00000000..0694ee28
--- /dev/null
+++ b/src/slang/physics/gravitational_radius.sl
@@ -0,0 +1,19 @@
+%!%+
+%\function{gravitational_radius}
+%\synopsis{calculates the gravitational radius in Meters defined as r_g=GM/c^,
+% for a mass given in units of M_sol.}
+%\usage{gravitational_radius(mass_in_solar)}
+%!%-
+define gravitational_radius() {
+   variable mass_in_solar;
+   
+   switch(_NARGS)
+   { case 1: mass_in_solar = ();}
+   { help(_function_name()); return; }
+
+   variable G = 6.673e-11;
+   variable c = 299792458.0;
+   variable Msol = 1.988e30;  % in kg
+   
+   return    G*Msol*mass_in_solar / (c^2);
+}
-- 
GitLab