internal Grib1WaveSpectra2DDirFreq(System.IO.Stream raf) { SupportClass.Skip(raf, 4); byte[] ver = GribNumbers.ReadBytes(raf, 4); System.Text.Encoding enc = System.Text.Encoding.ASCII; version = enc.GetString(ver); SupportClass.Skip(raf, 2); directionNumber = raf.ReadByte(); frequencyNumber = raf.ReadByte(); numberOfDirections = raf.ReadByte(); numberOfFrequencies = raf.ReadByte(); dirScaleFactor = GribNumbers.int4(raf); freqScaleFactor = GribNumbers.int4(raf); SupportClass.Skip(raf, 37); directions = new double[numberOfDirections]; frequencies = new double[numberOfFrequencies]; for (int i = 0; i < numberOfDirections; i++) { int val = GribNumbers.int4(raf); directions[i] = (double)val / (double)dirScaleFactor; } for (int i = 0; i < numberOfFrequencies; i++) { int val = GribNumbers.int4(raf); frequencies[i] = (double)val / (double)freqScaleFactor; } }
// *** constructors ******************************************************* /// <summary> Constructs a <tt>Grib2LocalUseSection</tt> object from a raf.</summary> /// <param name="raf"> /// </param> /// <throws> IOException if raf contains no valid GRIB product </throws> //UPGRADE_TODO: Class 'java.io.RandomAccessFile' was converted to 'System.IO.FileStream' which has a different behavior. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1073_javaioRandomAccessFile'" public Grib2LocalUseSection(System.IO.FileStream raf) { // octets 1-4 (Length of GDS) length = GribNumbers.int4(raf); //System.out.println( "LUS length=" + length ); section = raf.ReadByte(); // This is section 2 //System.out.println( "LUS Local Use is 2, section=" + section ); if (section != 2) { // no local use section SupportClass.Skip(raf, -5); return; } else { //SupportClass.Skip(raf, length - 5); bytes = new byte[length - 5]; int nb = raf.Read(bytes, 0, length - 5); if (nb != length - 5) { throw new NoValidGribException("Failed to read Local Use Section data"); } } } // end of Grib2LocalUseSection
private int Xlength; // length of the x axis // *** constructors ******************************************************* /// <summary> Constructor for a Grib2 Data Section.</summary> /// <param name="getData"> /// </param> /// <param name="raf"> /// </param> /// <param name="gds"> /// </param> /// <param name="drs"> /// </param> /// <param name="bms"> /// </param> /// <throws> IOException </throws> //UPGRADE_TODO: Class 'java.io.RandomAccessFile' was converted to 'System.IO.FileStream' which has a different behavior. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1073_javaioRandomAccessFile'" public Grib2DataSection(bool getData, System.IO.FileStream raf, Grib2GridDefinitionSection gds, Grib2DataRepresentationSection drs, Grib2BitMapSection bms) { //System.out.println( "raf.FilePointer=" + raf.FilePointer() ); // octets 1-4 (Length of DS) length = GribNumbers.int4(raf); //System.out.println( "DS length=" + length ); //System.out.println( "DS calculated end=" + ( raf.getFilePointer() + length -4 )); // octet 5 section 7 section = raf.ReadByte(); //System.out.println( "DS is 7, section=" + section ); if (!getData) { // skip data read //System.out.println( "raf.position before reposition="+raf.getFilePointer()); //System.out.println( "raf.length=" + raf.length() ); // sanity check for erronous ds length if (length > 0 && length < raf.Length) { SupportClass.Skip(raf, length - 5); //System.out.println( "raf.skipBytes = " + (length -5) ); } else { length = 5; // only read length and section } //System.out.println( "raf.position after skip=" + raf.getFilePointer() ); return; } int dtn = drs.DataTemplateNumber; //System.out.println( "DS dtn=" + dtn ); if (dtn == 0 || dtn == 1) { // 0: Grid point data - simple packing // 1: Matrix values - simple packing simpleUnpacking(raf, gds, drs, bms); } else if (dtn == 2) { // 2:Grid point data - complex packing complexUnpacking(raf, gds, drs); } else if (dtn == 3) { // 3: complex packing with spatial differencing complexUnpackingWithSpatial(raf, gds, drs); } else if (dtn == 40 || dtn == 40000) { // JPEG 2000 Stream Format jpeg2000Unpacking(raf, gds, drs, bms); } } // end Grib2DataSection
/// <summary> Creates an ensemble object for the product PDS.</summary> /// <param name="raf">RandomAccessFile. /// </param> /// <param name="parameterNumber"> /// </param> /// <throws> IOException </throws> //UPGRADE_TODO: Class 'java.io.RandomAccessFile' was converted to 'System.IO.FileStream' which has a different behavior. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1073_javaioRandomAccessFile'" internal Grib1Ensemble(System.IO.Stream raf, int parameterNumber) { // skip 12 bytes to start of ensemble SupportClass.Skip(raf, 12); // octet 41 id's ensemble int app = raf.ReadByte(); if (app != 1) { System.Console.Out.WriteLine("not ensemble product"); return; } // octet 42 Type eType = raf.ReadByte(); // octet 43 Identification number eId = raf.ReadByte(); // octet 44 Product Identifier eProd = raf.ReadByte(); // octet 45 Spatial Identifier eSpatial = raf.ReadByte(); if (parameterNumber == 191 || parameterNumber == 192) { // octet 46 Probability product definition epd = raf.ReadByte(); // octet 47 Probability type ept = raf.ReadByte(); // octet 48-51 Probability lower limit epll = GribNumbers.int4(raf); // octet 52-55 Probability upper limit epul = GribNumbers.int4(raf); // octet 56-60 reserved int reserved = GribNumbers.int4(raf); if (eType == 4 || eType == 5) { // octet 61 Ensemble size eSize = raf.ReadByte(); // octet 62 Cluster size eCSize = raf.ReadByte(); // octet 63 Number of clusters eCNumber = raf.ReadByte(); // octet 64 Clustering Method eCMethod = raf.ReadByte(); // octet 65-67 Northern latitude of clustering domain nLat = GribNumbers.int3(raf); // octet 68-70 Southern latitude of clustering domain sLat = GribNumbers.int3(raf); // octet 71-73 Eastern latitude of clustering domain eLat = GribNumbers.int3(raf); // octet 74-76 Western latitude of clustering domain wLat = GribNumbers.int3(raf); if (eType == 4) { // octets 77-86 Cluster Membership cm = new sbyte[10]; SupportClass.ReadInput(raf, cm, 0, cm.Length); } } } }
// *** constructors ******************************************************* /// <summary> Constructs a <tt>Grib1ProductDefinitionSection</tt> object from a raf. /// /// </summary> /// <param name="raf">with PDS content /// /// </param> /// <throws> NotSupportedException if raf contains no valid GRIB file </throws> //UPGRADE_TODO: Class 'java.io.RandomAccessFile' was converted to 'System.IO.FileStream' which has a different behavior. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1073_javaioRandomAccessFile'" public Grib1ProductDefinitionSection(System.IO.FileStream raf) { // octets 1-3 PDS length length = (int)GribNumbers.uint3(raf); //System.out.println( "PDS length = " + length ); // Paramter table octet 4 table_version = raf.ReadByte(); // Center octet 5 center_id = raf.ReadByte(); // octet 6 Generating Process - See Table A process_id = raf.ReadByte(); // octet 7 (id of grid type) - not supported yet grid_id = raf.ReadByte(); //octet 8 (flag for presence of GDS and BMS) int exists = raf.ReadByte(); //BKSystem.IO.BitStream s = new BKSystem.IO.BitStream(8); // s.WriteByte((byte)raf.ReadByte()); // s.Position = 0; // sbyte exists; // s.Read(out exists); // bms_exists = (exists & 64) == 64; gds_exists = (exists & 128) == 128; bms_exists = (exists & 64) == 64; // octet 9 (parameter and unit) parameterNumber = raf.ReadByte(); // octets 10-12 (level) int levelType = raf.ReadByte(); int levelValue1 = raf.ReadByte(); int levelValue2 = raf.ReadByte(); level = new GribPDSLevel(levelType, levelValue1, levelValue2); // octets 13-17 (base time for reference time) int year = raf.ReadByte(); int month = raf.ReadByte(); int day = raf.ReadByte(); int hour = raf.ReadByte(); int minute = raf.ReadByte(); // get info for forecast time // octet 18 Forecast time unit timeUnit = raf.ReadByte(); switch (timeUnit) { case 0: // minute tUnit = "minute"; break; case 1: // hours tUnit = "hour"; break; case 2: // day tUnit = "day"; break; case 3: // month tUnit = "month"; break; case 4: //1 year tUnit = "1year"; break; case 5: // decade tUnit = "decade"; break; case 6: // normal tUnit = "day"; break; case 7: // century tUnit = "century"; break; case 10: //3 hours tUnit = "3hours"; break; case 11: // 6 hours tUnit = "6hours"; break; case 12: // 12 hours tUnit = "12hours"; break; case 254: // second tUnit = "second"; break; default: System.Console.Error.WriteLine("PDS: Time Unit " + timeUnit + " is not yet supported"); break; } // octet 19 & 20 used to create Forecast time p1 = raf.ReadByte(); p2 = raf.ReadByte(); // octet 21 (time range indicator) timeRangeValue = raf.ReadByte(); // forecast time is always at the end of the range //System.out.println( "PDS timeRangeValue =" + timeRangeValue ); switch (timeRangeValue) { case 0: timeRange = "product valid at RT + P1"; forecastTime = p1; break; case 1: timeRange = "product valid for RT, P1=0"; forecastTime = 0; break; case 2: timeRange = "product valid from (RT + P1) to (RT + P2)"; forecastTime = p2; break; case 3: timeRange = "product is an average between (RT + P1) to (RT + P2)"; forecastTime = p2; break; case 4: timeRange = "product is an accumulation between (RT + P1) to (RT + P2)"; forecastTime = p2; break; case 5: timeRange = "product is the difference (RT + P2) - (RT + P1)"; forecastTime = p2; break; case 6: timeRange = "product is an average from (RT - P1) to (RT - P2)"; forecastTime = -p2; break; case 7: timeRange = "product is an average from (RT - P1) to (RT + P2)"; forecastTime = p2; break; case 10: timeRange = "product valid at RT + P1"; // p1 really consists of 2 bytes p1 and p2 forecastTime = p1 = GribNumbers.int2(p1, p2); p2 = 0; break; case 51: timeRange = "mean value from RT to (RT + P2)"; forecastTime = p2; break; default: System.Console.Error.WriteLine("PDS: Time Range Indicator " + timeRangeValue + " is not yet supported"); break; } // octet 22 & 23 int avgInclude = GribNumbers.int2(raf); // octet 24 int avgMissing = raf.ReadByte(); // octet 25 int century = raf.ReadByte() - 1; // octet 26, sub center subcenter_id = raf.ReadByte(); // octets 27-28 (decimal scale factor) decscale = GribNumbers.int2(raf); refTime = new DateTime(century * 100 + year, month, day, hour, minute, 0, DateTimeKind.Utc); baseTime = refTime; referenceTime = refTime.ToString(dateFormat); //System.out.println( "PDS referenceTime raw = " + //(century * 100 + year) +"-" + //month + "-" + day + "T" + hour +":" + minute +":00Z" ); //System.out.println( "PDS referenceTime = " + referenceTime ); // TODO /* * parameter_table = GribPDSParamTable.getParameterTable(center_id, subcenter_id, table_version); * parameter = parameter_table.getParameter(parameterNumber); */ parameter = new Parameter(); if (center_id == 7 && subcenter_id == 2) { // ensemble product epds = new Grib1Ensemble(raf, parameterNumber); } // Special handling of 2D Wave Spectra (single) else if (table_version == 140 && parameterNumber == 251) { SupportClass.Skip(raf, 12); int extDef = raf.ReadByte(); // Extension definition // Read ECMWF extension for 2D wave spectra single if (extDef == 13) { waveSpectra2DDirFreq = new Grib1WaveSpectra2DDirFreq(raf); } } else if (length != 28) { // check if all bytes read in section //lengthErr = true; int extra; for (int i = 29; i <= length; i++) { extra = raf.ReadByte(); } } } // end Grib1ProductDefinitionSection
// *** constructors ******************************************************* /// <summary> Constructs a Grib2ProductDefinitionSection object from a raf. /// /// </summary> /// <param name="raf">RandomAccessFile with PDS content /// /// </param> /// <throws> IOException if raf contains no valid GRIB file </throws> //UPGRADE_TODO: Class 'java.io.RandomAccessFile' was converted to 'System.IO.FileStream' which has a different behavior. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1073_javaioRandomAccessFile'" public Grib2ProductDefinitionSection(System.IO.FileStream raf) { // octets 1-4 (Length of PDS) length = GribNumbers.int4(raf); //System.out.println( "PDS length=" + length ); // octet 5 section = raf.ReadByte(); //System.out.println( "PDS is 4, section=" + section ); // octet 6-7 coordinates = GribNumbers.int2(raf); //System.out.println( "PDS coordinates=" + coordinates ); // octet 8-9 productDefinition = GribNumbers.int2(raf); //System.out.println( "PDS productDefinition=" + productDefinition ); switch (productDefinition) { // Analysis or forecast at a horizontal level or in a horizontal // layer at a point in time case 0: case 1: case 2: case 3: case 4: case 5: case 6: case 7: case 8: { // octet 10 parameterCategory = raf.ReadByte(); //System.out.println( "PDS parameterCategory=" + //parameterCategory ); // octet 11 parameterNumber = raf.ReadByte(); //System.out.println( "PDS parameterNumber=" + parameterNumber ); // octet 12 typeGenProcess = raf.ReadByte(); //System.out.println( "PDS typeGenProcess=" + typeGenProcess ); // octet 13 backGenProcess = raf.ReadByte(); //System.out.println( "PDS backGenProcess=" + backGenProcess ); // octet 14 analysisGenProcess = raf.ReadByte(); //System.out.println( "PDS analysisGenProcess=" + //analysisGenProcess ); // octet 15-16 hoursAfter = GribNumbers.int2(raf); //System.out.println( "PDS hoursAfter=" + hoursAfter ); // octet 17 minutesAfter = raf.ReadByte(); //System.out.println( "PDS minutesAfter=" + minutesAfter ); // octet 18 timeRangeUnit = raf.ReadByte(); //System.out.println( "PDS timeRangeUnit=" + timeRangeUnit ); // octet 19-22 forecastTime = GribNumbers.int4(raf); //System.out.println( "PDS forecastTime=" + forecastTime ); // octet 23 typeFirstFixedSurface = raf.ReadByte(); //System.out.println( "PDS typeFirstFixedSurface=" + // typeFirstFixedSurface ); // octet 24 int scaleFirstFixedSurface = raf.ReadByte(); //System.out.println( "PDS scaleFirstFixedSurface=" + // scaleFirstFixedSurface ); // octet 25-28 int valueFirstFixedSurface = GribNumbers.int4(raf); //System.out.println( "PDS valueFirstFixedSurface=" + // valueFirstFixedSurface ); //UPGRADE_WARNING: Data types in Visual C# might be different. Verify the accuracy of narrowing conversions. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1042'" FirstFixedSurfaceValue = (float)((scaleFirstFixedSurface == 0 || valueFirstFixedSurface == 0)?valueFirstFixedSurface:System.Math.Pow(valueFirstFixedSurface, -scaleFirstFixedSurface)); // octet 29 typeSecondFixedSurface = raf.ReadByte(); //System.out.println( "PDS typeSecondFixedSurface=" + //typeSecondFixedSurface ); // octet 30 int scaleSecondFixedSurface = raf.ReadByte(); //System.out.println( "PDS scaleSecondFixedSurface=" + //scaleSecondFixedSurface ); // octet 31-34 int valueSecondFixedSurface = GribNumbers.int4(raf); //System.out.println( "PDS valueSecondFixedSurface=" + //valueSecondFixedSurface ); //UPGRADE_WARNING: Data types in Visual C# might be different. Verify the accuracy of narrowing conversions. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1042'" SecondFixedSurfaceValue = (float)((scaleSecondFixedSurface == 0 || valueSecondFixedSurface == 0)?valueSecondFixedSurface:System.Math.Pow(valueSecondFixedSurface, -scaleSecondFixedSurface)); // Individual ensemble forecast, control and perturbed, at a // horizontal level or in a horizontal layer at a point in time if (productDefinition == 1) { System.Console.Out.WriteLine("PDS productDefinition == 1 not done"); //Derived forecast based on all ensemble members at a horizontal // level or in a horizontal layer at a point in time } else if (productDefinition == 2) { System.Console.Out.WriteLine("PDS productDefinition == 2 not done"); // Derived forecasts based on a cluster of ensemble members over // a rectangular area at a horizontal level or in a horizontal layer // at a point in time } else if (productDefinition == 3) { System.Console.Out.WriteLine("PDS productDefinition == 3 not done"); // Derived forecasts based on a cluster of ensemble members // over a circular area at a horizontal level or in a horizontal // layer at a point in time } else if (productDefinition == 4) { System.Console.Out.WriteLine("PDS productDefinition == 4 not done"); // Probability forecasts at a horizontal level or in a horizontal // layer at a point in time } else if (productDefinition == 5) { System.Console.Out.WriteLine("PDS productDefinition == 5 not done"); // Percentile forecasts at a horizontal level or in a horizontal layer // at a point in time } else if (productDefinition == 6) { System.Console.Out.WriteLine("PDS productDefinition == 6 not done"); // Analysis or forecast error at a horizontal level or in a horizontal // layer at a point in time } else if (productDefinition == 7) { System.Console.Out.WriteLine("PDS productDefinition == 7 not done"); // Average, accumulation, and/or extreme values at a horizontal // level or in a horizontal layer in a continuous or non-continuous // time interval } else if (productDefinition == 8) { //System.out.println( "PDS productDefinition == 8 " ); // 35-41 bytes int year = GribNumbers.int2(raf); int month = (raf.ReadByte()) - 1; int day = raf.ReadByte(); int hour = raf.ReadByte(); int minute = raf.ReadByte(); int second = raf.ReadByte(); //System.out.println( "PDS date:" + year +":" + month + //":" + day + ":" + hour +":" + minute +":" + second ); // 42 - 46 int timeRanges = raf.ReadByte(); //System.out.println( "PDS timeRanges=" + timeRanges ) ; int missingDataValues = GribNumbers.int4(raf); //System.out.println( "PDS missingDataValues=" + missingDataValues ) ; // 47 - 48 int outmostTimeRange = raf.ReadByte(); //System.out.println( "PDS outmostTimeRange=" + outmostTimeRange ) ; int missing = raf.ReadByte(); //System.out.println( "PDS missing=" + missing ) ; // 49 - 53 int statisticalProcess = raf.ReadByte(); //System.out.println( "PDS statisticalProcess=" + statisticalProcess ) ; int timeIncrement = GribNumbers.int4(raf); //System.out.println( "PDS timeIncrement=" + timeIncrement ) ; // 54 - 58 int indicatorTR = raf.ReadByte(); //System.out.println( "PDS indicatorTR=" + indicatorTR ) ; int lengthTR = GribNumbers.int4(raf); if (timeRanges == 1) { this.forecastTime += calculateIncrement(statisticalProcess, timeIncrement); } //System.out.println( "PDS lengthTR=" + lengthTR ) ; //int indicatorSF = raf.read(); //System.out.println( "PDS indicatorSF=" + indicatorSF ) ; //int incrementSF = GribNumbers.int4( raf ); //System.out.println( "PDS incrementSF=" + incrementSF ) ; } break; } // cases 0-8 // Radar product case 20: { parameterCategory = raf.ReadByte(); //System.out.println( "PDS parameterCategory=" + //parameterCategory ); parameterNumber = raf.ReadByte(); //System.out.println( "PDS parameterNumber=" + parameterNumber ); typeGenProcess = raf.ReadByte(); //System.out.println( "PDS typeGenProcess=" + typeGenProcess ); backGenProcess = raf.ReadByte(); //System.out.println( "PDS backGenProcess=" + backGenProcess ); hoursAfter = GribNumbers.int2(raf); //System.out.println( "PDS hoursAfter=" + hoursAfter ); minutesAfter = raf.ReadByte(); //System.out.println( "PDS minutesAfter=" + minutesAfter ); timeRangeUnit = raf.ReadByte(); //System.out.println( "PDS timeRangeUnit=" + timeRangeUnit ); forecastTime = GribNumbers.int4(raf); //System.out.println( "PDS forecastTime=" + forecastTime ); break; } // case 20 // Satellite Product case 30: { parameterCategory = raf.ReadByte(); //System.out.println( "PDS parameterCategory=" + parameterCategory ); parameterNumber = raf.ReadByte(); //System.out.println( "PDS parameterNumber=" + parameterNumber ); typeGenProcess = raf.ReadByte(); //System.out.println( "PDS typeGenProcess=" + typeGenProcess ); backGenProcess = raf.ReadByte(); //System.out.println( "PDS backGenProcess=" + backGenProcess ); nb = raf.ReadByte(); //System.out.println( "PDS nb =" + nb ); for (int j = 0; j < nb; j++) { SupportClass.Skip(raf, 10); } break; } // case 30 // CCITTIA5 character string case 254: { parameterCategory = raf.ReadByte(); //System.out.println( "PDS parameterCategory=" + //parameterCategory ); parameterNumber = raf.ReadByte(); //System.out.println( "PDS parameterNumber=" + parameterNumber ); //numberOfChars = GribNumbers.int4( raf ); //System.out.println( "PDS numberOfChars=" + //numberOfChars ); break; } // case 254 default: break; } // end switch }
/// <summary> Constructs a Grib1BinaryDataSection object from a raf. /// A bit map is defined. /// /// </summary> /// <param name="raf">raf with BDS content /// </param> /// <param name="decimalscale">the exponent of the decimal scale /// </param> /// <param name="bms">bit map section of GRIB record /// /// </param> /// <throws> NotSupportedException if stream contains no valid GRIB file </throws> //UPGRADE_TODO: Class 'java.io.RandomAccessFile' was converted to 'System.IO.FileStream' which has a different behavior. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1073_javaioRandomAccessFile'" public Grib1BinaryDataSection(System.IO.FileStream raf, int decimalscale, Grib1BitMapSection bms) { // octets 1-3 (section length) length = (int)GribNumbers.uint3(raf); //System.out.println( "BDS length = " + length ); // octet 4, 1st half (packing flag) int unusedbits = raf.ReadByte(); // TODO Check this!!! if ((unusedbits & 192) != 0) { throw new NotSupportedException("BDS: (octet 4, 1st half) not grid point data and simple packing "); } // octet 4, 2nd half (number of unused bits at end of this section) unusedbits = unusedbits & 15; //System.out.println( "BDS unusedbits = " + unusedbits ); // octets 5-6 (binary scale factor) int binscale = GribNumbers.int2(raf); // octets 7-10 (reference point = minimum value) float refvalue = GribNumbers.float4(raf); // octet 11 (number of bits per value) int numbits = raf.ReadByte(); //System.out.println( "BDS numbits = " + numbits ); if (numbits == 0) { isConstant = true; } //System.out.println( "BDS isConstant = " + isConstant ); // *** read values ******************************************************* //UPGRADE_WARNING: Data types in Visual C# might be different. Verify the accuracy of narrowing conversions. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1042'" float ref_Renamed = (float)(System.Math.Pow(10.0, -decimalscale) * refvalue); //UPGRADE_WARNING: Data types in Visual C# might be different. Verify the accuracy of narrowing conversions. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1042'" float scale = (float)(System.Math.Pow(10.0, -decimalscale) * System.Math.Pow(2.0, binscale)); if (bms != null) { bool[] bitmap = bms.Bitmap; values = new float[bitmap.Length]; for (int i = 0; i < bitmap.Length; i++) { if (bitmap[i]) { if (!isConstant) { values[i] = ref_Renamed + scale * bits2UInt(numbits, raf); } else { // rdg - added this to handle a constant valued parameter values[i] = ref_Renamed; } } else { values[i] = Grib1BinaryDataSection.UNDEFINED; } } } else { // bms is null if (!isConstant) { //System.out.println( "BDS values.size = " + //(((length - 11) * 8 - unusedbits) / numbits)); values = new float[((length - 11) * 8 - unusedbits) / numbits]; for (int i = 0; i < values.Length; i++) { values[i] = ref_Renamed + scale * bits2UInt(numbits, raf); } } else { // constant valued - same min and max int x = 0, y = 0; raf.Seek(raf.Position - 53, System.IO.SeekOrigin.Begin); // return to start of GDS length = (int)GribNumbers.uint3(raf); if (length == 42) { // Lambert/Mercator offset SupportClass.Skip(raf, 3); x = GribNumbers.int2(raf); y = GribNumbers.int2(raf); } else { SupportClass.Skip(raf, 7); length = (int)GribNumbers.uint3(raf); if (length == 32) { // Polar sterographic SupportClass.Skip(raf, 3); x = GribNumbers.int2(raf); y = GribNumbers.int2(raf); } else { x = y = 1; System.Console.Out.WriteLine("BDS constant value, can't determine array size"); } } values = new float[x * y]; for (int i = 0; i < values.Length; i++) { values[i] = ref_Renamed; } } } } // end Grib1BinaryDataSection
/// <summary> Constructs a <tt>Grib1GridDefinitionSection</tt> object from a raf. /// /// </summary> /// <param name="raf">RandomAccessFile with GDS content /// /// </param> /// <throws> NoValidGribException if raf contains no valid GRIB info </throws> //UPGRADE_TODO: Class 'java.io.RandomAccessFile' was converted to 'System.IO.FileStream' which has a different behavior. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1073_javaioRandomAccessFile'" public Grib1GridDefinitionSection(System.IO.FileStream raf) { int reserved; // used to read empty space // octets 1-3 (Length of GDS) length = (int)GribNumbers.uint3(raf); if (length == 0) { // there's a extra byte between PDS and GDS SupportClass.Skip(raf, -2); length = (int)GribNumbers.uint3(raf); } //System.out.println( "GDS length = " + length ); // TODO Fix Checksum stuff // get byte array for this gds, then reset raf to same position // calculate checksum for this gds via the byte array /* * long mark = raf.Position; * sbyte[] dst = new sbyte[length - 3]; * SupportClass.ReadInput(raf, dst, 0, dst.Length); * raf.Seek(mark, System.IO.SeekOrigin.Begin); * //UPGRADE_ISSUE: Class 'java.util.zip.CRC32' was not converted. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1000_javautilzipCRC32'" * //UPGRADE_ISSUE: Constructor 'java.util.zip.CRC32.CRC32' was not converted. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1000_javautilzipCRC32'" * CRC32 cs = new CRC32(); * //UPGRADE_ISSUE: Method 'java.util.zip.CRC32.update' was not converted. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1000_javautilzipCRC32'" * cs.update(dst); * //UPGRADE_ISSUE: Method 'java.util.zip.CRC32.getValue' was not converted. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1000_javautilzipCRC32'" * checksum = System.Convert.ToString(cs.getValue()); * //System.out.println( "GDS checksum =" + checksum ); */ // octets 4 NV int NV = raf.ReadByte(); //System.out.println( "GDS NV = " + NV ); // octet 5 PL the location (octet number) of the list of numbers of points in each row P_VorL = raf.ReadByte(); //System.out.println( "GDS PL = " + P_VorL ); // octet 6 (grid type) type = raf.ReadByte(); //System.out.println( "GDS grid type = " + type ); name = getName(type); if (type != 50) { // values same up to resolution // octets 7-8 (Nx - number of points along x-axis) nx = GribNumbers.int2(raf); nx = (nx == -1)?1:nx; // octets 9-10 (Ny - number of points along y-axis) ny = GribNumbers.int2(raf); ny = (ny == -1)?1:ny; // octets 11-13 (La1 - latitude of first grid point) lat1 = GribNumbers.int3(raf) / 1000.0; // octets 14-16 (Lo1 - longitude of first grid point) lon1 = GribNumbers.int3(raf) / 1000.0; // octet 17 (resolution and component flags). See Table 7 resolution = raf.ReadByte(); } switch (type) { // Latitude/Longitude grids , Arakawa semi-staggered e-grid rotated // Arakawa filled e-grid rotated case 0: case 4: case 40: case 201: case 202: // octets 18-20 (La2 - latitude of last grid point) lat2 = GribNumbers.int3(raf) / 1000.0; // octets 21-23 (Lo2 - longitude of last grid point) lon2 = GribNumbers.int3(raf) / 1000.0; // octets 24-25 (Dx - Longitudinal Direction Increment ) dx = GribNumbers.int2(raf) / 1000.0; // octets 26-27 (Dy - Latitudinal Direction Increment ) // Np - parallels between a pole and the equator if (type == 4) { np = GribNumbers.int2(raf); } else { dy = GribNumbers.int2(raf) / 1000.0; } // octet 28 (Scanning mode) See Table 8 scan = raf.ReadByte(); // octet 29-32 reserved reserved = GribNumbers.int4(raf); if (length > 32) { // getP_VorL(raf); // Vertical coordinates (NV) and thinned grids (PL) not supported - skip this SupportClass.Skip(raf, length - 32); } break; // end Latitude/Longitude grids case 1: // Mercator grids // octets 18-20 (La2 - latitude of last grid point) lat2 = GribNumbers.int3(raf) / 1000.0; // octets 21-23 (Lo2 - longitude of last grid point) lon2 = GribNumbers.int3(raf) / 1000.0; // octets 24-26 (Latin - latitude where cylinder intersects the earth latin1 = GribNumbers.int3(raf) / 1000.0; // octet 27 reserved reserved = raf.ReadByte(); // octet 28 (Scanning mode) See Table 8 scan = raf.ReadByte(); // octets 29-31 (Dx - Longitudinal Direction Increment ) dx = GribNumbers.int3(raf); // octets 32-34 (Dx - Longitudinal Direction Increment ) dy = GribNumbers.int3(raf); // octet 35-42 reserved reserved = GribNumbers.int4(raf); reserved = GribNumbers.int4(raf); if (length > 42) { // getP_VorL(raf); // Vertical coordinates (NV) and thinned grids (PL) not supported - skip this SupportClass.Skip(raf, length - 42); } break; // end Mercator grids case 3: // Lambert Conformal // octets 18-20 (Lov - Orientation of the grid - east lon parallel to y axis) lov = GribNumbers.int3(raf) / 1000.0; // octets 21-23 (Dx - the X-direction grid length) See Note 2 of Table D dx = GribNumbers.int3(raf); // octets 24-26 (Dy - the Y-direction grid length) See Note 2 of Table D dy = GribNumbers.int3(raf); // octets 27 (Projection Center flag) See Note 5 of Table D proj_center = raf.ReadByte(); // octet 28 (Scanning mode) See Table 8 scan = raf.ReadByte(); // octets 29-31 (Latin1 - first lat where secant cone cuts spherical earth latin1 = GribNumbers.int3(raf) / 1000.0; // octets 32-34 (Latin2 - second lat where secant cone cuts spherical earth) latin2 = GribNumbers.int3(raf) / 1000.0; // octets 35-37 (lat of southern pole) latsp = GribNumbers.int3(raf) / 1000.0; // octets 38-40 (lon of southern pole) lonsp = GribNumbers.int3(raf) / 1000.0; // octets 41-42 reserved = GribNumbers.int2(raf); if (length > 42) { // getP_VorL(raf); // Vertical coordinates (NV) and thinned grids (PL) not supported - skip this SupportClass.Skip(raf, length - 42); } break; // end Lambert Conformal case 5: // Polar Stereographic grids // octets 18-20 (Lov - Orientation of the grid - east lon parallel to y axis) lov = GribNumbers.int3(raf) / 1000.0; // octets 21-23 (Dx - Longitudinal Direction Increment ) dx = GribNumbers.int3(raf); // octets 24-26(Dy - Latitudinal Direction Increment ) dy = GribNumbers.int3(raf); // octets 27 (Projection Center flag) See Note 5 of Table D proj_center = raf.ReadByte(); // octet 28 (Scanning mode) See Table 8 scan = raf.ReadByte(); // octet 29-32 reserved reserved = GribNumbers.int4(raf); if (length > 32) { // getP_VorL(raf); // Vertical coordinates (NV) and thinned grids (PL) not supported - skip this SupportClass.Skip(raf, length - 32); } break; // end Polar Stereographic grids //rotated grid case 10: //18-20 this.lat2 = GribNumbers.int3(raf) / 1000.0; //21-23 this.lon2 = GribNumbers.int3(raf) / 1000.0; //24-25 var tmpDx = GribNumbers.int2(raf); this.dx = tmpDx == -1 ? -9999.0 : tmpDx / 1000.0; //26-27 var tmpDy = GribNumbers.int2(raf); this.dy = tmpDy == -1 ? -9999.0 : tmpDy / 1000.0; //28 scan = raf.ReadByte(); //29-32 GribNumbers.int4(raf); //33-35 this.latsp = GribNumbers.int3(raf) / 1000.0; //36-38 this.lonsp = GribNumbers.int3(raf) / 1000.0; //39-42 - angle GribNumbers.float4(raf); if (length > 42) { SupportClass.Skip(raf, length - 42); } break; default: System.Console.Out.WriteLine("Unknown Grid Type : " + type); break; } // end switch grid_type //System.out.println( "scan=" + scan ); if ((scan & 63) != 0) { throw new NoValidGribException("GDS: This scanning mode (" + scan + ") is not supported."); } } // end Grib1GridDefinitionSection( raf )