diff --git a/docs/zotero.bib b/docs/zotero.bib
index fe5dd37909..abaeea4c72 100644
--- a/docs/zotero.bib
+++ b/docs/zotero.bib
@@ -2873,8 +2873,9 @@ @article{Blewitt2003
 	month=feb
 }
 
-@article{Wang2012,
-	author={Wang, H., Xiang, L., Jia, L., Jiang, L., Wang, Z., Hu, B. and Gao, P.},
+@article{Wang2012-2,
+	author={Wang, H. and Xiang, L. and Jia, L. and Jiang, L. and Wang, Z. and Hu, B.
+	and Gao, P.},
 	year={2012},
 	title={Load Love numbers and Green's functions
 for elastic Earth models PREM, iasp91, ak135, and modified models with refined crustal structure from Crust 2.0},
@@ -2883,20 +2884,6 @@ @article{Wang2012
 	pages={190--199}
 }
 
-@article{Jansen2015,
-	title={Parameterization of eddy fluxes based on a mesoscale energy budget},
-	volume={92},
-	ISSN={1463-5003},
-	url={http://dx.doi.org/10.1016/j.ocemod.2015.05.007},
-	DOI={10.1016/j.ocemod.2015.05.007},
-	journal={Ocean Modelling},
-	publisher={Elsevier BV},
-	author={Jansen, Malte F. and Adcroft, Alistair J. and Hallberg, Robert and Held, Isaac M.},
-	year={2015},
-	month=aug,
-	pages={28–-41}
-}
-
 @inproceedings{Hallberg2003,
 	title={The ability of large-scale ocean models to accept parameterizations of boundary mixing, and a description of a refined bulk mixed-layer model},
 	author={Robert Hallberg},
diff --git a/src/parameterizations/lateral/MOM_MEKE.F90 b/src/parameterizations/lateral/MOM_MEKE.F90
index f3bd9e8934..c3859c3bd1 100644
--- a/src/parameterizations/lateral/MOM_MEKE.F90
+++ b/src/parameterizations/lateral/MOM_MEKE.F90
@@ -1892,7 +1892,7 @@ end subroutine MEKE_end
 !! \f$ U_b \f$ is a constant background bottom velocity scale and is
 !! typically not used (i.e. set to zero).
 !!
-!! Following Jansen et al., 2015, the projection of eddy energy on to the bottom
+!! Following \cite jansen2015, the projection of eddy energy on to the bottom
 !! is given by the ratio of bottom energy to column mean energy:
 !! \f[
 !! \gamma_b^2  = \frac{E_b}{E} = \gamma_{d0}
@@ -1924,12 +1924,12 @@ end subroutine MEKE_end
 !! \f[  \kappa_M = \gamma_\kappa \sqrt{ \gamma_t^2 U_e^2 A_\Delta } \f]
 !!
 !! where \f$ A_\Delta \f$ is the area of the grid cell.
-!! Following Jansen et al., 2015, we now use
+!! Following \cite jansen2015, we now use
 !!
 !! \f[  \kappa_M = \gamma_\kappa l_M \sqrt{ \gamma_t^2 U_e^2 } \f]
 !!
 !! where \f$ \gamma_\kappa \in [0,1] \f$ is a non-dimensional factor and,
-!! following Jansen et al., 2015, \f$\gamma_t^2\f$ is the ratio of barotropic
+!! following \cite jansen2015, \f$\gamma_t^2\f$ is the ratio of barotropic
 !! eddy energy to column mean eddy energy given by
 !! \f[
 !! \gamma_t^2  = \frac{E_t}{E} = \left( 1 + c_{t} \frac{L_d}{L_f} \right)^{-\frac{1}{4}}
diff --git a/src/parameterizations/lateral/MOM_hor_visc.F90 b/src/parameterizations/lateral/MOM_hor_visc.F90
index ea06795654..4a901beab0 100644
--- a/src/parameterizations/lateral/MOM_hor_visc.F90
+++ b/src/parameterizations/lateral/MOM_hor_visc.F90
@@ -3195,9 +3195,9 @@ end subroutine hor_visc_end
 !!
 !! \section section_horizontal_viscosity Horizontal viscosity in MOM
 !!
-!! This module contains the subroutine horizontal_viscosity() that calculates the
+!! This module contains the subroutine horizontal_viscosity that calculates the
 !! effects of horizontal viscosity, including parameterizations of the value of
-!! the viscosity itself. horizontal_viscosity() calculates the acceleration due to
+!! the viscosity itself. Subroutine horizontal_viscosity calculates the acceleration due to
 !! some combination of a biharmonic viscosity and a Laplacian viscosity. Either or
 !! both may use a coefficient that depends on the shear and strain of the flow.
 !! All metric terms are retained. The Laplacian is calculated as the divergence of
diff --git a/src/parameterizations/lateral/MOM_load_love_numbers.F90 b/src/parameterizations/lateral/MOM_load_love_numbers.F90
index c8bd8bbc3c..8faf3aafab 100644
--- a/src/parameterizations/lateral/MOM_load_love_numbers.F90
+++ b/src/parameterizations/lateral/MOM_load_love_numbers.F90
@@ -1461,7 +1461,7 @@ module MOM_load_love_numbers
 !! Variable Love_Data stores the Love numbers up to degree 1440. From left to right: degree, h, l, and k. Data in this
 !! module is imported from SAL calculation in Model for Prediction Across Scales (MPAS)-Ocean developed by Los Alamos
 !! National Laboratory and University of Michigan [\cite Barton2022 and \cite Brus2023]. The load Love numbers
-!! are from \cite Wang2012, which are in the center of mass of total Earth system reference frame (CM). When used,
+!! are from \cite Wang2012-2, which are in the center of mass of total Earth system reference frame (CM). When used,
 !! Love numbers with degree<2 should be converted to center of mass solid Earth reference frame (CF)
 !! [\cite Blewitt2003], as in subroutine calc_love_scaling in MOM_tidal_forcing module.
 !!