Remove MAPL_GenGridName from Base (#4835) (#4836)

Author: tclune

base/Base/Base_Base.F90

@@ -38,7 +38,6 @@ module MAPL_Base
   public MAPL_GetHorzIJIndex
   public MAPL_GetGlobalHorzIJIndex
   public MAPL_Reverse_Schmidt
-  public MAPL_GenGridName
    public MAPL_GridGetCorners
    public MAPL_GridGetInterior
 
@@ -133,13 +132,6 @@ module subroutine MAPL_Reverse_Schmidt(Grid, stretched, npts,lon,lat,lonR8,latR8
        integer,            optional, intent(out  ) :: rc  ! return code
      end subroutine MAPL_Reverse_Schmidt
 
-     module subroutine MAPL_GenGridName(im, jm, lon, lat, xyoffset, gridname, geos_style)
-       integer :: im, jm
-       character (len=*) :: gridname
-       real, optional    :: lon(:), lat(:)
-       integer, optional :: xyoffset
-       logical,  optional :: geos_style
-     end subroutine MAPL_GenGridName
 
        module function MAPL_StrUpCase(str) result(new)
          character(len=*), intent(IN) :: str

base/Base/Base_Base_implementation.F90

@@ -1340,91 +1340,6 @@ function grid_is_ok(grid) result(OK)
     end function
   end subroutine MAPL_GetGlobalHorzIJIndex
 
-  module subroutine MAPL_GenGridName(im, jm, lon, lat, xyoffset, gridname, geos_style)
-    integer :: im, jm
-    character (len=*) :: gridname
-    character(len=2)  :: dateline, pole
-    real, optional    :: lon(:), lat(:)
-    integer, optional :: xyoffset
-    logical,  optional :: geos_style
-
-    integer           :: I
-    real, parameter   :: eps=1.0e-4
-    character(len=16) :: imstr, jmstr
-    real              :: dlat
-
-    logical :: old_style
-    if (present(geos_style)) then
-       old_style = geos_style
-    else
-       old_style = .false.
-    end if
-
-    if (jm /= 6*im) then
-       ! Lat-Lon
-       dateline='UU' ! Undefined
-       pole='UU'     ! Undefined
-       if (present(LON) .and. present(LAT)) then
-
-          dlat = LAT(4) - LAT(3)
-          if ((abs(LAT(1) + 90.0) < eps).or.(abs(LAT(1) + 90.0 - 0.25*dLat) < eps)) then
-             pole='PC'
-          else if (abs(LAT(1) + 90.0 - 0.5*dLat) < eps) then
-             pole='PE'
-          end if
-          do I=0,1
-             if(abs(LON(1) + 180.0*I) < eps) then
-                dateline='DC'
-                exit
-             else if (abs(LON(1) + 180.0*I - 0.5*(LON(2)-LON(1))) < eps) then
-                dateline='DE'
-                exit
-             end if
-          end do
-
-       else if (present(xyoffset)) then
-          ! xyoffset Optional Flag for Grid Staggering (0:DcPc, 1:DePc, 2:DcPe, 3:DePe)
-          select case (xyoffset)
-          case (0)
-             dateline='DC'
-             pole='PC'
-          case (1)
-             dateline='DE'
-             pole='PC'
-          case (2)
-             dateline='DC'
-             pole='PE'
-          case (3)
-             dateline='DE'
-             pole='PE'
-          end select
-       endif
-
-       if (old_style) then
-          write(imstr,*) im
-          write(jmstr,*) jm
-          gridname =  pole // trim(adjustl(imstr))//'x'//&
-               trim(adjustl(jmstr))//'-'//dateline
-       else
-          write(gridname,'(a,i4.4,a,a,i4.4)') dateline,im,'x',pole,jm
-       end if
-    else
-       ! cubed-sphere
-       dateline='CF'
-       pole='6C'
-       if (old_style) then
-          pole='PE'
-          write(imstr,*) im
-          write(jmstr,*) jm
-          gridname =  pole // trim(adjustl(imstr))//'x'//&
-               trim(adjustl(jmstr))//'-CF'
-       else
-          write(gridname,'(a,i4.4,a,a)') dateline,im,'x',pole
-       end if
-    end if
-
-  end subroutine MAPL_GenGridName
-
   module subroutine MAPL_Reverse_Schmidt(Grid, stretched, npts, lon, lat, lonR8, latR8, lonRe, latRe, rc)
      type(ESMF_Grid), intent(inout) :: Grid
      logical, intent(out)           :: stretched