Add unit testing

src/aea.f

@@ -2,19 +2,6 @@
 C> @brief This subroutine converts ascii to ebcdic, or ebcdic to ascii
 C> @author desmarais @date 11-29-1982
 
-C>   Program history log:
-C>   - 11-29-1982  Desmarais
-C>   - 03-31-1988  R. E. Jones
-C>      - change logic so it works like a ibm370 translate instruction.
-C>   - 08-22-1988  R. E. Jones
-C>      - changes for microsoft fortran 4.10
-C>   - 09-04-1988  R. E. Jones
-C>      - change tables to 128 character set
-C>   - 01-31-1990  R. E. Jones
-C>      - convert to cray cft77 fortran cray does not allow char*1 to be set to hex
-C>   - 12-21-1998  Stephen Gilbert
-C>      - replaced function ichar with mova2i.
-C>
 C>   @param[in, out] IA character*1 array of ascii data  if nc < 0
 C>   @param[in, out] IE character*1 array of ebcdic data if nc > 0
 C>   @param[in] NC integer, contains character count to convert.

src/errexit.f

@@ -4,10 +4,6 @@
 
 C> Exit with a return code.
 C>
-C> Program history log:
-C> - 1998-06-04  Mark Iredell
-C> - 1999-01-26  Stephen Gilbert
-C>  - Changed to use XLF utility routine exit_(n) instead of exit(n).
 C> exit_(n) will return the proper value (n must be 4 byte int)
 C> to the sh/ksh shell status variable $? ($status for csh)
 C> on the IBM SP.

src/errmsg.f

@@ -4,9 +4,6 @@
 
 C> Write a message to stderr.
 C>
-C> Program history log:
-C>  - 1995-10-31  Mark Iredell
-C>
 C> @param[in] CMSG character*(*) message to write.
 C>
 C> @note This is a machine-dependent subprogram for Cray.

src/fparsei.f

@@ -5,9 +5,6 @@
 C> This subprogram extracts integers from a free-format
 C> character string. It is useful for parsing command arguments.
 C>
-C> Program history log:
-C> - 1998-09-03  Mark Iredell
-C>
 C> @param[in] CARG character*(*) string of ascii digits to parse.
 C> Integers may be separated by a comma or by blanks.
 C> @param[in] MARG integer maximum number of integers to parse.

src/fparser.f

@@ -5,9 +5,6 @@
 C> This subprogram extracts real numbers from a free-format
 C> character string. It is useful for parsing command arguments.
 C>
-C> Program history log:
-C> - 1998-09-03  Mark Iredell
-C>
 C> @param[in] CARG character*(*) string of ascii digits to parse.
 C> Real numbers may be separated by a comma or by blanks.
 C> @param[in] MARG integer maximum number of real numbers to parse.

src/gbyte.f

@@ -4,11 +4,6 @@
 
 C> This is the fortran version of gbyte
 C>
-C> Program history log:
-C> - Russell E. Jones 1991-03-DD
-C> Changes for SiliconGraphics IRIS-4D/25
-C> SiliconGraphics 3.3 FORTRAN 77.
-C>
 C> To unpack a byte into a target word. The unpacked byte is right-justified
 C> in the target word, and the remainder of the word is zero-filled.
 C>

src/gbytes.f

@@ -2,11 +2,6 @@
 C> @brief This is the fortran version of gbytes.
 C> @author Dr. Robert C. Gammill @date 1972-05
 
-C> Program history log:
-C> - Russell E. Jones 1991-03
-C> Changes for SiliconGraphics IRIS-4D/25
-C> SiliconGraphics 3.3 FORTRAN 77
-C>
 C> To unpack a series of bytes into a target
 C> array. Each unpacked byte is right-justified
 C> in its target word, and the remainder of the

src/getbit.f

@@ -8,9 +8,6 @@
 C> For particular binary and decimal scalings is computed.
 C> Grib bitmap masking for valid data is optionally used.
 C>
-C> Program history log:
-C> - Mark Iredell 1996-09-16
-C>
 C> @param[in] IBM Integer bitmap flag (=0 for no bitmap).
 C> @param[in] IBS Integer binary scaling (e.g. ibs=3 to round field
 C> to nearest eighth value).

src/getgb.f

@@ -16,11 +16,6 @@
 C> and the unpacked data. If the grib message is not found, then the
 C> return code will be nonzero.
 C>
-C> Program history log:
-C> - Mark Iredell 1994-04-01
-C> - Mark Iredell 1995-10-31 modularized portions of code into
-C> subprograms and allowed for unspecified index file
-C>
 C> @param[in] LUGB Integer unit of the unblocked grib data file
 C> @param[in] LUGI Integer unit of the unblocked grib index file
 C> (=0 to get index buffer from the grib file)

