From 01175479ecd58c73fef0d05721952f5de601270b Mon Sep 17 00:00:00 2001
From: Edward Hartnett <edward.hartnett@noaa.gov>
Date: Tue, 8 Aug 2023 03:02:06 -0600
Subject: [PATCH 1/3] converted getlocal() to F90

---
 src/CMakeLists.txt |   2 +-
 src/getlocal.F90   | 144 +++++++++++++++++++++++++++++++++++++++++++
 src/getlocal.f     | 150 ---------------------------------------------
 3 files changed, 145 insertions(+), 151 deletions(-)
 create mode 100644 src/getlocal.F90
 delete mode 100644 src/getlocal.f

diff --git a/src/CMakeLists.txt b/src/CMakeLists.txt
index 32b2e06f..c5abfff3 100644
--- a/src/CMakeLists.txt
+++ b/src/CMakeLists.txt
@@ -7,7 +7,7 @@
 set(fortran_src addfield.f addgrid.F90 addlocal.F90 cmplxpack.f compack.f
 comunpack.f drstemplates.F90 g2_gbytesc.F90 g2grids.F90 gb_info.F90
 getdim.f getfield.F90 getg2i.F90 getg2ir.F90 getgb2.F90 getgb2l.f getgb2p.F90
-getgb2r.f getgb2rp.f getgb2s.F90 getidx.F90 getlocal.f getpoly.f
+getgb2r.f getgb2rp.f getgb2s.F90 getidx.F90 getlocal.F90 getpoly.f
 gettemplates.F90 gf_free.F90 gf_getfld.F90 gf_unpack1.F90 gf_unpack2.F90
 gf_unpack3.F90 gf_unpack4.f gf_unpack5.f gf_unpack6.F90 gf_unpack7.f
 gribcreate.F90 gribend.F90 gribinfo.F90
diff --git a/src/getlocal.F90 b/src/getlocal.F90
new file mode 100644
index 00000000..eba61864
--- /dev/null
+++ b/src/getlocal.F90
@@ -0,0 +1,144 @@
+!> @file
+!> @brief This subroutine returns the contents of Section 2 from a
+!> GRIB2 message.
+!> @author Stephen Gilbert @date 2000-05-25
+
+!> This subroutine returns the contents of Section 2 from a GRIB2
+!> message.
+!>
+!> Since there can be multiple occurrences of Section 2 within a GRIB
+!> message, the calling routine indicates which occurrence is being
+!> requested with the localnum argument.
+!>
+!> @note Note that subroutine gb_info() can be used to first determine
+!> how many Local Use sections exist in a given GRIB message.
+!>
+!> @param[in] cgrib Character array that contains the GRIB2 message.
+!> @param[in] lcgrib Length (in bytes) of GRIB message array cgrib.
+!> @param[in] localnum The nth occurrence of Section 2 requested.
+!> @param[out] csec2 Character array containing information read from
+!> Section 2. The dimension of this array can be obtained in advance
+!> from argument maxlocal, which is returned from subroutine gb_info().
+!> @param[out] lcsec2 Number of bytes of character array csec2 read
+!> from Section 2.
+!> @param[out] ierr Error return code.
+!> - 0 no error.
+!> - 1 Beginning characters "GRIB" not found.
+!> - 2 GRIB message is not Edition 2.
+!> - 3 The data field request number was not positive.
+!> - 4 End string "7777" found, but not where expected.
+!> - 5 End string "7777" not found at end of message.
+!> - 6 GRIB message did not contain the requested number of data fields.
+!>
+!> @author Stephen Gilbert @date 2000-05-25
+subroutine getlocal(cgrib, lcgrib, localnum, csec2, lcsec2, ierr)
+  implicit none
+
+  character(len = 1), intent(in) :: cgrib(lcgrib)
+  integer, intent(in) :: lcgrib, localnum
+  character(len = 1), intent(out) :: csec2(*)
+  integer, intent(out) :: lcsec2, ierr
+
+  character(len = 4), parameter :: grib = 'GRIB', c7777 = '7777'
+  character(len = 4) :: ctemp
+  integer :: listsec0(2)
+  integer iofst, istart, numlocal
+  integer :: lengrib, lensec, lensec0, j, ipos, isecnum
+
+  ierr = 0
+  numlocal = 0
+
+  ! Check for valid request number.
+  if (localnum .le. 0) then
+     print *, 'getlocal: Request for local section must be positive.'
+     ierr = 3
+     return
+  endif
+
+  ! Check for beginning of GRIB message in the first 100 bytes
+  istart = 0
+  do j = 1, 100
+     ctemp = cgrib(j)//cgrib(j+1)//cgrib(j+2)//cgrib(j+3)
+     if (ctemp .eq. grib) then
+        istart = j
+        exit
+     endif
+  enddo
+  if (istart .eq. 0) then
+     print *, 'getlocal:  Beginning characters GRIB not found.'
+     ierr = 1
+     return
+  endif
+
+  ! Unpack Section 0 - Indicator Section
+  iofst = 8 * (istart + 5)
+  call g2_gbytec(cgrib, listsec0(1), iofst, 8)     ! Discipline
+  iofst = iofst + 8
+  call g2_gbytec(cgrib, listsec0(2), iofst, 8)     ! GRIB edition number
+  iofst = iofst + 8
+  iofst = iofst + 32
+  call g2_gbytec(cgrib, lengrib, iofst, 32)        ! Length of GRIB message
+  iofst = iofst + 32
+  lensec0 = 16
+  ipos = istart + lensec0
+
+  ! Currently handles only GRIB Edition 2.
+  if (listsec0(2) .ne. 2) then
+     print *, 'getlocal: can only decode GRIB edition 2.'
+     ierr = 2
+     return
+  endif
+
+  ! Loop through the remaining sections keeping track of the length of
+  ! each. Also check to see that if the current occurrence of Section
+  ! 2 is the same as the one requested.
+  do
+     ! Check to see if we are at end of GRIB message
+     ctemp = cgrib(ipos) // cgrib(ipos + 1) // cgrib(ipos + 2) // cgrib(ipos + 3)
+     if (ctemp .eq. c7777) then
+        ipos = ipos + 4
+        
+        ! If end of GRIB message not where expected, issue error
+        if (ipos .ne. (istart + lengrib)) then
+           print *, 'getlocal: "7777" found, but not where expected.'
+           ierr = 4
+           return
+        endif
+        exit
+     endif
+
+     ! Get length of Section and Section number
+     iofst = (ipos - 1) * 8
+     call g2_gbytec(cgrib, lensec, iofst, 32)        ! Get Length of Section
+     iofst = iofst + 32
+     call g2_gbytec(cgrib, isecnum, iofst, 8)         ! Get Section number
+     iofst = iofst + 8
+     
+     ! If found the requested occurrence of Section 2, 
+     ! return the section contents.
+     if (isecnum .eq. 2) then
+        numlocal = numlocal + 1
+        if (numlocal.eq.localnum) then
+           lcsec2 = lensec - 5
+           csec2(1:lcsec2) = cgrib(ipos + 5:ipos + lensec - 1)
+           return
+        endif
+     endif
+     
+     ! Check to see if we read pass the end of the GRIB
+     ! message and missed the terminator string '7777'.
+     ipos = ipos + lensec                 ! Update beginning of section pointer
+     if (ipos .gt. (istart + lengrib)) then
+        print *, 'getlocal: "7777"  not found at end of GRIB message.'
+        ierr = 5
+        return
+     endif
+  enddo
+
+  ! If exited from above loop, the end of the GRIB message was reached
+  ! before the requested occurrence of section 2 was found.
+  print *, 'getlocal: GRIB message contained ', numlocal, ' local sections.'
+  print *, 'getlocal: The request was for the ', localnum,  ' occurrence.'
+  ierr = 6
+
+end subroutine getlocal
diff --git a/src/getlocal.f b/src/getlocal.f
deleted file mode 100644
index 4094f46f..00000000
--- a/src/getlocal.f
+++ /dev/null
@@ -1,150 +0,0 @@
-!> @file
-!> @brief This subroutine returns the contents of Section 2 from a
-!> GRIB2 message.
-!> @author Stephen Gilbert @date 2000-05-25
-
-!> This subroutine returns the contents of Section 2 from a GRIB2
-!> message.
-!>
-!> Since there can be multiple occurrences of Section 2 within a GRIB
-!> message, the calling routine indicates which occurrence is being
-!> requested with the localnum argument.
-!>
-!> @param[in] cgrib Character array that contains the GRIB2 message.
-!> @param[in] lcgrib Length (in bytes) of GRIB message array cgrib.
-!> @param[in] localnum The nth occurrence of Section 2 requested.
-!> @param[out] csec2 Character array containing information read from
-!> Section 2. The dimension of this array can be obtained in advance
-!> from argument maxlocal, which is returned from subroutine gb_info().
-!> @param[out] lcsec2 Number of bytes of character array csec2 read
-!> from Section 2.
-!> @param[out] ierr Error return code.
-!> - 0 no error.
-!> - 1 Beginning characters "GRIB" not found.
-!> - 2 GRIB message is not Edition 2.
-!> - 3 The data field request number was not positive.
-!> - 4 End string "7777" found, but not where expected.
-!> - 5 End string "7777" not found at end of message.
-!> - 6 GRIB message did not contain the requested number of data fields.
-!>
-!> @note Note that subroutine gb_info() can be used to first determine
-!> how many Local Use sections exist in a given GRIB message.
-!>
-!> @author Stephen Gilbert @date 2000-05-25
-      subroutine getlocal(cgrib,lcgrib,localnum,csec2,lcsec2,ierr)
-      implicit none
-
-      character(len=1),intent(in) :: cgrib(lcgrib)
-      integer,intent(in) :: lcgrib,localnum
-      character(len=1),intent(out) :: csec2(*)
-      integer,intent(out) :: lcsec2,ierr
-      
-      character(len=4),parameter :: grib='GRIB',c7777='7777'
-      character(len=4) :: ctemp
-      integer :: listsec0(2)
-      integer iofst,istart,numlocal
-      integer :: lengrib, lensec, lensec0, j, ipos, isecnum
-
-      ierr=0
-      numlocal=0
-
-!     Check for valid request number
-      if (localnum.le.0) then
-        print *,'getlocal: Request for local section must be positive.'
-        ierr=3
-        return
-      endif
-
-!     Check for beginning of GRIB message in the first 100 bytes
-      istart=0
-      do j=1,100
-        ctemp=cgrib(j)//cgrib(j+1)//cgrib(j+2)//cgrib(j+3)
-        if (ctemp.eq.grib ) then
-          istart=j
-          exit
-        endif
-      enddo
-      if (istart.eq.0) then
-        print *,'getlocal:  Beginning characters GRIB not found.'
-        ierr=1
-        return
-      endif
-
-!     Unpack Section 0 - Indicator Section 
-      iofst=8*(istart+5)
-      call g2_gbytec(cgrib,listsec0(1),iofst,8)     ! Discipline
-      iofst=iofst+8
-      call g2_gbytec(cgrib,listsec0(2),iofst,8)     ! GRIB edition number
-      iofst=iofst+8
-      iofst=iofst+32
-      call g2_gbytec(cgrib,lengrib,iofst,32)        ! Length of GRIB message
-      iofst=iofst+32
-      lensec0=16
-      ipos=istart+lensec0
-
-!     Currently handles only GRIB Edition 2.
-      if (listsec0(2).ne.2) then
-        print *,'getlocal: can only decode GRIB edition 2.'
-        ierr=2
-        return
-      endif
-
-!     Loop through the remaining sections keeping track of the length of
-!     each. Also check to see that if the current occurrence of Section
-!     2 is the same as the one requested.
-      do
-        !    Check to see if we are at end of GRIB message
-        ctemp=cgrib(ipos)//cgrib(ipos+1)//cgrib(ipos+2)//cgrib(ipos+3)
-        if (ctemp.eq.c7777 ) then
-          ipos=ipos+4
-          !    If end of GRIB message not where expected, issue error
-          if (ipos.ne.(istart+lengrib)) then
-            print *,'getlocal: "7777" found, but not where expected.'
-            ierr=4
-            return
-          endif
-          exit
-        endif
-        !     Get length of Section and Section number
-        iofst=(ipos-1)*8
-        call g2_gbytec(cgrib,lensec,iofst,32)        ! Get Length of Section
-        iofst=iofst+32
-        call g2_gbytec(cgrib,isecnum,iofst,8)         ! Get Section number
-        iofst=iofst+8
-        !   If found the requested occurrence of Section 2,
-        !   return the section contents.
-        if (isecnum.eq.2) then
-          numlocal=numlocal+1
-          if (numlocal.eq.localnum) then
-            lcsec2=lensec-5
-            csec2(1:lcsec2)=cgrib(ipos+5:ipos+lensec-1)
-            return
-          endif
-        endif
-        !   Check to see if we read pass the end of the GRIB
-        !   message and missed the terminator string '7777'.
-        ipos=ipos+lensec                 ! Update beginning of section pointer
-        if (ipos.gt.(istart+lengrib)) then
-          print *,'getlocal: "7777"  not found at end of GRIB message.'
-          ierr=5
-          return
-        endif
-        
-      enddo
-
-!     If exited from above loop, the end of the GRIB message was reached
-!     before the requested occurrence of section 2 was found.
-      print *,'getlocal: GRIB message contained ',numlocal,
-     &        ' local sections.'
-      print *,'getlocal: The request was for the ',localnum,
-     &        ' occurrence.'
-      ierr=6
-
-      end
-
-
-
-
-
-
-

From 3cddfddf60fe2b1432c8685560ef10a937a67269 Mon Sep 17 00:00:00 2001
From: Edward Hartnett <edward.hartnett@noaa.gov>
Date: Tue, 8 Aug 2023 03:10:34 -0600
Subject: [PATCH 2/3] converted gridtemplates to F90

---
 src/CMakeLists.txt    |   2 +-
 src/gridtemplates.F90 | 415 +++++++++++++++++++++++++++++++++++++++++
 src/gridtemplates.f   | 419 ------------------------------------------
 3 files changed, 416 insertions(+), 420 deletions(-)
 create mode 100644 src/gridtemplates.F90
 delete mode 100644 src/gridtemplates.f

diff --git a/src/CMakeLists.txt b/src/CMakeLists.txt
index c5abfff3..82852616 100644
--- a/src/CMakeLists.txt
+++ b/src/CMakeLists.txt
@@ -11,7 +11,7 @@ getgb2r.f getgb2rp.f getgb2s.F90 getidx.F90 getlocal.F90 getpoly.f
 gettemplates.F90 gf_free.F90 gf_getfld.F90 gf_unpack1.F90 gf_unpack2.F90
 gf_unpack3.F90 gf_unpack4.f gf_unpack5.f gf_unpack6.F90 gf_unpack7.f
 gribcreate.F90 gribend.F90 gribinfo.F90
-${CMAKE_CURRENT_BINARY_DIR}/gribmod.F90 gridtemplates.f intmath.f
+${CMAKE_CURRENT_BINARY_DIR}/gribmod.F90 gridtemplates.F90 intmath.f
 ixgb2.f jpcpack.F90 jpcunpack.F90 misspack.f mkieee.f pack_gp.f
 params_ecmwf.F90 params.F90 pdstemplates.F90 pngpack.F90 pngunpack.F90
 putgb2.F90 rdieee.f realloc.f reduce.f simpack.f simunpack.F90 skgb.F90
diff --git a/src/gridtemplates.F90 b/src/gridtemplates.F90
new file mode 100644
index 00000000..13a34c25
--- /dev/null
+++ b/src/gridtemplates.F90
@@ -0,0 +1,415 @@
+!> @file
+!> @brief This Fortran module contains info on all the available
+!> GRIB2 Grid Definition Templates used in Section 3 - the Grid
+!> Definition Section (GDS).
+!> @author Stephen Gilbert @date 2000-05-09
+
+!> This Fortran module contains info on all the available GRIB2 Grid
+!> Definition Templates used in [Section 3 - the Grid Definition
+!> Section
+!> (GDS)](https://www.nco.ncep.noaa.gov/pmb/docs/grib2/grib2_doc/grib2_sect3.shtml).
+!>
+!> Each Template has three parts:
+!> 1. The number of entries in the template (mapgridlen);
+!> 2. A map of the template (mapgrid), which contains the number of
+!> octets in which to pack each of the template values;
+!> 3. A logical value (needext) that indicates whether the Template
+!> needs to be extended. In some cases the number of entries in a
+!> template can vary depending upon values specified in the "static"
+!> part of the template. (See Template 3.120 as an example).
+!>
+!> This module also contains two subroutines:
+!> * getgridtemplate() returns the octet map for a specified Template
+!> number
+!> * extgridtemplate() calculates the extended octet map of a template
+!> that needs extension.
+!>
+!> Array mapgrid contains the number of bytes in which the
+!> corresponding template values will be stored. A negative value in
+!> mapgrid is used to indicate that the corresponding template entry
+!> can contain negative values. This information is used later when
+!> packing/unpacking the template data values.
+!>
+!> Negative data values in GRIB are stored with the left most bit set to
+!> one, and a negative number of bytes value in mapgrid indicates that
+!> this possibility should be considered. The number of bytes used to
+!> store the data value in this case would be the absolute value of the
+!> negative value in mapgrid.
+!>
+!> @author Stephen Gilbert @date 2000-05-09
+module gridtemplates
+
+  integer, parameter :: MAXLEN = 200 !< maximum number of octets in mapgrid
+  integer, parameter :: MAXTEMP = 31 !< maximum number of entries in the template
+
+  type gridtemplate
+     integer :: template_num
+     integer :: mapgridlen
+     integer, dimension(MAXLEN) :: mapgrid
+     logical :: needext
+  end type gridtemplate
+
+  type(gridtemplate), dimension(MAXTEMP) :: templates !< a gridtemplate type
+
+  data templates(1)%template_num /0/     !  Lat/Lon
+  data templates(1)%mapgridlen /19/
+  data templates(1)%needext /.false./
+  data (templates(1)%mapgrid(j), j = 1, 19)
+  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1/
+
+  data templates(2)%template_num /1/     !  Rotated Lat/Lon
+  data templates(2)%mapgridlen /22/
+  data templates(2)%needext /.false./
+  data (templates(2)%mapgrid(j), j = 1, 22)
+  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1, -4, 4, 4/
+
+  data templates(3)%template_num /2/     !  Stretched Lat/Lon
+  data templates(3)%mapgridlen /22/
+  data templates(3)%needext /.false./
+  data (templates(3)%mapgrid(j), j = 1, 22)
+  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1, -4, 4, -4/
+
+  data templates(4)%template_num /3/     !  Stretched & Rotated Lat/Lon
+  data templates(4)%mapgridlen /25/
+  data templates(4)%needext /.false./
+  data (templates(4)%mapgrid(j), j = 1, 25)
+  &       /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1, -4, 4, 4, -4, 4, -4/
+
+  data templates(5)%template_num /10/     !  Mercator
+  data templates(5)%mapgridlen /19/
+  data templates(5)%needext /.false./
+  data (templates(5)%mapgrid(j), j = 1, 19)
+  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, -4, 4, 1, -4, -4, 4, 1, 4, 4, 4/
+
+  data templates(6)%template_num /20/     !  Polar Stereographic
+  data templates(6)%mapgridlen /18/
+  data templates(6)%needext /.false./
+  data (templates(6)%mapgrid(j), j = 1, 18)
+  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1, 1/
+
+  data templates(7)%template_num /30/     !  Lambert Conformal
+  data templates(7)%mapgridlen /22/
+  data templates(7)%needext /.false./
+  data (templates(7)%mapgrid(j), j = 1, 22)
+  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1, 1, -4, -4, -4, 4/
+
+  data templates(8)%template_num /40/     !  Gaussian Lat/Lon
+  data templates(8)%mapgridlen /19/
+  data templates(8)%needext /.false./
+  data (templates(8)%mapgrid(j), j = 1, 19)
+  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1/
+
+  data templates(9)%template_num /41/     !  Rotated Gaussian Lat/Lon
+  data templates(9)%mapgridlen /22/
+  data templates(9)%needext /.false./
+  data (templates(9)%mapgrid(j), j = 1, 22)
+  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1, -4, 4, 4/
+
+  data templates(10)%template_num /42/     !  Stretched Gaussian Lat/Lon
+  data templates(10)%mapgridlen /22/
+  data templates(10)%needext /.false./
+  data (templates(10)%mapgrid(j), j = 1, 22)
+  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1, -4, 4, -4/
+
+  data templates(11)%template_num /43/     !  Strtchd and Rot'd Gaus Lat/Lon
+  data templates(11)%mapgridlen /25/
+  data templates(11)%needext /.false./
+  data (templates(11)%mapgrid(j), j = 1, 25)
+  &          /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1, -4, 4, 4, -4, 4, -4/
+
+  data templates(12)%template_num /50/    !  Spherical Harmonic Coefficients
+  data templates(12)%mapgridlen /5/
+  data templates(12)%needext /.false./
+  data (templates(12)%mapgrid(j), j = 1, 5) /4, 4, 4, 1, 1/
+
+  data templates(13)%template_num /51/   !  Rotated Spherical Harmonic Coeff
+  data templates(13)%mapgridlen /8/
+  data templates(13)%needext /.false./
+  data (templates(13)%mapgrid(j), j = 1, 8) /4, 4, 4, 1, 1, -4, 4, 4/
+
+  data templates(14)%template_num /52/   !  Stretch Spherical Harmonic Coeff
+  data templates(14)%mapgridlen /8/
+  data templates(14)%needext /.false./
+  data (templates(14)%mapgrid(j), j = 1, 8) /4, 4, 4, 1, 1, -4, 4, -4/
+
+  data templates(15)%template_num /53/   !  Strch and Rot Spher Harm Coeffs
+  data templates(15)%mapgridlen /11/
+  data templates(15)%needext /.false./
+  data (templates(15)%mapgrid(j), j = 1, 11) /4, 4, 4, 1, 1, -4, 4, 4, -4, 4, -4/
+
+  data templates(16)%template_num /90/     !  Space view Perspective
+  data templates(16)%mapgridlen /21/
+  data templates(16)%needext /.false./
+  data (templates(16)%mapgrid(j), j = 1, 21)
+  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, -4, 4, 1, 4, 4, 4, 4, 1, 4, 4, 4, 4/
+
+  data templates(17)%template_num /100/    !  Triangular grid (icosahedron)
+  data templates(17)%mapgridlen /11/
+  data templates(17)%needext /.false./
+  data (templates(17)%mapgrid(j), j = 1, 11) /1, 1, 2, 1, -4, 4, 4, 1, 1, 1, 4/
+
+  data templates(18)%template_num /110/ !  Equatorial Azimuthal equidistant
+  data templates(18)%mapgridlen /16/
+  data templates(18)%needext /.false./
+  data (templates(18)%mapgrid(j), j = 1, 16)
+  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, -4, 4, 1, 4, 4, 1, 1/
+
+  data templates(19)%template_num /120/     !  Azimuth-range
+  data templates(19)%mapgridlen /7/
+  data templates(19)%needext /.true./
+  data (templates(19)%mapgrid(j), j = 1, 7) /4, 4, -4, 4, 4, 4, 1/
+
+  data templates(20)%template_num /1000/     !  Cross Section Grid
+  data templates(20)%mapgridlen /20/
+  data templates(20)%needext /.true./
+  data (templates(20)%mapgrid(j), j = 1, 20)
+  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, -4, 4, 1, 4, 4, 1, 2, 1, 1, 2/
+
+  data templates(21)%template_num /1100/     !  Hovmoller Diagram Grid
+  data templates(21)%mapgridlen /28/
+  data templates(21)%needext /.false./
+  data (templates(21)%mapgrid(j), j = 1, 28)
+  &    /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, -4, 4, 1, -4, 4, 1, 4, 1, -4, 1, 1, -4, 2, 1, 1, 1, 1, 1/
+
+  data templates(22)%template_num /1200/     !  Time Section Grid
+  data templates(22)%mapgridlen /16/
+  data templates(22)%needext /.true./
+  data (templates(22)%mapgrid(j), j = 1, 16)
+  &              /4, 1, -4, 1, 1, -4, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2/
+
+  data templates(23)%template_num /31/     !  Albers Equal Area
+  data templates(23)%mapgridlen /22/
+  data templates(23)%needext /.false./
+  data (templates(23)%mapgrid(j), j = 1, 22)
+  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1, 1, -4, -4, -4, 4/
+
+  data templates(24)%template_num /204/     !  Curilinear Orthogonal Grids
+  data templates(24)%mapgridlen /19/
+  data templates(24)%needext /.false./
+  data (templates(24)%mapgrid(j), j = 1, 19)
+  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1/
+
+  data templates(25)%template_num /32768/   ! Rotate Lat/Lon E-grid
+  data templates(25)%mapgridlen /19/
+  data templates(25)%needext /.false./
+  data (templates(25)%mapgrid(j), j = 1, 19)
+  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1/
+
+  data templates(26)%template_num /32769/   ! Rotate Lat/Lon Non-E Stagger grid
+  data templates(26)%mapgridlen /21/
+  data templates(26)%needext /.false./
+  data (templates(26)%mapgrid(j), j = 1, 21)
+  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1, 4, 4/
+  !
+  !     GDT 3.4 Added (08/05/2013)
+  !
+  data templates(27)%template_num /4/     !  Variable resolution Latitude/Longitude
+  data templates(27)%mapgridlen /13/
+  data templates(27)%needext /.true./
+  data (templates(27)%mapgrid(j), j = 1, 13)
+  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, 1, 1/
+  !
+  !     GDT 3.5 Added (08/05/2013)
+  !
+  data templates(28)%template_num /5/     !  Variable resolution rotate Latitude/Longitude
+  data templates(28)%mapgridlen /16/
+  data templates(28)%needext /.true./
+  data (templates(28)%mapgrid(j), j = 1, 16)
+  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, 1, 1, -4, 4, 4/
+  !
+  !     GDT 3.12 Added (08/05/2013)
+  !
+  data templates(29)%template_num /12/     !  Transverse Mercator
+  data templates(29)%mapgridlen /22/
+  data templates(29)%needext /.false./
+  data (templates(29)%mapgrid(j), j = 1, 22)
+  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, -4, 4, 1, -4, 4, 4, 1, 4, 4, -4, -4, -4, -4/
+  !
+  !     GDT 3.101 Added (08/05/2013)
+  !
+  data templates(30)%template_num /101/     !  General unstructured grid
+  data templates(30)%mapgridlen /4/
+  data templates(30)%needext /.false./
+  data (templates(30)%mapgrid(j), j = 1, 4)
+  &              /1, 4, 1, -4/
+  !
+  !     GDT 3.140 Added (08/05/2013)
+  !
+  data templates(31)%template_num /140/     !  Lambert Azimuthal Equal Area Projection
+  data templates(31)%mapgridlen /17/
+  data templates(31)%needext /.false./
+  data (templates(31)%mapgrid(j), j = 1, 17)
+  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, -4, 4, 4, 4, 1, 4, 4, 1/
+
+contains
+
+  !> Return the index of specified Grid Definition Template.
+  !>
+  !> @param[in] number NN, indicating the number of the Grid
+  !> Definition Template 3.NN that is being requested.
+  !>
+  !> @return Index of the grid template in array templates, if
+  !> template exists.  -1, otherwise.
+  !>
+  !> @author Stephen Gilbert @date 2001-06-28
+  integer function getgridindex(number)
+    implicit none
+    integer, intent(in) :: number
+    integer :: j
+
+    getgridindex = -1
+
+    do j = 1, MAXTEMP
+       if (number.eq.templates(j)%template_num) then
+          getgridindex = j
+          return
+       endif
+    enddo
+
+  end function getgridindex
+
+  !> Get the grid template information for a specified Grid Definition
+  !> Template.
+  !>
+  !> The number of entries in the template is returned along with a
+  !> map of the number of octets occupied by each entry. Also, a flag
+  !> is returned to indicate whether the template would need to be
+  !> extended.
+  !>
+  !> @param[in] number NN, indicating the number of the Grid
+  !> Definition Template that is being requested.
+  !> @param[out] nummap Number of entries in the Template.
+  !> @param[out] map An array containing the number of octets that
+  !> each template entry occupies when packed up into the GDS.
+  !> @param[out] needext Logical variable indicating whether the Grid
+  !> Defintion Template has to be extended.
+  !> @param[out] iret Error return code.
+  !> - 0 no error.
+  !> - 1 Undefine Grid Template number.
+  !>
+  !> @author Stephen Gilbert @date 2000-05-09
+  subroutine getgridtemplate(number, nummap, map, needext, iret)
+    implicit none
+
+    integer, intent(in) :: number
+    integer, intent(out) :: nummap, map(*), iret
+    logical, intent(out) :: needext
+    integer :: index
+
+    iret = 0
+
+    index = getgridindex(number)
+
+    if (index.ne.-1) then
+       nummap = templates(index)%mapgridlen
+       needext = templates(index)%needext
+       map(1:nummap) = templates(index)%mapgrid(1:nummap)
+    else
+       nummap = 0
+       needext = .false.
+       print *, 'getgridtemplate: Grid Template ', number, 
+       &               ' not defined.'
+       iret = 1
+    endif
+  end subroutine getgridtemplate
+
+  !> Generate the remaining octet map for a given Grid Definition
+  !> Template, if required.
+  !>
+  !> Some Templates can vary depending on data values given in an
+  !> earlier part of the Template, and it is necessary to know some of
+  !> the earlier entry values to generate the full octet map of the
+  !> Template.
+  !>
+  !> @param[in] number NN, indicating the number of the Grid
+  !> Definition Template 3.NN that is being requested.
+  !> @param[in] list The list of values for each entry in the Grid
+  !> Definition Template.
+  !> @param[out] nummap Number of entries in the Template.
+  !> @param[out] map An array containing the number of octets that
+  !> each template entry occupies when packed up into the GDS.
+  !>
+  !>    @author Stephen Gilbert @date 2000-05-09
+  subroutine extgridtemplate(number, list, nummap, map)
+    implicit none
+
+    integer, intent(in) :: number, list(*)
+    integer, intent(out) :: nummap, map(*)
+    integer i,  index,  n,  ni,  nj
+
+    index = getgridindex(number)
+    if (index.eq.-1) return
+
+    if (.not. templates(index)%needext) return
+    nummap = templates(index)%mapgridlen
+    map(1:nummap) = templates(index)%mapgrid(1:nummap)
+
+    if (number.eq.120) then
+       N = list(2)
+       do i = 1, N
+          map(nummap + 1) = 2
+          map(nummap + 2) = -2
+          nummap = nummap + 2
+       enddo
+    elseif (number.eq.4) then
+       NI = list(8)
+       do i = 1, NI
+          map(nummap + 1) = 4
+          nummap = nummap + 1
+       enddo
+       NJ = list(9)
+       do i = 1, NJ
+          map(nummap + 1) = -4
+          nummap = nummap + 1
+       enddo
+    elseif (number.eq.5) then
+       NI = list(8)
+       do i = 1, NI
+          map(nummap + 1) = 4
+          nummap = nummap + 1
+       enddo
+       NJ = list(9)
+       do i = 1, NJ
+          map(nummap + 1) = -4
+          nummap = nummap + 1
+       enddo
+    elseif (number.eq.1000) then
+       N = list(20)
+       do i = 1, N
+          map(nummap + 1) = 4
+          nummap = nummap + 1
+       enddo
+    elseif (number.eq.1200 ) then
+       N = list(16)
+       do i = 1, N
+          map(nummap + 1) = 4
+          nummap = nummap + 1
+       enddo
+    endif
+  end subroutine extgridtemplate
+
+  !> This function returns the initial length (number of entries) in
+  !> the static part of specified Grid Definition Template.
+  !>
+  !> @param[in] number NN,  indicating the number of the Grid
+  !> Definition Template that is being requested.
+  !>
+  !> @return Number of entries in the static part of the grid
+  !> definition template, or 0, if requested template is not found.
+  !>
+  !> @author Stephen Gilbert @date 2004-05-11
+  integer function getgdtlen(number)
+    implicit none
+
+    integer, intent(in) :: number
+    integer :: index
+
+    getgdtlen = 0
+
+    index = getgridindex(number)
+
+    if (index .ne. -1) then
+       getgdtlen = templates(index)%mapgridlen
+    endif
+  end function getgdtlen
+end module gridtemplates
diff --git a/src/gridtemplates.f b/src/gridtemplates.f
deleted file mode 100644
index 90a7a08c..00000000
--- a/src/gridtemplates.f
+++ /dev/null
@@ -1,419 +0,0 @@
-!> @file
-!> @brief This Fortran module contains info on all the available
-!> GRIB2 Grid Definition Templates used in Section 3 - the Grid
-!> Definition Section (GDS).
-!> @author Stephen Gilbert @date 2000-05-09
-
-!> This Fortran module contains info on all the available GRIB2 Grid
-!> Definition Templates used in [Section 3 - the Grid Definition
-!> Section
-!> (GDS)](https://www.nco.ncep.noaa.gov/pmb/docs/grib2/grib2_doc/grib2_sect3.shtml).
-!>
-!> Each Template has three parts:
-!> 1. The number of entries in the template (mapgridlen);
-!> 2. A map of the template (mapgrid), which contains the number of
-!> octets in which to pack each of the template values;
-!> 3. A logical value (needext) that indicates whether the Template
-!> needs to be extended. In some cases the number of entries in a
-!> template can vary depending upon values specified in the "static"
-!> part of the template. (See Template 3.120 as an example).
-!>
-!> This module also contains two subroutines:
-!> * getgridtemplate() returns the octet map for a specified Template
-!> number
-!> * extgridtemplate() calculates the extended octet map of a template
-!> that needs extension.
-!>
-!> Array mapgrid contains the number of bytes in which the
-!> corresponding template values will be stored. A negative value in
-!> mapgrid is used to indicate that the corresponding template entry
-!> can contain negative values. This information is used later when
-!> packing/unpacking the template data values.
-!>
-!> Negative data values in GRIB are stored with the left most bit set to
-!> one, and a negative number of bytes value in mapgrid indicates that
-!> this possibility should be considered. The number of bytes used to
-!> store the data value in this case would be the absolute value of the
-!> negative value in mapgrid.
-!>
-!> @author Stephen Gilbert @date 2000-05-09
-      module gridtemplates
-
-      integer,parameter :: MAXLEN=200 !< maximum number of octets in mapgrid
-      integer,parameter :: MAXTEMP=31 !< maximum number of entries in the template
-
-      type gridtemplate
-          integer :: template_num
-          integer :: mapgridlen
-          integer,dimension(MAXLEN) :: mapgrid
-          logical :: needext
-      end type gridtemplate
-
-      type(gridtemplate),dimension(MAXTEMP) :: templates !< a gridtemplate type
-
-      data templates(1)%template_num /0/     !  Lat/Lon 
-      data templates(1)%mapgridlen /19/
-      data templates(1)%needext /.false./
-      data (templates(1)%mapgrid(j),j=1,19) 
-     &              /1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1/
-
-      data templates(2)%template_num /1/     !  Rotated Lat/Lon 
-      data templates(2)%mapgridlen /22/
-      data templates(2)%needext /.false./
-      data (templates(2)%mapgrid(j),j=1,22) 
-     &              /1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1,-4,4,4/
-
-      data templates(3)%template_num /2/     !  Stretched Lat/Lon 
-      data templates(3)%mapgridlen /22/
-      data templates(3)%needext /.false./
-      data (templates(3)%mapgrid(j),j=1,22) 
-     &              /1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1,-4,4,-4/
-
-      data templates(4)%template_num /3/     !  Stretched & Rotated Lat/Lon 
-      data templates(4)%mapgridlen /25/
-      data templates(4)%needext /.false./
-      data (templates(4)%mapgrid(j),j=1,25) 
-     &       /1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1,-4,4,4,-4,4,-4/
-
-      data templates(5)%template_num /10/     !  Mercator
-      data templates(5)%mapgridlen /19/
-      data templates(5)%needext /.false./
-      data (templates(5)%mapgrid(j),j=1,19)
-     &              /1,1,4,1,4,1,4,4,4,-4,4,1,-4,-4,4,1,4,4,4/
-
-      data templates(6)%template_num /20/     !  Polar Stereographic
-      data templates(6)%mapgridlen /18/
-      data templates(6)%needext /.false./
-      data (templates(6)%mapgrid(j),j=1,18) 
-     &              /1,1,4,1,4,1,4,4,4,-4,4,1,-4,4,4,4,1,1/
-
-      data templates(7)%template_num /30/     !  Lambert Conformal
-      data templates(7)%mapgridlen /22/
-      data templates(7)%needext /.false./
-      data (templates(7)%mapgrid(j),j=1,22) 
-     &              /1,1,4,1,4,1,4,4,4,-4,4,1,-4,4,4,4,1,1,-4,-4,-4,4/
-
-      data templates(8)%template_num /40/     !  Gaussian Lat/Lon
-      data templates(8)%mapgridlen /19/
-      data templates(8)%needext /.false./
-      data (templates(8)%mapgrid(j),j=1,19) 
-     &              /1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1/
-
-      data templates(9)%template_num /41/     !  Rotated Gaussian Lat/Lon
-      data templates(9)%mapgridlen /22/
-      data templates(9)%needext /.false./
-      data (templates(9)%mapgrid(j),j=1,22) 
-     &              /1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1,-4,4,4/
-
-      data templates(10)%template_num /42/     !  Stretched Gaussian Lat/Lon
-      data templates(10)%mapgridlen /22/
-      data templates(10)%needext /.false./
-      data (templates(10)%mapgrid(j),j=1,22) 
-     &              /1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1,-4,4,-4/
-
-      data templates(11)%template_num /43/     !  Strtchd and Rot'd Gaus Lat/Lon
-      data templates(11)%mapgridlen /25/
-      data templates(11)%needext /.false./
-      data (templates(11)%mapgrid(j),j=1,25) 
-     &          /1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1,-4,4,4,-4,4,-4/
-
-      data templates(12)%template_num /50/    !  Spherical Harmonic Coefficients
-      data templates(12)%mapgridlen /5/
-      data templates(12)%needext /.false./
-      data (templates(12)%mapgrid(j),j=1,5) /4,4,4,1,1/
-
-      data templates(13)%template_num /51/   !  Rotated Spherical Harmonic Coeff
-      data templates(13)%mapgridlen /8/
-      data templates(13)%needext /.false./
-      data (templates(13)%mapgrid(j),j=1,8) /4,4,4,1,1,-4,4,4/
-
-      data templates(14)%template_num /52/   !  Stretch Spherical Harmonic Coeff
-      data templates(14)%mapgridlen /8/
-      data templates(14)%needext /.false./
-      data (templates(14)%mapgrid(j),j=1,8) /4,4,4,1,1,-4,4,-4/
-
-      data templates(15)%template_num /53/   !  Strch and Rot Spher Harm Coeffs
-      data templates(15)%mapgridlen /11/
-      data templates(15)%needext /.false./
-      data (templates(15)%mapgrid(j),j=1,11) /4,4,4,1,1,-4,4,4,-4,4,-4/
-
-      data templates(16)%template_num /90/     !  Space view Perspective
-      data templates(16)%mapgridlen /21/
-      data templates(16)%needext /.false./
-      data (templates(16)%mapgrid(j),j=1,21) 
-     &              /1,1,4,1,4,1,4,4,4,-4,4,1,4,4,4,4,1,4,4,4,4/
-
-      data templates(17)%template_num /100/    !  Triangular grid (icosahedron)
-      data templates(17)%mapgridlen /11/
-      data templates(17)%needext /.false./
-      data (templates(17)%mapgrid(j),j=1,11) /1,1,2,1,-4,4,4,1,1,1,4/
-
-      data templates(18)%template_num /110/ !  Equatorial Azimuthal equidistant
-      data templates(18)%mapgridlen /16/
-      data templates(18)%needext /.false./
-      data (templates(18)%mapgrid(j),j=1,16) 
-     &              /1,1,4,1,4,1,4,4,4,-4,4,1,4,4,1,1/
-
-       data templates(19)%template_num /120/     !  Azimuth-range 
-       data templates(19)%mapgridlen /7/
-       data templates(19)%needext /.true./
-       data (templates(19)%mapgrid(j),j=1,7) /4,4,-4,4,4,4,1/
-
-       data templates(20)%template_num /1000/     !  Cross Section Grid 
-       data templates(20)%mapgridlen /20/
-       data templates(20)%needext /.true./
-       data (templates(20)%mapgrid(j),j=1,20) 
-     &              /1,1,4,1,4,1,4,4,4,4,-4,4,1,4,4,1,2,1,1,2/
-
-       data templates(21)%template_num /1100/     !  Hovmoller Diagram Grid 
-       data templates(21)%mapgridlen /28/
-       data templates(21)%needext /.false./
-       data (templates(21)%mapgrid(j),j=1,28) 
-     &    /1,1,4,1,4,1,4,4,4,4,-4,4,1,-4,4,1,4,1,-4,1,1,-4,2,1,1,1,1,1/
-
-       data templates(22)%template_num /1200/     !  Time Section Grid 
-       data templates(22)%mapgridlen /16/
-       data templates(22)%needext /.true./
-       data (templates(22)%mapgrid(j),j=1,16) 
-     &              /4,1,-4,1,1,-4,2,1,1,1,1,1,2,1,1,2/
-
-      data templates(23)%template_num /31/     !  Albers Equal Area
-      data templates(23)%mapgridlen /22/
-      data templates(23)%needext /.false./
-      data (templates(23)%mapgrid(j),j=1,22) 
-     &              /1,1,4,1,4,1,4,4,4,-4,4,1,-4,4,4,4,1,1,-4,-4,-4,4/
-
-      data templates(24)%template_num /204/     !  Curilinear Orthogonal Grids
-      data templates(24)%mapgridlen /19/
-      data templates(24)%needext /.false./
-      data (templates(24)%mapgrid(j),j=1,19)
-     &              /1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1/
-
-      data templates(25)%template_num /32768/   ! Rotate Lat/Lon E-grid
-      data templates(25)%mapgridlen /19/
-      data templates(25)%needext /.false./
-      data (templates(25)%mapgrid(j),j=1,19)
-     &              /1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1/
-
-      data templates(26)%template_num /32769/   ! Rotate Lat/Lon Non-E Stagger grid
-      data templates(26)%mapgridlen /21/
-      data templates(26)%needext /.false./
-      data (templates(26)%mapgrid(j),j=1,21)
-     &              /1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1,4,4/
-!
-!     GDT 3.4 Added (08/05/2013)
-!
-      data templates(27)%template_num /4/     !  Variable resolution Latitude/Longitude
-      data templates(27)%mapgridlen /13/
-      data templates(27)%needext /.true./
-      data (templates(27)%mapgrid(j),j=1,13)
-     &              /1,1,4,1,4,1,4,4,4,4,4,1,1/
-!
-!     GDT 3.5 Added (08/05/2013)
-!
-      data templates(28)%template_num /5/     !  Variable resolution rotate Latitude/Longitude
-      data templates(28)%mapgridlen /16/
-      data templates(28)%needext /.true./
-      data (templates(28)%mapgrid(j),j=1,16)
-     &              /1,1,4,1,4,1,4,4,4,4,4,1,1,-4,4,4/
-!
-!     GDT 3.12 Added (08/05/2013)
-!
-      data templates(29)%template_num /12/     !  Transverse Mercator
-      data templates(29)%mapgridlen /22/
-      data templates(29)%needext /.false./
-      data (templates(29)%mapgrid(j),j=1,22)
-     &              /1,1,4,1,4,1,4,4,4,-4,4,1,-4,4,4,1,4,4,-4,-4,-4,-4/
-!
-!     GDT 3.101 Added (08/05/2013)
-!
-      data templates(30)%template_num /101/     !  General unstructured grid
-      data templates(30)%mapgridlen /4/
-      data templates(30)%needext /.false./
-      data (templates(30)%mapgrid(j),j=1,4)
-     &              /1,4,1,-4/
-!
-!     GDT 3.140 Added (08/05/2013)
-!
-      data templates(31)%template_num /140/     !  Lambert Azimuthal Equal Area Projection
-      data templates(31)%mapgridlen /17/
-      data templates(31)%needext /.false./
-      data (templates(31)%mapgrid(j),j=1,17)
-     &              /1,1,4,1,4,1,4,4,4,-4,4,4,4,1,4,4,1/
-
-      contains
-
-!>    Return the index of specified Grid Definition Template.
-!>
-!>    @param[in] number NN, indicating the number of the Grid Definition
-!>    Template 3.NN that is being requested.
-!>
-!>    @return Index of the grid template in array templates, if template exists.
-!>    -1, otherwise.
-!>
-!>    @author Stephen Gilbert @date 2001-06-28
-         integer function getgridindex(number)
-           implicit none
-           integer,intent(in) :: number
-           integer :: j
-
-           getgridindex=-1
-
-           do j=1,MAXTEMP
-              if (number.eq.templates(j)%template_num) then
-                 getgridindex=j
-                 return
-              endif
-           enddo
-
-         end function
-
-!>    Get the grid template information for a specified Grid 
-!>    Definition Template.
-!>
-!>    The number of entries in the template is returned along with a map
-!>    of the number of octets occupied by each entry. Also, a flag is
-!>    returned to indicate whether the template would need to be
-!>    extended.
-!>
-!>    @param[in] number NN, indicating the number of the Grid Definition
-!>    Template that is being requested.
-!>    @param[out] nummap Number of entries in the Template.
-!>    @param[out] map An array containing the number of octets that each
-!>    template entry occupies when packed up into the GDS.
-!>    @param[out] needext Logical variable indicating whether the Grid
-!>    Defintion Template has to be extended.
-!>    @param[out] iret Error return code.
-!>    - 0 no error.
-!>    - 1 Undefine Grid Template number.
-!>
-!>    @author Stephen Gilbert @date 2000-05-09
-         subroutine getgridtemplate(number,nummap,map,needext,iret)
-           implicit none
-
-           integer,intent(in) :: number
-           integer,intent(out) :: nummap,map(*),iret
-           logical,intent(out) :: needext
-           integer :: index
-
-           iret=0
-
-           index=getgridindex(number)
-
-           if (index.ne.-1) then
-              nummap=templates(index)%mapgridlen
-              needext=templates(index)%needext
-              map(1:nummap)=templates(index)%mapgrid(1:nummap)
-           else
-             nummap=0
-             needext=.false.
-             print *,'getgridtemplate: Grid Template ',number,
-     &               ' not defined.'
-             iret=1
-           endif
-
-         end subroutine
-
-!>    Generate the remaining octet map for a given Grid
-!>    Definition Template, if required.
-!>    
-!>    Some Templates can vary
-!>    depending on data values given in an earlier part of the Template,
-!>    and it is necessary to know some of the earlier entry values to
-!>    generate the full octet map of the Template.
-!>    
-!>    @param[in] number NN, indicating the number of the Grid Definition
-!>    Template 3.NN that is being requested.
-!>    @param[in] list The list of values for each entry in the Grid
-!>    Definition Template.
-!>    @param[out] nummap Number of entries in the Template.
-!>    @param[out] map An array containing the number of octets that each
-!>    template entry occupies when packed up into the GDS.
-!>
-!>    @author Stephen Gilbert @date 2000-05-09
-         subroutine extgridtemplate(number,list,nummap,map)
-           implicit none
-
-           integer,intent(in) :: number,list(*)
-           integer,intent(out) :: nummap,map(*)
-           integer i, index, n, ni, nj
-
-           index=getgridindex(number)
-           if (index.eq.-1) return
-
-           if ( .not. templates(index)%needext ) return
-           nummap=templates(index)%mapgridlen
-           map(1:nummap)=templates(index)%mapgrid(1:nummap)
-
-           if ( number.eq.120 ) then
-              N=list(2)
-              do i=1,N
-                map(nummap+1)=2
-                map(nummap+2)=-2
-                nummap=nummap+2
-              enddo
-           elseif ( number.eq.4 ) then
-              NI=list(8)
-              do i=1,NI
-                map(nummap+1)=4
-                nummap=nummap+1
-              enddo
-              NJ=list(9)
-              do i=1,NJ
-                map(nummap+1)=-4
-                nummap=nummap+1
-              enddo
-           elseif ( number.eq.5 ) then
-              NI=list(8)
-              do i=1,NI
-                map(nummap+1)=4
-                nummap=nummap+1
-              enddo
-              NJ=list(9)
-              do i=1,NJ
-                map(nummap+1)=-4
-                nummap=nummap+1
-              enddo
-           elseif ( number.eq.1000 ) then
-              N=list(20)
-              do i=1,N
-                map(nummap+1)=4
-                nummap=nummap+1
-              enddo
-           elseif ( number.eq.1200 ) then
-              N=list(16)
-              do i=1,N
-                map(nummap+1)=4
-                nummap=nummap+1
-              enddo
-           endif
-
-         end subroutine
-
-!>    This function returns the initial length (number of entries) in
-!>    the static part of specified Grid Definition Template.
-!>
-!>    @param[in] number NN, indicating the number of the Grid Definition
-!>    Template that is being requested.
-!>
-!>    @return Number of entries in the static part of the grid definition template,
-!>    or 0, if requested template is not found.
-!>
-!>    @author Stephen Gilbert @date 2004-05-11
-         integer function getgdtlen(number)
-           implicit none
-
-           integer,intent(in) :: number
-           integer :: index
-
-           getgdtlen=0
-
-           index=getgridindex(number)
-
-           if (index.ne.-1) then
-              getgdtlen=templates(index)%mapgridlen
-           endif
-
-           end function
-      end
-

From 4112d2636ee93596876412f48a332c40fc035ce6 Mon Sep 17 00:00:00 2001
From: Edward Hartnett <edward.hartnett@noaa.gov>
Date: Tue, 8 Aug 2023 03:14:22 -0600
Subject: [PATCH 3/3] converted gridtemplates to F90

---
 src/gridtemplates.F90 | 105 +++++++++++++++++++++---------------------
 1 file changed, 52 insertions(+), 53 deletions(-)

diff --git a/src/gridtemplates.F90 b/src/gridtemplates.F90
index 13a34c25..8e94cb61 100644
--- a/src/gridtemplates.F90
+++ b/src/gridtemplates.F90
@@ -54,68 +54,68 @@ module gridtemplates
   data templates(1)%template_num /0/     !  Lat/Lon
   data templates(1)%mapgridlen /19/
   data templates(1)%needext /.false./
-  data (templates(1)%mapgrid(j), j = 1, 19)
-  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1/
+  data (templates(1)%mapgrid(j), j = 1, 19) &
+       /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1/
 
   data templates(2)%template_num /1/     !  Rotated Lat/Lon
   data templates(2)%mapgridlen /22/
   data templates(2)%needext /.false./
-  data (templates(2)%mapgrid(j), j = 1, 22)
-  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1, -4, 4, 4/
+  data (templates(2)%mapgrid(j), j = 1, 22) &
+       /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1, -4, 4, 4/
 
   data templates(3)%template_num /2/     !  Stretched Lat/Lon
   data templates(3)%mapgridlen /22/
   data templates(3)%needext /.false./
-  data (templates(3)%mapgrid(j), j = 1, 22)
-  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1, -4, 4, -4/
+  data (templates(3)%mapgrid(j), j = 1, 22) &
+       /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1, -4, 4, -4/
 
-  data templates(4)%template_num /3/     !  Stretched & Rotated Lat/Lon
+  data templates(4)%template_num /3/     !  Stretch and Rotated Lat/Lon
   data templates(4)%mapgridlen /25/
   data templates(4)%needext /.false./
-  data (templates(4)%mapgrid(j), j = 1, 25)
-  &       /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1, -4, 4, 4, -4, 4, -4/
+  data (templates(4)%mapgrid(j), j = 1, 25) &
+       /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1, -4, 4, 4, -4, 4, -4/
 
   data templates(5)%template_num /10/     !  Mercator
   data templates(5)%mapgridlen /19/
   data templates(5)%needext /.false./
-  data (templates(5)%mapgrid(j), j = 1, 19)
-  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, -4, 4, 1, -4, -4, 4, 1, 4, 4, 4/
+  data (templates(5)%mapgrid(j), j = 1, 19) &
+       /1, 1, 4, 1, 4, 1, 4, 4, 4, -4, 4, 1, -4, -4, 4, 1, 4, 4, 4/
 
   data templates(6)%template_num /20/     !  Polar Stereographic
   data templates(6)%mapgridlen /18/
   data templates(6)%needext /.false./
-  data (templates(6)%mapgrid(j), j = 1, 18)
-  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1, 1/
+  data (templates(6)%mapgrid(j), j = 1, 18) &
+       /1, 1, 4, 1, 4, 1, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1, 1/
 
   data templates(7)%template_num /30/     !  Lambert Conformal
   data templates(7)%mapgridlen /22/
   data templates(7)%needext /.false./
-  data (templates(7)%mapgrid(j), j = 1, 22)
-  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1, 1, -4, -4, -4, 4/
+  data (templates(7)%mapgrid(j), j = 1, 22) &
+       /1, 1, 4, 1, 4, 1, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1, 1, -4, -4, -4, 4/
 
   data templates(8)%template_num /40/     !  Gaussian Lat/Lon
   data templates(8)%mapgridlen /19/
   data templates(8)%needext /.false./
-  data (templates(8)%mapgrid(j), j = 1, 19)
-  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1/
+  data (templates(8)%mapgrid(j), j = 1, 19) &
+       /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1/
 
   data templates(9)%template_num /41/     !  Rotated Gaussian Lat/Lon
   data templates(9)%mapgridlen /22/
   data templates(9)%needext /.false./
-  data (templates(9)%mapgrid(j), j = 1, 22)
-  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1, -4, 4, 4/
+  data (templates(9)%mapgrid(j), j = 1, 22) &
+       /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1, -4, 4, 4/
 
   data templates(10)%template_num /42/     !  Stretched Gaussian Lat/Lon
   data templates(10)%mapgridlen /22/
   data templates(10)%needext /.false./
-  data (templates(10)%mapgrid(j), j = 1, 22)
-  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1, -4, 4, -4/
+  data (templates(10)%mapgrid(j), j = 1, 22) &
+       /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1, -4, 4, -4/
 
   data templates(11)%template_num /43/     !  Strtchd and Rot'd Gaus Lat/Lon
   data templates(11)%mapgridlen /25/
   data templates(11)%needext /.false./
-  data (templates(11)%mapgrid(j), j = 1, 25)
-  &          /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1, -4, 4, 4, -4, 4, -4/
+  data (templates(11)%mapgrid(j), j = 1, 25) &
+       /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1, -4, 4, 4, -4, 4, -4/
 
   data templates(12)%template_num /50/    !  Spherical Harmonic Coefficients
   data templates(12)%mapgridlen /5/
@@ -140,8 +140,8 @@ module gridtemplates
   data templates(16)%template_num /90/     !  Space view Perspective
   data templates(16)%mapgridlen /21/
   data templates(16)%needext /.false./
-  data (templates(16)%mapgrid(j), j = 1, 21)
-  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, -4, 4, 1, 4, 4, 4, 4, 1, 4, 4, 4, 4/
+  data (templates(16)%mapgrid(j), j = 1, 21) &
+       /1, 1, 4, 1, 4, 1, 4, 4, 4, -4, 4, 1, 4, 4, 4, 4, 1, 4, 4, 4, 4/
 
   data templates(17)%template_num /100/    !  Triangular grid (icosahedron)
   data templates(17)%mapgridlen /11/
@@ -151,8 +151,8 @@ module gridtemplates
   data templates(18)%template_num /110/ !  Equatorial Azimuthal equidistant
   data templates(18)%mapgridlen /16/
   data templates(18)%needext /.false./
-  data (templates(18)%mapgrid(j), j = 1, 16)
-  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, -4, 4, 1, 4, 4, 1, 1/
+  data (templates(18)%mapgrid(j), j = 1, 16) &
+       /1, 1, 4, 1, 4, 1, 4, 4, 4, -4, 4, 1, 4, 4, 1, 1/
 
   data templates(19)%template_num /120/     !  Azimuth-range
   data templates(19)%mapgridlen /7/
@@ -162,84 +162,84 @@ module gridtemplates
   data templates(20)%template_num /1000/     !  Cross Section Grid
   data templates(20)%mapgridlen /20/
   data templates(20)%needext /.true./
-  data (templates(20)%mapgrid(j), j = 1, 20)
-  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, -4, 4, 1, 4, 4, 1, 2, 1, 1, 2/
+  data (templates(20)%mapgrid(j), j = 1, 20) &
+       /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, -4, 4, 1, 4, 4, 1, 2, 1, 1, 2/
 
   data templates(21)%template_num /1100/     !  Hovmoller Diagram Grid
   data templates(21)%mapgridlen /28/
   data templates(21)%needext /.false./
-  data (templates(21)%mapgrid(j), j = 1, 28)
-  &    /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, -4, 4, 1, -4, 4, 1, 4, 1, -4, 1, 1, -4, 2, 1, 1, 1, 1, 1/
+  data (templates(21)%mapgrid(j), j = 1, 28) &
+       /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, -4, 4, 1, -4, 4, 1, 4, 1, -4, 1, 1, -4, 2, 1, 1, 1, 1, 1/
 
   data templates(22)%template_num /1200/     !  Time Section Grid
   data templates(22)%mapgridlen /16/
   data templates(22)%needext /.true./
-  data (templates(22)%mapgrid(j), j = 1, 16)
-  &              /4, 1, -4, 1, 1, -4, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2/
+  data (templates(22)%mapgrid(j), j = 1, 16) &
+       /4, 1, -4, 1, 1, -4, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2/
 
   data templates(23)%template_num /31/     !  Albers Equal Area
   data templates(23)%mapgridlen /22/
   data templates(23)%needext /.false./
-  data (templates(23)%mapgrid(j), j = 1, 22)
-  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1, 1, -4, -4, -4, 4/
+  data (templates(23)%mapgrid(j), j = 1, 22) &
+       /1, 1, 4, 1, 4, 1, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1, 1, -4, -4, -4, 4/
 
   data templates(24)%template_num /204/     !  Curilinear Orthogonal Grids
   data templates(24)%mapgridlen /19/
   data templates(24)%needext /.false./
-  data (templates(24)%mapgrid(j), j = 1, 19)
-  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1/
+  data (templates(24)%mapgrid(j), j = 1, 19) &
+       /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1/
 
   data templates(25)%template_num /32768/   ! Rotate Lat/Lon E-grid
   data templates(25)%mapgridlen /19/
   data templates(25)%needext /.false./
-  data (templates(25)%mapgrid(j), j = 1, 19)
-  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1/
+  data (templates(25)%mapgrid(j), j = 1, 19) &
+       /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1/
 
   data templates(26)%template_num /32769/   ! Rotate Lat/Lon Non-E Stagger grid
   data templates(26)%mapgridlen /21/
   data templates(26)%needext /.false./
-  data (templates(26)%mapgrid(j), j = 1, 21)
-  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1, 4, 4/
+  data (templates(26)%mapgrid(j), j = 1, 21) &
+       /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, -4, 4, 1, -4, 4, 4, 4, 1, 4, 4/
   !
   !     GDT 3.4 Added (08/05/2013)
   !
   data templates(27)%template_num /4/     !  Variable resolution Latitude/Longitude
   data templates(27)%mapgridlen /13/
   data templates(27)%needext /.true./
-  data (templates(27)%mapgrid(j), j = 1, 13)
-  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, 1, 1/
+  data (templates(27)%mapgrid(j), j = 1, 13) &
+       /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, 1, 1/
   !
   !     GDT 3.5 Added (08/05/2013)
   !
   data templates(28)%template_num /5/     !  Variable resolution rotate Latitude/Longitude
   data templates(28)%mapgridlen /16/
   data templates(28)%needext /.true./
-  data (templates(28)%mapgrid(j), j = 1, 16)
-  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, 1, 1, -4, 4, 4/
+  data (templates(28)%mapgrid(j), j = 1, 16) &
+       /1, 1, 4, 1, 4, 1, 4, 4, 4, 4, 4, 1, 1, -4, 4, 4/
   !
   !     GDT 3.12 Added (08/05/2013)
   !
   data templates(29)%template_num /12/     !  Transverse Mercator
   data templates(29)%mapgridlen /22/
   data templates(29)%needext /.false./
-  data (templates(29)%mapgrid(j), j = 1, 22)
-  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, -4, 4, 1, -4, 4, 4, 1, 4, 4, -4, -4, -4, -4/
+  data (templates(29)%mapgrid(j), j = 1, 22) &
+       /1, 1, 4, 1, 4, 1, 4, 4, 4, -4, 4, 1, -4, 4, 4, 1, 4, 4, -4, -4, -4, -4/
   !
   !     GDT 3.101 Added (08/05/2013)
   !
   data templates(30)%template_num /101/     !  General unstructured grid
   data templates(30)%mapgridlen /4/
   data templates(30)%needext /.false./
-  data (templates(30)%mapgrid(j), j = 1, 4)
-  &              /1, 4, 1, -4/
+  data (templates(30)%mapgrid(j), j = 1, 4) &
+       /1, 4, 1, -4/
   !
   !     GDT 3.140 Added (08/05/2013)
   !
   data templates(31)%template_num /140/     !  Lambert Azimuthal Equal Area Projection
   data templates(31)%mapgridlen /17/
   data templates(31)%needext /.false./
-  data (templates(31)%mapgrid(j), j = 1, 17)
-  &              /1, 1, 4, 1, 4, 1, 4, 4, 4, -4, 4, 4, 4, 1, 4, 4, 1/
+  data (templates(31)%mapgrid(j), j = 1, 17) &
+       /1, 1, 4, 1, 4, 1, 4, 4, 4, -4, 4, 4, 4, 1, 4, 4, 1/
 
 contains
 
@@ -307,8 +307,7 @@ subroutine getgridtemplate(number, nummap, map, needext, iret)
     else
        nummap = 0
        needext = .false.
-       print *, 'getgridtemplate: Grid Template ', number, 
-       &               ' not defined.'
+       print *, 'getgridtemplate: Grid Template ', number, ' not defined.'
        iret = 1
     endif
   end subroutine getgridtemplate