g2cprod.c

Go to the documentation of this file.

 
#include "grib2_int.h"
 
extern G2C_FILE_INFO_T g2c_file[G2C_MAX_FILES + 1];
 
int
g2c_get_prod(int g2cid, int msg_num, int prod_num, int *num_data_points, float *data)
{
    G2C_MESSAGE_INFO_T *msg;
    G2C_SECTION_INFO_T *sec3, *sec4, *sec5, *sec7;
    G2C_SECTION3_INFO_T *sec3_info;
    /* G2C_SECTION4_INFO_T *sec4_info; */
    G2C_SECTION5_INFO_T *sec5_info;
    unsigned char *buf;
    size_t bytes_read;
    int ret = G2C_NOERROR;
        
    /* Check inputs. */
    if (g2cid < 0 || g2cid > G2C_MAX_FILES)
    return G2C_EBADID;
    if (msg_num < 0 || prod_num < 0)
    return G2C_EINVAL;
    if (g2c_file[g2cid].g2cid != g2cid)
    return G2C_EBADID;
 
    /* Find the message. */
    for (msg = g2c_file[g2cid].msg; msg; msg = msg->next)
    if (msg->msg_num == msg_num)
        break;
    if (!msg)
    return G2C_ENOMSG;
 
    /* Find the product. After this, sec4 will point to the
     * appropropriate section 4 G2C_SECTION_INFO_T. */
    for (sec4 = msg->sec; sec4; sec4 = sec4->next)
    if (sec4->sec_num == 4 && ((G2C_SECTION4_INFO_T *)sec4->sec_info)->field_num == prod_num)
        break;
    if (!sec4)
    return G2C_ENOPRODUCT;
    /* sec4_info = (G2C_SECTION4_INFO_T *)sec4->sec_info; */
 
    /* Find the grid definiton section, section 3. It will come
     * earlier in the list. */
    for (sec3 = sec4; sec3; sec3 = sec3->prev)
    if (sec3->sec_num == 3)
        break;
    if (!sec3)
    return G2C_ENOSECTION;
    sec3_info = (G2C_SECTION3_INFO_T *)sec3->sec_info;
 
    /* Find the section 5, data representation section, to learn how
     * this product is compressed. Section 5 is after section 4 in the
     * list. */
    for (sec5 = sec4; sec5; sec5 = sec5->next)
        if (sec5->sec_num == 5)
            break;
    if (!sec5)
    return G2C_ENOSECTION;
    sec5_info = (G2C_SECTION5_INFO_T *)sec5->sec_info;
 
    /* Find the section 7, data section. */
    for (sec7 = sec5; sec7; sec7 = sec7->next)
        if (sec7->sec_num == 7)
            break;
    if (!sec7)
    return G2C_ENOSECTION;
    
    /* Give the caller number of data points, if desired. */
    if (num_data_points)
        *num_data_points = sec5_info->num_data_points;
 
    /* If user doesn't want data, we need go no further. */
    if (!data)
        return G2C_NOERROR;
 
    /* Allocate a char buffer to hold the packed data. */
    if (!(buf = malloc(sizeof(char) * sec7->sec_len)))
        return G2C_ENOMEM;
    
    /* Jump to this section in the file. */
    if (fseek(g2c_file[g2cid].f, sec7->bytes_to_sec + sec7->msg->bytes_to_msg, SEEK_SET))
    return G2C_ERROR;
 
    /* Read the product into a char buffer. */
    if ((bytes_read = fread(buf, 1, sec7->sec_len, g2c_file[g2cid].f)) != sec7->sec_len)
    return G2C_EFILE;
 
    /* Unpack the char buffer into a float array, which must be
     * allocated by the caller. */
    ret = g2c_unpack7(buf, sec3_info->grid_def, sec3->template_len, sec3->template,
              sec5_info->data_def, sec5->template_len, sec5->template,
              sec5_info->num_data_points, data);
 
    /* Free the char buffer. */
    free(buf);
 
    return ret;
}