public void Should_be_able_to_decode_file_when_all_it_is_correct_v2(string correctFileResource)
{
var input = _fileIndex [correctFileResource];
var stream = new MemoryStream (input);
_testedv2 = new Grib2Input (stream);
_testedv2.scan (true, true);
Assert.IsTrue (_testedv2.ProductsCount > 0);
}
static void GribTwo()
{
#region Grib 2 Code
Grib2Input input = new Grib2Input(RandomAccessFile);
if (!input.scan(false, false))
{
Console.WriteLine("Failed to successfully scan grib file");
return;
}
Grib2Data data = new Grib2Data(RandomAccessFile);
var records = input.Records;
foreach (Grib2Record record in records)
{
IGrib2IndicatorSection iis = record.Is;
IGrib2IdentificationSection id = record.ID;
IGrib2ProductDefinitionSection pdsv = record.PDS;
IGrib2GridDefinitionSection gdsv = record.GDS;
long time = id.RefTime.AddTicks(record.PDS.ForecastTime * 3600000).Ticks;
Console.WriteLine("Record description at " + " forecast " + new DateTime(time) + ": " + string.Format("{0} {1} {2}", iis.Discipline, pdsv.ParameterCategory, pdsv.ParameterNumber));
float[] values = data.getData(record.getGdsOffset(), record.getPdsOffset());
if ((iis.Discipline == 0) && (pdsv.ParameterCategory == 2) && (pdsv.ParameterNumber == 2))
{
// U-component_of_wind
int c = 0;
for (double lat = gdsv.La1; lat >= gdsv.La2; lat = lat - gdsv.Dy)
{
for (double lon = gdsv.Lo1; lon <= gdsv.Lo2; lon = lon + gdsv.Dx)
{
Console.WriteLine("U-Wind " + lat + "\t" + lon + "\t" + values[c]);
c++;
}
}
}
else if ((iis.Discipline == 0) && (pdsv.ParameterCategory == 2) && (pdsv.ParameterNumber == 3))
{
// V-component_of_wind
int c = 0;
for (double lat = gdsv.La1; lat >= gdsv.La2; lat = lat - gdsv.Dy)
{
for (double lon = gdsv.Lo1; lon <= gdsv.Lo2; lon = lon + gdsv.Dx)
{
Console.WriteLine("V-Wind " + lat + "\t" + lon + "\t" + values[c]);
c++;
}
}
}
else if ((iis.Discipline == 0) && (pdsv.ParameterCategory == 1) && (pdsv.ParameterNumber == 1))
{
// RH
int c = 0;
for (double lat = gdsv.La1; lat >= gdsv.La2; lat = lat - gdsv.Dy)
{
for (double lon = gdsv.Lo1; lon <= gdsv.Lo2; lon = lon + gdsv.Dx)
{
Console.WriteLine("RH " + lat + "\t" + lon + "\t" + values[c]);
c++;
}
}
}
}
#endregion
}
