diff --git a/ndsl/constants.py b/ndsl/constants.py
index b7e91b18..0416c29b 100644
--- a/ndsl/constants.py
+++ b/ndsl/constants.py
@@ -78,86 +78,130 @@ class ConstantVersions(Enum):
 # Physical constants
 #####################
 if CONST_VERSION == ConstantVersions.GEOS:
-    RADIUS = Float(6.371e6)  # Radius of the Earth [m]
+    RADIUS = Float(6.371e6)
+    """Radius of the Earth [m]"""
     PI_8 = np.float64(3.14159265358979323846)
     PI = Float(PI_8)
-    OMEGA = Float(2.0) * PI / Float(86164.0)  # Rotation of the earth
-    GRAV = Float(9.80665)  # Acceleration due to gravity [m/s^2].04
-    RGRAV = Float(1.0) / GRAV  # Inverse of gravitational acceleration
-    RDGAS = Float(8314.47) / Float(
-        28.965
-    )  # Gas constant for dry air [J/kg/deg] ~287.04
-    RVGAS = Float(8314.47) / Float(18.015)  # Gas constant for water vapor [J/kg/deg]
-    HLV = Float(2.4665e6)  # Latent heat of evaporation [J/kg]
-    HLF = Float(3.3370e5)  # Latent heat of fusion [J/kg]  ~3.34e5
-    KAPPA = RDGAS / (Float(3.5) * RDGAS)  # Specific heat capacity of dry air at
+    OMEGA = Float(2.0) * PI / Float(86164.0)
+    """Rotation of the earth"""
+    GRAV = Float(9.80665)
+    """Acceleration due to gravity [m/s^2].04"""
+    RGRAV = Float(1.0) / GRAV
+    """Inverse of gravitational acceleration"""
+    RDGAS = Float(8314.47) / Float(28.965)
+    """Gas constant for dry air [J/kg/deg] ~287.04"""
+    RVGAS = Float(8314.47) / Float(18.015)
+    """Gas constant for water vapor [J/kg/deg]"""
+    HLV = Float(2.4665e6)
+    """Latent heat of evaporation [J/kg]"""
+    HLF = Float(3.3370e5)
+    """Latent heat of fusion [J/kg]  ~3.34e5"""
+    KAPPA = RDGAS / (Float(3.5) * RDGAS)
+    """Specific heat capacity of dry air at"""
     CP_AIR = RDGAS / KAPPA
-    TFREEZE = Float(273.16)  # Freezing temperature of fresh water [K]
+    TFREEZE = Float(273.16)
+    """Freezing temperature of fresh water [K]"""
     SAT_ADJUST_THRESHOLD = Float(1.0e-6)
 elif CONST_VERSION == ConstantVersions.UFS:
-    RADIUS = Float(6.3712e6)  # Radius of the Earth [m]
+    RADIUS = Float(6.3712e6)
+    """Radius of the Earth [m]"""
     PI = Float(3.1415926535897931)
-    OMEGA = Float(7.2921e-5)  # Rotation of the earth
-    GRAV = Float(9.80665)  # Acceleration due to gravity [m/s^2].04
-    RGRAV = Float(1.0 / GRAV)  # Inverse of gravitational acceleration
-    RDGAS = Float(287.05)  # Gas constant for dry air [J/kg/deg] ~287.04
-    RVGAS = Float(461.50)  # Gas constant for water vapor [J/kg/deg]
-    HLV = Float(2.5e6)  # Latent heat of evaporation [J/kg]
-    HLF = Float(3.3358e5)  # Latent heat of fusion [J/kg]  ~3.34e5
+    OMEGA = Float(7.2921e-5)
+    """Rotation of the earth"""
+    GRAV = Float(9.80665)
+    """Acceleration due to gravity [m/s^2].04"""
+    RGRAV = Float(1.0 / GRAV)
+    """Inverse of gravitational acceleration"""
+    RDGAS = Float(287.05)
+    """Gas constant for dry air [J/kg/deg] ~287.04"""
+    RVGAS = Float(461.50)
+    """Gas constant for water vapor [J/kg/deg]"""
+    HLV = Float(2.5e6)
+    """Latent heat of evaporation [J/kg]"""
+    HLF = Float(3.3358e5)
+    """Latent heat of fusion [J/kg]  ~3.34e5"""
     CP_AIR = Float(1004.6)
