Added dOxyGen comments for calc_isoneutral_slopes

src/core/MOM_isopycnal_slopes.F90

@@ -24,13 +24,14 @@ subroutine calc_isoneutral_slopes(G, GV, h, e, tv, dt_kappa_smooth, &
   real, dimension(SZI_(G),SZJ_(G),SZK_(G)),    intent(in)    :: h    !< Layer thicknesses, in H (usually m or kg m-2)
   real, dimension(SZI_(G),SZJ_(G),SZK_(G)+1),  intent(in)    :: e    !< Interface heights (m)
   type(thermo_var_ptrs),                       intent(in)    :: tv   !< A structure pointing to various thermodynamic variables
-  real,                                        intent(in)    :: dt_kappa_smooth
+  real,                                        intent(in)    :: dt_kappa_smooth !< A vertical diffusive smoothing timescale, in s.
   real, dimension(SZIB_(G),SZJ_(G),SZK_(G)+1), intent(inout) :: slope_x !< Isopycnal slope in i-direction (nondim)
   real, dimension(SZI_(G),SZJB_(G),SZK_(G)+1), intent(inout) :: slope_y !< Isopycnal slope in j-direction (nondim)
-  real, dimension(SZIB_(G),SZJ_(G),SZK_(G)+1), intent(inout) :: N2_u !< Brunt-Vaisala frequency squared at u-points (s-2)
-  real, dimension(SZI_(G),SZJB_(G),SZK_(G)+1), intent(inout) :: N2_v !< Brunt-Vaisala frequency squared at u-points (s-2)
-  optional                                                   :: N2_u, N2_v
-  integer, optional,                           intent(in)    :: halo !< Halo width over which to compute
+  real, dimension(SZIB_(G),SZJ_(G),SZK_(G)+1), &
+                                     optional, intent(inout) :: N2_u !< Brunt-Vaisala frequency squared at u-points (s-2)
+  real, dimension(SZI_(G),SZJB_(G),SZK_(G)+1), &
+                                     optional, intent(inout) :: N2_v !< Brunt-Vaisala frequency squared at u-points (s-2)
+  integer,                           optional, intent(in)    :: halo !< Halo width over which to compute
   ! Local variables
   real, dimension(SZI_(G), SZJ_(G), SZK_(G)) :: &
     T, &          ! The temperature (or density) in C, with the values in