aboutsummaryrefslogtreecommitdiffstats
diff options
context:
space:
mode:
authorVolker Hoffmann <volker@cheleb.net>2015-08-17 11:13:18 +0200
committerVolker Hoffmann <volker@cheleb.net>2015-08-17 11:13:18 +0200
commit70137eb4c78662d8c51042348fe3e7acf2b046c5 (patch)
treefbde29701ee2f012b0b8e945ea7bd96ab3b142b4
parent486d7a54b8a919908f832e6c98270d8c86277b25 (diff)
docs: add units
-rw-r--r--Helpers/impact_helpers.py4
1 files changed, 2 insertions, 2 deletions
diff --git a/Helpers/impact_helpers.py b/Helpers/impact_helpers.py
index 69a1cde..cc90715 100644
--- a/Helpers/impact_helpers.py
+++ b/Helpers/impact_helpers.py
@@ -144,7 +144,7 @@ def impactor_mass(D, v_i, theta, rho_t=3.34, rho_i=2.00, grav=1.622):
# Fix Units
grav /= 1000.0 # m/s2 => km/s2
- # Compute Impactor Diameter
+ # Compute Impactor Diameter (km)
L = D
L /= 1.161 * (rho_i/rho_t)**(1.0/3.0) * \
v_i**0.44 * \
@@ -152,7 +152,7 @@ def impactor_mass(D, v_i, theta, rho_t=3.34, rho_i=2.00, grav=1.622):
np.sin(theta)**(1.0/3.0)
L = L**(1.0/0.78)
- # Compute (Spherical) Impactor Mass
+ # Compute (Spherical) Impactor Mass (kg)
M = 4.0/3.0 * np.pi * (L/2.0)**3.0 * rho_i/1000.0*(1000.0*100.0)**(3.0)
# Return [Float or Numpy Float Array; Same as D]