src/getgb1.f

@@ -14,11 +14,6 @@
 C> and the unpacked data. If the grib message is not found, then the
 C> return code will be nonzero.
 C>
-C> Program history log:
-C>  - Mark Iredell 1994-04-01
-C>  - Ralph Jones 1995-05-10  Add one more parameter to getgb and
-C> change name to getgb1.
-C>
 C> @param[in] lugb logical unit of the unblocked grib data file.
 C> @param[in] lugi logical unit of the unblocked grib index file.
 C> @param[in] jf integer maximum number of data points to unpack.

src/getgb1r.f

@@ -2,10 +2,6 @@
 C> @brief Reads and unpacks a grib message.
 C> @author Mark Iredell @date 1995-10-31
 
-C> Program history log:
-C> - Mark Iredell 1995-10-31
-C> - Chuang 2004-07-22 Add packing bit number nbitss in the argument
-C> list because eta grib files need it to repack grib file.
 C> @param[in] LUGB Integer unit of the unblocked grib data file.
 C> @param[in] LSKIP Integer number of bytes to skip.
 C> @param[in] LGRIB Integer number of bytes to read.

src/getgb1re.f

@@ -4,10 +4,6 @@
 
 C> Reads and unpacks a grib message.
 C>
-C> PROGRAM HISTORY LOG:
-C>   Mark Iredell 1995-10-31
-C>   Y. Zhu 1997-02-11 Included probability and cluster arguments.
-C>
 C> @param[in] LUGB Integer unit of the unblocked grib data file.
 C> @param[in] LSKIP Integer number of bytes to skip.
 C> @param[in] LGRIB Integer number of bytes to read.

src/getgb1s.f

@@ -8,10 +8,6 @@
 c> and the unpacked pds and gds parameters. (A requested parameter
 c> of -1 means to allow any value of this parameter to be found.)
 C>
-C> Program history log:
-C> - Mark Iredell 1995-10-31
-C> - Mark Iredell 2001-06-05 Apply linux port by ebisuzaki.
-C>
 C> @param[in] CBUF Character*1 (nlen*nnum) buffer containing index data.
 C> @param[in] NLEN Integer length of each index record in bytes.
 C> @param[in] NNUM Integer number of index records.

src/getgbe.f

@@ -16,11 +16,6 @@
 C> and the unpacked data. If the grib message is not found, then the
 C> return code will be nonzero.
 C>
-C> Program history log:
-C> - Mark Iredell 1994-04-01
-C> - Mark Iredell 1995-10-31 Modularized portions of code into subprograms
-C> and allowed for unspecified index file.
-C>
 C> @param[in] lugb Integer unit of the unblocked grib data file.
 C> @param[in] lugi Integer unit of the unblocked grib index file
 C> (=0 to get index buffer from the grib file).

src/getgbeh.f

@@ -13,10 +13,6 @@
 C> returned along with the unpacked pds and gds parameters. If the
 C> grib message is not found, then the return code will be nonzero.
 C>
-C> Program History:
-C> - 1995-10-31  Mark Iredell Modularized portions of code into subprograms
-C> and allowed for unspecified index file.
-C>
 C> @param[in] LUGB Integer unit of the unblocked grib data file.
 C> (only used if lugi=0)
 C> @param[in] LUGI Integer unit of the unblocked grib index file.

src/getgbem.f

@@ -15,11 +15,6 @@
 C> and the unpacked data. if the grib message is not found, then the
 C> return code will be nonzero.
 C>
-C> Program history log:
-C> - Mark Iredell 1994-04-01
-C> - Mark Iredell 1995-10-31 Modularized portions of code into subprograms
-C> and allowed for unspecified index file.
-C>
 C> @param[in] LUGB Integer unit of the unblocked grib data file.
 C> @param[in] LUGI Integer unit of the unblocked grib index file
 C> (=0 to get index buffer from the grib file).

src/getgbemh.f

@@ -12,11 +12,6 @@
 C> returned along with the unpacked pds and gds parameters. If the
 C> grib message is not found, then the return code will be nonzero.
 C>
-C> Program history log:
-C> - Mark Iredell 1994-04-01
-C> - Mark Iredell 1995-10-31 Modularized portions of code into subprograms
-C> and allowed for unspecified index file
-C>
 C> @param[in] lugb integer unit of the unblocked grib data file
 C> (only used if lugi=0).
 C> @param[in] lugi integer unit of the unblocked grib index file

src/getgbemn.f

@@ -15,13 +15,6 @@
 C> and the unpacked data. If the grib message is not found, then the
 C> return code will be nonzero.
 C>
-C> Program history log:
-C>  - Mark Iredell 1994-04-01
-C>  - Mark Iredell 1995-10-31  Modularized portions of code into subprograms
-C>  and allowed for unspecified index file
-C>  - Chuang 2004-07-22 add packing bit number nbitss in the argument
-C>  list because eta grib files need it to repack grib file
-C>
 C> @param[in] lugb integer unit of the unblocked grib data file.
 C> @param[in] lugi integer unit of the unblocked grib index file
 C> (=0 to get index buffer from the grib file).

src/getgbemp.f

@@ -13,11 +13,6 @@
 C> pds and gds parameters and the packed grib message. If the grib
 C> message is not found, then the return code will be nonzero.
 C>
-C> Program history log:
-C>  - Mark Iredell 1994-04-01
-C>  - Mark Iredell 1995-10-31 Modularized portions of code into subprograms
-C>  and allowed for unspecified index file.
-C>
 C> @param[in] lugb integer unit of the unblocked grib data file.
 C> @param[in] lugi integer unit of the unblocked grib index file
 C> (=0 to get index buffer from the grib file).

src/getgbens.f

@@ -16,11 +16,6 @@
 C> return code will be nonzero.
 C> This obsolescent version has been replaced by getgbe.
 C>
-C> Program history log:
-C>  - Mark Iredell 1994-04-01
-C>  - Mark Iredell 1995-10-31 Modularized portions of code into subprograms
-C>  and allowed for unspecified index file.
-C>
 C> @param[in] lugb integer unit of the unblocked grib data file.
 C> @param[in] lugi integer unit of the unblocked grib index file
 C> (=0 to get index buffer from the grib file).

src/getgbep.f

@@ -14,11 +14,6 @@
 C> pds and gds parameters and the packed grib message. If the grib
 C> message is not found, then the return code will be nonzero.
 C>
-C> Program history log:
-C> - Mark Iredell 1994-04-01
-C> - Mark Iredell 1995-10-31 Modularized portions of code into subprograms
-C> and allowed for unspecified index file.
-C>
 C> @param[in] lugb integer unit of the unblocked grib data file.
 C> @param[in] lugi integer unit of the unblocked grib index file
 C> (=0 to get index buffer from the grib file).

src/getgbex.f

@@ -15,12 +15,6 @@
 C> and the unpacked data. If the grib message is not found, then the
 C> return code will be nonzero.
 C>
-C> Program history log:
-C> - Mark Iredell 1994-04-01
-C> - Mark Iredell 1995-10-31 Modularized portions of code into subprograms
-C> and allowed for unspecified index file.
-C> - Y. Zhu 1997-02-11 Included probability and cluster arguments.
-C>
 C> @param[in] lugb integer unit of the unblocked grib data file.
 C> @param[in] lugi integer unit of the unblocked grib index file
 C> (=0 to get index buffer from the grib file).

src/getgbexm.f

@@ -14,12 +14,6 @@
 C> and the unpacked data. If the grib message is not found, then the
 C> return code will be nonzero.
 C>
-C> Program history log:
-C> - Mark Iredell 1994-04-01
-C> - Mark Iredell 1995-10-31 Modularized portions of code into subprograms
-C> and allowed for unspecified index file.
-C> - Y. Zhu 1997-02-11 Included probability and cluster arguments.
-C>
 C> @param[in] lugb integer unit of the unblocked grib data file.
 C> @param[in] lugi integer unit of the unblocked grib index file
 C> (=0 to get index buffer from the grib file).

src/getgbh.f

@@ -13,11 +13,6 @@
 C> returned along with the unpacked pds and gds parameters. If the
 C> grib message is not found, then the return code will be nonzero.
 C>
-C> Program history log:
-C> - Mark Iredell 1994-04-01
-C> - Mark Iredell 1995-10-31 Modularized portions of code into subprograms
-C> and allowed for unspecified index file.
-C>
 C> @param[in] lugb integer unit of the unblocked grib data file
 C> (only used if lugi=0).
 C> @param[in] lugi integer unit of the unblocked grib index file

src/getgbm.f

@@ -14,14 +14,6 @@
 C> and the unpacked data. If the grib message is not found, then the
 C> return code will be nonzero.
 C>
-C> Program history log:
-C> - Mark Iredell 1994-04-01
-C> - Mark Iredell 1995-10-31  Modularized portions of code into subprograms
-C> and allowed for unspecified index file.
-C> - Chuang 2004-07-22  Add nbitss to the argument list of getgb1r that
-C> is called in this subroutine.
-C> - Wang 2010-03-02  wang Increase msk1 to 256000000 for nemsio files.
-C>
 C> @param[in] lugb integer unit of the unblocked grib data file.
 C> @param[in] lugi integer unit of the unblocked grib index file
 C> (=0 to get index buffer from the grib file).

