diff --git a/hack/lundquist_number.py b/hack/lundquist_number.py
new file mode 100644
index 0000000..acabd4f
--- /dev/null
+++ b/hack/lundquist_number.py
@@ -0,0 +1,13 @@
+import astropy.units as u
+
+from plasmapy.formulary import dimensionless, speeds
+
+
+def Lundquist_Num(
+        L: u.m,
+        B: u.T,
+        density: (u.m**-3, u.kg / u.m**3),
+        sigma: u.S / u.m
+) -> u.dimensionless_unscaled:
+    alfven = speeds.Alfven_speed(B, density)
+    return dimensionless.Mag_Reynolds(alfven, L, sigma)