Exemple #1
0
        public void Section_Check()
        {
            using var stream = File.OpenRead(GribFileSamples.WmoOneDataSetMessage);
            var reader = new Grib2Reader(stream);

            var messages = reader.ReadMessages();

            Check.That(messages).HasOneElementOnly();
            Check.That(messages.Single().DataSets).HasOneElementOnly();

            var record = messages.Single().DataSets.Single();

            CheckIndicatorSection(record);

            CheckIdentificationSection(record);

            CheckGridDefinitionSection(record);

            CheckProductDefinitionSection(record);

            CheckDataRepresentationSection(record);

            CheckBitmapSection(record);

            CheckDataSection(record);
        }
Exemple #2
0
        public void Read_PerturbationFile_Test()
        {
            using var stream = File.OpenRead(GribFileSamples.NcepGefsPerturbationFile);
            var reader = new Grib2Reader(stream);

            var datasets = reader.ReadAllDataSets().ToArray();

            const int expectedEnsembleForecastsNumber = 20;
            const int expectedPerturbationNumber      = 8;

            int perturbationNumber;
            int ensembleForecastsNumber;

            var tmpDataSet = datasets.Single(datasets => datasets.Parameter.Equals(Parameter.Temperature));

            Check.That(tmpDataSet.ProductDefinitionSection.ProductDefinition.TryGet(ProductDefinitionContent.EnsembleForecastsNumber, out ensembleForecastsNumber)).IsTrue();
            Check.That(ensembleForecastsNumber).IsEqualTo(expectedEnsembleForecastsNumber);

            Check.That(tmpDataSet.ProductDefinitionSection.TryGet(ProductDefinitionContent.PerturbationNumber, out perturbationNumber)).IsTrue();
            Check.That(perturbationNumber).IsEqualTo(expectedPerturbationNumber);

            var apcpDataset = datasets.Single(d => d.Parameter.Equals(Parameter.TotalPrecipitation));

            Check.That(apcpDataset.ProductDefinitionSection.ProductDefinition.TryGet(ProductDefinitionContent.EnsembleForecastsNumber, out ensembleForecastsNumber)).IsTrue();
            Check.That(ensembleForecastsNumber).IsEqualTo(expectedEnsembleForecastsNumber);

            Check.That(apcpDataset.ProductDefinitionSection.TryGet(ProductDefinitionContent.PerturbationNumber, out perturbationNumber)).IsTrue();
            Check.That(perturbationNumber).IsEqualTo(expectedPerturbationNumber);
        }
Exemple #3
0
        public void Read_AvgFile_Test()
        {
            using var stream = File.OpenRead(GribFileSamples.NcepGefsAvgFile);
            var reader = new Grib2Reader(stream);

            var datasets = reader.ReadAllDataSets().ToArray();

            var apcpDataset = datasets.Single(d => d.Parameter.Equals(Parameter.TotalPrecipitation));

            Check.That(apcpDataset.ProductDefinitionSection.ProductDefinition.TryGet(ProductDefinitionContent.EnsembleForecastsNumber, out var ensembleForecastsNumber)).IsTrue();
            Check.That(ensembleForecastsNumber).IsEqualTo(20);
        }
        public void Test()
        {
            var reader = new Grib2Reader(GribFileSamples.MfAromeTemperatureFile);

            var dsTemperature = reader.ReadAllDataSets().Single();

            var items = new Dictionary <Coordinate, float?>(reader.ReadDataSetValues(dsTemperature));

            double resolution = 0.025, north = 53, west = -8, east = 12;
            double lat = north;
            double lon = west;

            foreach (var item in expectedAromeValues)
            {
                var shouldChangeLatitude = lon >= east;

                Check.That(items[(lat, lon)]).IsEqualTo(item);
Exemple #5
0
        public void Test()
        {
            using var stream = File.OpenRead(GribFileSamples.NcepGfsTmpApcpFile);
            var reader = new Grib2Reader(stream);

            var datasets = reader.ReadAllDataSets().ToArray();

            Check.That(datasets.Select(d => d.Parameter)).ContainsExactly(
                Parameter.Temperature,
                Parameter.TotalPrecipitation,                // APCP 6 hours accumulation
                Parameter.TotalPrecipitation                 // APCP 96 hours accumulation
                );

            var temperatureDs = datasets.Single(d => d.Parameter.Equals(Parameter.Temperature));

            var data = reader.ReadDataSetValues(temperatureDs).ToList();

            // Excepted values read using Panoply
            Check.That(data).ContainsExactly(new Dictionary <Coordinate, float?>
            {
                { (-21.25, 55), 301.405579f },
        public void Test()
        {
            var reader = new Grib2Reader(GribFileSamples.DwdCosmoTotalPrecipitationFile);

            var dsTemperature = reader.ReadAllDataSets().Single();

            var items = new Dictionary <Coordinate, float?>(reader.ReadDataSetValues(dsTemperature));

            double resolution = 0.02, south = 43, west = 0, east = 17;
            double lat = south;
            double lon = west;

            // Expected values from Panoply
            // ordered from West to East then South to North
            var expectedCosmoValues = TestUtils.ReadCsvValues("references/DWD-COSMO-Total_precipitation_rate_surface_12_Hour_Accumulation.csv");

            foreach (var item in expectedCosmoValues)
            {
                var shouldChangeLatitude = lon >= east - 1e-8;

                Check.That(items[(lat, lon)]).IsEqualTo(item);
Exemple #7
0
        public void Data_Check()
        {
            // Expected values retrieved from Panoply
            var expectedValues = new[]
            {
                5407.5f, 5416.5f, 5360.6f, 5470.8f, 5512.7f,
                5449.3f, 5376.0f, 5503.3f, 5515.2f, 5448.4f,
                5405.9f, 5467.4f, 5440.7f, 5448.5f, 5403.4f,
                5348.1f, 5378.2f, 5417.6f, 5388.2f, 5520.9f,
                5367.6f, 5465.7f, 5357.2f, 5475.4f, 5407.3f
            };

            using var stream = File.OpenRead(GribFileSamples.WmoOneDataSetMessage);
            var reader = new Grib2Reader(stream);

            var dataSets = reader.ReadAllDataSets().ToArray();

            Check.That(dataSets).HasOneElementOnly();

            var data = reader.ReadDataSetRawData(dataSets.Single()).ToArray();

            Check.That(data).ContainsExactly(expectedValues);
        }