import numpy as np

__all__ == ['lsun', 'pc', 'lightspeed', 'jansky_mks', 'to_cgs', 'loge']


# Useful constants
lsun = 3.846e33
pc = 3.085677581467192e18  # in cm
lightspeed = 2.998e18  # AA/s
jansky_mks = 1e-26
# value to go from L_sun/AA to erg/s/cm^2/AA at 10pc
to_cgs = lsun/(4.0 * np.pi * (pc*10)**2)

# change base
loge = np.log10(np.e)

@@ -58,28 +58,82 @@
 ', '
-jansky_mks', 'to_cgs
+ckms',%0A            'jansky_mks', 'jansky_cgs',%0A            'to_cgs_at_10pc
 ', '
@@ -245,26 +245,68 @@
 A/s%0A
-jansky_mks = 1e-26
+ckms = 2.998e5  # km/s%0Ajansky_mks = 1e-26%0Ajansky_cgs = 1e-23
 %0A# v
@@ -365,16 +365,26 @@
 _cgs
+_at_10pc
  = lsun
-/
+ / 
 (4.0

clean up constants module