#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <stddef.h>
#include <math.h>
#include <float.h>
#define UNDEFINED 9.999e20
#define UNDEFINED_LOW 9.9989e20
#define UNDEFINED_HIGH 9.9991e20
#define UNDEF_ERR 0.00001
#define UNDEFINED_VAL(x) ((x) >= UNDEFINED_LOW && (x) <= UNDEFINED_HIGH)
unsigned char *rd_GDS(unsigned char *pds, char *filename, int grid_type);
unsigned int size_24section(unsigned char *section);
int wrt_24section(unsigned char *section, char *filename);
int append_24section(unsigned char *section, FILE *output);
unsigned char *rd_24section(char *filename);
unsigned char *cpGRIBsec(unsigned char *section);
void wrt_grib_msg(FILE *output, unsigned char *pds, unsigned char *gds, unsigned
char *bms, unsigned char *bds);
void set_Center(unsigned char *pds, int center);
int get_Center(unsigned char *pds);
void set_DecScale(unsigned char *pds, int dec_scale);
int get_DecScale(unsigned char *pds);
void set_PDSGridType(unsigned char *pds, int n);
int get_PDSGridType(unsigned char *pds);
void set_HasGDS(unsigned char *pds);
void clr_HasGDS(unsigned char *pds);
int get_HasGDS(unsigned char *pds);
void set_HasBMS(unsigned char *pds);
void clr_HasBMS(unsigned char *pds);
int get_HasBMS(unsigned char *pds);
void set_InitTime(unsigned char *pds, unsigned int time);
unsigned int get_InitTime(unsigned char *pds);
void set_InitDate(unsigned char *pds, int year, int month, int day,
int hour, int minute);
void get_InitDate(unsigned char *pds, int *year, int *month, int *day,
int *hour, int *minute);
unsigned int get_NextYear(unsigned int time);
unsigned int get_NextMonth(unsigned int time);
unsigned int get_NextDay(unsigned int time);
#define get_Hour(time) ((time) % 100)
#define get_Day(time) (((time)/100) % 100)
#define get_Month(time) (((time)/10000) % 100)
#define get_Year(time) ((time)/1000000)
void set_NextYear(unsigned char *pds);
void set_NextMonth(unsigned char *pds);
void set_NextDay(unsigned char *pds);
void set_BDSMaxBits(int n);
void set_BDSMinBits(int n);
int get_BDSMaxBits();
int get_BDSMinBits();
void list2bitstream(int *list, unsigned char *bitstream, int ndata, int nbits);
void flist2bitstream(float *list, unsigned char *bitstream, int ndata, int nbits);
void set_NMCparm(unsigned char *pds, char *name);
char *get_NMCparm(unsigned char *pds);
void set_NMCnparm(unsigned char *pds, int parm);
int get_NMCnparm(unsigned char *pds);
unsigned char *seek_grib(FILE *file, long *pos, long *len_grib,
unsigned char *buffer, unsigned int buf_len);
int read_grib(FILE *file, long pos, long len_grib, unsigned char *buffer);
double int_power(double x, int y);
int flt2ibm(float x, unsigned char *ibm);
double ibm2flt(unsigned char *ibm);
/* gio.c */
float *rd_f77data(int *n, FILE *in);
void rd_n_f77data(float *a, int n, FILE *in);
void wrt_n_f77data(float *a, int n, FILE *out);
unsigned char *mk_void_BMS(unsigned char *pds);
unsigned char *mk_BMS(unsigned char *pds, float *bindata, int *n, float undef_low, float undef_hi);
unsigned char *mk_BDS(unsigned char *pds, float *bindata, int n);
int get_nxny(unsigned char *pds, unsigned char *gds, unsigned char *bms, unsigned char *bds);
int rd_grib_msg(FILE *input, long int pos, unsigned char **pds,
unsigned char **gds, unsigned char **bms, unsigned char **bds);
float *get_unpk_bds(unsigned char *pds, unsigned char *gds, unsigned char *bms,
unsigned char *bds);
int unpk_bds(float *array, unsigned char *pds, unsigned char *gds,
unsigned char *bms, unsigned char *bds, int array_size);
int *get_bitmap(unsigned char *pds, unsigned char *gds, unsigned char *bms,
unsigned char *bds);
void unpk_bms(int *bitmap, unsigned char *bms, int nxny);
unsigned char *gds_grid2();
unsigned char *NCEP_GDS(unsigned char *pds, int grid_type);
/* private */
void set_int3(unsigned char *string, int n);
void set_int2(unsigned char *string, int n);
void set_PDSlevel(unsigned char *pds, int type, int value);
void set_TimeRange(unsigned char *pds, int time_range, int p1, int p2, \
int units, int nave, int nmissing);
void get_TimeRange(unsigned char *pds, int *time_range, int *p1, int *p2, \
int *units, int *nave, int *nmissing);
#define _LEN24(pds) ((int) ((pds[0]<<16)+(pds[1]<<8)+pds[2]))
#define __LEN24(pds) ((pds) == NULL ? 0 : (int) ((pds[0]<<16)+(pds[1]<<8)+pds[2]))
/* version 3.4 of grib headers w. ebisuzaki */
/* this version is incomplete */
#ifndef INT2
#define INT2(a,b) ((1-(int) ((unsigned) (a & 0x80) >> 6)) * (int) (((a & 0x7f) << 8) + b))
#endif
#define PDS_Len1(pds) (pds[0])
#define PDS_Len2(pds) (pds[1])
#define PDS_Len3(pds) (pds[2])
#define PDS_LEN(pds) ((int) ((pds[0]<<16)+(pds[1]<<8)+pds[2]))
#define PDS_Vsn(pds) (pds[3])
#define PDS_Center(pds) (pds[4])
#define PDS_Model(pds) (pds[5])
#define PDS_Grid(pds) (pds[6])
#define PDS_HAS_GDS(pds) ((pds[7] & 128) != 0)
#define PDS_HAS_BMS(pds) ((pds[7] & 64) != 0)
#define PDS_PARAM(pds) (pds[8])
#define PDS_L_TYPE(pds) (pds[9])
#define PDS_LEVEL1(pds) (pds[10])
#define PDS_LEVEL2(pds) (pds[11])
#define PDS_KPDS5(pds) (pds[8])
#define PDS_KPDS6(pds) (pds[9])
#define PDS_KPDS7(pds) ((int) ((pds[10]<<8) + pds[11]))
/* this requires a 32-bit default integer machine */
#define PDS_Field(pds) ((pds[8]<<24)+(pds[9]<<16)+(pds[10]<<8)+pds[11])
#define PDS_Year(pds) (pds[12])
#define PDS_Month(pds) (pds[13])
#define PDS_Day(pds) (pds[14])
#define PDS_Hour(pds) (pds[15])
#define PDS_Minute(pds) (pds[16])
#define PDS_ForecastTimeUnit(pds) (pds[17])
#define PDS_P1(pds) (pds[18])
#define PDS_P2(pds) (pds[19])
#define PDS_TimeRange(pds) (pds[20])
#define PDS_NumAve(pds) ((int) ((pds[21]<<8)+pds[22]))
#define PDS_NumMissing(pds) (pds[23])
#define PDS_Century(pds) (pds[24])
#define PDS_Subcenter(pds) (pds[25])
#define PDS_DecimalScale(pds) INT2(pds[26],pds[27])
#define PDS_Year4(pds) (pds[12] + 100*(pds[24] - (pds[12] != 0)))
/* various centers */
#define NMC 7
#define ECMWF 98
/* ECMWF Extensions */
#define PDS_EcLocalId(pds) (PDS_LEN(pds) >= 41 ? (pds[40]) : 0)
#define PDS_EcClass(pds) (PDS_LEN(pds) >= 42 ? (pds[41]) : 0)
#define PDS_EcType(pds) (PDS_LEN(pds) >= 43 ? (pds[42]) : 0)
#define PDS_EcStream(pds) (PDS_LEN(pds) >= 45 ? (INT2(pds[43], pds[44])) : 0)
/* v0.1 w. ebisuzaki */
/* need to put in all values */
/* used by TimeRange */
#define MINUTE 0
#define HOUR 1
#define DAY 2
#define MONTH 3
#define YEAR 4
#define DECADE 5
#define NORMAL 6
#define CENTURY 7
#define SECOND 254
#define set_lev_sfc(pds) set_PDSlevel(pds,1,0)
#define set_lev_toa(pds) set_PDSlevel(pds,8,0)
#define set_lev_mb(pds,lev) set_PDSlevel(pds,100,lev)
#define set_lev_hp(pds,lev) set_PDSlevel(pds,100,lev)
#define set_lev_mb2mb(pds,top,bottom) set_PDSlevel(pds,101,((int) top/10)*256+ (int) bottom/10)
#define set_lev_hp2hp(pds,top,bottom) set_PDSlevel(pds,101,((int) top/10)*256+ (int) bottom/10)
#define set_lev_msl set_PDSlevel(pds,102,0)
#define set_lev_m_MSL(pds,meters) set_PDSlevel(pds,103,meters)
#define set_lev_m2m_MSL(pds,top,bot) set_PDSlevel(pds,104,((int) top/10)*256+(int) bot/10)
#define set_lev_m_gnd(pds,meters) set_PDSlevel(pds,105,meters)
#define set_lev_m2m_gnd(pds,top,bot) set_PDSlevel(pds,106,((int) top/10)*256+(int) bot/10)
#define set_lev_sigma(pds,sigma) set_PDSlevel(pds,107,(int) (sigma*10000))
#define set_lev_cm_ugnd(pds,cm) set_PDSlevel(pds,111,(int) cm)
#define set_lev_cm2cm_ugnd(pds,top,bot) set_PDSlevel(pds,112,(int) top * 256 + (int) bot)
#define set_lev_K(pds,pot_tmp) set_PDSlevel(pds,113, pot_tmp)
#define set_lev_clm(pds) set_PDSlevel(pds,200,0)
#define set_lev_ocean(pds) set_PDSlevel(pds,201,0)
/* tr p1 p2 unit nave nmissing */
#define set_Analysis(pds) set_TimeRange(pds,0,0,0,1,0,0)
#define set_UninitAnalysis(pds) set_TimeRange(pds,0,0,0,1,0,0)
#define set_InitAnalysis(pds) set_TimeRange(pds,0,0,0,1,0,0)
#define set_Average(pds,p1,p2,unit,nave,nmissing) set_TimeRange(pds,3,p1,p2,unit,nave,nmissing)
#define set_Accumulation(pds,p1,p2,unit,nave,nmissing) set_TimeRange(pds,4,p1,p2,unit,nave,nmissing)
#define set_Forecast(pds,time,unit) set_TimeRange(pds,10,time>>8,time&255,unit,0,0)
#define set_Climo(pds,dt,unit,nyear) set_TimeRange(pds,51,0,dt,unit,nyear,0)
#define set_DiurnalClimo(pds,dt,unit,nyear) set_TimeRange(pds,51,1,dt,unit,nyear,0)
#define set_AveForecast(pds,p1,p2,unit,nave,nmissing) set_TimeRange(pds,113,p1,p2,unit,nave,nmissing)
#define set_AveAnalysis(pds,p2,unit,nave,nmissing) set_TimeRange(pds,113,0,p2,unit,nave,nmissing)
#define set_AccumFcst(pds,p1,p2,unit,nave,nmissing) set_TimeRange(pds,113,p1,p2,unit,nave,nmissing)
#define set_AccumAnalysis(pds,p2,unit,nave,nmissing) set_TimeRange(pds,113,0,p2,unit,nave,nmissing)
#define set_Variance(pds,p1,dt,nave,nmissing) set_TimeRange(pds,118,p1,dt,unit,nave,nmissing)
/*
* Takes a gds-less file and adds a gds
* -- only for NCEP files --
*
* Wesley Ebisuzaki
*
*/
int main(int argc, char **argv) {
long int len_grib, pos = 0;
unsigned char *pds, *gds, *bms, *bds, *gds0;
FILE *input, *output;
int nrec, ngds;
/* preliminaries .. open up all files */
if (argc != 3) {
fprintf(stderr, "%s [in gribfile] [out gribfile]\n", argv[0]);
fprintf(stderr, "adds GDS to NCEP files\n");
exit(8);
}
if ((input = fopen(argv[1],"rb")) == NULL) {
fprintf(stderr,"could not open file: %s\n", argv[1]);
exit(7);
}
if ((output = fopen(argv[2],"wb")) == NULL) {
fprintf(stderr,"could not open file: %s\n", argv[2]);
exit(7);
}
nrec = ngds = pos = 0;
for(;;) { /* loop over all records */
len_grib = rd_grib_msg(input, pos, &pds, &gds, &bms, &bds);
if (len_grib <= 0) break;
nrec++;
if (PDS_HAS_GDS(pds) || PDS_Center(pds) != NMC) {
wrt_grib_msg(output, pds, gds, bms, bds);
}
else {
gds0 = NCEP_GDS(pds, PDS_Grid(pds));
wrt_grib_msg(output, pds, gds0, bms, bds);
if (gds0) free(gds0);
ngds++;
}
pos += len_grib;
}
fclose(output);
fclose(input);
printf("%d GDS sections added out of %d records\n", ngds, nrec);
return 0;
}
/*
* rd_grib_msg.c
*
* Wesley Ebisuzaki
*
* returns pointers to various sections of a grib file
* NOTE: must save contents of pds .... bds before calling
* rd_grib_msg again.
*/
#define BUFF_ALLOC0 40000
#define MSEEK 1024
static unsigned char *buffer = NULL;
static int buffer_size = 0;
int rd_grib_msg(FILE *input, long int pos, unsigned char **pds,
unsigned char **gds, unsigned char **bms, unsigned char **bds) {
unsigned char *msg, *pointer, *ppds;
long int len_grib;
/* setup grib buffer */
if (buffer == NULL) {
if ((buffer = (unsigned char *) malloc(BUFF_ALLOC0)) == NULL) {
fprintf(stderr,"not enough memory: rd_grib_msg\n");
buffer_size = BUFF_ALLOC0;
}
}
/* find start and length of message */
msg = seek_grib(input, &pos, &len_grib, buffer, MSEEK);
if (msg == NULL) {
return -1;
}
/* read all whole grib record */
if (len_grib + msg - buffer > buffer_size) {
buffer_size = len_grib + msg - buffer + 1000;
buffer = (unsigned char *) realloc((void *) buffer, buffer_size);
if (buffer == NULL) {
fprintf(stderr,"ran out of memory\n");
exit(8);
}
}
read_grib(input, pos, len_grib, buffer);
/* parse grib message */
msg = buffer;
*pds = ppds = (msg + 8);
pointer = ppds + PDS_LEN(ppds);
if (PDS_HAS_GDS(ppds)) {
*gds = pointer;
pointer += _LEN24(pointer);
}
else {
*gds = NULL;
}
if (PDS_HAS_BMS(ppds)) {
*bms = pointer;
pointer += _LEN24(pointer);
}
else {
*bms = NULL;
}
*bds = pointer;
pointer += _LEN24(pointer);
/* end section - "7777" in ascii */
if (pointer[0] != 0x37 || pointer[1] != 0x37 ||
pointer[2] != 0x37 || pointer[3] != 0x37) {
fprintf(stderr,"\n\n missing end section\n");
fprintf(stderr, "%2x %2x %2x %2x\n", pointer[0], pointer[1],
pointer[2], pointer[3]);
exit(8);
}
return len_grib;
}
/* v1.0 w. ebisuzaki
*
* takes the sections (pds, gds, bms and bds) and writes a grib message
*/
static unsigned char header[8] = {'G', 'R', 'I', 'B', ' ', ' ', ' ', '\1'},
trailer[4] = {'7', '7', '7', '7'};
void wrt_grib_msg(FILE *output, unsigned char *pds, unsigned char *gds,
unsigned char *bms, unsigned char *bds) {
unsigned int size;
/* header */
size = __LEN24(pds) + __LEN24(gds) + __LEN24(bms) + __LEN24(bds) + 12;
header[4] = (size >> 16) & 255;
header[5] = (size >> 8) & 255;
header[6] = (size ) & 255;
if (fwrite(header, sizeof (unsigned char), 8, output) != 8) {
fprintf(stderr,"error writing output\n");
exit(8);
}
append_24section(pds, output);
append_24section(gds, output);
append_24section(bms, output);
append_24section(bds, output);
if (fwrite(trailer, sizeof (unsigned char), 4, output) != 4) {
fprintf(stderr,"error writing output\n");
exit(8);
}
}
/*
* find next grib header
*
* file = what do you think?
* pos = initial position to start looking at ( = 0 for 1st call)
* returns with position of next grib header (units=bytes)
* len_grib = length of the grib record (bytes)
* buffer[buf_len] = buffer for reading/writing
*
* returns (char *) to start of GRIB header+PDS
* NULL if not found
*
* adapted from SKGB (Mark Iredell)
*
* v1.1 9/94 Wesley Ebisuzaki
* v1.2 3/96 Wesley Ebisuzaki handles short records at end of file
*
*/
#ifndef min
#define min(a,b) ((a) < (b) ? (a) : (b))
#endif
#define NTRY 100
/* #define LEN_HEADER_PDS (28+42+100) */
#define LEN_HEADER_PDS (28+8)
unsigned char *seek_grib(FILE *file, long *pos, long *len_grib,
unsigned char *buffer, unsigned int buf_len) {
int i, j, len;
for (j = 0; j < NTRY; j++) {
if (fseek(file, *pos, SEEK_SET) == -1) {
*len_grib = 0;
return (unsigned char *) NULL;
}
i = fread(buffer, sizeof (unsigned char), buf_len, file);
len = i - LEN_HEADER_PDS;
for (i = 0; i < len; i++) {
if (buffer[i] == 'G' && buffer[i+1] == 'R' && buffer[i+2] == 'I'
&& buffer[i+3] == 'B' && buffer[i+7] == 1) {
*pos = i + *pos;
*len_grib = (buffer[i+4] << 16) + (buffer[i+5] << 8) +
buffer[i+6];
return (buffer+i);
}
}
*pos = *pos + (buf_len - LEN_HEADER_PDS);
}
*len_grib = 0;
return (unsigned char *) NULL;
}
/*
* read_grib.c
*
* v1.0 9/94 Wesley Ebisuzaki
*
*/
int read_grib(FILE *file, long pos, long len_grib, unsigned char *buffer) {
int i;
if (fseek(file, pos, SEEK_SET) == -1) {
return 0;
}
i = fread(buffer, sizeof (unsigned char), len_grib, file);
return (i == len_grib);
}
static unsigned char gds1[ 42] = {
0, 0, 42, 0, 255, 1, 0, 73, 0, 23,
128, 187, 218, 0, 0, 0, 128, 0, 187, 218,
0, 0, 0, 0, 87, 228, 0, 64, 7, 214,
133, 7, 214, 133, 0, 0, 0, 0, 0, 0,
0, 0,
};
static unsigned char gds2[ 32] = {
0, 0, 32, 0, 255, 0, 0, 144, 0, 73,
1, 95, 144, 0, 0, 0, 128, 129, 95, 144,
128, 9, 196, 9, 196, 9, 196, 0, 0, 0,
0, 0,
};
static unsigned char gds3[ 32] = {
0, 0, 32, 0, 255, 0, 1, 104, 0, 181,
1, 95, 144, 0, 0, 0, 128, 129, 95, 144,
128, 3, 232, 3, 232, 3, 232, 0, 0, 0,
0, 0,
};
static unsigned char gds4[ 32] = {
0, 0, 32, 0, 255, 0, 2, 208, 1, 105,
1, 95, 144, 0, 0, 0, 128, 129, 95, 144,
128, 1, 244, 1, 244, 1, 244, 0, 0, 0,
0, 0,
};
static unsigned char gds5[ 32] = {
0, 0, 32, 0, 255, 5, 0, 53, 0, 57,
0, 29, 223, 130, 9, 67, 8, 129, 154, 40,
2, 232, 36, 2, 232, 36, 0, 64, 0, 0,
0, 0,
};
static unsigned char gds6[ 32] = {
0, 0, 32, 0, 255, 5, 0, 53, 0, 45,
0, 29, 223, 130, 9, 67, 8, 129, 154, 40,
2, 232, 36, 2, 232, 36, 0, 64, 0, 0,
0, 0,
};
static unsigned char gds21[ 32] = {
0, 0, 32, 0, 255, 0, 0, 37, 0, 37,
0, 0, 0, 0, 0, 0, 128, 1, 95, 144,
2, 191, 32, 19, 136, 9, 196, 64, 0, 0,
0, 0,
};
static unsigned char gds22[ 32] = {
0, 0, 32, 0, 255, 0, 0, 37, 0, 37,
0, 0, 0, 130, 191, 32, 128, 1, 95, 144,
0, 0, 0, 19, 136, 9, 196, 64, 0, 0,
0, 0,
};
static unsigned char gds23[ 32] = {
0, 0, 32, 0, 255, 0, 0, 37, 0, 37,
129, 95, 144, 0, 0, 0, 128, 0, 0, 0,
2, 191, 32, 19, 136, 9, 196, 64, 0, 0,
0, 0,
};
static unsigned char gds24[ 32] = {
0, 0, 32, 0, 255, 0, 0, 37, 0, 37,
129, 95, 144, 130, 191, 32, 128, 0, 0, 0,
0, 0, 0, 19, 136, 9, 196, 64, 0, 0,
0, 0,
};
static unsigned char gds25[ 32] = {
0, 0, 32, 0, 255, 0, 0, 72, 0, 19,
0, 0, 0, 0, 0, 0, 128, 1, 95, 144,
5, 106, 184, 19, 136, 19, 136, 64, 0, 0,
0, 0,
};
static unsigned char gds26[ 32] = {
0, 0, 32, 0, 255, 0, 0, 72, 0, 19,
129, 95, 144, 0, 0, 0, 128, 0, 0, 0,
5, 106, 184, 19, 136, 19, 136, 64, 0, 0,
0, 0,
};
static unsigned char gds27[ 32] = {
0, 0, 32, 0, 255, 5, 0, 65, 0, 65,
128, 79, 62, 129, 232, 72, 8, 129, 56, 128,
5, 208, 72, 5, 208, 72, 0, 64, 0, 0,
0, 0,
};
static unsigned char gds28[ 32] = {
0, 0, 32, 0, 255, 5, 0, 65, 0, 65,
0, 79, 62, 2, 54, 104, 8, 1, 134, 160,
5, 208, 72, 5, 208, 72, 128, 64, 0, 0,
0, 0,
};
static unsigned char gds29[ 32] = {
0, 0, 32, 0, 255, 0, 0, 145, 0, 37,
0, 0, 0, 0, 0, 0, 128, 1, 95, 144,
5, 126, 64, 9, 196, 9, 196, 64, 0, 0,
0, 0,
};
static unsigned char gds30[ 32] = {
0, 0, 32, 0, 255, 0, 0, 145, 0, 37,
0, 0, 0, 129, 95, 144, 128, 0, 0, 0,
5, 126, 64, 9, 196, 9, 196, 64, 0, 0,
0, 0,
};
static unsigned char gds33[ 32] = {
0, 0, 32, 0, 255, 0, 0, 181, 0, 46,
0, 0, 0, 0, 0, 0, 128, 1, 95, 144,
5, 126, 64, 7, 208, 7, 208, 64, 0, 0,
0, 0,
};
static unsigned char gds34[ 32] = {
0, 0, 32, 0, 255, 0, 0, 181, 0, 46,
129, 95, 144, 0, 0, 0, 128, 0, 0, 0,
5, 126, 64, 7, 208, 7, 208, 64, 0, 0,
0, 0,
};
static unsigned char gds37[178] = {
0, 0, 178, 0, 33, 0, 255, 255, 0, 73,
0, 0, 0, 128, 117, 48, 128, 1, 95, 144,
0, 234, 96, 255, 255, 4, 226, 64, 0, 0,
0, 0, 0, 73, 0, 73, 0, 73, 0, 73,
0, 73, 0, 73, 0, 73, 0, 73, 0, 72,
0, 72, 0, 72, 0, 71, 0, 71, 0, 71,
0, 70, 0, 70, 0, 69, 0, 69, 0, 68,
0, 67, 0, 67, 0, 66, 0, 65, 0, 65,
0, 64, 0, 63, 0, 62, 0, 61, 0, 60,
0, 60, 0, 59, 0, 58, 0, 57, 0, 56,
0, 55, 0, 54, 0, 52, 0, 51, 0, 50,
0, 49, 0, 48, 0, 47, 0, 45, 0, 44,
0, 43, 0, 42, 0, 40, 0, 39, 0, 38,
0, 36, 0, 35, 0, 33, 0, 32, 0, 30,
0, 29, 0, 28, 0, 26, 0, 25, 0, 23,
0, 22, 0, 20, 0, 19, 0, 17, 0, 16,
0, 14, 0, 12, 0, 11, 0, 9, 0, 8,
0, 6, 0, 5, 0, 3, 0, 2,
};
static unsigned char gds38[178] = {
0, 0, 178, 0, 33, 0, 255, 255, 0, 73,
0, 0, 0, 0, 234, 96, 128, 1, 95, 144,
2, 73, 240, 255, 255, 4, 226, 64, 0, 0,
0, 0, 0, 73, 0, 73, 0, 73, 0, 73,
0, 73, 0, 73, 0, 73, 0, 73, 0, 72,
0, 72, 0, 72, 0, 71, 0, 71, 0, 71,
0, 70, 0, 70, 0, 69, 0, 69, 0, 68,
0, 67, 0, 67, 0, 66, 0, 65, 0, 65,
0, 64, 0, 63, 0, 62, 0, 61, 0, 60,
0, 60, 0, 59, 0, 58, 0, 57, 0, 56,
0, 55, 0, 54, 0, 52, 0, 51, 0, 50,
0, 49, 0, 48, 0, 47, 0, 45, 0, 44,
0, 43, 0, 42, 0, 40, 0, 39, 0, 38,
0, 36, 0, 35, 0, 33, 0, 32, 0, 30,
0, 29, 0, 28, 0, 26, 0, 25, 0, 23,
0, 22, 0, 20, 0, 19, 0, 17, 0, 16,
0, 14, 0, 12, 0, 11, 0, 9, 0, 8,
0, 6, 0, 5, 0, 3, 0, 2,
};
static unsigned char gds39[178] = {
0, 0, 178, 0, 33, 0, 255, 255, 0, 73,
0, 0, 0, 2, 73, 240, 128, 1, 95, 144,
129, 212, 192, 255, 255, 4, 226, 64, 0, 0,
0, 0, 0, 73, 0, 73, 0, 73, 0, 73,
0, 73, 0, 73, 0, 73, 0, 73, 0, 72,
0, 72, 0, 72, 0, 71, 0, 71, 0, 71,
0, 70, 0, 70, 0, 69, 0, 69, 0, 68,
0, 67, 0, 67, 0, 66, 0, 65, 0, 65,
0, 64, 0, 63, 0, 62, 0, 61, 0, 60,
0, 60, 0, 59, 0, 58, 0, 57, 0, 56,
0, 55, 0, 54, 0, 52, 0, 51, 0, 50,
0, 49, 0, 48, 0, 47, 0, 45, 0, 44,
0, 43, 0, 42, 0, 40, 0, 39, 0, 38,
0, 36, 0, 35, 0, 33, 0, 32, 0, 30,
0, 29, 0, 28, 0, 26, 0, 25, 0, 23,
0, 22, 0, 20, 0, 19, 0, 17, 0, 16,
0, 14, 0, 12, 0, 11, 0, 9, 0, 8,
0, 6, 0, 5, 0, 3, 0, 2,
};
static unsigned char gds40[178] = {
0, 0, 178, 0, 33, 0, 255, 255, 0, 73,
0, 0, 0, 129, 212, 192, 128, 1, 95, 144,
128, 117, 48, 255, 255, 4, 226, 64, 0, 0,
0, 0, 0, 73, 0, 73, 0, 73, 0, 73,
0, 73, 0, 73, 0, 73, 0, 73, 0, 72,
0, 72, 0, 72, 0, 71, 0, 71, 0, 71,
0, 70, 0, 70, 0, 69, 0, 69, 0, 68,
0, 67, 0, 67, 0, 66, 0, 65, 0, 65,
0, 64, 0, 63, 0, 62, 0, 61, 0, 60,
0, 60, 0, 59, 0, 58, 0, 57, 0, 56,
0, 55, 0, 54, 0, 52, 0, 51, 0, 50,
0, 49, 0, 48, 0, 47, 0, 45, 0, 44,
0, 43, 0, 42, 0, 40, 0, 39, 0, 38,
0, 36, 0, 35, 0, 33, 0, 32, 0, 30,
0, 29, 0, 28, 0, 26, 0, 25, 0, 23,
0, 22, 0, 20, 0, 19, 0, 17, 0, 16,
0, 14, 0, 12, 0, 11, 0, 9, 0, 8,
0, 6, 0, 5, 0, 3, 0, 2,
};
static unsigned char gds41[178] = {
0, 0, 178, 0, 33, 0, 255, 255, 0, 73,
129, 95, 144, 128, 117, 48, 128, 0, 0, 0,
0, 234, 96, 255, 255, 4, 226, 64, 0, 0,
0, 0, 0, 2, 0, 3, 0, 5, 0, 6,
0, 8, 0, 9, 0, 11, 0, 12, 0, 14,
0, 16, 0, 17, 0, 19, 0, 20, 0, 22,
0, 23, 0, 25, 0, 26, 0, 28, 0, 29,
0, 30, 0, 32, 0, 33, 0, 35, 0, 36,
0, 38, 0, 39, 0, 40, 0, 42, 0, 43,
0, 44, 0, 45, 0, 47, 0, 48, 0, 49,
0, 50, 0, 51, 0, 52, 0, 54, 0, 55,
0, 56, 0, 57, 0, 58, 0, 59, 0, 60,
0, 60, 0, 61, 0, 62, 0, 63, 0, 64,
0, 65, 0, 65, 0, 66, 0, 67, 0, 67,
0, 68, 0, 69, 0, 69, 0, 70, 0, 70,
0, 71, 0, 71, 0, 71, 0, 72, 0, 72,
0, 72, 0, 73, 0, 73, 0, 73, 0, 73,
0, 73, 0, 73, 0, 73, 0, 73,
};
static unsigned char gds42[178] = {
0, 0, 178, 0, 33, 0, 255, 255, 0, 73,
129, 95, 144, 0, 234, 96, 128, 0, 0, 0,
2, 73, 240, 255, 255, 4, 226, 64, 0, 0,
0, 0, 0, 2, 0, 3, 0, 5, 0, 6,
0, 8, 0, 9, 0, 11, 0, 12, 0, 14,
0, 16, 0, 17, 0, 19, 0, 20, 0, 22,
0, 23, 0, 25, 0, 26, 0, 28, 0, 29,
0, 30, 0, 32, 0, 33, 0, 35, 0, 36,
0, 38, 0, 39, 0, 40, 0, 42, 0, 43,
0, 44, 0, 45, 0, 47, 0, 48, 0, 49,
0, 50, 0, 51, 0, 52, 0, 54, 0, 55,
0, 56, 0, 57, 0, 58, 0, 59, 0, 60,
0, 60, 0, 61, 0, 62, 0, 63, 0, 64,
0, 65, 0, 65, 0, 66, 0, 67, 0, 67,
0, 68, 0, 69, 0, 69, 0, 70, 0, 70,
0, 71, 0, 71, 0, 71, 0, 72, 0, 72,
0, 72, 0, 73, 0, 73, 0, 73, 0, 73,
0, 73, 0, 73, 0, 73, 0, 73,
};
static unsigned char gds43[178] = {
0, 0, 178, 0, 33, 0, 255, 255, 0, 73,
129, 95, 144, 2, 73, 240, 128, 0, 0, 0,
129, 212, 192, 255, 255, 4, 226, 64, 0, 0,
0, 0, 0, 2, 0, 3, 0, 5, 0, 6,
0, 8, 0, 9, 0, 11, 0, 12, 0, 14,
0, 16, 0, 17, 0, 19, 0, 20, 0, 22,
0, 23, 0, 25, 0, 26, 0, 28, 0, 29,
0, 30, 0, 32, 0, 33, 0, 35, 0, 36,
0, 38, 0, 39, 0, 40, 0, 42, 0, 43,
0, 44, 0, 45, 0, 47, 0, 48, 0, 49,
0, 50, 0, 51, 0, 52, 0, 54, 0, 55,
0, 56, 0, 57, 0, 58, 0, 59, 0, 60,
0, 60, 0, 61, 0, 62, 0, 63, 0, 64,
0, 65, 0, 65, 0, 66, 0, 67, 0, 67,
0, 68, 0, 69, 0, 69, 0, 70, 0, 70,
0, 71, 0, 71, 0, 71, 0, 72, 0, 72,
0, 72, 0, 73, 0, 73, 0, 73, 0, 73,
0, 73, 0, 73, 0, 73, 0, 73,
};
static unsigned char gds44[178] = {
0, 0, 178, 0, 33, 0, 255, 255, 0, 73,
129, 95, 144, 129, 212, 192, 128, 0, 0, 0,
128, 117, 48, 255, 255, 4, 226, 64, 0, 0,
0, 0, 0, 2, 0, 3, 0, 5, 0, 6,
0, 8, 0, 9, 0, 11, 0, 12, 0, 14,
0, 16, 0, 17, 0, 19, 0, 20, 0, 22,
0, 23, 0, 25, 0, 26, 0, 28, 0, 29,
0, 30, 0, 32, 0, 33, 0, 35, 0, 36,
0, 38, 0, 39, 0, 40, 0, 42, 0, 43,
0, 44, 0, 45, 0, 47, 0, 48, 0, 49,
0, 50, 0, 51, 0, 52, 0, 54, 0, 55,
0, 56, 0, 57, 0, 58, 0, 59, 0, 60,
0, 60, 0, 61, 0, 62, 0, 63, 0, 64,
0, 65, 0, 65, 0, 66, 0, 67, 0, 67,
0, 68, 0, 69, 0, 69, 0, 70, 0, 70,
0, 71, 0, 71, 0, 71, 0, 72, 0, 72,
0, 72, 0, 73, 0, 73, 0, 73, 0, 73,
0, 73, 0, 73, 0, 73, 0, 73,
};
static unsigned char gds45[ 32] = {
0, 0, 32, 0, 255, 0, 1, 32, 0, 145,
1, 95, 144, 0, 0, 0, 128, 129, 95, 144,
128, 4, 226, 4, 226, 4, 226, 0, 0, 0,
0, 0,
};
static unsigned char gds55[ 32] = {
0, 0, 32, 0, 255, 5, 0, 87, 0, 71,
128, 42, 195, 130, 90, 177, 8, 129, 154, 40,
3, 224, 48, 3, 224, 48, 0, 64, 0, 0,
0, 0,
};
static unsigned char gds56[ 32] = {
0, 0, 32, 0, 255, 5, 0, 87, 0, 71,
0, 29, 223, 130, 9, 67, 8, 129, 154, 40,
1, 240, 24, 1, 240, 24, 0, 64, 0, 0,
0, 0,
};
static unsigned char gds61[ 32] = {
0, 0, 32, 0, 255, 0, 0, 91, 0, 46,
0, 0, 0, 0, 0, 0, 128, 1, 95, 144,
2, 191, 32, 7, 208, 7, 208, 64, 0, 0,
0, 0,
};
static unsigned char gds62[ 32] = {
0, 0, 32, 0, 255, 0, 0, 91, 0, 46,
0, 0, 0, 130, 191, 32, 128, 1, 95, 144,
0, 0, 0, 7, 208, 7, 208, 64, 0, 0,
0, 0,
};
static unsigned char gds63[ 32] = {
0, 0, 32, 0, 255, 0, 0, 91, 0, 46,
0, 0, 0, 129, 95, 144, 128, 0, 0, 0,
2, 191, 32, 7, 208, 7, 208, 64, 0, 0,
0, 0,
};
static unsigned char gds64[ 32] = {
0, 0, 32, 0, 255, 0, 0, 91, 0, 46,
129, 95, 144, 130, 191, 32, 128, 0, 0, 0,
0, 0, 0, 7, 208, 7, 208, 64, 0, 0,
0, 0,
};
static unsigned char gds85[ 32] = {
0, 0, 32, 0, 255, 0, 1, 104, 0, 90,
0, 1, 244, 0, 1, 244, 128, 1, 93, 156,
5, 124, 76, 3, 232, 3, 232, 64, 0, 0,
0, 0,
};
static unsigned char gds86[ 32] = {
0, 0, 32, 0, 255, 0, 1, 104, 0, 90,
129, 93, 156, 0, 1, 244, 128, 128, 1, 244,
5, 124, 76, 3, 232, 3, 232, 64, 0, 0,
0, 0,
};
static unsigned char gds87[ 32] = {
0, 0, 32, 0, 255, 5, 0, 81, 0, 62,
0, 89, 92, 129, 214, 171, 8, 129, 154, 40,
1, 10, 57, 1, 10, 57, 0, 64, 0, 0,
0, 0,
};
static unsigned char gds90[ 32] = {
0, 0, 32, 0, 255, 201, 50, 102, 0, 1,
0, 0, 182, 130, 73, 127, 136, 0, 0, 92,
0, 0, 141, 2, 65, 2, 26, 64, 0, 0,
0, 0,
};
static unsigned char gds91[ 32] = {
0, 0, 32, 0, 255, 202, 100, 203, 0, 1,
0, 0, 182, 130, 73, 127, 136, 0, 0, 183,
0, 0, 141, 2, 65, 2, 26, 64, 0, 0,
0, 0,
};
static unsigned char gds92[ 32] = {
0, 0, 32, 0, 255, 201, 94, 98, 0, 1,
0, 37, 206, 129, 247, 58, 136, 0, 0, 127,
0, 0, 191, 1, 22, 1, 7, 64, 0, 0,
0, 0,
};
static unsigned char gds93[ 32] = {
0, 0, 32, 0, 255, 202, 188, 195, 0, 1,
0, 37, 206, 129, 247, 58, 136, 0, 0, 253,
0, 0, 191, 1, 22, 1, 7, 64, 0, 0,
0, 0,
};
static unsigned char gds94[ 32] = {
0, 0, 32, 0, 255, 201, 191, 20, 0, 1,
0, 37, 206, 129, 247, 58, 136, 0, 0, 181,
0, 1, 15, 0, 194, 0, 185, 64, 0, 0,
0, 0,
};
static unsigned char gds95[ 32] = {
0, 0, 32, 0, 255, 202, 126, 39, 0, 1,
0, 37, 206, 129, 247, 58, 136, 0, 1, 105,
0, 1, 15, 0, 194, 0, 185, 64, 0, 0,
0, 0,
};
static unsigned char gds96[ 32] = {
0, 0, 32, 0, 255, 201, 162, 158, 0, 1,
128, 13, 113, 130, 69, 63, 136, 0, 0, 160,
0, 1, 5, 1, 77, 1, 52, 64, 0, 0,
0, 0,
};
static unsigned char gds97[ 32] = {
0, 0, 32, 0, 255, 202, 69, 59, 0, 1,
128, 13, 113, 130, 69, 63, 136, 0, 1, 63,
0, 1, 5, 1, 77, 1, 52, 64, 0, 0,
0, 0,
};
static unsigned char gds98[ 32] = {
0, 0, 32, 0, 255, 4, 0, 192, 0, 94,
1, 89, 222, 0, 0, 0, 128, 129, 89, 222,
128, 3, 170, 7, 83, 0, 47, 0, 0, 0,
0, 0,
};
static unsigned char gds100[ 32] = {
0, 0, 32, 0, 255, 5, 0, 83, 0, 83,
0, 66, 212, 129, 249, 16, 8, 129, 154, 40,
1, 101, 60, 1, 101, 60, 0, 64, 0, 0,
0, 0,
};
static unsigned char gds101[ 32] = {
0, 0, 32, 0, 255, 5, 0, 113, 0, 91,
0, 41, 32, 130, 23, 186, 8, 129, 154, 40,
1, 101, 60, 1, 101, 60, 0, 64, 0, 0,
0, 0,
};
static unsigned char gds103[ 32] = {
0, 0, 32, 0, 255, 5, 0, 65, 0, 56,
0, 87, 133, 129, 218, 8, 8, 129, 154, 40,
1, 101, 60, 1, 101, 60, 0, 64, 0, 0,
0, 0,
};
static unsigned char gds104[ 32] = {
0, 0, 32, 0, 255, 5, 0, 147, 0, 110,
128, 1, 12, 130, 32, 211, 8, 129, 154, 40,
1, 98, 131, 1, 98, 131, 0, 64, 0, 0,
0, 0,
};
static unsigned char gds105[ 32] = {
0, 0, 32, 0, 255, 5, 0, 83, 0, 83,
0, 68, 121, 129, 249, 16, 8, 129, 154, 40,
1, 98, 131, 1, 98, 131, 0, 64, 0, 0,
0, 0,
};
static unsigned char gds106[ 32] = {
0, 0, 32, 0, 255, 5, 0, 165, 0, 117,
0, 68, 125, 129, 249, 16, 8, 129, 154, 40,
0, 177, 61, 0, 177, 61, 0, 64, 0, 0,
0, 0,
};
static unsigned char gds107[ 32] = {
0, 0, 32, 0, 255, 5, 0, 120, 0, 92,
0, 91, 142, 129, 213, 104, 8, 129, 154, 40,
0, 177, 61, 0, 177, 61, 0, 64, 0, 0,
0, 0,
};
static unsigned char gds126[ 32] = {
0, 0, 32, 0, 255, 4, 1, 128, 0, 190,
1, 92, 189, 0, 0, 0, 128, 129, 92, 189,
128, 3, 170, 3, 170, 0, 95, 0, 0, 0,
0, 0,
};
static unsigned char gds192[ 32] = {
0, 0, 32, 0, 255, 203, 0, 237, 1, 131, 128, 13,
113, 130, 69, 63, 136, 0, 195, 80, 129, 177, 152,
0, 225, 0, 207, 64, 0, 0, 0, 0
};
static unsigned char gds201[ 32] = {
0, 0, 32, 0, 255, 5, 0, 65, 0, 65,
128, 81, 90, 130, 73, 240, 8, 129, 154, 40,
5, 208, 72, 5, 208, 72, 0, 64, 0, 0,
0, 0,
};
static unsigned char gds202[ 32] = {
0, 0, 32, 0, 255, 5, 0, 65, 0, 43,
0, 30, 158, 130, 38, 228, 8, 129, 154, 40,
2, 232, 36, 2, 232, 36, 0, 64, 0, 0,
0, 0,
};
static unsigned char gds203[ 32] = {
0, 0, 32, 0, 255, 5, 0, 45, 0, 39,
0, 74, 188, 130, 213, 237, 8, 130, 73, 240,
2, 232, 36, 2, 232, 36, 0, 64, 0, 0,
0, 0,
};
static unsigned char gds204[ 42] = {
0, 0, 42, 0, 255, 1, 0, 93, 0, 68,
128, 97, 168, 1, 173, 176, 128, 0, 236, 228,
129, 170, 73, 0, 78, 32, 0, 64, 2, 113,
0, 2, 113, 0, 0, 0, 0, 0, 0, 0,
0, 0,
};
static unsigned char gds205[ 32] = {
0, 0, 32, 0, 255, 5, 0, 45, 0, 39,
0, 2, 104, 129, 75, 168, 8, 128, 234, 96,
2, 232, 36, 2, 232, 36, 0, 64, 0, 0,
0, 0,
};
static unsigned char gds206[ 42] = {
0, 0, 42, 0, 255, 3, 0, 51, 0, 41,
0, 87, 17, 129, 204, 231, 8, 129, 115, 24,
1, 61, 119, 1, 61, 119, 0, 64, 0, 97,
168, 0, 97, 168, 0, 0, 0, 0, 0, 0,
0, 0,
};
static unsigned char gds207[ 32] = {
0, 0, 32, 0, 255, 5, 0, 49, 0, 35,
0, 164, 101, 130, 174, 25, 8, 130, 73, 240,
1, 116, 18, 1, 116, 18, 0, 64, 0, 0,
0, 0,
};
static unsigned char gds208[ 42] = {
0, 0, 42, 0, 255, 1, 0, 29, 0, 27,
0, 36, 127, 130, 141, 147, 128, 0, 109, 188,
130, 57, 214, 0, 78, 32, 0, 64, 1, 56,
128, 1, 56, 128, 0, 0, 0, 0, 0, 0,
0, 0,
};
static unsigned char gds209[ 42] = {
0, 0, 42, 0, 255, 3, 0, 101, 0, 81,
0, 87, 17, 129, 204, 231, 8, 129, 115, 24,
0, 158, 187, 0, 158, 187, 0, 64, 0, 97,
168, 0, 97, 168, 0, 0, 0, 0, 0, 0,
0, 0,
};
static unsigned char gds210[ 42] = {
0, 0, 42, 0, 255, 1, 0, 25, 0, 25,
0, 35, 40, 129, 44, 200, 128, 0, 103, 54,
128, 229, 1, 0, 78, 32, 0, 64, 1, 56,
128, 1, 56, 128, 0, 0, 0, 0, 0, 0,
0, 0,
};
static unsigned char gds211[ 42] = {
0, 0, 42, 0, 255, 3, 0, 93, 0, 65,
0, 47, 158, 130, 9, 83, 8, 129, 115, 24,
1, 61, 119, 1, 61, 119, 0, 64, 0, 97,
168, 0, 97, 168, 0, 0, 0, 0, 0, 0,
0, 0,
};
static unsigned char gds212[ 42] = {
0, 0, 42, 0, 255, 3, 0, 185, 0, 129,
0, 47, 158, 130, 9, 83, 8, 129, 115, 24,
0, 158, 187, 0, 158, 187, 0, 64, 0, 97,
168, 0, 97, 168, 0, 0, 0, 0, 0, 0,
0, 0,
};
static unsigned char gds213[ 32] = {
0, 0, 32, 0, 255, 5, 0, 129, 0, 85,
0, 30, 158, 130, 38, 228, 8, 129, 154, 40,
1, 116, 18, 1, 116, 18, 0, 64, 0, 0,
0, 0,
};
static unsigned char gds214[ 32] = {
0, 0, 32, 0, 255, 5, 0, 97, 0, 69,
0, 164, 101, 130, 174, 25, 8, 130, 73, 240,
0, 186, 9, 0, 186, 9, 0, 64, 0, 0,
0, 0,
};
static unsigned char gds215[ 42] = {
0, 0, 42, 0, 255, 3, 1, 113, 1, 1,
0, 47, 158, 130, 9, 83, 8, 129, 115, 24,
0, 79, 94, 0, 79, 94, 0, 64, 0, 97,
168, 0, 97, 168, 0, 0, 0, 0, 0, 0,
0, 0,
};
static unsigned char *std_ncep_gds[256] = {
NULL, /* 0 */
gds1,
gds2,
gds3,
gds4,
gds5,
gds6,
NULL, /* 7 */
NULL, /* 8 */
NULL, /* 9 */
NULL, /* 10 */
NULL, /* 11 */
NULL, /* 12 */
NULL, /* 13 */
NULL, /* 14 */
NULL, /* 15 */
NULL, /* 16 */
NULL, /* 17 */
NULL, /* 18 */
NULL, /* 19 */
NULL, /* 20 */
gds21,
gds22,
gds23,
gds24,
gds25,
gds26,
gds27,
gds28,
gds29,
gds30,
NULL, /* 31 */
NULL, /* 32 */
gds33,
gds34,
NULL, /* 35 */
NULL, /* 36 */
gds37,
gds38,
gds39,
gds40,
gds41,
gds42,
gds43,
gds44,
gds45,
NULL, /* 46 */
NULL, /* 47 */
NULL, /* 48 */
NULL, /* 49 */
NULL, /* 50 */
NULL, /* 51 */
NULL, /* 52 */
NULL, /* 53 */
NULL, /* 54 */
gds55,
gds56,
NULL, /* 57 */
NULL, /* 58 */
NULL, /* 59 */
NULL, /* 60 */
gds61,
gds62,
gds63,
gds64,
NULL, /* 65 */
NULL, /* 66 */
NULL, /* 67 */
NULL, /* 68 */
NULL, /* 69 */
NULL, /* 70 */
NULL, /* 71 */
NULL, /* 72 */
NULL, /* 73 */
NULL, /* 74 */
NULL, /* 75 */
NULL, /* 76 */
NULL, /* 77 */
NULL, /* 78 */
NULL, /* 79 */
NULL, /* 80 */
NULL, /* 81 */
NULL, /* 82 */
NULL, /* 83 */
NULL, /* 84 */
gds85,
gds86,
gds87,
NULL, /* 88 */
NULL, /* 89 */
gds90,
gds91,
gds92,
gds93,
gds94,
gds95,
gds96,
gds97,
gds98,
NULL, /* 99 */
gds100,
gds101,
NULL, /* 102 */
gds103,
gds104,
gds105,
gds106,
gds107,
NULL, /* 108 */
NULL, /* 109 */
NULL, /* 110 */
NULL, /* 111 */
NULL, /* 112 */
NULL, /* 113 */
NULL, /* 114 */
NULL, /* 115 */
NULL, /* 116 */
NULL, /* 117 */
NULL, /* 118 */
NULL, /* 119 */
NULL, /* 120 */
NULL, /* 121 */
NULL, /* 122 */
NULL, /* 123 */
NULL, /* 124 */
NULL, /* 125 */
gds126,
NULL, /* 127 */
NULL, /* 128 */
NULL, /* 129 */
NULL, /* 130 */
NULL, /* 131 */
NULL, /* 132 */
NULL, /* 133 */
NULL, /* 134 */
NULL, /* 135 */
NULL, /* 136 */
NULL, /* 137 */
NULL, /* 138 */
NULL, /* 139 */
NULL, /* 140 */
NULL, /* 141 */
NULL, /* 142 */
NULL, /* 143 */
NULL, /* 144 */
NULL, /* 145 */
NULL, /* 146 */
NULL, /* 147 */
NULL, /* 148 */
NULL, /* 149 */
NULL, /* 150 */
NULL, /* 151 */
NULL, /* 152 */
NULL, /* 153 */
NULL, /* 154 */
NULL, /* 155 */
NULL, /* 156 */
NULL, /* 157 */
NULL, /* 158 */
NULL, /* 159 */
NULL, /* 160 */
NULL, /* 161 */
NULL, /* 162 */
NULL, /* 163 */
NULL, /* 164 */
NULL, /* 165 */
NULL, /* 166 */
NULL, /* 167 */
NULL, /* 168 */
NULL, /* 169 */
NULL, /* 170 */
NULL, /* 171 */
NULL, /* 172 */
NULL, /* 173 */
NULL, /* 174 */
NULL, /* 175 */
NULL, /* 176 */
NULL, /* 177 */
NULL, /* 178 */
NULL, /* 179 */
NULL, /* 180 */
NULL, /* 181 */
NULL, /* 182 */
NULL, /* 183 */
NULL, /* 184 */
NULL, /* 185 */
NULL, /* 186 */
NULL, /* 187 */
NULL, /* 188 */
NULL, /* 189 */
NULL, /* 190 */
NULL, /* 191 */
gds192,
NULL, /* 193 */
NULL, /* 194 */
NULL, /* 195 */
NULL, /* 196 */
NULL, /* 197 */
NULL, /* 198 */
NULL, /* 199 */
NULL, /* 200 */
gds201,
gds202,
gds203,
gds204,
gds205,
gds206,
gds207,
gds208,
gds209,
gds210,
gds211,
gds212,
gds213,
gds214,
gds215,
NULL, /* 216 */
NULL, /* 217 */
NULL, /* 218 */
NULL, /* 219 */
NULL, /* 220 */
NULL, /* 221 */
NULL, /* 222 */
NULL, /* 223 */
NULL, /* 224 */
NULL, /* 225 */
NULL, /* 226 */
NULL, /* 227 */
NULL, /* 228 */
NULL, /* 229 */
NULL, /* 230 */
NULL, /* 231 */
NULL, /* 232 */
NULL, /* 233 */
NULL, /* 234 */
NULL, /* 235 */
NULL, /* 236 */
NULL, /* 237 */
NULL, /* 238 */
NULL, /* 239 */
NULL, /* 240 */
NULL, /* 241 */
NULL, /* 242 */
NULL, /* 243 */
NULL, /* 244 */
NULL, /* 245 */
NULL, /* 246 */
NULL, /* 247 */
NULL, /* 248 */
NULL, /* 249 */
NULL, /* 250 */
NULL, /* 251 */
NULL, /* 252 */
NULL, /* 253 */
NULL, /* 254 */
NULL, /* 255 */
};
unsigned char *NCEP_GDS(unsigned char *pds, int grid_type) {
if (std_ncep_gds[grid_type] == NULL) {
set_PDSGridType(pds, grid_type);
clr_HasGDS(pds);
return NULL;
}
set_PDSGridType(pds, grid_type);
set_HasGDS(pds);
return (cpGRIBsec(std_ncep_gds[grid_type]));
}
/* w. ebisuzaki
*
* creates a new copy of a 24section
*
*/
unsigned char *cpGRIBsec(unsigned char *section) {
int len;
unsigned char *newsection;
if (section == NULL) return NULL;
len = _LEN24(section);
if ((newsection = (unsigned char *) malloc(len)) == NULL) {
fprintf(stderr,"ran out of memory in cp_section\n");
exit(8);
}
memcpy(newsection, section, len);
return newsection;
}
/*
* 1996 Wesley Ebisuzaki
*
* wrt_24section
* writes a section (bds, gds, bms, etc)
* (section must start with 3 byte size)
*
*/
int append_24section(unsigned char *section, FILE *output) {
int i;
i = __LEN24(section);
if (fwrite(section, sizeof (unsigned char), i, output) != i) {
fprintf(stderr,"error writing output\n");
return -1;
}
return 0;
}
void set_Center(unsigned char *pds, int center) {
pds[4] = center;
}
int get_Center(unsigned char *pds) {
return pds[4];
}
void set_DecScale(unsigned char *pds, int dec_scale) {
if (dec_scale >= 0 && dec_scale < 32768) {
pds[26] = ((unsigned int) dec_scale) >> 8;
pds[27] = dec_scale & 255;
}
else if (dec_scale < 0 && dec_scale > -32768) {
dec_scale = - dec_scale;
pds[26] = (((unsigned int) dec_scale) >> 8) | 128;
pds[27] = dec_scale & 255;
}
else {
fprintf(stderr,"Bad dec_scale %d\n", dec_scale);
exit(8);
}
}
int get_DecScale(unsigned char *pds) {
return (int) (PDS_DecimalScale(pds));
}
void set_PDSGridType(unsigned char *pds, int n) {
pds[6] = n;
}
int get_PDSGridTYPE (unsigned char *pds) {
return pds[6];
}
void set_HasGDS(unsigned char *pds) {
pds[7] |= 128;
}
void clr_HasGDS(unsigned char *pds) {
pds[7] &= 127;
}
int get_HasGDS(unsigned char *pds) {
return (pds[7] & 128) != 0;
}
void set_HasBMS(unsigned char *pds) {
pds[7] |= 0x40;
}
void clr_HasBMS(unsigned char *pds) {
pds[7] &= 0xbf;
}
int get_HasBMS(unsigned char *pds) {
return (pds[7] & 0x40) != 0;
}