示例#1
0
文件: gb2param.c 项目: Unidata/LDM
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;

}
示例#2
0
文件: gb2param.c 项目: khallock/LDM
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;

}