diff --git a/gsw/_wrapped_ufuncs.py b/gsw/_wrapped_ufuncs.py
index 2759e49..e9eeeac 100644
--- a/gsw/_wrapped_ufuncs.py
+++ b/gsw/_wrapped_ufuncs.py
@@ -5135,7 +5135,7 @@ def rho(SA, CT, p):
 
     Returns
     -------
-    rho : array-like, kg/m
+    rho : array-like, kg/m^3
         in-situ density
 
 
@@ -5194,7 +5194,7 @@ def rho_alpha_beta(SA, CT, p):
 
     Returns
     -------
-    rho : array-like, kg/m
+    rho : array-like, kg/m^3
         in-situ density
     alpha : array-like, 1/K
         thermal expansion coefficient
diff --git a/tools/fix_wrapped_ufunc_typos.py b/tools/fix_wrapped_ufunc_typos.py
index 1f195cd..f83b8de 100644
--- a/tools/fix_wrapped_ufunc_typos.py
+++ b/tools/fix_wrapped_ufunc_typos.py
@@ -36,6 +36,7 @@
     ('equlibrium', 'equilibrium'),
     (' apendix ', ' appendix '),
     (' slighty ', ' slightly '),
+    ('rho : array-like, kg/m', 'rho : array-like, kg/m^3'),
 ]
 
 with open(wrapmod) as f: