Technical validation for GRIB Edition 2 product definition templates 67 and 68

IPET-DRMM-PFC-II_distribution.docx

pdt67.grib2

pdt68.grib2

Validation.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Technical validation: GRIB product definition templates 67 and 68\n",
    "\n",
    "This validation uses a modified version of wgrib2, available [here](https://github.com/erget/wgrib2/commit/06df9d798ae85649d12b1f9856fc80d37205ad45), and tables for ecCodes which have been updated as found [here](https://github.com/erget/eccodes/commit/ceb7d9800df2cbe1c13748456d7954dba87b13d8). Pull requests will be sent to the repository maintainers upon acceptance of the new templates.\n",
    "\n",
    "Only the GRIB tables, not the ecCodes version itself, was updated - no changes to the source codes were necessary."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\r\n",
      "eccodes Version 0.13.0\r\n",
      "\r\n",
      "Definition files path from environment variable ECCODES_DEFINITION_PATH=/home/dlee/.local/share/eccodes/definitions.edzw:/home/dlee/.local/share/eccodes/definitions\r\n",
      "\r\n",
      "Default SAMPLES path is used: /home/dlee/.local/share/eccodes/samples\r\n",
      "SAMPLES path can be changed setting ECCODES_SAMPLES_PATH environment variable\r\n",
      "\r\n"
     ]
    }
   ],
   "source": [
    "codes_info"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "/home/dlee/.local/share/eccodes/definitions.edzw:/home/dlee/.local/share/eccodes/definitions\n",
      "/home/dlee/.local/share/eccodes/definitions"
     ]
    }
   ],
   "source": [
    "codes_info -d\n",
    "defs=$(codes_info -d | cut -d \":\" -f 2)\n",
    "echo $defs"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "pdt67.grib2 pdt68.grib2\r\n"
     ]
    }
   ],
   "source": [
    "validation_files=pdt6*.grib2\n",
    "echo $validation_files"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "ecCodes is able to dump the contents of both GRIBs."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "***** FILE: pdt67.grib2 \r\n",
      "#==============   MESSAGE 1 ( length=209 )                 ==============\r\n",
      "1-4       identifier = GRIB\r\n",
      "5-6       reserved = MISSING\r\n",
      "7         discipline = 0 [Meteorological products  (grib2/tables/19/0.0.table) ]\r\n",
      "8         editionNumber = 2\r\n",
      "9-16      totalLength = 209\r\n",
      "======================   SECTION_1 ( length=21, padding=0 )    ======================\r\n",
      "1-4       section1Length = 21\r\n",
      "5         numberOfSection = 1\r\n",
      "6-7       centre = 78 [Offenbach  (RSMC)  (grib1/0.table) ]\r\n",
      "8-9       subCentre = 255\r\n",
      "10        tablesVersion = 19 [Unknown code table entry (grib2/tables/1.0.table) ]\r\n",
      "11        localTablesVersion = 1 [Unknown code table entry (grib2/tables/19/1.1.table) ]\r\n",
      "12        significanceOfReferenceTime = 1 [Start of forecast  (grib2/tables/19/1.2.table) ]\r\n",
      "13-14     year = 2014\r\n",
      "15        month = 4\r\n",
      "16        day = 3\r\n",
      "17        hour = 0\r\n",
      "18        minute = 0\r\n",
      "19        second = 0\r\n",
      "20        productionStatusOfProcessedData = 2 [Research products  (grib2/tables/19/1.3.table) ]\r\n",
      "21        typeOfProcessedData = 1 [Forecast products  (grib2/tables/19/1.4.table , grib2/tables/local/edzw/1/1.4.table) ]\r\n",
      "======================   SECTION_3 ( length=35, padding=0 )    ======================\r\n",
      "1-4       section3Length = 35\r\n",
      "5         numberOfSection = 3\r\n",
      "6         sourceOfGridDefinition = 0 [Specified in Code table 3.1  (grib2/tables/19/3.0.table) ]\r\n",
      "7-10      numberOfDataPoints = 1\r\n",
      "11        numberOfOctectsForNumberOfPoints = 0\r\n",
      "12        interpretationOfNumberOfPoints = 0 [There is no appended list  (grib2/tables/19/3.11.table) ]\r\n",
      "13-14     gridDefinitionTemplateNumber = 101 [General unstructured grid  (grib2/tables/19/3.1.table) ]\r\n",
      "15        shapeOfTheEarth = 6 [Earth assumed spherical with radius of 6 371 229.0 m  (grib2/tables/19/3.2.table) ]\r\n",
      "16-18     numberOfGridUsed = 26\r\n",
      "19        numberOfGridInReference = 1\r\n",
      "20-35     uuidOfHGrid = 16 {\r\n",
      "                     a2, 7b, 8d, e6, 18, c4, 11, e4, 82, 0a, b5, b0, 98, c6, a5, c0\r\n",
      "                  } # bytes uuidOfHGrid \r\n",
      "======================   SECTION_4 ( length=101, padding=24 )   ======================\r\n",
      "1-4       section4Length = 101\r\n",
      "5         numberOfSection = 4\r\n",
      "6-7       NV = 0\r\n",
      "8-9       productDefinitionTemplateNumber = 67 [Average, accumulation and/or extreme values or other statistically processed values at a horizontal level or in a horizontal layer in a continuous or non-continuous time interval for atmospheric chemical constituents based on a distribution function (grib2/tables/19/4.0.table , grib2/tables/local/edzw/1/4.0.table) ]\r\n",
      "10        parameterCategory = 20 [Atmospheric chemical constituents  (grib2/tables/19/4.1.0.table , grib2/tables/local/edzw/1/4.1.0.table) ]\r\n",
      "11        parameterNumber = 6 [Surface dry deposition mass flux  (kg m-2 s-1)  (grib2/tables/19/4.2.0.20.table) ]\r\n",
      "12-13     constituentType = 62001 [Dust dry (grib2/tables/19/4.230.table) ]\r\n",
      "14-15     numberOfModeOfDistribution = 3\r\n",
      "16-17     modeNumber = 1\r\n",
      "18-19     typeOfDistributionFunction = 8 [No distribution function. The encoded variable is derived from variables characterized by type of distribution function of type no. 7 with fixed variance (p1)  (grib2/tables/19/4.240.table) ]\r\n",
      "20        numberOfDistributionFunctionParameters = 2\r\n",
      "21        scaleFactorOfDistributionFunctionParameter = 3\r\n",
      "22-25     scaledValueOfDistributionFunctionParameter = 1700\r\n",
      "26        scaleFactorOfDistributionFunctionParameter = 0\r\n",
      "27-30     scaledValueOfDistributionFunctionParameter = 2650\r\n",
      "31        typeOfGeneratingProcess = 2 [Forecast  (grib2/tables/19/4.3.table , grib2/tables/local/edzw/1/4.3.table) ]\r\n",
      "32        backgroundProcess = 0\r\n",
      "33        generatingProcessIdentifier = 1\r\n",
      "34-35     hoursAfterDataCutoff = 0\r\n",
      "36        minutesAfterDataCutoff = 0\r\n",
      "37        indicatorOfUnitOfTimeRange = 1 [Hour  (grib2/tables/19/4.4.table) ]\r\n",
      "38-41     forecastTime = 0\r\n",
      "42        typeOfFirstFixedSurface = 1 [Ground or water surface  (-)  (grib2/tables/19/4.5.table , grib2/tables/local/edzw/1/4.5.table) ]\r\n",
      "43        scaleFactorOfFirstFixedSurface = MISSING\r\n",
      "44-47     scaledValueOfFirstFixedSurface = MISSING\r\n",
      "48        typeOfSecondFixedSurface = 255 [Missing  (grib2/tables/19/4.5.table , grib2/tables/local/edzw/1/4.5.table) ]\r\n",
      "49        scaleFactorOfSecondFixedSurface = MISSING\r\n",
      "50-53     scaledValueOfSecondFixedSurface = MISSING\r\n",
      "54-55     yearOfEndOfOverallTimeInterval = 2014\r\n",
      "56        monthOfEndOfOverallTimeInterval = 4\r\n",
      "57        dayOfEndOfOverallTimeInterval = 5\r\n",
      "58        hourOfEndOfOverallTimeInterval = 0\r\n",
      "59        minuteOfEndOfOverallTimeInterval = 0\r\n",
      "60        secondOfEndOfOverallTimeInterval = 0\r\n",
      "61        numberOfTimeRange = 1\r\n",
      "62-65     numberOfMissingInStatisticalProcess = 0\r\n",
      "66        typeOfStatisticalProcessing = 1 [Accumulation  (grib2/tables/19/4.10.table) ]\r\n",
      "67        typeOfTimeIncrement = 2 [Successive times processed have same start time of forecast, forecast time is incremented  (grib2/tables/19/4.11.table , grib2/tables/local/edzw/1/4.11.table) ]\r\n",
      "68        indicatorOfUnitForTimeRange = 1 [Hour  (grib2/tables/19/4.4.table) ]\r\n",
      "69-72     lengthOfTimeRange = 48\r\n",
      "73        indicatorOfUnitForTimeIncrement = 255 [Missing  (grib2/tables/19/4.4.table) ]\r\n",
      "74-77     timeIncrement = 0\r\n",
      "======================   SECTION_5 ( length=21, padding=0 )    ======================\r\n",
      "1-4       section5Length = 21\r\n",
      "5         numberOfSection = 5\r\n",
      "6-9       numberOfValues = 1\r\n",
      "10-11     dataRepresentationTemplateNumber = 0 [Grid point data - simple packing  (grib2/tables/19/5.0.table) ]\r\n",
      "12-15     referenceValue = 999\r\n",
      "16-17     binaryScaleFactor = -73\r\n",
      "18-19     decimalScaleFactor = 0\r\n",
      "20        bitsPerValue = 0\r\n",
      "21        typeOfOriginalFieldValues = 0 [Floating point  (grib2/tables/19/5.1.table) ]\r\n",
      "======================   SECTION_6 ( length=6, padding=0 )     ======================\r\n",
      "1-4       section6Length = 6\r\n",
      "5         numberOfSection = 6\r\n",
      "6         bitMapIndicator = 255 [A bit map does not apply to this product  (grib2/tables/19/6.0.table) ]\r\n",
      "======================   SECTION_7 ( length=5, padding=0 )     ======================\r\n",
      "1-4       section7Length = 5\r\n",
      "5         numberOfSection = 7\r\n",
      "======================   SECTION_8 ( length=4, padding=0 )     ======================\r\n",
      "1-4       7777 = 7777\r\n",
      "#==============   MESSAGE 2 ( length=209 )                 ==============\r\n",
      "1-4       identifier = GRIB\r\n",
      "5-6       reserved = MISSING\r\n",
      "7         discipline = 0 [Meteorological products  (grib2/tables/19/0.0.table) ]\r\n",
      "8         editionNumber = 2\r\n",
      "9-16      totalLength = 209\r\n",
      "======================   SECTION_1 ( length=21, padding=0 )    ======================\r\n",
      "1-4       section1Length = 21\r\n",
      "5         numberOfSection = 1\r\n",
      "6-7       centre = 78 [Offenbach  (RSMC)  (grib1/0.table) ]\r\n",
      "8-9       subCentre = 255\r\n",
      "10        tablesVersion = 19 [Unknown code table entry (grib2/tables/1.0.table) ]\r\n",
      "11        localTablesVersion = 1 [Unknown code table entry (grib2/tables/19/1.1.table) ]\r\n",
      "12        significanceOfReferenceTime = 1 [Start of forecast  (grib2/tables/19/1.2.table) ]\r\n",
      "13-14     year = 2014\r\n",
      "15        month = 4\r\n",
      "16        day = 3\r\n",
      "17        hour = 0\r\n",
      "18        minute = 0\r\n",
      "19        second = 0\r\n",
      "20        productionStatusOfProcessedData = 2 [Research products  (grib2/tables/19/1.3.table) ]\r\n",
      "21        typeOfProcessedData = 1 [Forecast products  (grib2/tables/19/1.4.table , grib2/tables/local/edzw/1/1.4.table) ]\r\n",
      "======================   SECTION_3 ( length=35, padding=0 )    ======================\r\n",
      "1-4       section3Length = 35\r\n",
      "5         numberOfSection = 3\r\n",
      "6         sourceOfGridDefinition = 0 [Specified in Code table 3.1  (grib2/tables/19/3.0.table) ]\r\n",
      "7-10      numberOfDataPoints = 1\r\n",
      "11        numberOfOctectsForNumberOfPoints = 0\r\n",
      "12        interpretationOfNumberOfPoints = 0 [There is no appended list  (grib2/tables/19/3.11.table) ]\r\n",
      "13-14     gridDefinitionTemplateNumber = 101 [General unstructured grid  (grib2/tables/19/3.1.table) ]\r\n",
      "15        shapeOfTheEarth = 6 [Earth assumed spherical with radius of 6 371 229.0 m  (grib2/tables/19/3.2.table) ]\r\n",
      "16-18     numberOfGridUsed = 26\r\n",
      "19        numberOfGridInReference = 1\r\n",
      "20-35     uuidOfHGrid = 16 {\r\n",
      "                     a2, 7b, 8d, e6, 18, c4, 11, e4, 82, 0a, b5, b0, 98, c6, a5, c0\r\n",
      "                  } # bytes uuidOfHGrid \r\n",
      "======================   SECTION_4 ( length=101, padding=24 )   ======================\r\n",
      "1-4       section4Length = 101\r\n",
      "5         numberOfSection = 4\r\n",
      "6-7       NV = 0\r\n",
      "8-9       productDefinitionTemplateNumber = 67 [Average, accumulation and/or extreme values or other statistically processed values at a horizontal level or in a horizontal layer in a continuous or non-continuous time interval for atmospheric chemical constituents based on a distribution function (grib2/tables/19/4.0.table , grib2/tables/local/edzw/1/4.0.table) ]\r\n",
      "10        parameterCategory = 20 [Atmospheric chemical constituents  (grib2/tables/19/4.1.0.table , grib2/tables/local/edzw/1/4.1.0.table) ]\r\n",
      "11        parameterNumber = 6 [Surface dry deposition mass flux  (kg m-2 s-1)  (grib2/tables/19/4.2.0.20.table) ]\r\n",
      "12-13     constituentType = 62001 [Dust dry (grib2/tables/19/4.230.table) ]\r\n",
      "14-15     numberOfModeOfDistribution = 3\r\n",
      "16-17     modeNumber = 2\r\n",
      "18-19     typeOfDistributionFunction = 8 [No distribution function. The encoded variable is derived from variables characterized by type of distribution function of type no. 7 with fixed variance (p1)  (grib2/tables/19/4.240.table) ]\r\n",
      "20        numberOfDistributionFunctionParameters = 2\r\n",
      "21        scaleFactorOfDistributionFunctionParameter = 3\r\n",
      "22-25     scaledValueOfDistributionFunctionParameter = 1600\r\n",
      "26        scaleFactorOfDistributionFunctionParameter = 0\r\n",
      "27-30     scaledValueOfDistributionFunctionParameter = 2650\r\n",
      "31        typeOfGeneratingProcess = 2 [Forecast  (grib2/tables/19/4.3.table , grib2/tables/local/edzw/1/4.3.table) ]\r\n",
      "32        backgroundProcess = 0\r\n",
      "33        generatingProcessIdentifier = 1\r\n",
      "34-35     hoursAfterDataCutoff = 0\r\n",
      "36        minutesAfterDataCutoff = 0\r\n",
      "37        indicatorOfUnitOfTimeRange = 1 [Hour  (grib2/tables/19/4.4.table) ]\r\n",
      "38-41     forecastTime = 0\r\n",
      "42        typeOfFirstFixedSurface = 1 [Ground or water surface  (-)  (grib2/tables/19/4.5.table , grib2/tables/local/edzw/1/4.5.table) ]\r\n",
      "43        scaleFactorOfFirstFixedSurface = MISSING\r\n",
      "44-47     scaledValueOfFirstFixedSurface = MISSING\r\n",
      "48        typeOfSecondFixedSurface = 255 [Missing  (grib2/tables/19/4.5.table , grib2/tables/local/edzw/1/4.5.table) ]\r\n",
      "49        scaleFactorOfSecondFixedSurface = MISSING\r\n",
      "50-53     scaledValueOfSecondFixedSurface = MISSING\r\n",
      "54-55     yearOfEndOfOverallTimeInterval = 2014\r\n",
      "56        monthOfEndOfOverallTimeInterval = 4\r\n",
      "57        dayOfEndOfOverallTimeInterval = 5\r\n",
      "58        hourOfEndOfOverallTimeInterval = 0\r\n",
      "59        minuteOfEndOfOverallTimeInterval = 0\r\n",
      "60        secondOfEndOfOverallTimeInterval = 0\r\n",
      "61        numberOfTimeRange = 1\r\n",
      "62-65     numberOfMissingInStatisticalProcess = 0\r\n",
      "66        typeOfStatisticalProcessing = 1 [Accumulation  (grib2/tables/19/4.10.table) ]\r\n",
      "67        typeOfTimeIncrement = 2 [Successive times processed have same start time of forecast, forecast time is incremented  (grib2/tables/19/4.11.table , grib2/tables/local/edzw/1/4.11.table) ]\r\n",
      "68        indicatorOfUnitForTimeRange = 1 [Hour  (grib2/tables/19/4.4.table) ]\r\n",
      "69-72     lengthOfTimeRange = 48\r\n",
      "73        indicatorOfUnitForTimeIncrement = 255 [Missing  (grib2/tables/19/4.4.table) ]\r\n",
      "74-77     timeIncrement = 0\r\n",
      "======================   SECTION_5 ( length=21, padding=0 )    ======================\r\n",
      "1-4       section5Length = 21\r\n",
      "5         numberOfSection = 5\r\n",
      "6-9       numberOfValues = 1\r\n",
      "10-11     dataRepresentationTemplateNumber = 0 [Grid point data - simple packing  (grib2/tables/19/5.0.table) ]\r\n",
      "12-15     referenceValue = 999\r\n",
      "16-17     binaryScaleFactor = -76\r\n",
      "18-19     decimalScaleFactor = 0\r\n",
      "20        bitsPerValue = 0\r\n",
      "21        typeOfOriginalFieldValues = 0 [Floating point  (grib2/tables/19/5.1.table) ]\r\n",
      "======================   SECTION_6 ( length=6, padding=0 )     ======================\r\n",
      "1-4       section6Length = 6\r\n",
      "5         numberOfSection = 6\r\n",
      "6         bitMapIndicator = 255 [A bit map does not apply to this product  (grib2/tables/19/6.0.table) ]\r\n",
      "======================   SECTION_7 ( length=5, padding=0 )     ======================\r\n",
      "1-4       section7Length = 5\r\n",
      "5         numberOfSection = 7\r\n",
      "======================   SECTION_8 ( length=4, padding=0 )     ======================\r\n",
      "1-4       7777 = 7777\r\n",
      "#==============   MESSAGE 3 ( length=209 )                 ==============\r\n",
      "1-4       identifier = GRIB\r\n",
      "5-6       reserved = MISSING\r\n",
      "7         discipline = 0 [Meteorological products  (grib2/tables/19/0.0.table) ]\r\n",
      "8         editionNumber = 2\r\n",
      "9-16      totalLength = 209\r\n",
      "======================   SECTION_1 ( length=21, padding=0 )    ======================\r\n",
      "1-4       section1Length = 21\r\n",
      "5         numberOfSection = 1\r\n",
      "6-7       centre = 78 [Offenbach  (RSMC)  (grib1/0.table) ]\r\n",
      "8-9       subCentre = 255\r\n",
      "10        tablesVersion = 19 [Unknown code table entry (grib2/tables/1.0.table) ]\r\n",
      "11        localTablesVersion = 1 [Unknown code table entry (grib2/tables/19/1.1.table) ]\r\n",
      "12        significanceOfReferenceTime = 1 [Start of forecast  (grib2/tables/19/1.2.table) ]\r\n",
      "13-14     year = 2014\r\n",
      "15        month = 4\r\n",
      "16        day = 3\r\n",
      "17        hour = 0\r\n",
      "18        minute = 0\r\n",
      "19        second = 0\r\n",
      "20        productionStatusOfProcessedData = 2 [Research products  (grib2/tables/19/1.3.table) ]\r\n",
      "21        typeOfProcessedData = 1 [Forecast products  (grib2/tables/19/1.4.table , grib2/tables/local/edzw/1/1.4.table) ]\r\n",
      "======================   SECTION_3 ( length=35, padding=0 )    ======================\r\n",
      "1-4       section3Length = 35\r\n",
      "5         numberOfSection = 3\r\n",
      "6         sourceOfGridDefinition = 0 [Specified in Code table 3.1  (grib2/tables/19/3.0.table) ]\r\n",
      "7-10      numberOfDataPoints = 1\r\n",
      "11        numberOfOctectsForNumberOfPoints = 0\r\n",
      "12        interpretationOfNumberOfPoints = 0 [There is no appended list  (grib2/tables/19/3.11.table) ]\r\n",
      "13-14     gridDefinitionTemplateNumber = 101 [General unstructured grid  (grib2/tables/19/3.1.table) ]\r\n",
      "15        shapeOfTheEarth = 6 [Earth assumed spherical with radius of 6 371 229.0 m  (grib2/tables/19/3.2.table) ]\r\n",
      "16-18     numberOfGridUsed = 26\r\n",
      "19        numberOfGridInReference = 1\r\n",
      "20-35     uuidOfHGrid = 16 {\r\n",
      "                     a2, 7b, 8d, e6, 18, c4, 11, e4, 82, 0a, b5, b0, 98, c6, a5, c0\r\n",
      "                  } # bytes uuidOfHGrid \r\n",
      "======================   SECTION_4 ( length=101, padding=24 )   ======================\r\n",
      "1-4       section4Length = 101\r\n",
      "5         numberOfSection = 4\r\n",
      "6-7       NV = 0\r\n",
      "8-9       productDefinitionTemplateNumber = 67 [Average, accumulation and/or extreme values or other statistically processed values at a horizontal level or in a horizontal layer in a continuous or non-continuous time interval for atmospheric chemical constituents based on a distribution function (grib2/tables/19/4.0.table , grib2/tables/local/edzw/1/4.0.table) ]\r\n",
      "10        parameterCategory = 20 [Atmospheric chemical constituents  (grib2/tables/19/4.1.0.table , grib2/tables/local/edzw/1/4.1.0.table) ]\r\n",
      "11        parameterNumber = 6 [Surface dry deposition mass flux  (kg m-2 s-1)  (grib2/tables/19/4.2.0.20.table) ]\r\n",
      "12-13     constituentType = 62001 [Dust dry (grib2/tables/19/4.230.table) ]\r\n",
      "14-15     numberOfModeOfDistribution = 3\r\n",
      "16-17     modeNumber = 3\r\n",
      "18-19     typeOfDistributionFunction = 8 [No distribution function. The encoded variable is derived from variables characterized by type of distribution function of type no. 7 with fixed variance (p1)  (grib2/tables/19/4.240.table) ]\r\n",
      "20        numberOfDistributionFunctionParameters = 2\r\n",
      "21        scaleFactorOfDistributionFunctionParameter = 3\r\n",
      "22-25     scaledValueOfDistributionFunctionParameter = 1500\r\n",
      "26        scaleFactorOfDistributionFunctionParameter = 0\r\n",
      "27-30     scaledValueOfDistributionFunctionParameter = 2650\r\n",
      "31        typeOfGeneratingProcess = 2 [Forecast  (grib2/tables/19/4.3.table , grib2/tables/local/edzw/1/4.3.table) ]\r\n",
      "32        backgroundProcess = 0\r\n",
      "33        generatingProcessIdentifier = 1\r\n",
      "34-35     hoursAfterDataCutoff = 0\r\n",
      "36        minutesAfterDataCutoff = 0\r\n",
      "37        indicatorOfUnitOfTimeRange = 1 [Hour  (grib2/tables/19/4.4.table) ]\r\n",
      "38-41     forecastTime = 0\r\n",
      "42        typeOfFirstFixedSurface = 1 [Ground or water surface  (-)  (grib2/tables/19/4.5.table , grib2/tables/local/edzw/1/4.5.table) ]\r\n",
      "43        scaleFactorOfFirstFixedSurface = MISSING\r\n",
      "44-47     scaledValueOfFirstFixedSurface = MISSING\r\n",
      "48        typeOfSecondFixedSurface = 255 [Missing  (grib2/tables/19/4.5.table , grib2/tables/local/edzw/1/4.5.table) ]\r\n",
      "49        scaleFactorOfSecondFixedSurface = MISSING\r\n",
      "50-53     scaledValueOfSecondFixedSurface = MISSING\r\n",
      "54-55     yearOfEndOfOverallTimeInterval = 2014\r\n",
      "56        monthOfEndOfOverallTimeInterval = 4\r\n",
      "57        dayOfEndOfOverallTimeInterval = 5\r\n",
      "58        hourOfEndOfOverallTimeInterval = 0\r\n",
      "59        minuteOfEndOfOverallTimeInterval = 0\r\n",
      "60        secondOfEndOfOverallTimeInterval = 0\r\n",
      "61        numberOfTimeRange = 1\r\n",
      "62-65     numberOfMissingInStatisticalProcess = 0\r\n",
      "66        typeOfStatisticalProcessing = 1 [Accumulation  (grib2/tables/19/4.10.table) ]\r\n",
      "67        typeOfTimeIncrement = 2 [Successive times processed have same start time of forecast, forecast time is incremented  (grib2/tables/19/4.11.table , grib2/tables/local/edzw/1/4.11.table) ]\r\n",
      "68        indicatorOfUnitForTimeRange = 1 [Hour  (grib2/tables/19/4.4.table) ]\r\n",
      "69-72     lengthOfTimeRange = 48\r\n",
      "73        indicatorOfUnitForTimeIncrement = 255 [Missing  (grib2/tables/19/4.4.table) ]\r\n",
      "74-77     timeIncrement = 0\r\n",
      "======================   SECTION_5 ( length=21, padding=0 )    ======================\r\n",
      "1-4       section5Length = 21\r\n",
      "5         numberOfSection = 5\r\n",
      "6-9       numberOfValues = 1\r\n",
      "10-11     dataRepresentationTemplateNumber = 0 [Grid point data - simple packing  (grib2/tables/19/5.0.table) ]\r\n",
      "12-15     referenceValue = 999\r\n",
      "16-17     binaryScaleFactor = -10\r\n",
      "18-19     decimalScaleFactor = 0\r\n",
      "20        bitsPerValue = 0\r\n",
      "21        typeOfOriginalFieldValues = 0 [Floating point  (grib2/tables/19/5.1.table) ]\r\n",
      "======================   SECTION_6 ( length=6, padding=0 )     ======================\r\n",
      "1-4       section6Length = 6\r\n",
      "5         numberOfSection = 6\r\n",
      "6         bitMapIndicator = 255 [A bit map does not apply to this product  (grib2/tables/19/6.0.table) ]\r\n",
      "======================   SECTION_7 ( length=5, padding=0 )     ======================\r\n",
      "1-4       section7Length = 5\r\n",
      "5         numberOfSection = 7\r\n",
      "======================   SECTION_8 ( length=4, padding=0 )     ======================\r\n",
      "1-4       7777 = 7777\r\n",
      "***** FILE: pdt68.grib2 \r\n",
      "#==============   MESSAGE 4 ( length=212 )                 ==============\r\n",
      "1-4       identifier = GRIB\r\n",
      "5-6       reserved = MISSING\r\n",
      "7         discipline = 0 [Meteorological products  (grib2/tables/19/0.0.table) ]\r\n",
      "8         editionNumber = 2\r\n",
      "9-16      totalLength = 212\r\n",
      "======================   SECTION_1 ( length=21, padding=0 )    ======================\r\n",
      "1-4       section1Length = 21\r\n",
      "5         numberOfSection = 1\r\n",
      "6-7       centre = 78 [Offenbach  (RSMC)  (grib1/0.table) ]\r\n",
      "8-9       subCentre = 255\r\n",
      "10        tablesVersion = 19 [Unknown code table entry (grib2/tables/1.0.table) ]\r\n",
      "11        localTablesVersion = 1 [Unknown code table entry (grib2/tables/19/1.1.table) ]\r\n",
      "12        significanceOfReferenceTime = 1 [Start of forecast  (grib2/tables/19/1.2.table) ]\r\n",
      "13-14     year = 2014\r\n",
      "15        month = 4\r\n",
      "16        day = 3\r\n",
      "17        hour = 0\r\n",
      "18        minute = 0\r\n",
      "19        second = 0\r\n",
      "20        productionStatusOfProcessedData = 2 [Research products  (grib2/tables/19/1.3.table) ]\r\n",
      "21        typeOfProcessedData = 1 [Forecast products  (grib2/tables/19/1.4.table , grib2/tables/local/edzw/1/1.4.table) ]\r\n",
      "======================   SECTION_3 ( length=35, padding=0 )    ======================\r\n",
      "1-4       section3Length = 35\r\n",
      "5         numberOfSection = 3\r\n",
      "6         sourceOfGridDefinition = 0 [Specified in Code table 3.1  (grib2/tables/19/3.0.table) ]\r\n",
      "7-10      numberOfDataPoints = 1\r\n",
      "11        numberOfOctectsForNumberOfPoints = 0\r\n",
      "12        interpretationOfNumberOfPoints = 0 [There is no appended list  (grib2/tables/19/3.11.table) ]\r\n",
      "13-14     gridDefinitionTemplateNumber = 101 [General unstructured grid  (grib2/tables/19/3.1.table) ]\r\n",
      "15        shapeOfTheEarth = 6 [Earth assumed spherical with radius of 6 371 229.0 m  (grib2/tables/19/3.2.table) ]\r\n",
      "16-18     numberOfGridUsed = 26\r\n",
      "19        numberOfGridInReference = 1\r\n",
      "20-35     uuidOfHGrid = 16 {\r\n",
      "                     a2, 7b, 8d, e6, 18, c4, 11, e4, 82, 0a, b5, b0, 98, c6, a5, c0\r\n",
      "                  } # bytes uuidOfHGrid \r\n",
      "======================   SECTION_4 ( length=104, padding=24 )   ======================\r\n",
      "1-4       section4Length = 104\r\n",
      "5         numberOfSection = 4\r\n",
      "6-7       NV = 0\r\n",
      "8-9       productDefinitionTemplateNumber = 68 [Individual ensemble forecast, control and perturbed, at a horizontal level or oin a horizontal layer in a continuous or non-continuous time interval for atmospheric chemical constituents based on a distribution function (grib2/tables/19/4.0.table , grib2/tables/local/edzw/1/4.0.table) ]\r\n",
      "10        parameterCategory = 20 [Atmospheric chemical constituents  (grib2/tables/19/4.1.0.table , grib2/tables/local/edzw/1/4.1.0.table) ]\r\n",
      "11        parameterNumber = 6 [Surface dry deposition mass flux  (kg m-2 s-1)  (grib2/tables/19/4.2.0.20.table) ]\r\n",
      "12-13     constituentType = 62001 [Dust dry (grib2/tables/19/4.230.table) ]\r\n",
      "14-15     numberOfModeOfDistribution = 3\r\n",
      "16-17     modeNumber = 1\r\n",
      "18-19     typeOfDistributionFunction = 8 [No distribution function. The encoded variable is derived from variables characterized by type of distribution function of type no. 7 with fixed variance (p1)  (grib2/tables/19/4.240.table) ]\r\n",
      "20        numberOfDistributionFunctionParameters = 2\r\n",
      "21        scaleFactorOfDistributionFunctionParameter = 3\r\n",
      "22-25     scaledValueOfDistributionFunctionParameter = 1700\r\n",
      "26        scaleFactorOfDistributionFunctionParameter = 0\r\n",
      "27-30     scaledValueOfDistributionFunctionParameter = 2650\r\n",
      "31        typeOfGeneratingProcess = 2 [Forecast  (grib2/tables/19/4.3.table , grib2/tables/local/edzw/1/4.3.table) ]\r\n",
      "32        backgroundProcess = 0\r\n",
      "33        generatingProcessIdentifier = 1\r\n",
      "34-35     hoursAfterDataCutoff = 0\r\n",
      "36        minutesAfterDataCutoff = 0\r\n",
      "37        indicatorOfUnitOfTimeRange = 1 [Hour  (grib2/tables/19/4.4.table) ]\r\n",
      "38-41     forecastTime = 0\r\n",
      "42        typeOfFirstFixedSurface = 1 [Ground or water surface  (-)  (grib2/tables/19/4.5.table , grib2/tables/local/edzw/1/4.5.table) ]\r\n",
      "43        scaleFactorOfFirstFixedSurface = MISSING\r\n",
      "44-47     scaledValueOfFirstFixedSurface = MISSING\r\n",
      "48        typeOfSecondFixedSurface = 255 [Missing  (grib2/tables/19/4.5.table , grib2/tables/local/edzw/1/4.5.table) ]\r\n",
      "49        scaleFactorOfSecondFixedSurface = MISSING\r\n",
      "50-53     scaledValueOfSecondFixedSurface = MISSING\r\n",
      "54        typeOfEnsembleForecast = 192 [other types of ensemble forecasts (grib2/tables/19/4.6.table , grib2/tables/local/edzw/1/4.6.table) ]\r\n",
      "55        perturbationNumber = 13\r\n",
      "56        numberOfForecastsInEnsemble = 20\r\n",
      "57-58     yearOfEndOfOverallTimeInterval = 2014\r\n",
      "59        monthOfEndOfOverallTimeInterval = 4\r\n",
      "60        dayOfEndOfOverallTimeInterval = 5\r\n",
      "61        hourOfEndOfOverallTimeInterval = 0\r\n",
      "62        minuteOfEndOfOverallTimeInterval = 0\r\n",
      "63        secondOfEndOfOverallTimeInterval = 0\r\n",
      "64        numberOfTimeRange = 1\r\n",
      "65-68     numberOfMissingInStatisticalProcess = 0\r\n",
      "69        typeOfStatisticalProcessing = 1 [Accumulation  (grib2/tables/19/4.10.table) ]\r\n",
      "70        typeOfTimeIncrement = 2 [Successive times processed have same start time of forecast, forecast time is incremented  (grib2/tables/19/4.11.table , grib2/tables/local/edzw/1/4.11.table) ]\r\n",
      "71        indicatorOfUnitForTimeRange = 1 [Hour  (grib2/tables/19/4.4.table) ]\r\n",
      "72-75     lengthOfTimeRange = 48\r\n",
      "76        indicatorOfUnitForTimeIncrement = 255 [Missing  (grib2/tables/19/4.4.table) ]\r\n",
      "77-80     timeIncrement = 0\r\n",
      "======================   SECTION_5 ( length=21, padding=0 )    ======================\r\n",
      "1-4       section5Length = 21\r\n",
      "5         numberOfSection = 5\r\n",
      "6-9       numberOfValues = 1\r\n",
      "10-11     dataRepresentationTemplateNumber = 0 [Grid point data - simple packing  (grib2/tables/19/5.0.table) ]\r\n",
      "12-15     referenceValue = 999\r\n",
      "16-17     binaryScaleFactor = -73\r\n",
      "18-19     decimalScaleFactor = 0\r\n",
      "20        bitsPerValue = 0\r\n",
      "21        typeOfOriginalFieldValues = 0 [Floating point  (grib2/tables/19/5.1.table) ]\r\n",
      "======================   SECTION_6 ( length=6, padding=0 )     ======================\r\n",
      "1-4       section6Length = 6\r\n",
      "5         numberOfSection = 6\r\n",
      "6         bitMapIndicator = 255 [A bit map does not apply to this product  (grib2/tables/19/6.0.table) ]\r\n",
      "======================   SECTION_7 ( length=5, padding=0 )     ======================\r\n",
      "1-4       section7Length = 5\r\n",
      "5         numberOfSection = 7\r\n",
      "======================   SECTION_8 ( length=4, padding=0 )     ======================\r\n",
      "1-4       7777 = 7777\r\n",
      "#==============   MESSAGE 5 ( length=212 )                 ==============\r\n",
      "1-4       identifier = GRIB\r\n",
      "5-6       reserved = MISSING\r\n",
      "7         discipline = 0 [Meteorological products  (grib2/tables/19/0.0.table) ]\r\n",
      "8         editionNumber = 2\r\n",
      "9-16      totalLength = 212\r\n",
      "======================   SECTION_1 ( length=21, padding=0 )    ======================\r\n",
      "1-4       section1Length = 21\r\n",
      "5         numberOfSection = 1\r\n",
      "6-7       centre = 78 [Offenbach  (RSMC)  (grib1/0.table) ]\r\n",
      "8-9       subCentre = 255\r\n",
      "10        tablesVersion = 19 [Unknown code table entry (grib2/tables/1.0.table) ]\r\n",
      "11        localTablesVersion = 1 [Unknown code table entry (grib2/tables/19/1.1.table) ]\r\n",
      "12        significanceOfReferenceTime = 1 [Start of forecast  (grib2/tables/19/1.2.table) ]\r\n",
      "13-14     year = 2014\r\n",
      "15        month = 4\r\n",
      "16        day = 3\r\n",
      "17        hour = 0\r\n",
      "18        minute = 0\r\n",
      "19        second = 0\r\n",
      "20        productionStatusOfProcessedData = 2 [Research products  (grib2/tables/19/1.3.table) ]\r\n",
      "21        typeOfProcessedData = 1 [Forecast products  (grib2/tables/19/1.4.table , grib2/tables/local/edzw/1/1.4.table) ]\r\n",
      "======================   SECTION_3 ( length=35, padding=0 )    ======================\r\n",
      "1-4       section3Length = 35\r\n",
      "5         numberOfSection = 3\r\n",
      "6         sourceOfGridDefinition = 0 [Specified in Code table 3.1  (grib2/tables/19/3.0.table) ]\r\n",
      "7-10      numberOfDataPoints = 1\r\n",
      "11        numberOfOctectsForNumberOfPoints = 0\r\n",
      "12        interpretationOfNumberOfPoints = 0 [There is no appended list  (grib2/tables/19/3.11.table) ]\r\n",
      "13-14     gridDefinitionTemplateNumber = 101 [General unstructured grid  (grib2/tables/19/3.1.table) ]\r\n",
      "15        shapeOfTheEarth = 6 [Earth assumed spherical with radius of 6 371 229.0 m  (grib2/tables/19/3.2.table) ]\r\n",
      "16-18     numberOfGridUsed = 26\r\n",
      "19        numberOfGridInReference = 1\r\n",
      "20-35     uuidOfHGrid = 16 {\r\n",
      "                     a2, 7b, 8d, e6, 18, c4, 11, e4, 82, 0a, b5, b0, 98, c6, a5, c0\r\n",
      "                  } # bytes uuidOfHGrid \r\n",
      "======================   SECTION_4 ( length=104, padding=24 )   ======================\r\n",
      "1-4       section4Length = 104\r\n",
      "5         numberOfSection = 4\r\n",
      "6-7       NV = 0\r\n",
      "8-9       productDefinitionTemplateNumber = 68 [Individual ensemble forecast, control and perturbed, at a horizontal level or oin a horizontal layer in a continuous or non-continuous time interval for atmospheric chemical constituents based on a distribution function (grib2/tables/19/4.0.table , grib2/tables/local/edzw/1/4.0.table) ]\r\n",
      "10        parameterCategory = 20 [Atmospheric chemical constituents  (grib2/tables/19/4.1.0.table , grib2/tables/local/edzw/1/4.1.0.table) ]\r\n",
      "11        parameterNumber = 6 [Surface dry deposition mass flux  (kg m-2 s-1)  (grib2/tables/19/4.2.0.20.table) ]\r\n",
      "12-13     constituentType = 62001 [Dust dry (grib2/tables/19/4.230.table) ]\r\n",
      "14-15     numberOfModeOfDistribution = 3\r\n",
      "16-17     modeNumber = 2\r\n",
      "18-19     typeOfDistributionFunction = 8 [No distribution function. The encoded variable is derived from variables characterized by type of distribution function of type no. 7 with fixed variance (p1)  (grib2/tables/19/4.240.table) ]\r\n",
      "20        numberOfDistributionFunctionParameters = 2\r\n",
      "21        scaleFactorOfDistributionFunctionParameter = 3\r\n",
      "22-25     scaledValueOfDistributionFunctionParameter = 1600\r\n",
      "26        scaleFactorOfDistributionFunctionParameter = 0\r\n",
      "27-30     scaledValueOfDistributionFunctionParameter = 2650\r\n",
      "31        typeOfGeneratingProcess = 2 [Forecast  (grib2/tables/19/4.3.table , grib2/tables/local/edzw/1/4.3.table) ]\r\n",
      "32        backgroundProcess = 0\r\n",
      "33        generatingProcessIdentifier = 1\r\n",
      "34-35     hoursAfterDataCutoff = 0\r\n",
      "36        minutesAfterDataCutoff = 0\r\n",
      "37        indicatorOfUnitOfTimeRange = 1 [Hour  (grib2/tables/19/4.4.table) ]\r\n",
      "38-41     forecastTime = 0\r\n",
      "42        typeOfFirstFixedSurface = 1 [Ground or water surface  (-)  (grib2/tables/19/4.5.table , grib2/tables/local/edzw/1/4.5.table) ]\r\n",
      "43        scaleFactorOfFirstFixedSurface = MISSING\r\n",
      "44-47     scaledValueOfFirstFixedSurface = MISSING\r\n",
      "48        typeOfSecondFixedSurface = 255 [Missing  (grib2/tables/19/4.5.table , grib2/tables/local/edzw/1/4.5.table) ]\r\n",
      "49        scaleFactorOfSecondFixedSurface = MISSING\r\n",
      "50-53     scaledValueOfSecondFixedSurface = MISSING\r\n",
      "54        typeOfEnsembleForecast = 192 [other types of ensemble forecasts (grib2/tables/19/4.6.table , grib2/tables/local/edzw/1/4.6.table) ]\r\n",
      "55        perturbationNumber = 13\r\n",
      "56        numberOfForecastsInEnsemble = 20\r\n",
      "57-58     yearOfEndOfOverallTimeInterval = 2014\r\n",
      "59        monthOfEndOfOverallTimeInterval = 4\r\n",
      "60        dayOfEndOfOverallTimeInterval = 5\r\n",
      "61        hourOfEndOfOverallTimeInterval = 0\r\n",
      "62        minuteOfEndOfOverallTimeInterval = 0\r\n",
      "63        secondOfEndOfOverallTimeInterval = 0\r\n",
      "64        numberOfTimeRange = 1\r\n",
      "65-68     numberOfMissingInStatisticalProcess = 0\r\n",
      "69        typeOfStatisticalProcessing = 1 [Accumulation  (grib2/tables/19/4.10.table) ]\r\n",
      "70        typeOfTimeIncrement = 2 [Successive times processed have same start time of forecast, forecast time is incremented  (grib2/tables/19/4.11.table , grib2/tables/local/edzw/1/4.11.table) ]\r\n",
      "71        indicatorOfUnitForTimeRange = 1 [Hour  (grib2/tables/19/4.4.table) ]\r\n",
      "72-75     lengthOfTimeRange = 48\r\n",
      "76        indicatorOfUnitForTimeIncrement = 255 [Missing  (grib2/tables/19/4.4.table) ]\r\n",
      "77-80     timeIncrement = 0\r\n",
      "======================   SECTION_5 ( length=21, padding=0 )    ======================\r\n",
      "1-4       section5Length = 21\r\n",
      "5         numberOfSection = 5\r\n",
      "6-9       numberOfValues = 1\r\n",
      "10-11     dataRepresentationTemplateNumber = 0 [Grid point data - simple packing  (grib2/tables/19/5.0.table) ]\r\n",
      "12-15     referenceValue = 999\r\n",
      "16-17     binaryScaleFactor = -76\r\n",
      "18-19     decimalScaleFactor = 0\r\n",
      "20        bitsPerValue = 0\r\n",
      "21        typeOfOriginalFieldValues = 0 [Floating point  (grib2/tables/19/5.1.table) ]\r\n",
      "======================   SECTION_6 ( length=6, padding=0 )     ======================\r\n",
      "1-4       section6Length = 6\r\n",
      "5         numberOfSection = 6\r\n",
      "6         bitMapIndicator = 255 [A bit map does not apply to this product  (grib2/tables/19/6.0.table) ]\r\n",
      "======================   SECTION_7 ( length=5, padding=0 )     ======================\r\n",
      "1-4       section7Length = 5\r\n",
      "5         numberOfSection = 7\r\n",
      "======================   SECTION_8 ( length=4, padding=0 )     ======================\r\n",
      "1-4       7777 = 7777\r\n",
      "#==============   MESSAGE 6 ( length=212 )                 ==============\r\n",
      "1-4       identifier = GRIB\r\n",
      "5-6       reserved = MISSING\r\n",
      "7         discipline = 0 [Meteorological products  (grib2/tables/19/0.0.table) ]\r\n",
      "8         editionNumber = 2\r\n",
      "9-16      totalLength = 212\r\n",
      "======================   SECTION_1 ( length=21, padding=0 )    ======================\r\n",
      "1-4       section1Length = 21\r\n",
      "5         numberOfSection = 1\r\n",
      "6-7       centre = 78 [Offenbach  (RSMC)  (grib1/0.table) ]\r\n",
      "8-9       subCentre = 255\r\n",
      "10        tablesVersion = 19 [Unknown code table entry (grib2/tables/1.0.table) ]\r\n",
      "11        localTablesVersion = 1 [Unknown code table entry (grib2/tables/19/1.1.table) ]\r\n",
      "12        significanceOfReferenceTime = 1 [Start of forecast  (grib2/tables/19/1.2.table) ]\r\n",
      "13-14     year = 2014\r\n",
      "15        month = 4\r\n",
      "16        day = 3\r\n",
      "17        hour = 0\r\n",
      "18        minute = 0\r\n",
      "19        second = 0\r\n",
      "20        productionStatusOfProcessedData = 2 [Research products  (grib2/tables/19/1.3.table) ]\r\n",
      "21        typeOfProcessedData = 1 [Forecast products  (grib2/tables/19/1.4.table , grib2/tables/local/edzw/1/1.4.table) ]\r\n",
      "======================   SECTION_3 ( length=35, padding=0 )    ======================\r\n",
      "1-4       section3Length = 35\r\n",
      "5         numberOfSection = 3\r\n",
      "6         sourceOfGridDefinition = 0 [Specified in Code table 3.1  (grib2/tables/19/3.0.table) ]\r\n",
      "7-10      numberOfDataPoints = 1\r\n",
      "11        numberOfOctectsForNumberOfPoints = 0\r\n",
      "12        interpretationOfNumberOfPoints = 0 [There is no appended list  (grib2/tables/19/3.11.table) ]\r\n",
      "13-14     gridDefinitionTemplateNumber = 101 [General unstructured grid  (grib2/tables/19/3.1.table) ]\r\n",
      "15        shapeOfTheEarth = 6 [Earth assumed spherical with radius of 6 371 229.0 m  (grib2/tables/19/3.2.table) ]\r\n",
      "16-18     numberOfGridUsed = 26\r\n",
      "19        numberOfGridInReference = 1\r\n",
      "20-35     uuidOfHGrid = 16 {\r\n",
      "                     a2, 7b, 8d, e6, 18, c4, 11, e4, 82, 0a, b5, b0, 98, c6, a5, c0\r\n",
      "                  } # bytes uuidOfHGrid \r\n",
      "======================   SECTION_4 ( length=104, padding=24 )   ======================\r\n",
      "1-4       section4Length = 104\r\n",
      "5         numberOfSection = 4\r\n",
      "6-7       NV = 0\r\n",
      "8-9       productDefinitionTemplateNumber = 68 [Individual ensemble forecast, control and perturbed, at a horizontal level or oin a horizontal layer in a continuous or non-continuous time interval for atmospheric chemical constituents based on a distribution function (grib2/tables/19/4.0.table , grib2/tables/local/edzw/1/4.0.table) ]\r\n",
      "10        parameterCategory = 20 [Atmospheric chemical constituents  (grib2/tables/19/4.1.0.table , grib2/tables/local/edzw/1/4.1.0.table) ]\r\n",
      "11        parameterNumber = 6 [Surface dry deposition mass flux  (kg m-2 s-1)  (grib2/tables/19/4.2.0.20.table) ]\r\n",
      "12-13     constituentType = 62001 [Dust dry (grib2/tables/19/4.230.table) ]\r\n",
      "14-15     numberOfModeOfDistribution = 3\r\n",
      "16-17     modeNumber = 3\r\n",
      "18-19     typeOfDistributionFunction = 8 [No distribution function. The encoded variable is derived from variables characterized by type of distribution function of type no. 7 with fixed variance (p1)  (grib2/tables/19/4.240.table) ]\r\n",
      "20        numberOfDistributionFunctionParameters = 2\r\n",
      "21        scaleFactorOfDistributionFunctionParameter = 3\r\n",
      "22-25     scaledValueOfDistributionFunctionParameter = 1500\r\n",
      "26        scaleFactorOfDistributionFunctionParameter = 0\r\n",
      "27-30     scaledValueOfDistributionFunctionParameter = 2650\r\n",
      "31        typeOfGeneratingProcess = 2 [Forecast  (grib2/tables/19/4.3.table , grib2/tables/local/edzw/1/4.3.table) ]\r\n",
      "32        backgroundProcess = 0\r\n",
      "33        generatingProcessIdentifier = 1\r\n",
      "34-35     hoursAfterDataCutoff = 0\r\n",
      "36        minutesAfterDataCutoff = 0\r\n",
      "37        indicatorOfUnitOfTimeRange = 1 [Hour  (grib2/tables/19/4.4.table) ]\r\n",
      "38-41     forecastTime = 0\r\n",
      "42        typeOfFirstFixedSurface = 1 [Ground or water surface  (-)  (grib2/tables/19/4.5.table , grib2/tables/local/edzw/1/4.5.table) ]\r\n",
      "43        scaleFactorOfFirstFixedSurface = MISSING\r\n",
      "44-47     scaledValueOfFirstFixedSurface = MISSING\r\n",
      "48        typeOfSecondFixedSurface = 255 [Missing  (grib2/tables/19/4.5.table , grib2/tables/local/edzw/1/4.5.table) ]\r\n",
      "49        scaleFactorOfSecondFixedSurface = MISSING\r\n",
      "50-53     scaledValueOfSecondFixedSurface = MISSING\r\n",
      "54        typeOfEnsembleForecast = 192 [other types of ensemble forecasts (grib2/tables/19/4.6.table , grib2/tables/local/edzw/1/4.6.table) ]\r\n",
      "55        perturbationNumber = 13\r\n",
      "56        numberOfForecastsInEnsemble = 20\r\n",
      "57-58     yearOfEndOfOverallTimeInterval = 2014\r\n",
      "59        monthOfEndOfOverallTimeInterval = 4\r\n",
      "60        dayOfEndOfOverallTimeInterval = 5\r\n",
      "61        hourOfEndOfOverallTimeInterval = 0\r\n",
      "62        minuteOfEndOfOverallTimeInterval = 0\r\n",
      "63        secondOfEndOfOverallTimeInterval = 0\r\n",
      "64        numberOfTimeRange = 1\r\n",
      "65-68     numberOfMissingInStatisticalProcess = 0\r\n",
      "69        typeOfStatisticalProcessing = 1 [Accumulation  (grib2/tables/19/4.10.table) ]\r\n",
      "70        typeOfTimeIncrement = 2 [Successive times processed have same start time of forecast, forecast time is incremented  (grib2/tables/19/4.11.table , grib2/tables/local/edzw/1/4.11.table) ]\r\n",
      "71        indicatorOfUnitForTimeRange = 1 [Hour  (grib2/tables/19/4.4.table) ]\r\n",
      "72-75     lengthOfTimeRange = 48\r\n",
      "76        indicatorOfUnitForTimeIncrement = 255 [Missing  (grib2/tables/19/4.4.table) ]\r\n",
      "77-80     timeIncrement = 0\r\n",
      "======================   SECTION_5 ( length=21, padding=0 )    ======================\r\n",
      "1-4       section5Length = 21\r\n",
      "5         numberOfSection = 5\r\n",
      "6-9       numberOfValues = 1\r\n",
      "10-11     dataRepresentationTemplateNumber = 0 [Grid point data - simple packing  (grib2/tables/19/5.0.table) ]\r\n",
      "12-15     referenceValue = 999\r\n",
      "16-17     binaryScaleFactor = -10\r\n",
      "18-19     decimalScaleFactor = 0\r\n",
      "20        bitsPerValue = 0\r\n",
      "21        typeOfOriginalFieldValues = 0 [Floating point  (grib2/tables/19/5.1.table) ]\r\n",
      "======================   SECTION_6 ( length=6, padding=0 )     ======================\r\n",
      "1-4       section6Length = 6\r\n",
      "5         numberOfSection = 6\r\n",
      "6         bitMapIndicator = 255 [A bit map does not apply to this product  (grib2/tables/19/6.0.table) ]\r\n",
      "======================   SECTION_7 ( length=5, padding=0 )     ======================\r\n",
      "1-4       section7Length = 5\r\n",
      "5         numberOfSection = 7\r\n",
      "======================   SECTION_8 ( length=4, padding=0 )     ======================\r\n",
      "1-4       7777 = 7777\r\n"
     ]
    }
   ],
   "source": [
    "grib_dump $validation_files -O"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "wgrib2 can also show the GRIBs' contents."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "wgrib2 output for pdt67.grib2:\r\n",
      "1:0:00Z03apr2014:SDDMFLX Surface Dry Deposition Mass Flux [kg/m^2/s]:lvl1=(1,missing) lvl2=(255,missing):surface:anl:chemical=Dust Dry\r\n",
      "2:209:00Z03apr2014:SDDMFLX Surface Dry Deposition Mass Flux [kg/m^2/s]:lvl1=(1,missing) lvl2=(255,missing):surface:anl:chemical=Dust Dry\r\n",
      "3:418:00Z03apr2014:SDDMFLX Surface Dry Deposition Mass Flux [kg/m^2/s]:lvl1=(1,missing) lvl2=(255,missing):surface:anl:chemical=Dust Dry\r\n",
      "wgrib2 output for pdt68.grib2:\r\n",
      "1:0:00Z03apr2014:SDDMFLX Surface Dry Deposition Mass Flux [kg/m^2/s]:lvl1=(1,missing) lvl2=(255,missing):surface:anl:ENS=? table4.6=192 pert=13:chemical=Dust Dry\r\n",
      "2:212:00Z03apr2014:SDDMFLX Surface Dry Deposition Mass Flux [kg/m^2/s]:lvl1=(1,missing) lvl2=(255,missing):surface:anl:ENS=? table4.6=192 pert=13:chemical=Dust Dry\r\n",
      "3:424:00Z03apr2014:SDDMFLX Surface Dry Deposition Mass Flux [kg/m^2/s]:lvl1=(1,missing) lvl2=(255,missing):surface:anl:ENS=? table4.6=192 pert=13:chemical=Dust Dry\r\n"
     ]
    }
   ],
   "source": [
    "for grib in $validation_files\n",
    "do\n",
    "    echo wgrib2 output for $grib:\n",
    "    # Using wgrib2 installed from sources linked to above\n",
    "    wgrib2 -v2 $grib\n",
    "done"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Bash",
   "language": "bash",
   "name": "bash"
  },
  "language_info": {
   "codemirror_mode": "shell",
   "file_extension": ".sh",
   "mimetype": "text/x-sh",
   "name": "bash"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}

VREPORT_StatisticalProcess4Distribution.docx