diff --git a/gsw/stability.py b/gsw/stability.py
index 76189a2..c4acbb0 100644
--- a/gsw/stability.py
+++ b/gsw/stability.py
@@ -47,9 +47,10 @@ def Nsquared(SA, CT, p, lat=None, axis=0):
     Returns
     -------
     N2 : array
-        Buoyancy frequency-squared at pressure midpoints, 1/s.
+        Buoyancy frequency-squared at pressure midpoints, 1/s^2.
         The shape along the pressure axis dimension is one
         less than that of the inputs.
+        (Frequency N is in radians per second.)
     p_mid : array
         Pressure at midpoints of p, dbar.
         The array shape matches N2.