g2_unpack6.c

Go to the documentation of this file.

 
#include <stdio.h>
#include <stdlib.h>
#include "grib2.h"
 
g2int
g2_unpack6(unsigned char *cgrib, g2int *iofst, g2int ngpts, g2int *ibmap,
           g2int **bmap)
{
    g2int j, ierr = 0, isecnum;
    g2int *lbmap = 0;
    g2int *intbmap;
 
    *bmap = NULL;
 
    *iofst = *iofst + 32;    /* skip Length of Section */
    gbit(cgrib, &isecnum, *iofst, 8);         /* Get Section Number */
    *iofst = *iofst + 8;
 
    if (isecnum != 6)
    {
        ierr = 2;
        fprintf(stderr, "g2_unpack6: Not Section 6 data.\n");
        return(ierr);
    }
 
    gbit(cgrib, ibmap, *iofst, 8);    /* Get bit-map indicator */
    *iofst = *iofst + 8;
 
    if (*ibmap == 0)
    {               /* Unpack bitmap */
        if (ngpts > 0)
            lbmap = calloc(ngpts, sizeof(g2int));
        if (lbmap == 0)
        {
            ierr = 6;
            return(ierr);
        }
        else
        {
            *bmap = lbmap;
        }
        intbmap = calloc(ngpts, sizeof(g2int));
        gbits(cgrib, intbmap, *iofst, 1, 0, ngpts);
        *iofst = *iofst + ngpts;
        for (j = 0; j < ngpts; j++)
        {
            lbmap[j] = (g2int)intbmap[j];
        }
        free(intbmap);
        /*
          else if (*ibmap.eq.254)               ! Use previous bitmap
          return(ierr);
          else if (*ibmap.eq.255)               ! No bitmap in message
          bmap(1:ngpts)=.true.
          else {
          print *,'gf_unpack6: Predefined bitmap ',*ibmap,' not recognized.'
          ierr=4;
        */
    }
 
    return(ierr);    /* End of Section 6 processing */
 
}