gridtemplates.c

Go to the documentation of this file.

 
#include <stdlib.h>
#include "grib2.h"
#include "gridtemplates.h"
 
g2int
getgridindex(g2int number)
{
    g2int j, getgridindex = -1;
 
    for (j = 0; j < MAXGRIDTEMP; j++)
    {
        if (number == templatesgrid[j].template_num)
        {
            getgridindex = j;
            return(getgridindex);
        }
    }
 
    return(getgridindex);
}
 
gtemplate *
getgridtemplate(g2int number)
{
    g2int index;
    gtemplate *new;
 
    index = getgridindex(number);
 
    if (index != -1)
    {
        new = malloc(sizeof(gtemplate));
        new->type = 3;
        new->num = templatesgrid[index].template_num;
        new->maplen = templatesgrid[index].mapgridlen;
        new->needext = templatesgrid[index].needext;
        new->map = (g2int *)templatesgrid[index].mapgrid;
        new->extlen = 0;
        new->ext = NULL;
        return(new);
    }
    else
    {
        printf("getgridtemplate: GDT Template 3.%d not defined.\n", (int)number);
        return(NULL);
    }
 
    return(NULL);
}
 
gtemplate *
extgridtemplate(g2int number, g2int *list)
{
    gtemplate *new;
    g2int index, i;
 
    index = getgridindex(number);
    if (index == -1)
        return(0);
 
    new = getgridtemplate(number);
 
    if (!new->needext)
        return(new);
 
    if (number == 120)
    {
        new->extlen = list[1] * 2;
        new->ext = malloc(sizeof(g2int) * new->extlen);
        for (i = 0; i < new->extlen; i++)
        {
            if (i % 2 == 0)
                new->ext[i] = 2;
            else
                new->ext[i] = -2;
        }
    }
    else if (number == 4)
    {
        new->extlen = list[7];
        new->ext = malloc(sizeof(g2int) * new->extlen);
        for (i = 0; i < new->extlen; i++)
        {
            new->ext[i] = 4;
        }
        new->extlen = list[8];
        new->ext = malloc(sizeof(g2int) * new->extlen);
        for (i = 0; i < new->extlen; i++)
        {
            new->ext[i] = -4;
        }
    }
    else if (number == 5)
    {
        new->extlen = list[7];
        new->ext = malloc(sizeof(g2int) * new->extlen);
        for (i = 0; i < new->extlen; i++)
        {
            new->ext[i] = 4;
        }
        new->extlen = list[8];
        new->ext = malloc(sizeof(g2int) * new->extlen);
        for (i = 0; i < new->extlen; i++)
        {
            new->ext[i] = -4;
        }
    }
    else if (number == 1000)
    {
        new->extlen = list[19];
        new->ext = malloc(sizeof(g2int) * new->extlen);
        for (i = 0; i < new->extlen; i++)
        {
            new->ext[i] = 4;
        }
    }
    else if (number == 1200)
    {
        new->extlen = list[15];
        new->ext = malloc(sizeof(g2int) * new->extlen);
        for (i = 0; i < new->extlen; i++)
        {
            new->ext[i] = 4;
        }
    }
 
    return(new);
}