diff --git a/docs/src/UniversalFunctions.md b/docs/src/UniversalFunctions.md
index 5db34e3d..354d3b2e 100644
--- a/docs/src/UniversalFunctions.md
+++ b/docs/src/UniversalFunctions.md
@@ -1,6 +1,6 @@
 # Universal Functions
 
-UniversalFunctions.jl provides universal functions for SurfaceFluxes.jl. Here, we reproduce some plots from literature, specifically from Gryanik et al. 2020, and Businger.
+UniversalFunctions.jl provides universal functions for SurfaceFluxes.jl. The functions are defined in `src/UniversalFunctions.jl` and the plots are generated in `docs/src/plot_universal_functions.jl`. Here, we reproduce some plots from literature, specifically from Gryanik et al. 2020, Businger 1971, and Bonan 2019. Note that Bonan uses the forms of $\phi$ and $\psi$ from Dyer and Hicks 1970; Dyer 1974; Brutsaert 1982, pp. 68–71; Garratt 1992, pp. 52–54. 
 
 ```@example
 include("plot_universal_functions.jl")
@@ -24,3 +24,10 @@ include("plot_universal_functions.jl")
 ![](Businger_phi_h.svg)
 ![](Businger_phi_m.svg)
 
+
+## Figs 1,2 (Bonan)
+
+![](Bonan_phi_h.svg)
+![](Bonan_phi_m.svg)
+![](Bonan_psi_h.svg)
+![](Bonan_psi_m.svg)
diff --git a/docs/src/plot_universal_functions.jl b/docs/src/plot_universal_functions.jl
index 85e535d8..5173898a 100644
--- a/docs/src/plot_universal_functions.jl
+++ b/docs/src/plot_universal_functions.jl
@@ -59,10 +59,10 @@ function save_ψ_figs(ufts, ζ; ylims = nothing, fig_prefix = "", xaxis = :ident
 end
 
 
+# Gryanik Plots
 save_ϕ_figs(ufts, FT(0):FT(0.01):FT(15); ylims = (0, 30), fig_prefix = "Gryanik12")
 save_ψ_figs(ufts, FT(0):FT(0.01):FT(15); ylims = (-25, 0), fig_prefix = "Gryanik12")
 
-
 save_ϕ_figs(
     ufts,
     10 .^ (FT(-3):0.1:FT(2));
@@ -73,4 +73,9 @@ save_ϕ_figs(
 )
 
 
+# Businger Plots
 save_ϕ_figs(ufts, FT(-2.5):FT(0.01):FT(2); ylims = (-1, 8), fig_prefix = "Businger")
+
+# Bonan Plots
+save_ϕ_figs(ufts, FT(-2):FT(0.01):FT(1); ylims = (0, 6), fig_prefix = "Bonan")
+save_ψ_figs(ufts, FT(-2):FT(0.01):FT(1); ylims = (-5, 4), fig_prefix = "Bonan")