GRIBEX - History of changes

000370 April 30, 2009

• grchk1.F: added class 17 la LACE ALADIN,18 yt YOTC

000360 September 29, 2008

• getsetValues.c: rename copyName to copyNameLoc due to multiple definition in shared librar
  y
• gribtemplates: added localDefinitionTemplate_098_000_031
• grchk1.F: added stream 1040 Ensemble Forecast Hindcast Statistics (EFHS)
• grchk1.F: added stream 1240 Eurosip Monthly Means, 1241 Eurosip Hindcast Monthly Means
• dmesec2 emesec2: Bug fix decode/encode Mercator data

000350 February 12, 2008

• calcop.F: packcf.F unpackcf.F : Increased maximum Truncation to be handled to T2047
• grchk1.F: added class 16 dt (Data Targeting System)
• gribtemplate:
  • added - localDefinitionTemplate_214_098_245
  • updated - localDefinitionTemplate_098_235_015, localDefinitionTemplate_098_235_004, localDefinitionTemplate_098_235_016

  PBIO routines

  • pbio.c: Bug fix - setvbuf to be called every time new file is open
  • jopnggsm.c: Bug fix - calculation of coefficient file size
  • gbyte_alpha.c: Bug fix for 64 bits packing

000340 September 03, 2007

• Added header with LGPL license rules in each program

000320 April 05, 2007

• grchk1.F: added streams
  • 1032 efho Ensemble Forecast Hindcast Overlap
  • 1033 enfh Ensemble Forecast Hindcasts
  • 1078 ewho Ensemble Forecast Wave Hindcast Overlap
  • 1079 enwh Ensemble Forecast Wave Hindcasts
• grchk1.F: added satellite identifiers 56,171
• grchk1.F: added class 15 SREPS Short-Range Ensemble Prediction System
• gribtemplates/localDefinitionTemplate_214_098_244:Definition for SREPS from Spain
• fortranInterface.c: Support for A8 (8 character string in local definitions)
• handleLocalDefinitions.c: Support for A8 (8 character string in local definitions)
• pbio/pbio.c: New version of this routine

000310 November 22, 2006

• Updated Grib tables
• Added Local Definitions 30 - Forecasting Systems with Variable Resolution

000300 March 14, 2006

• Fix in Gribex.F, computation of scale factor with a minimum equal to the reference value.
• grchk1.F: Add STREAM = 1075 Monthly means of Daily forecast accumulations
• grchk1.F: Add Time range indicator = 120 Average of N Forecast
• grchk1.F: Add satellite identifier 199
• grchk1.F:Add CLASS = 12 TIGGE
• 211 level type has been replaced with 209 because 211 is in use by NCEP
• gdecode.c : Bug fix for Section 4 length
• gribex.F: to read correctly encoded big gribs (where first bit of section 0 length does not mean size*120), anyway, gribex will pack it with the "*120" workaround"
• grchk1.F: Fix check for satellites images
• d2ordr.F: Fix to avoid bad check for large grids(second order packing).
• grchk1.F: new streams have been added 1030 - Ensemble Data Assimilation,
  1088 - Ensemble Wave Data Assimilation
  

000281 Jun 20, 2005

• Modified the behavior for the accessing the local definition files
  it was : "if -file dont exist- go to defaultDir if you have not found the file within $LOCAL_DEFINITION_TEMPLATES"
  It is now : "if file dont exist try to reconstruct the file name with center = 098(ECMWF) and subcenter = 0" (findLocalDefinitionFile.c)
• Changed EMOSLIB version number to 000281
• Fix to check that the real GRIB length is consistent with the lenght value in the header file for very large grib products (ie >8Mb)
  This might not be the case because for very large grib products, the value in section 0 is multiplied by 120. If this is not consistent, KSEC0(1) is updated with the real value (gribex.F).
• Fix in dsect4a.F (2nd order pack of spectral coefficients) to leave the file pointer at the end of the data when very large product is found. The former procedure used the value in the gribfile, section 4 to position the pointer at the end when the function returns. However, because it is not possible to code a value bigger than 8388607, therefore the restored position was restored from bad value. The function now restore the position from the point reached after the last coefficient is read (only if big product is found).
• Added : STREAM - 1024 daily archive hindcast , 1034 ensemble forecast overlap, 1039 ensemble forecast hindcasts (changed name), 1084 ensemble wave hindcast (changed name), 1085 wave hindcast, 1086 wave ensemble forecast overlap, 1087. TYPE - 29 Probability Boundaries, 36 Probability Distribution, 37 Cluster information. CLASS - 10 cosmo leps, CLASS 11 Ensemles, TYPE of LEVEL - 211 Ocean wave data, 212 - Wave spectra (grchk1.F)
• Fix in handleLocalDefinitions.c to add support for arrays of signed integers.

000263 September 22, 2004

• Correction to encoding/decoding of satellite image GRIBs (dsvsec2.F, esvsec2.F, gribex.F).
  
