Corrected index capitialization

  Corrected the capitalization of the indices in several arrays, to follow the
MOM6 convention of lowercase indices at tracer points and capital indices at
vorticity points.  All answers are bitwise identical.

src/core/MOM_barotropic.F90

@@ -1422,7 +1422,7 @@ subroutine btstep(U_in, V_in, eta_in, dt, bc_accel_u, bc_accel_v, &
                       (CS%ice_strength_length**2 * dtbt)
 
       ! Units of dyn_coef: m2 s-2 H-1
-      dyn_coef_eta(I,j) = min(dyn_coef_max, ice_strength * GV%H_to_m)
+      dyn_coef_eta(i,j) = min(dyn_coef_max, ice_strength * GV%H_to_m)
     enddo ; enddo ; endif
   endif
 
@@ -1627,7 +1627,7 @@ subroutine btstep(U_in, V_in, eta_in, dt, bc_accel_u, bc_accel_v, &
       if (CS%dynamic_psurf) then
 !GOMP do
         do j=jsv-1,jev+1 ; do i=isv-1,iev+1
-          p_surf_dyn(i,j) = dyn_coef_eta(I,j) * (eta_pred(i,j) - eta(i,j))
+          p_surf_dyn(i,j) = dyn_coef_eta(i,j) * (eta_pred(i,j) - eta(i,j))
         enddo ; enddo
       endif
     endif

src/parameterizations/vertical/MOM_set_viscosity.F90

@@ -1617,12 +1617,12 @@ subroutine set_viscous_ML(u, v, h, tv, forces, visc, dt, G, GV, CS, symmetrize)
 
     do_any_shelf = .false.
     if (associated(forces%frac_shelf_v)) then
-      do I=Is,Ie
+      do i=is,ie
         if (forces%frac_shelf_v(i,J)*G%mask2dCv(i,J) == 0.0) then
-          do_i(I) = .false.
+          do_i(i) = .false.
           visc%tbl_thick_shelf_v(i,J) = 0.0 ; visc%kv_tbl_shelf_v(i,J) = 0.0
         else
-          do_i(I) = .true. ; do_any_shelf = .true.
+          do_i(i) = .true. ; do_any_shelf = .true.
         endif
       enddo
     endif