Change "end if" to endif and "end do" to enddo

  Replaced "end if" with "endif" in 35 places in 10 files and "end do" with
"enddo" in 12 places in 4 files to comply with the MOM6 style guide.  In
addition, in 23 places in 7 files, code like `if(test)` were replaced with
`if (test)`, also to align with the MOM6 style guide.  All answers are
bitwise identical.

Author: Hallberg-NOAA

config_src/drivers/STALE_mct_cap/mom_ocean_model_mct.F90

@@ -585,7 +585,7 @@ subroutine update_ocean_model(Ice_ocean_boundary, OS, Ocean_sfc, &
   call disable_averaging(OS%diag)
   Master_time = OS%Time ; Time1 = OS%Time
 
-  if(OS%offline_tracer_mode) then
+  if (OS%offline_tracer_mode) then
     call step_offline(OS%forces, OS%fluxes, OS%sfc_state, Time1, dt_coupling, OS%MOM_CSp)
 
   elseif ((.not.do_thermo) .or. (.not.do_dyn)) then
@@ -808,7 +808,7 @@ subroutine initialize_ocean_public_type(input_domain, Ocean_sfc, diag, maskmap,
 
   call mpp_get_layout(input_domain,layout)
   call mpp_get_global_domain(input_domain, xsize=xsz, ysize=ysz)
-  if(PRESENT(maskmap)) then
+  if (PRESENT(maskmap)) then
     call mpp_define_domains((/1,xsz,1,ysz/),layout,Ocean_sfc%Domain, maskmap=maskmap)
   else
     call mpp_define_domains((/1,xsz,1,ysz/),layout,Ocean_sfc%Domain)

config_src/drivers/STALE_mct_cap/mom_surface_forcing_mct.F90

@@ -439,14 +439,14 @@ subroutine convert_IOB_to_fluxes(IOB, fluxes, index_bounds, Time, valid_time, G,
       fluxes%lrunoff(i,j) = kg_m2_s_conversion * IOB%rofl_flux(i-i0,j-j0) * G%mask2dT(i,j)
     else if (associated(IOB%runoff)) then
       fluxes%lrunoff(i,j) = kg_m2_s_conversion * IOB%runoff(i-i0,j-j0) * G%mask2dT(i,j)
-    end if
+    endif
 
     ! ice runoff flux
     if (associated(IOB%rofi_flux)) then
       fluxes%frunoff(i,j) = kg_m2_s_conversion * IOB%rofi_flux(i-i0,j-j0) * G%mask2dT(i,j)
     else if (associated(IOB%calving)) then
       fluxes%frunoff(i,j) = kg_m2_s_conversion * IOB%calving(i-i0,j-j0) * G%mask2dT(i,j)
-    end if
+    endif
 
     if (associated(IOB%ustar_berg)) &
       fluxes%ustar_berg(i,j) = US%m_to_Z*US%T_to_s * IOB%ustar_berg(i-i0,j-j0) * G%mask2dT(i,j)

config_src/drivers/STALE_mct_cap/ocn_comp_mct.F90

@@ -201,7 +201,7 @@ subroutine ocn_init_mct( EClock, cdata_o, x2o_o, o2x_o, NLFilename )
 
     !  set the shr log io unit number
     call shr_file_setLogUnit(stdout)
-  end if
+  endif
 
   call set_calendar_type(NOLEAP)  !TODO: confirm this
 
@@ -316,13 +316,13 @@ subroutine ocn_init_mct( EClock, cdata_o, x2o_o, o2x_o, NLFilename )
     close(nu)
     if (is_root_pe()) then
       write(stdout,*) 'Reading restart file(s): ',trim(restartfiles)
-    end if
+    endif
     call shr_file_freeUnit(nu)
     call ocean_model_init(glb%ocn_public, glb%ocn_state, time0, time_start, input_restart_file=trim(restartfiles))
   endif
   if (is_root_pe()) then
     write(stdout,'(/12x,a/)') '======== COMPLETED MOM INITIALIZATION ========'
-  end if
+  endif
 
   ! Initialize ocn_state%sfc_state out of sight
   call ocean_model_init_sfc(glb%ocn_state, glb%ocn_public)
@@ -384,14 +384,14 @@ subroutine ocn_init_mct( EClock, cdata_o, x2o_o, o2x_o, NLFilename )
   if (mom_cpl_dt /= ocn_cpl_dt) then
     write(stdout,*) 'ERROR mom_cpl_dt and ocn_cpl_dt must be identical'
     call exit(0)
-  end if
+  endif
 
   ! send initial state to driver
 
   !TODO:
   ! if ( lsend_precip_fact )  then
   !    call seq_infodata_PutData( infodata, precip_fact=precip_fact)
-  ! end if
+  ! endif
 
   if (debug .and. root_pe().eq.pe_here()) print *, "calling ocn_export"
   call ocn_export(glb%ind, glb%ocn_public, glb%grid, o2x_o%rattr, mom_cpl_dt, ncouple_per_day)
@@ -412,7 +412,7 @@ subroutine ocn_init_mct( EClock, cdata_o, x2o_o, o2x_o, NLFilename )
   if (is_root_pe()) then
     call shr_file_setLogUnit (shrlogunit)
     call shr_file_setLogLevel(shrloglev)
-  end if
+  endif
 
 end subroutine ocn_init_mct
 
@@ -488,10 +488,10 @@ subroutine ocn_run_mct( EClock, cdata_o, x2o_o, o2x_o)
       time_start = time_start-coupling_timestep
       ! double the first coupling interval (to account for the missing coupling interval to due to lag)
       coupling_timestep = coupling_timestep*2
-    end if
+    endif
 
     firstCall = .false.
-  end if
+  endif
 
   ! Debugging clocks
   if (debug .and. is_root_pe()) then
@@ -526,7 +526,7 @@ subroutine ocn_run_mct( EClock, cdata_o, x2o_o, o2x_o)
           c1=glb%c1, c2=glb%c2, c3=glb%c3, c4=glb%c4)
   else
     call ocn_import(x2o_o%rattr, glb%ind,  glb%grid, Ice_ocean_boundary, glb%ocn_public, stdout, Eclock )
-  end if
+  endif
 
   ! Update internal ocean
   call update_ocean_model(ice_ocean_boundary, glb%ocn_state, glb%ocn_public, time_start, coupling_timestep)
@@ -760,7 +760,7 @@ end subroutine ocn_domain_mct
   else
     write(stdout,*) 'ocn_comp_mct ERROR: unknown starttype'
     call exit(0)
-  end if
+  endif
   return
 
 end function

config_src/drivers/nuopc_cap/mom_cap.F90

@@ -431,7 +431,7 @@ subroutine InitializeP0(gcomp, importState, exportState, clock, rc)
   if (ChkErr(rc,__LINE__,u_FILE_u)) return
   if (isPresent .and. isSet) then
      if (trim(value) .eq. '.true.') restart_eor = .true.
-  end if
+  endif
 
   if (localPet == 0) call cap_profiling("mom", "InitializeP0", "E")
 
@@ -506,7 +506,7 @@ subroutine InitializeAdvertise(gcomp, importState, exportState, clock, rc)
 
   if (localPet == 0) call cap_profiling("mom", "InitializeAdvertise", "B")
 
-  if(write_runtimelog) timeiads = MPI_Wtime()
+  if (write_runtimelog) timeiads = MPI_Wtime()
 
   call ESMF_LogWrite(subname//' enter', ESMF_LOGMSG_INFO)
 
@@ -571,7 +571,7 @@ subroutine InitializeAdvertise(gcomp, importState, exportState, clock, rc)
 
       if (cesm_coupled) then
         ! Multiinstance logfile name needs a correction
-        if(len_trim(inst_suffix) > 0) then
+        if (len_trim(inst_suffix) > 0) then
           n = index(logfile, '.')
           logfile = logfile(1:n-1)//trim(inst_suffix)//logfile(n:)
         endif
@@ -660,7 +660,7 @@ subroutine InitializeAdvertise(gcomp, importState, exportState, clock, rc)
     if (is_root_pe()) then
         write(stdout,*) 'ice_ncat = ', Ice_ocean_boundary%ice_ncat
     endif
-  end if
+  endif
 
   if (is_root_pe()) then
     write(stdout,*) subname//'start time: y,m,d-',year,month,day,'h,m,s=',hour,minute,second
@@ -945,7 +945,7 @@ subroutine InitializeAdvertise(gcomp, importState, exportState, clock, rc)
     call NUOPC_Advertise(exportState, standardName=fldsFrOcn(n)%stdname, name=fldsFrOcn(n)%shortname, rc=rc)
     if (ChkErr(rc,__LINE__,u_FILE_u)) return
   enddo
-  if(write_runtimelog .and. is_root_pe()) write(stdout,*) 'In ',trim(subname),' time ', MPI_Wtime()-timeiads
+  if (write_runtimelog .and. is_root_pe()) write(stdout,*) 'In ',trim(subname),' time ', MPI_Wtime()-timeiads
 
   if (localPet == 0) call cap_profiling("mom", "InitializeAdvertise", "E")
 
@@ -1044,7 +1044,7 @@ subroutine InitializeRealize(gcomp, importState, exportState, clock, rc)
 
   if (localPet == 0) call cap_profiling("mom", "InitializeRealize", "B")
 
-  if(write_runtimelog) timeirls = MPI_Wtime()
+  if (write_runtimelog) timeirls = MPI_Wtime()
 
   call shr_log_setLogUnit (stdout)
 
@@ -1253,7 +1253,7 @@ subroutine InitializeRealize(gcomp, importState, exportState, clock, rc)
     do n = 1,numOwnedElements
       lonMesh(n) = ownedElemCoords(2*n-1)
       latMesh(n) = ownedElemCoords(2*n)
-    end do
+    enddo
 
     elemMaskArray = ESMF_ArrayCreate(Distgrid, maskMesh, rc=rc)
     if (ChkErr(rc,__LINE__,u_FILE_u)) return
@@ -1270,8 +1270,8 @@ subroutine InitializeRealize(gcomp, importState, exportState, clock, rc)
         mask(n) = ocean_grid%mask2dT(ig,jg)
         lon(n)  = ocean_grid%geolonT(ig,jg)
         lat(n)  = ocean_grid%geolatT(ig,jg)
-      end do
-    end do
+      enddo
+    enddo
 
     eps_omesh = get_eps_omesh(ocean_state)
     do n = 1,lsize
@@ -1297,7 +1297,7 @@ subroutine InitializeRealize(gcomp, importState, exportState, clock, rc)
         write(err_msg, frmt)n,maskMesh(n),mask(n)
         call MOM_error(FATAL, err_msg)
       endif
-    end do
+    enddo
 
     ! realize the import and export fields using the mesh
     call MOM_RealizeFields(importState, fldsToOcn_num, fldsToOcn, "Ocn import", &
@@ -1342,8 +1342,8 @@ subroutine InitializeRealize(gcomp, importState, exportState, clock, rc)
           mod2med_areacor(k) = model_areas(k) / mesh_areas(k)
           med2mod_areacor(k) = mesh_areas(k) / model_areas(k)
         endif
-      end do
-    end do
+      enddo
+    enddo
     deallocate(mesh_areas)
     deallocate(model_areas)
 
@@ -1621,7 +1621,7 @@ subroutine InitializeRealize(gcomp, importState, exportState, clock, rc)
   if (use_cdeps_inline) then
      call mom_inline_init(gcomp, clock, eMesh, localPet, rc=rc)
      if (ChkErr(rc,__LINE__,u_FILE_u)) return
-  end if
+  endif
 
   !---------------------------------
   ! write out diagnostics
@@ -1632,7 +1632,7 @@ subroutine InitializeRealize(gcomp, importState, exportState, clock, rc)
   !if (ChkErr(rc,__LINE__,u_FILE_u)) return
 
   timere = 0.
-  if(write_runtimelog .and. is_root_pe()) write(stdout,*) 'In ',trim(subname),' time ', MPI_Wtime()-timeirls
+  if (write_runtimelog .and. is_root_pe()) write(stdout,*) 'In ',trim(subname),' time ', MPI_Wtime()-timeirls
 
   if (localPet == 0) call cap_profiling("mom", "InitializeRealize", "E")
 
@@ -1669,7 +1669,7 @@ subroutine DataInitialize(gcomp, rc)
 
   if (localPet == 0) call cap_profiling("mom", "DataInitialize", "B")
 
-  if(write_runtimelog) timedis = MPI_Wtime()
+  if (write_runtimelog) timedis = MPI_Wtime()
 
   ! query the Component for its clock, importState and exportState
   call ESMF_GridCompGet(gcomp, clock=clock, importState=importState, exportState=exportState, rc=rc)
@@ -1731,7 +1731,7 @@ subroutine DataInitialize(gcomp, rc)
     enddo
   endif
 
-  if(write_runtimelog .and. is_root_pe()) write(stdout,*) 'In ',trim(subname),' time ', MPI_Wtime()-timedis
+  if (write_runtimelog .and. is_root_pe()) write(stdout,*) 'In ',trim(subname),' time ', MPI_Wtime()-timedis
 
   if (localPet == 0) call cap_profiling("mom", "DataInitialize", "E")
 
@@ -1793,10 +1793,10 @@ subroutine ModelAdvance(gcomp, rc)
 
   if (localPet == 0) call cap_profiling("mom", "ModelAdvance", "B")
 
-  if(profile_memory) call ESMF_VMLogMemInfo("Entering MOM Model_ADVANCE: ")
-  if(write_runtimelog) then
+  if (profile_memory) call ESMF_VMLogMemInfo("Entering MOM Model_ADVANCE: ")
+  if (write_runtimelog) then
      timers = MPI_Wtime()
-     if(timere>0. .and. is_root_pe()) write(stdout,*) 'In ',trim(subname),' time since last time step ',timers-timere
+     if (timere>0. .and. is_root_pe()) write(stdout,*) 'In ',trim(subname),' time since last time step ',timers-timere
   endif
 
   call shr_log_setLogUnit (stdout)
@@ -1911,7 +1911,7 @@ subroutine ModelAdvance(gcomp, rc)
     if (use_cdeps_inline) then
       call mom_inline_run(clock, ocean_public, ocean_grid, ice_ocean_boundary, dbug, rc=rc)
       if (ChkErr(rc,__LINE__,u_FILE_u)) return
-    end if
+    endif
 
     !---------------
     ! Update MOM6
@@ -1956,7 +1956,7 @@ subroutine ModelAdvance(gcomp, rc)
       ! turn off the alarm
       call ESMF_AlarmRingerOff(restart_alarm, rc=rc )
       if (ChkErr(rc,__LINE__,u_FILE_u)) return
-    end if
+    endif
 
     write_restart_eor = .false.
     if (restart_eor) then
@@ -1966,8 +1966,8 @@ subroutine ModelAdvance(gcomp, rc)
         ! turn off the alarm
         call ESMF_AlarmRingerOff(stop_alarm, rc=rc )
         if (ChkErr(rc,__LINE__,u_FILE_u)) return
-      end if
-    end if
+      endif
+    endif
 
 #ifndef CESMCOUPLED
     call is_restart_fh(clock, restartfh_info, write_restartfh)
@@ -2071,12 +2071,12 @@ subroutine ModelAdvance(gcomp, rc)
     enddo
   endif
 
-  if(write_runtimelog) then
+  if (write_runtimelog) then
     timere = MPI_Wtime()
-    if(is_root_pe()) write(stdout,*) 'In ',trim(subname),' time ', timere-timers
+    if (is_root_pe()) write(stdout,*) 'In ',trim(subname),' time ', timere-timers
   endif
 
-  if(profile_memory) call ESMF_VMLogMemInfo("Leaving MOM Model_ADVANCE: ")
+  if (profile_memory) call ESMF_VMLogMemInfo("Leaving MOM Model_ADVANCE: ")
 
   if (localPet == 0) call cap_profiling("mom", "ModelAdvance", "E")
 
@@ -2299,7 +2299,7 @@ subroutine ocean_model_finalize(gcomp, rc)
     write(stdout,*) 'MOM: --- finalize called ---'
   endif
   rc = ESMF_SUCCESS
-  if(write_runtimelog) timefs = MPI_Wtime()
+  if (write_runtimelog) timefs = MPI_Wtime()
 
   call ESMF_GridCompGetInternalState(gcomp, ocean_internalstate, rc)
   if (ChkErr(rc,__LINE__,u_FILE_u)) return
@@ -2334,7 +2334,7 @@ subroutine ocean_model_finalize(gcomp, rc)
   call outputlog_run(clock, .true., rc=rc)
   if (ChkErr(rc,__LINE__,u_FILE_u)) return
 
-  if(write_runtimelog .and. is_root_pe()) write(stdout,*) 'In ',trim(subname),' time ', MPI_Wtime()-timefs
+  if (write_runtimelog .and. is_root_pe()) write(stdout,*) 'In ',trim(subname),' time ', MPI_Wtime()-timefs
 
   if (localPet == 0) call cap_profiling("mom", "ocean_model_finalize", "E")