• Added handling for Mercator GRIBs (dmesec2.F, emesec2.F, gribex.F).
• Corrected display of parameter information when centre is not ECMWF but sub-centre is ECMWF (chktab2.F).
• Corrected display of step when time unit = 2 (days) (gdecode.c).
• Add satellite identifier 55 to table of allowed ids (grchk1.F).
• Adjust the tolerance check when calculating integer scaling (gribex.F).
  

000240 September 1, 2003

• Fix handling of the post-auxiliary array flag during encoding/decoding of ECMWF local definition 4 (ocean model data).
  
• Add ECMWF local definition 21 (sensitive area products) for TOST.
  

000230 August 1, 2003

• Add ECMWF local definition 20 (4D variational increments)
  
• Add class 8 for ELDAS
  
• Add types 34 (gridded observations) and 90 - 99 (Mars development)
  
• Add streams 1061 - 1069 (Mars development) and 1100 and 1101 (ELDAS)
  
• Add ECMWF local table 200 (Differences between first guess and analysis)
  

000221 May 22, 2003

• Add ECMWF local definition 20 (4D variational increments)
  
• Add class 8 for ELDAS
  
• Add types 34 (gridded observations) and 90 - 99 (Mars development)
  
• Add streams 1061 - 1069 (Mars development) and 1100 and 1101 (ELDAS)
  
• Add ECMWF local table 200 (Differences between first guess and analysis)
  

000220 March 20, 2003

• Use INTEGER_8 to set number of bits per integer (JPBITS).
  
• Use E_DEF_ and D_DEF_ for local definitions (gdecode).
  
• Add types 41 and 42
  
• Add local definition 190
  
• In GRIBEX
• Force the rounding calculation for very large products whose length won't fit into 23 bits, eg lat/lon 0.01 x 0.01 degree grid.
  
• Modify handling of length of local definition 14.
• Add EFI order to local definition 19.
  
• Display data bytes for local definitions 190 and 191
  
• Fix display of expver on linux and 64-bit hp
  
• Add total number of ensembles in definition 15

000210 February 28, 2003

• Add ECMWF local definition 191 (free format data descriptor) and local definition 19 (extreme forecast index).
  
• Add generic Member State class 121.
  
• Return the number of values (bytes) in ksec4(21) when the 'B' option is used (GRIBEX).
  
• use C routines for encoding/decoding ECMWF local section 1 extensions

000200 October 11, 2002

• Allow compilation without using POINTERs (USE_NO_POINTERS).
  
• Add facility to redirect output of GRIBEX diagnostics or GRPRS routines print of GRIB headers to fortran file using environment variable GRPRS_STREAM (default is stdout).
  
• Read all environment variables in one routine (grsdef).
  
• Handle little-endian ordering of letters in the experiment version number in ECMWF local extensions to GRIB section 1.
  
• Modify coordinate printout in ECMWF local definition 4 (ocean): I9 -> I10.
  
• Add ECMWF GRIB local definition 50 (Member State data)
  
• Add stream 1074, MSDA, monthly/daily standard deviation and covariance

000180 March 26, 2002

• Remove emosnum.F from gribex routines.

000170 March 26, 2002

• Modify print of parameter description to include name and units; and fix diplay of experiment version in section 1 on linux (grprs1).
  
• Handle encoding when calculated scaling factor clashes with number of bits per value (happens for fields whose range is a power of 2 and looses the max value).
  
• Speedup decoding of pseudo-grib (decops2 and offset2).
  
• Adjust ECMWF definition 18 (multi-analysis ensemble data) for analyses.
  
• Ensure 'R' decode option uses 'nearest neighbour' interpolation for vegetation parameters (table 128, numbers 27/28/29/30).
  

000150 March 7, 2001

• Allow 'R' option to work on rotated/stretched gaussian and lat/long grids.
  
• Fix 'J' option to report correctly the numbers of missing and non-missing values in bitmapped fields.
  
• Add local definition 18 for multi-analysis ensembles.
  
• Add streams 1093, 1094, 1095 and 1096.
  
• Fix calculation of length of section 1 in local definition 16.
  
• Change GRIBEX version number to 13.034.
  

000140 February 8, 2001

• Modify C versions of gbytes/gbyte/sbytes/sbyte to accept void pointers to avoid compiler diagnostics when REAL arrays used instead of INTEGER.
  

000120 October 31, 2000

• Add diagnostics when calculation of reference value in GRIB floating point format does not find value lower than original number.
  
• Expand parameter description character string lengths (to 81).
  
• Change GRIBEX version number to 13.033.
  
• Refine selection of debug output via GRIBEX_DEBUG (= ON/OFF, 1, 2).

000110 September 4, 2000

• Extend limit checks to allow upto T1023, N480.
  
• Make loop calculating P power coefficents in complex packing scalar on VPPS to avoid differences on 700 and 5000.

000100 July 5, 2000

• Replace calls to (fortran) FREE by JFREE.
  
• Remove redundant routines (search, dswmrs, ldc1cr, delsp, ccf1cr, analsw, parval).
  
• Add checking switch via environment variable GRIBEX_CHECK (ON/OFF).
  
• Add streams 1026 (SCDA) and 1027 (SCWV).

000090 May 16, 2000