-    KAPPA = RDGAS / CP_AIR  # Specific heat capacity of dry air at
-    TFREEZE = Float(273.15)  # Freezing temperature of fresh water [K]
+    KAPPA = RDGAS / CP_AIR
+    """Specific heat capacity of dry air at"""
+    TFREEZE = Float(273.15)
+    """Freezing temperature of fresh water [K]"""
     SAT_ADJUST_THRESHOLD = Float(1.0e-8)
 elif CONST_VERSION == ConstantVersions.GFDL:
-    RADIUS = Float(6371.0e3)  # Radius of the Earth [m] #6371.0e3
-    PI = Float(3.14159265358979323846)  # 3.14159265358979323846
-    OMEGA = Float(7.292e-5)  # Rotation of the earth  # 7.292e-5
-    GRAV = Float(9.80)  # Acceleration due to gravity [m/s^2].04
-    RGRAV = Float(1.0) / GRAV  # Inverse of gravitational acceleration
-    RDGAS = Float(287.04)  # Gas constant for dry air [J/kg/deg] ~287.04
-    RVGAS = Float(461.50)  # Gas constant for water vapor [J/kg/deg]
-    HLV = Float(2.500e6)  # Latent heat of evaporation [J/kg]
-    HLF = Float(3.34e5)  # Latent heat of fusion [J/kg]  ~3.34e5
+    RADIUS = Float(6371.0e3)
+    """Radius of the Earth [m] #6371.0e3"""
+    PI = Float(3.14159265358979323846)
+    """3.14159265358979323846"""
+    OMEGA = Float(7.292e-5)
+    """Rotation of the earth  # 7.292e-5"""
+    GRAV = Float(9.80)
+    """Acceleration due to gravity [m/s^2].04"""
+    RGRAV = Float(1.0) / GRAV
+    """Inverse of gravitational acceleration"""
+    RDGAS = Float(287.04)
+    """Gas constant for dry air [J/kg/deg] ~287.04"""
+    RVGAS = Float(461.50)
+    """Gas constant for water vapor [J/kg/deg]"""
+    HLV = Float(2.500e6)
+    """Latent heat of evaporation [J/kg]"""
+    HLF = Float(3.34e5)
+    """Latent heat of fusion [J/kg]  ~3.34e5"""
     KAPPA = Float(2.0) / Float(7.0)
-    CP_AIR = RDGAS / KAPPA  # Specific heat capacity of dry air at
-    TFREEZE = Float(273.16)  # Freezing temperature of fresh water [K]
+    CP_AIR = RDGAS / KAPPA
+    """Specific heat capacity of dry air at"""
+    TFREEZE = Float(273.16)
+    """Freezing temperature of fresh water [K]"""
     SAT_ADJUST_THRESHOLD = Float(1.0e-8)
 else:
     raise RuntimeError("Constant selector failed, bad code.")
 
 SECONDS_PER_DAY = Float(86400.0)
 DZ_MIN = Float(2.0)
-CV_AIR = CP_AIR - RDGAS  # Heat capacity of dry air at constant volume
+CV_AIR = CP_AIR - RDGAS
+"""Heat capacity of dry air at constant volume"""
 RDG = -RDGAS / GRAV
 CNST_0P20 = Float(0.2)
 K1K = RDGAS / CV_AIR
 CNST_0P20 = Float(0.2)
