print_data.c

print_data.c How to print all the data from a grib message.

00001 
00010 /*
00011  * C Implementation: print_data
00012  *
00013  * Description: prints all the data contained in a grib file
00014  *
00015  * Author: Enrico Fucile <enrico.fucile@ecmwf.int>
00016  *
00017  *
00018  */
00019 #include <stdio.h>
00020 #include <stdlib.h>
00021 
00022 #include "grib_api.h"
00023 
00024 void usage(char* prog) {
00025   printf("usage: %s filename\n",prog);
00026   exit(1);
00027 }
00028 
00029 int main(int argc, char** argv) {
00030   int err = 0,i;
00031   double *values = NULL;
00032   double max,min,average;
00033   size_t values_len= 0;
00034 
00035   FILE* in = NULL;
00036   char* filename ;
00037   grib_handle *h = NULL;
00038 
00039   if (argc<2) usage(argv[0]);
00040   filename=argv[1];
00041 
00042   in = fopen(filename,"r");
00043   if(!in) {
00044     printf("ERROR: unable to open file %s\n",filename);
00045     return 1;
00046   }
00047 
00048   /* create new handle from a message in a file*/
00049   h = grib_handle_new_from_file(0,in,&err);
00050   if (h == NULL) {
00051     printf("Error46   }