• Print 'Year of reference time of data' in GRIB section 1 as YY and as (YYYY). (grprs1)
• Add ECMWF GRIB section 1 local definition 16 (seasonal monthly mean forecasts).
• Add ECMWF GRIB section 1 local definition 15 (ocean ensemble seasonal forecasts).
  
• Allow 4-digit counts in printout of 'no. of pts along latitude' for reduced gaussian grid GRIB header section 2.
  
• Allow stream 1090 to use post-auxiliary array in ECMWF local definition 4. In this definition, ensemble forecast number occupies 2 bytes; GRIB section 1 printout adjusted accordingly.
  
• Allow stream 1026, short cut-off daily archive (boundary conditions).
  
• Use 'old-style' calculation of reference values (remove use of LEAST value when calculation of reference fails).
  
• Add CLASS=6 for DEMETER and CLASS=7 for PROVOST.
  
• Add system and method number capability to ECMWF local definition 13.
  
• Add stream number 1082 for wave seasonal forecast.

000080 January 25, 2000

• Add C routines for encoding/decoding sign bit in GRIB octets (csgnbt.c and dsgnbt.c).
• Correct printing of GRIB section 1 to use century and year for displayed 4-digit year (grprs1).
• Initialise all local variables.
• Add debug switch via environment variable GRIBEX_DEBUG (ON/OFF).
• Allow 2-byte ensemble number in ECMWF ensemble seasonal forecasts (stream 1090) in ECMWF GRIB section 1 local definitions 1 and 4.

000070 November 9, 1999

• Check parameter number and local table 2 number against values in WMO standard versions of table 2 (parameters), and local tables from other centres. The directory of tables may be specified through a preprocessing directive or via an environment variable.
• Changes to the initialisation mechanism for global variables.
• KSEC4(22:24) used to give access to extra GRIBEX variables (as KSEC4(21) for non-missing points of bit-masked field, namely bit-map pointer to primary bit-map (if any), full length and unused bit count of section 4.
• Allow values 1 to 255 for identification of centre (code table 0).
• Allow for degenerate case when 'number of bits per packed value' is zero; use 8 bits per value instead.
• Modify binary scale factor computation when 'number of bits per packed value' is close to number of bits per computer integer.
• Decimal scale factor taken into account for 'X' option.
• Use of pre-computed bit widths for some optimisations.
• Enable simple packing as a backup method for HOPER='C' when no 2nd-order method may be applied.
• Reintroduction of run-time switches used to keep a control on extended methods and sub-methods (these switches being now set initially to on for consistency with GRIBEX 13.023) in exhaustive mode (kept to off).

000060 October 10, 1999

• Fix parameter table checking when invalid table number given more than once.
• Allow values 1 to 255 for identification of centre (code table 0).
• Add HP double precision compilation options.
• Add linux compilation options.
  

000050 May 25, 1999

• Add post-auxiliary array in ocean products.
• Add post-auxiliary array in ocean products.
• Allow type 72 (ocean flux) in GRIB section 1.
• Allow comment lines at start of local table 2 description files. (Comment lines subsequently deleted).
• Check parameter number and local table 2 number against values in local definitions files.
• Remove extend2o preprocessing directive(s); 2nd order packing now standard.

000040 March 16, 1999

• Modify calculation of exponent and mantissa for GRIB floating point values (eg. reference value).
• Add SWHP (229) and MWPP (232) to local table 2 version 131 (probabilities).
• Add ECMWF ocean data type (=192) with new GRIB section 2 definition; also add
• MARS types FR/AR (70/71).
• Handle GRIB section 2 encoding/decoding in separate subroutines.
• Check parameter number and local table 2 number against values in local definitions files. (Removed 1999/03/29).

000020 November 20, 1998

• Pseudo GRIB encoding/decoding to store/recover dates internally held with reference to 1900 as baseline date.
• ECMWF local definition 6 (SSTs used by analysis) modified to hold dates internally with reference to 19000000 as baseline date.
• Compile gbyte with optimization (and no debug).
• 2nd order packing routines optimized.
  
• Decoding of complex packed spectral data in GRIB edition 0 added.
  
• GRIBEX version number changed to 13.021.
  
• Update local_table_2_versions_128/129/140/160.
  
• Add ECMWF GRIB section 1 local definition 12 (means, averages, etc).
  
• Add EMOS version number in local table 2 definition pathname.

000010 October 1, 1998

• add print of non-missing values for field using a bitmap
• optimize vectorization on VPP of bitmap handling routines
• fix logic for encoding level types 102/200/201
• make all parameter numbers 1 to 255 available for tables 128/129
• allow Member State class numbers 100 to 120
• add class 5 for ERA 40
• add streams 1072 and 1073 for ERA monthly data
• make debug printout from bit handling routine (inxbit) a compile time option
• add ECMWF local definition 14 for brightness temperatures
• gribex version number set to 13.011
• expand internal array size for handling 2nd order packing
• display missing year as 255 in GRIB header printout
• change comments in gsbite.F header to emphasise that KMASK argument is an array and must be suitably dimensioned
• add a Fortran-callable routine (yyyy2cy) to convert YYYY year to century and year for gribex call
  

000000 June 23, 1998