Removed spaces between subroutine and arguments

  Changed 'call sub (...)' to 'call sub(...)' in several places.  All answers
are bitwise identical, and only white-space is changed.

config_src/coupled_driver/MOM_surface_forcing.F90

@@ -1009,7 +1009,7 @@ subroutine surface_forcing_init(Time, G, param_file, diag, CS, restore_salt, res
 
   CS%diag => diag
 
-  call write_version_number (version)
+  call write_version_number(version)
   ! Read all relevant parameters and write them to the model log.
   call log_version(param_file, mdl, version, "")
 

src/core/MOM_PressureForce_analytic_FV.F90

@@ -240,9 +240,9 @@ subroutine PressureForce_AFV_nonBouss(h, tv, PFu, PFv, G, GV, CS, ALE_CSp, p_atm
   ! of freedeom needed to know the linear profile).
   if ( use_ALE ) then
     if ( CS%Recon_Scheme == 1 ) then
-      call pressure_gradient_plm (ALE_CSp, S_t, S_b, T_t, T_b, G, GV, tv, h, CS%boundary_extrap)
+      call pressure_gradient_plm(ALE_CSp, S_t, S_b, T_t, T_b, G, GV, tv, h, CS%boundary_extrap)
     elseif ( CS%Recon_Scheme == 2) then
-      call pressure_gradient_ppm (ALE_CSp, S_t, S_b, T_t, T_b, G, GV, tv, h, CS%boundary_extrap)
+      call pressure_gradient_ppm(ALE_CSp, S_t, S_b, T_t, T_b, G, GV, tv, h, CS%boundary_extrap)
     endif
   endif
 
@@ -253,7 +253,7 @@ subroutine PressureForce_AFV_nonBouss(h, tv, PFu, PFv, G, GV, CS, ALE_CSp, p_atm
     if (use_EOS) then
       if ( use_ALE ) then
         if ( CS%Recon_Scheme == 1 ) then
-          call int_spec_vol_dp_generic_plm ( T_t(:,:,k), T_b(:,:,k), &
+          call int_spec_vol_dp_generic_plm( T_t(:,:,k), T_b(:,:,k), &
                     S_t(:,:,k), S_b(:,:,k), p(:,:,K), p(:,:,K+1), &
                     alpha_ref, dp_neglect, p(:,:,nz+1), G%HI, &
                     tv%eqn_of_state, dza(:,:,k), intp_dza(:,:,k), &
@@ -665,14 +665,14 @@ subroutine PressureForce_AFV_Bouss(h, tv, PFu, PFv, G, GV, CS, ALE_CSp, p_atm, p
       ! where the layers are located.
       if ( use_ALE ) then
         if ( CS%Recon_Scheme == 1 ) then
-          call int_density_dz_generic_plm ( T_t(:,:,k), T_b(:,:,k), &
+          call int_density_dz_generic_plm( T_t(:,:,k), T_b(:,:,k), &
                     S_t(:,:,k), S_b(:,:,k), e(:,:,K), e(:,:,K+1), &
                     rho_ref, CS%Rho0, GV%g_Earth,    &
                     dz_neglect, G%bathyT, G%HI, G%HI, &
                     tv%eqn_of_state, dpa, intz_dpa, intx_dpa, inty_dpa, &
                     useMassWghtInterp = CS%useMassWghtInterp)
         elseif ( CS%Recon_Scheme == 2 ) then
-          call int_density_dz_generic_ppm ( tv%T(:,:,k), T_t(:,:,k), T_b(:,:,k), &
+          call int_density_dz_generic_ppm( tv%T(:,:,k), T_t(:,:,k), T_b(:,:,k), &
                     tv%S(:,:,k), S_t(:,:,k), S_b(:,:,k), e(:,:,K), e(:,:,K+1), &
                     rho_ref, CS%Rho0, GV%g_Earth, &
                     G%HI, G%HI, tv%eqn_of_state, dpa, intz_dpa,    &

src/core/MOM_PressureForce_blocked_AFV.F90

@@ -663,14 +663,14 @@ subroutine PressureForce_blk_AFV_Bouss(h, tv, PFu, PFv, G, GV, CS, ALE_CSp, p_at
         ! where the layers are located.
         if ( use_ALE ) then
           if ( CS%Recon_Scheme == 1 ) then
-            call int_density_dz_generic_plm ( T_t(:,:,k), T_b(:,:,k), &
+            call int_density_dz_generic_plm( T_t(:,:,k), T_b(:,:,k), &
                       S_t(:,:,k), S_b(:,:,k), e(:,:,K), e(:,:,K+1), &
                       rho_ref, CS%Rho0, GV%g_Earth,    &
                       dz_neglect, G%bathyT, G%HI, G%Block(n), &
                       tv%eqn_of_state, dpa_bk, intz_dpa_bk, intx_dpa_bk, inty_dpa_bk, &
                       useMassWghtInterp = CS%useMassWghtInterp)
           elseif ( CS%Recon_Scheme == 2 ) then
-            call int_density_dz_generic_ppm ( tv%T(:,:,k), T_t(:,:,k), T_b(:,:,k), &
+            call int_density_dz_generic_ppm( tv%T(:,:,k), T_t(:,:,k), T_b(:,:,k), &
                       tv%S(:,:,k), S_t(:,:,k), S_b(:,:,k), e(:,:,K), e(:,:,K+1), &
                       rho_ref, CS%Rho0, GV%g_Earth, &
                       G%HI, G%Block(n), tv%eqn_of_state, dpa_bk, intz_dpa_bk,    &

src/diagnostics/MOM_sum_output.F90

@@ -223,7 +223,7 @@ subroutine MOM_sum_output_init(G, param_file, directory, ntrnc, &
                  "The maximum velocity allowed before the velocity \n"//&
                  "components are truncated.", units="m s-1", default=3.0e8)
     CS%max_Energy = 10.0 * maxvel**2
-    call log_param (param_file, mdl, "MAX_ENERGY as used", CS%max_Energy)
+    call log_param(param_file, mdl, "MAX_ENERGY as used", CS%max_Energy)
   endif
 
   call get_param(param_file, mdl, "ENERGYFILE", energyfile, &

src/equation_of_state/MOM_EOS.F90

@@ -1680,10 +1680,10 @@ subroutine int_density_dz_generic_ppm (T, T_t, T_b, S, S_t, S_b, &
     S_node(9) = 0.5 * ( S_node(6) + S_node(8) )
     S_node(7) = 0.5 * ( S_node(3) + S_node(4) )
 
-    call calculate_density ( T_node, S_node, p_node, r_node, 1, 9, EOS )
+    call calculate_density( T_node, S_node, p_node, r_node, 1, 9, EOS )
     r_node = r_node - rho_ref
 
-    call compute_integral_quadratic ( x, y, r_node, intx_dpa(i-ioff,j-joff) )
+    call compute_integral_quadratic( x, y, r_node, intx_dpa(i-ioff,j-joff) )
 
     intx_dpa(i-ioff,j-joff) = intx_dpa(i-ioff,j-joff) * G_e
 
@@ -1755,8 +1755,8 @@ subroutine compute_integral_quadratic ( x, y, f, integral )
   do k = 1,9
 
     ! Evaluate shape functions and gradients for isomorphism
-    call evaluate_shape_bilinear ( xi(k), eta(k), phiiso, &
-                                   dphiisodxi, dphiisodeta )
+    call evaluate_shape_bilinear( xi(k), eta(k), phiiso, &
+                                  dphiisodxi, dphiisodeta )
 
     ! Determine gradient of global coordinate at integration point
     dxdxi  = 0.0
@@ -1775,7 +1775,7 @@ subroutine compute_integral_quadratic ( x, y, f, integral )
     jacobian_k = dxdxi*dydeta - dydxi*dxdeta
 
     ! Evaluate shape functions for interpolation
-    call evaluate_shape_quadratic ( xi(k), eta(k), phi, dphidxi, dphideta )
+    call evaluate_shape_quadratic( xi(k), eta(k), phi, dphidxi, dphideta )
 
     ! Evaluate function at integration point
     f_k = 0.0

src/framework/MOM_horizontal_regridding.F90

@@ -497,7 +497,7 @@ subroutine horiz_interp_and_extrap_tracer_record(filename, varnam,  conversion,
 
     call mpp_sync()
     call mpp_broadcast(tr_inp,id*jdp,root_PE())
-    call mpp_sync_self ()
+    call mpp_sync_self()
 
     mask_in=0.0
 
@@ -770,7 +770,7 @@ subroutine horiz_interp_and_extrap_tracer_fms_id(fms_id,  Time, conversion, G, t
 
     call mpp_sync()
     call mpp_broadcast(tr_inp,id*jdp,root_PE())
-    call mpp_sync_self ()
+    call mpp_sync_self()
 
     mask_in=0.0
 

src/tracer/MOM_offline_aux.F90

@@ -578,7 +578,7 @@ subroutine offline_add_diurnal_SW(fluxes, G, Time_start, Time_end)
 
     call diurnal_solar(G%geoLatT(i,j)*rad, G%geoLonT(i,j)*rad, Time_start, cosz=cosz_dt, &
                        fracday=fracday_dt, rrsun=rrsun_dt, dt_time=dt_here)
-    call daily_mean_solar (G%geoLatT(i,j)*rad, time_since_ae, cosz_day, fracday_day, rrsun_day)
+    call daily_mean_solar(G%geoLatT(i,j)*rad, time_since_ae, cosz_day, fracday_day, rrsun_day)
     diurnal_factor = cosz_dt*fracday_dt*rrsun_dt / &
                      max(1e-30, cosz_day*fracday_day*rrsun_day)