void gb2_param ( char *wmovartbl, char *lclvartbl, Gribmsg *cmsg,
char *param, int *scal, float *msng, int *iret )
/************************************************************************
* gb2_param *
* *
* This routine gets the parameter values *
* from the GRIB2 PDS, and obtains the GEMPAK parameter information *
* from the appropriate GRIB2 parameter table. *
* *
* If either wmovartbl or lclvartbl are NULL, the default tables are *
* read. *
* *
* gb2_param ( wmovartbl, lclvartbl, cmsg, param, scal, msng, iret ) *
* *
* Input parameters: *
* *wmovartbl char WMO parameter table *
* *lclvartbl char Local parameter table *
* *cmsg struct Gribmsg GRIB2 message structure *
* *
* Output parameters: *
* *param char 12 character parameter name. *
* *scal int scale factor associated with *
* current parameter *
* *msng float missing value associated with *
* current parameter *
* *iret int return code *
* 1 = No gempak param name *
* defined for this grid *
** *
* Log: *
* S. Gilbert/NCEP 12/04 *
* S. Gilbert/NCEP 10/05 Fix null character location *
* S. Gilbert/NCEP 10/05 Use new routines to read tables *
***********************************************************************/
{
int ier, disc, cat, id, pdtn, iver, lclver, len;
const char* filename;
G2Vinfo g2var;
G2vars_t *g2vartbl;
/*---------------------------------------------------------------------*/
*iret = 0;
strncpy( param, "UNKNOWN", 12);
/*
* Get Parameter information from Parameter table(s).
*/
iver=cmsg->gfld->idsect[2];
lclver=cmsg->gfld->idsect[3];
disc=cmsg->gfld->discipline;
cat=cmsg->gfld->ipdtmpl[0];
id=cmsg->gfld->ipdtmpl[1];
pdtn=cmsg->gfld->ipdtnum;
/*
* According to the GRIB2 documentation,
* <http://www.wmo.int/pages/prog/www/WMOCodes/Guides/GRIB/GRIB2_062006.pdf>,
* all the following conditions are true for a GRIB2 message that uses the
* GRIB Master Table, which is maintained by the WMO Secretariat:
* - Master Table version number isn't missing (255)
* - Discipline number isn't reserved for local use (192-254)
* - Category number isn't reserved for local use (192-254)
* - Parameter number isn't reserved for local use (192-254)
* - Product definition template number isn't reserved for local use
* (32768-65534)
* - Local Table version number is zero
*
* We've seen many GRIB2 messages from NCEP that violate these conditions.
*
* Steve Emmerson 2018-07-26
*/
if ((iver != 255) &&
((disc < 192 || disc == 255 ) &&
(cat < 192 || cat == 255 ) &&
(id < 192 || id == 255 ) &&
(pdtn < 32768 || pdtn == 65535) &&
(lclver == 0 || lclver == 255))) { // Allow missing value
/*
* Get WMO Parameter table.
*/
gb2_gtvartbl(wmovartbl, "wmo", iver, &g2vartbl, &filename, &ier);
}
else {
/*
* Get local Parameter table.
*/
gb2_gtvartbl(lclvartbl, cmsg->origcntr, lclver, &g2vartbl, &filename,
&ier);
}
if (ier != 0) {
log_add("Couldn't get parameter table: iver=%d, disc=%d, cat=%d, "
"id=%d, pdtn=%d, center=%.*s, lclver=%d", iver, disc, cat, id,
pdtn, (int)sizeof(cmsg->origcntr), cmsg->origcntr, lclver);
*iret = 1;
return;
}
/*
* Get parameter information from table.
*/
gb2_skvar(disc, cat, id, pdtn, g2vartbl, &g2var, &ier);
if (ier == -1) {
log_warning("Couldn't get parameter info: iver=%d, disc=%d, cat=%d, "
"id=%d, pdtn=%d, center=%.*s, lclver=%d, file=%s",
iver, disc, cat, id, pdtn, (int)sizeof(cmsg->origcntr),
cmsg->origcntr, lclver, filename);
*iret = 1;
return;
}
if (ier) {
log_warning("Using parameter with different PDTN: "
"iver=%d, disc=%d, cat=%d, id=%d, desired pdtn=%d, "
"used pdtn=%d, center=%.*s, lclver=%d, file=%s",
iver, disc, cat, id, pdtn, g2var.pdtnmbr,
(int)sizeof(cmsg->origcntr), cmsg->origcntr, lclver, filename);
}
/*
* Insert time range period in parameter abbreviation,
* if necessary.
*/
gb2_ctim ( cmsg->tmrange, g2var.gemname );
/*
* Adjust ensemble information in parameter abbreviation,
* if necessary.
* NOT DESIRED AT THIS TIME
gb2_ens ( cmsg->gfld, g2var.gemname );
*/
/*
* Adjust probability information in parameter abbreviation,
* if necessary.
*/
gb2_prob ( cmsg->gfld, g2var.gemname );
/*
* Add generating process information in parameter abbreviation
* if necessary.
*/
gb2_proc ( cmsg->gfld, g2var.gemname );
len = strlen(g2var.gemname);
strncpy( param, g2var.gemname, 12);
if ( len > 12 ) param[12]='\0';
if ( len < 12 ) { /* pad gempak parameter name with blanks */
char blanks[13]=" ";
strncat( param, blanks, 12-len );
}
*scal = g2var.scale;
*msng = g2var.missing;
}
void gb2_param ( char *wmovartbl, char *lclvartbl, Gribmsg *cmsg,
char *param, int *scal, float *msng, int *iret )
/************************************************************************
* gb2_param *
* *
* This routine gets the parameter values *
* from the GRIB2 PDS, and obtains the GEMPAK parameter information *
* from the appropriate GRIB2 parameter table. *
* *
* If either wmovartbl or lclvartbl are NULL, the default tables are *
* read. *
* *
* gb2_param ( wmovartbl, lclvartbl, cmsg, param, scal, msng, iret ) *
* *
* Input parameters: *
* *wmovartbl char WMO parameter table *
* *lclvartbl char Local parameter table *
* *cmsg struct Gribmsg GRIB2 message structure *
* *
* Output parameters: *
* *param char 12 character parameter name. *
* *scal int scale factor associated with *
* current parameter *
* *msng float missing value associated with *
* current parameter *
* *iret int return code *
* 1 = No gempak param name *
* defined for this grid *
** *
* Log: *
* S. Gilbert/NCEP 12/04 *
* S. Gilbert/NCEP 10/05 Fix null character location *
* S. Gilbert/NCEP 10/05 Use new routines to read tables *
***********************************************************************/
{
int ier, disc, cat, id, pdtn, iver, lclver, len;
const char* filename;
G2Vinfo g2var;
G2vars_t *g2vartbl;
/*---------------------------------------------------------------------*/
*iret = 0;
strncpy( param, "UNKNOWN", 12);
/*
* Get Parameter information from Parameter table(s).
*/
iver=cmsg->gfld->idsect[2];
lclver=cmsg->gfld->idsect[3];
disc=cmsg->gfld->discipline;
cat=cmsg->gfld->ipdtmpl[0];
id=cmsg->gfld->ipdtmpl[1];
pdtn=cmsg->gfld->ipdtnum;
#if 0
/*
* Some GRIB2 messages from NCEP *don't* have the Master Table Version
* number set to 255 (which is wrong). Consequently, the following hack sets
* things right.
*/
if (iver != 255 && strcmp(cmsg->origcntr, "wmo")) {
unotice("Setting Master Table Version to 255 for non-WMO originating center: "
"iver=%d, disc=%d, cat=%d, id=%d, pdtn=%d, "
"center=%.*s, lclver=%d",
iver, disc, cat, id, pdtn,
(int)sizeof(cmsg->origcntr), cmsg->origcntr, lclver);
cmsg->gfld->idsect[2] = iver = 255;
}
#endif
if ((iver != 255) &&
((disc < 192 || disc == 255 ) &&
(cat < 192 || cat == 255 ) &&
(id < 192 || id == 255 ) &&
(pdtn < 32768 || pdtn == 65535))) {
/*
* Get WMO Parameter table.
*/
gb2_gtvartbl(wmovartbl, "wmo", iver, &g2vartbl, &filename, &ier);
}
else {
/*
* Get Local Parameter table.
*/
gb2_gtvartbl(lclvartbl, cmsg->origcntr, lclver, &g2vartbl, &filename,
&ier);
}
if (ier != 0) {
unotice("Couldn't get parameter table: iver=%d, disc=%d, cat=%d, "
"id=%d, pdtn=%d, center=%.*s, lclver=%d", iver, disc, cat, id,
pdtn, (int)sizeof(cmsg->origcntr), cmsg->origcntr, lclver);
*iret = 1;
return;
}
/*
* Get parameter information from table.
*/
gb2_skvar(disc, cat, id, pdtn, g2vartbl, &g2var, &ier);
if (ier != 0) {
unotice("Couldn't get parameter info: iver=%d, disc=%d, cat=%d, id=%d, "
"pdtn=%d, center=%.*s, lclver=%d, file=%s", iver, disc, cat, id,
pdtn, (int)sizeof(cmsg->origcntr), cmsg->origcntr, lclver,
filename);
*iret = 1;
return;
}
/*
* Insert time range period in parameter abbreviation,
* if necessary.
*/
gb2_ctim ( cmsg->tmrange, g2var.gemname );
/*
* Adjust ensemble information in parameter abbreviation,
* if necessary.
* NOT DESIRED AT THIS TIME
gb2_ens ( cmsg->gfld, g2var.gemname );
*/
/*
* Adjust probability information in parameter abbreviation,
* if necessary.
*/
gb2_prob ( cmsg->gfld, g2var.gemname );
/*
* Add generating process information in parameter abbreviation
* if necessary.
*/
gb2_proc ( cmsg->gfld, g2var.gemname );
len = strlen(g2var.gemname);
strncpy( param, g2var.gemname, 12);
if ( len > 12 ) param[12]='\0';
if ( len < 12 ) { /* pad gempak parameter name with blanks */
char blanks[13]=" ";
strncat( param, blanks, 12-len );
}
*scal = g2var.scale;
*msng = g2var.missing;
}