gridtemplates.c

Go to the documentation of this file.

 
#include <stdlib.h>
#include "grib2_int.h"
#define MAXGRIDTEMP 31 
#define MAXGRIDMAPLEN 200 
struct gridtemplate
{
    g2int template_num; 
    g2int mapgridlen; 
    g2int needext; 
    g2int mapgrid[MAXGRIDMAPLEN]; 
};
 
static const struct gridtemplate templatesgrid[MAXGRIDTEMP] =
{
    /* 3.0: Lat/Lon grid */
    { 0, 19, 0, {1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1} },
    /* 3.1: Rotated Lat/Lon grid */
    { 1, 22, 0, {1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1,-4,4,4} },
    /* 3.2: Stretched Lat/Lon grid */
    { 2, 22, 0, {1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1,-4,4,-4} },
    /* 3.3: Stretched & Rotated Lat/Lon grid */
    { 3, 25, 0, {1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1,-4,4,4,-4,4,-4} },
    /* Added GDT 3.4,3.5    (08/05/2013) */
    /* 3.4: Variable resolution Latitude/Longitude */
    { 4, 13, 1, {1,1,4,1,4,1,4,4,4,4,4,1,1} },
    /* 3.5: Variable resolution rotate Latitude/Longitude */
    { 5, 16, 1, {1,1,4,1,4,1,4,4,4,4,4,1,1,-4,4,4} },
    /* 3.12: Transverse Mercator */
    {12, 22, 0, {1,1,4,1,4,1,4,4,4,-4,4,1,-4,4,4,1,4,4,-4,-4,-4,-4} },
    /* 3.101: General unstructured grid */
    {101, 4, 0, {1,4,1,-4} },
    /* 3.140: Lambert Azimuthal Equal Area Projection */
    {140, 17, 0, {1,1,4,1,4,1,4,4,4,-4,4,4,4,1,4,4,1} },
 
    /* 3.10: Mercator */
    {10, 19, 0, {1,1,4,1,4,1,4,4,4,-4,4,1,-4,-4,4,1,4,4,4} },
    /* 3.20: Polar Stereographic Projection */
    {20, 18, 0, {1,1,4,1,4,1,4,4,4,-4,4,1,-4,4,4,4,1,1} },
    /* 3.30: Lambert Conformal */
    {30, 22, 0, {1,1,4,1,4,1,4,4,4,-4,4,1,-4,4,4,4,1,1,-4,-4,-4,4} },
    /* 3.31: Albers equal area */
    {31, 22, 0, {1,1,4,1,4,1,4,4,4,-4,4,1,-4,4,4,4,1,1,-4,-4,-4,4} },
    /* 3.40: Guassian Lat/Lon */
    {40, 19, 0, {1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1} },
    /* 3.41: Rotated Gaussian Lat/Lon */
    {41, 22, 0, {1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1,-4,4,4} },
    /* 3.42: Stretched Gaussian Lat/Lon */
    {42, 22, 0, {1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1,-4,4,-4} },
    /* 3.43: Stretched and Rotated Gaussian Lat/Lon */
    {43, 25, 0, {1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1,-4,4,4,-4,4,-4} },
    /* 3.50: Spherical Harmonic Coefficients */
    {50, 5, 0, {4,4,4,1,1} },
    /* 3.51: Rotated Spherical Harmonic Coefficients */
    {51, 8, 0, {4,4,4,1,1,-4,4,4} },
    /* 3.52: Stretched Spherical Harmonic Coefficients */
    {52, 8, 0, {4,4,4,1,1,-4,4,-4} },
    /* 3.53: Stretched and Rotated Spherical Harmonic Coefficients */
    {53, 11, 0, {4,4,4,1,1,-4,4,4,-4,4,-4} },
    /* 3.90: Space View Perspective or orthographic */
    {90, 21, 0, {1,1,4,1,4,1,4,4,4,-4,4,1,4,4,4,4,1,4,4,4,4} },
    /* 3.100: Triangular grid based on an icosahedron */
    {100, 11, 0, {1,1,2,1,-4,4,4,1,1,1,4} },
    /* 3.110: Equatorial Azimuthal equidistant */
    {110, 16, 0, {1,1,4,1,4,1,4,4,4,-4,4,1,4,4,1,1} },
    /* 3.120: Azimuth-range projection */
    {120, 7, 1, {4,4,-4,4,4,4,1} },
    /* 3.204: Curvilinear Orthogonal Grid */
    {204, 19, 0, {1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1} },
    /* 3.32768: Rot Lat/Lon E-grid (Arakawa) */
    {32768, 19, 0, {1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1} },
    /* 3.32769: Rot Lat/Lon Non-E Staggered grid (Arakawa) */
    {32769, 21, 0, {1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1,4,4} },
    /* 3.1000: Cross Section Grid */
    {1000, 20, 1, {1,1,4,1,4,1,4,4,4,4,-4,4,1,4,4,1,2,1,1,2} },
    /* 3.1100: Hovmoller Diagram Grid */
    {1100, 28, 0, {1,1,4,1,4,1,4,4,4,4,-4,4,1,-4,4,1,4,1,-4,1,1,-4,2,1,1,1,1,1} },
    /* 3.1200: Time Section Grid */
    {1200, 16, 1, {4,1,-4,1,1,-4,2,1,1,1,1,1,2,1,1,2} }
} ;
 
static 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);
}