src/getgbmh.f

@@ -12,11 +12,6 @@
 C> returned along with the unpacked pds and gds parameters. If the
 C> grib message is not found, then the return code will be nonzero.
 C>
-C> Program history log:
-C> - Mark Iredell 1994-04-01
-C> - Mark Iredell 1995-10-31 Modularized portions of code into subprograms
-C> and allowed for unspecified index file.
-C>
 C> @param[in] lugb integer unit of the unblocked grib data file
 C> (only used if lugi=0).
 C> @param[in] lugi integer unit of the unblocked grib index file

src/getgbmp.f

@@ -13,11 +13,6 @@
 C> pds and gds parameters and the packed grib message. If the grib
 C> message is not found, then the return code will be nonzero.
 C>
-C> Program history log:
-C> - Mark Iredell 1994-04-01  iredell
-C> - Mark Iredell 1995-10-31  iredell Modularized portions of code into subprograms
-C> and allowed for unspecified index file.
-C>
 C> @param[in] lugb integer unit of the unblocked grib data file.
 C> @param[in] lugi integer unit of the unblocked grib index file
 C> (=0 to get index buffer from the grib file).

src/getgbp.f

@@ -14,11 +14,6 @@
 C> pds and gds parameters and the packed grib message. If the grib
 C> message is not found, then the return code will be nonzero.
 C>
-C> Program history log:
-C> - Mark Iredell 1994-04-01
-C> - Mark Iredell 1995-10-31 Modularized portions of code into subprograms
-C> and allowed for unspecified index file.
-C>
 C> @param[in] lugb integer unit of the unblocked grib data file.
 C> @param[in] lugi integer unit of the unblocked grib index file
 C> (=0 to get index buffer from the grib file).

src/getgi.f

@@ -25,10 +25,6 @@
 C> - byte 173-184: (optional) bytes 29-40 of the pds.
 C> - byte 185-320: (optional) bytes 43-178 of the gds.
 C>
-C> Program history log:
-C> - Mark Iredell  1995-10-31
-C> - Mark Iredell  1996-10-31 Augmented optional definitions to byte 320.
-C>
 C> @param[in] lugi integer unit of the unblocked grib index file.
 C> @param[in] mnum integer number of index records to skip (usually 0).
 C> @param[in] mbuf integer length of cbuf in bytes.

src/getgir.f

@@ -19,10 +19,6 @@
 C> - byte 173-184: (optional) bytes 29-40 of the pds.
 C> - byte 185-320: (optional) bytes 43-178 of the gds.
 C>
-C> Program history log:
-C> - Mark Iredell 1995-10-31
-C> - Mark Iredell 1996-10-31 Augmented optional definitions to byte 320.
-C>
 C> @param[in] lugb integer unit of the unblocked grib file.
 C> @param[in] msk1 integer number of bytes to search for first message.
 C> @param[in] msk2 integer number of bytes to search for other messages.

src/gtbits.f

@@ -8,9 +8,6 @@
 C> The minimum and maximum rounded field values are also returned.
 C> Grib bitmap masking for valid data is optionally used.
 C>
-C> Program history log:
-C> - Mark Iredell 1992-10-31
-C>
 C> @param[in] ibm integer bitmap flag (=0 for no bitmap).
 c> @param[in] ids integer decimal scaling
 c> (e.g. ids=3 to round field to nearest milli-value).

src/idsdef.f

@@ -5,9 +5,6 @@
 C> Sets decimal scalings defaults for various parameters.
 C> A decimal scaling of -3 means data is packed in kilo-si units.
 C>
-C> Program history log:
-C> - Mark Iredell 1992-10-31
-C>
 C> @param[in] IPTV parameter table version (only 1 or 2 is recognized).
 C> @param[out] IDS integer (255) decimal scalings
 C> (unknown decimal scalings will not be set).

src/instrument.f

@@ -7,12 +7,6 @@
 !> of a program is invoked as well as the minimum, maximum
 !> and total wall-clock time spent in the given section.
 !>
-!> Program history log:
-!> - Mark Iredell 1998-07-16
-!> - Frimel and Kalina 2019-09-17
-!> Decompose return statistcs if statement
-!> - Boi Vuong 2020-04-02 Check for ka > 0 and modifiled ifblock statement
-!> into two separate ifblock statements.
 !> @param[in] K Integer positive section number
 !> or maximum section number in the first invocation
 !> or zero to reset all wall-clock statistics