#include #include #include #include #include #include #include "gribw.h" #include "pds4.h" /* * weasd 2 snoc * * Wesley Ebisuzaki * */ #define N 100000 #define SNOC 238 #define WEASD 65 void main(int argc, char **argv) { long int len_grib, pos = 0; unsigned char *pds, *gds, *bms, *bds; FILE *input, *output; int i, nxny; float *array; int array_size; /* preliminaries .. open up all files */ if (argc != 3) { fprintf(stderr, "%s [in gribfile] [out gribfile]\n", argv[0]); 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); } array = (float *) malloc(N * sizeof(float)); array_size = N; if (array == NULL) { fprintf(stderr,"memeory allocation failed\n"); exit(8); } pos = 0; for(;;) { /* loop over all records */ len_grib = rd_grib_msg(input, pos, &pds, &gds, &bms, &bds); if (len_grib <= 0) break; if (PDS_PARAM(pds) == WEASD) { nxny = get_nxny(pds, gds, bms, bds); if (nxny > array_size) { free(array); array_size = nxny; array = (float *) malloc(array_size * sizeof(float)); if (array == NULL) { fprintf(stderr,"memory allocation failed\n"); exit(8); } } unpk_bds(array, pds, gds, bms, bds, nxny); for (i = 0; i < nxny; i++) { if (! UNDEFINED_VAL(array[i])) { array[i] = array[i] > 0.0 ? 100.0 : 0.0; } } pds[8] = SNOC; set_DecScale(pds, -2); bms = mk_BMS(pds, array, &nxny, UNDEFINED_LOW, UNDEFINED_HIGH); bds = mk_BDS(pds, array, nxny); wrt_grib_msg(output, pds, gds, bms, bds); free(bds); if (bms) free(bms); } else { wrt_grib_msg(output, pds, gds, bms, bds); } pos += len_grib; } exit(0); }