-CV_VAP = Float(3.0) * RVGAS  # Heat capacity of water vapor at constant volume
-ZVIR = RVGAS / RDGAS - Float(1)  # con_fvirt in Fortran physics
-C_ICE = Float(1972.0)  # Heat capacity of ice at -15 degrees Celsius
-C_ICE_0 = Float(2106.0)  # Heat capacity of ice at 0 degrees Celsius
-C_LIQ = Float(4.1855e3)  # Heat capacity of water at 15 degrees Celsius
-CP_VAP = Float(4.0) * RVGAS  # Heat capacity of water vapor at constant pressure
-TICE = Float(273.16)  # Freezing temperature
-DC_ICE = C_LIQ - C_ICE  # Isobaric heating / cooling
-DC_VAP = CP_VAP - C_LIQ  # Isobaric heating / cooling
-D2ICE = DC_VAP + DC_ICE  # Isobaric heating / cooling
+CV_VAP = Float(3.0) * RVGAS
+"""Heat capacity of water vapor at constant volume"""
+ZVIR = RVGAS / RDGAS - Float(1)
+"""con_fvirt in Fortran physics"""
+C_ICE = Float(1972.0)
+"""Heat capacity of ice at -15 degrees Celsius"""
+C_ICE_0 = Float(2106.0)
+"""Heat capacity of ice at 0 degrees Celsius"""
+C_LIQ = Float(4.1855e3)
+"""Heat capacity of water at 15 degrees Celsius"""
+CP_VAP = Float(4.0) * RVGAS
+"""Heat capacity of water vapor at constant pressure"""
+TICE = Float(273.16)
+"""Freezing temperature"""
+DC_ICE = C_LIQ - C_ICE
+"""Isobaric heating / cooling"""
+DC_VAP = CP_VAP - C_LIQ
+"""Isobaric heating / cooling"""
+D2ICE = DC_VAP + DC_ICE
+"""Isobaric heating / cooling"""
 LI0 = HLF - DC_ICE * TICE
 EPS = RDGAS / RVGAS
 EPSM1 = EPS - Float(1.0)
-LV0 = (
-    HLV - DC_VAP * TICE
-)  # 3.13905782e6, evaporation latent heat coefficient at 0 degrees Kelvin
-LI00 = (
-    HLF - DC_ICE * TICE
-)  # -2.7105966e5, fusion latent heat coefficient at 0 degrees Kelvin
-LI2 = (
-    LV0 + LI00
-)  # 2.86799816e6, sublimation latent heat coefficient at 0 degrees Kelvin
-E00 = Float(611.21)  # Saturation vapor pressure at 0 degrees Celsius (Pa)
-PSAT = Float(610.78)  # Saturation vapor pressure at H2O 3pt (Pa)
-T_WFR = TICE - Float(40.0)  # homogeneous freezing temperature
+LV0 = HLV - DC_VAP * TICE
+"""3.13905782e6, evaporation latent heat coefficient at 0 degrees Kelvin"""
+LI00 = HLF - DC_ICE * TICE
+"""-2.7105966e5, fusion latent heat coefficient at 0 degrees Kelvin"""
+LI2 = LV0 + LI00
+"""2.86799816e6, sublimation latent heat coefficient at 0 degrees Kelvin"""
+E00 = Float(611.21)
+"""Saturation vapor pressure at 0 degrees Celsius (Pa)"""
+PSAT = Float(610.78)
+"""Saturation vapor pressure at H2O 3pt (Pa)"""
+T_WFR = TICE - Float(40.0)
+"""homogeneous freezing temperature"""
 TICE0 = TICE - Float(0.01)
-T_MIN = Float(178.0)  # Minimum temperature to freeze-dry all water vapor
+T_MIN = Float(178.0)
+"""Minimum temperature to freeze-dry all water vapor"""
 T_SAT_MIN = TICE - Float(160.0)
-LAT2 = np.power((HLV + HLF), 2, dtype=Float)  # used in bigg mechanism
+LAT2 = np.power((HLV + HLF), 2, dtype=Float)
+"""Used in bigg mechanism"""
+TTP = 2.7316e2
+"""Temperature of H2O triple point"""