getpoly.c

Go to the documentation of this file.

 
#include <stdio.h>
#include <stdlib.h>
#include "grib2.h"
 
g2int
getpoly(unsigned char *csec3, g2int *jj, g2int *kk, g2int *mm)
{
 
    g2int *igdstmpl, *list_opt;
    g2int *igds;
    g2int iofst, igdtlen, num_opt;
 
    iofst = 0;       /* set offset to beginning of section */
    if (!g2_unpack3(csec3, &iofst, &igds, &igdstmpl, &igdtlen,
            &list_opt, &num_opt))
    {
        switch (igds[4])     /*  Template number */
        {
        case 50:     /* Spherical harmonic coefficients */
        case 51:
        case 52:
        case 53:
        {
            *jj = igdstmpl[0];
            *kk = igdstmpl[1];
            *mm = igdstmpl[2];
            break;
        }
        default:
        {
            *jj = 0;
            *kk = 0;
            *mm = 0;
            break;
        }
        }     /* end switch */
    }
    else
    {
        *jj = 0;
        *kk = 0;
        *mm = 0;
    }
 
    if (igds)
    free(igds);
    if (igdstmpl)
    free(igdstmpl);
    if (list_opt)
    free(list_opt);
 
    return 0;
}