static struct gribtab_s *search_gribtab(struct gribtab_s *p, unsigned char **sec){
int discipline, center, mastertab, localtab, parmcat, parmnum;
discipline = GB2_Discipline(sec);
center = GB2_Center(sec);
mastertab = GB2_MasterTable(sec);
localtab = GB2_LocalTable(sec);
parmcat = GB2_ParmCat(sec);
parmnum = GB2_ParmNum(sec);
if (mastertab == 0) mastertab = 1;
for (; p->disc >= 0; p++) {
if (parmnum < 192) {
if (discipline == p->disc && (mastertab == p->mtab || p->mtab == -1) &&
parmcat == p->pcat && parmnum == p->pnum) {
return p;
}
}
else {
if (discipline == p->disc && center == p->cntr &&
localtab == p->ltab &&
parmcat == p->pcat && parmnum == p->pnum) {
return p;
break;
}
}
}
return NULL;
}
unsigned char *flag_table_3_3_location(unsigned char **sec) {
int grid_template;
unsigned char *gds;
grid_template = code_table_3_1(sec);
gds = sec[3];
switch (grid_template) {
case 0:
case 1:
case 2:
case 3:
case 40:
case 41:
case 42:
case 43:
case 140:
case 204:
return gds+54; break;
case 4:
case 5:
case 10:
case 12:
case 20:
case 30:
case 31:
case 90:
case 110:
return gds+46; break;
case 32768:
if (GB2_Center(sec) == NCEP) return gds+54;
return NULL;
case 32769:
if (GB2_Center(sec) == NCEP) return gds+54;
return NULL;
break;
default: return NULL; break;
}
return NULL;
}
int f_center(ARG0) {
int ctr;
char tmp[20];
const char *string;
if (mode >= 0) {
ctr = GB2_Center(sec);
switch (ctr) {
#include "code_table0.dat"
}
sprintf(inv_out,"center=%s", string);
}
return 0;
}
/*
* HEADER:-1:code_table_4.3:inv:0:code table 4.3 Type of Generating Process
*/
int f_code_table_4_3(ARG0) {
int val, center;
const char *string;
if (mode >= 0) {
val = code_table_4_3(sec);
center = GB2_Center(sec);
string = NULL;
switch(val) {
#include "CodeTable_4.3.dat"
}
if (string == NULL) sprintf(inv_out,"code table 4.3=%d", val);
else sprintf(inv_out,"code table 4.3=%d %s", val, string);
}
return 0;
}
static struct gribtable_s *search_gribtable(struct gribtable_s *p, unsigned char **sec){
int discipline, center, mastertab, localtab, parmcat, parmnum;
int use_local_table;
static int count = 0;
if (p == NULL) return NULL;
discipline = GB2_Discipline(sec);
center = GB2_Center(sec);
mastertab = GB2_MasterTable(sec);
localtab = GB2_LocalTable(sec);
parmcat = GB2_ParmCat(sec);
parmnum = GB2_ParmNum(sec);
// // if (mastertab == 0) mastertab = 1;
// if (mastertab >= 0 && mastertab <= 10) mastertab = 1;
use_local_table = (mastertab == 255) ? 1 : 0;
if (parmnum >= 192 || parmcat >= 192) use_local_table = 1;
if (use_local_table == 1 && localtab == 0) {
if (count++ < 6 ) fprintf(stderr,"**** ERROR: local table = 0 is not allowed, set to 1 ***\n");
localtab = 1;
}
if (use_local_table == 1 && localtab == 255) {
fatal_error("local gribtable is undefined (255)","");
}
if (! use_local_table) {
for (; p->disc >= 0; p++) {
if (discipline == p->disc && (mastertab >= p->mtab_low) && (mastertab <= p->mtab_high) &&
parmcat == p->pcat && parmnum == p->pnum) {
return p;
}
}
}
else {
// printf(">> cname local find: disc %d center %d localtab %d pcat %d pnum %d\n", discipline, center, localtab, parmcat, parmnum);
for (; p->disc >= 0; p++) {
if (discipline == p->disc && center == p->cntr && localtab == p->ltab &&
parmcat == p->pcat && parmnum == p->pnum) {
return p;
}
}
}
return NULL;
}
int getName(unsigned char **sec, int mode, char *inv_out, char *name, char *desc, char *unit) {
int discipline, center, mastertab, localtab, parmcat, parmnum;
struct gribtab_s *p;
p = NULL;
#ifdef USE_TIGGE
if (tigge) p = search_gribtab(tigge_gribtab, sec);
#endif
if (p == NULL) p = search_gribtab(gribtab, sec);
if (p) {
if (name) strcpy(name, p->name);
if (desc) strcpy(desc, p->desc);
if (unit) strcpy(unit, p->unit);
if (inv_out) {
sprintf(inv_out, "%s", p->name);
inv_out += strlen(inv_out);
if (mode) sprintf(inv_out," %s [%s]", p->desc, p->unit);
}
}
else {
if (name) *name = 0;
if (desc) *desc = 0;
if (unit) *unit = 0;
discipline = GB2_Discipline(sec);
center = GB2_Center(sec);
mastertab = GB2_MasterTable(sec);
localtab = GB2_LocalTable(sec);
parmcat = GB2_ParmCat(sec);
parmnum = GB2_ParmNum(sec);
if (parmnum < 192) {
sprintf(inv_out,"var discipline=%d master_table=%d parmcat=%d parm=%d",
discipline, mastertab, parmcat, parmnum);
}
else {
sprintf(inv_out,"var discipline=%d center=%d local_table=%d parmcat=%d parm=%d",
discipline, center, localtab, parmcat, parmnum);
}
}
return 0;
}
int f_subcenter(ARG0) {
int ctr, subctr;
char tmp[20];
const char *string;
if (mode >= 0) {
ctr = GB2_Center(sec);
subctr = GB2_Subcenter(sec);
string = NULL;
if (ctr == 7) {
switch (subctr) {
#include "ncep_tableC.dat"
}
}
if (string == NULL) sprintf(inv_out,"subcenter=%d", subctr);
else sprintf(inv_out,"subcenter=%s", string);
}
return 0;
}
int fix_ncep_3(unsigned char **sec) {
int j;
if (GB2_Center(sec) != NCEP) return 0; // only for NCEP files
j = code_table_5_0(sec);
if (j == 0 || j == 1 || j == 40 || j == 41 || j == 40000 || j == 40010) { // simple, jpeg, png
if (sec[5][19] == 0) {
sec[5][17] = sec[5][18] = 0;
}
}
else if (j == 2 || j == 3) { // complex packing
if (uint4(sec[5]+31) == 0) { // number of groups
sec[5][17] = sec[5][18] = 0;
}
}
else if (j == 50 || j == 51) { // spectral
if (sec[5][19] == 0) {
sec[5][17] = sec[5][18] = 0;
}
}
return 0;
}
/*
* HEADER:-1:code_table_4.230:inv:0:code table 4.230 chemical constituent type
*/
int f_code_table_4_230(ARG0) {
int val;
const char *string;
static int error_count = 0;
if (mode >= 0) {
val = code_table_4_230(sec);
if (val >= 0) {
string = NULL;
if (GB2_MasterTable(sec) <= 4) {
if (error_count++ <= 10) fprintf(stderr,
"Warning: if file made with ECMWF API, may be using incompatible chemistry table\n");
}
switch(val) {
#include "CodeTable_4.230.dat"
}
if (GB2_MasterTable(sec) <= 4 && GB2_Center(sec) == ECMWF) string = NULL;
if (string == NULL) sprintf(inv_out,"code table 4.230=%d", val);
else sprintf(inv_out,"code table 4.230=%d %s", val, string);
}
}
return 0;
}
unsigned char *code_table_3_2_location(unsigned char **sec) {
int grid_def;
grid_def = code_table_3_1(sec);
if (grid_def < 50) return sec[3]+14;
switch (grid_def) {
case 90:
case 110:
case 130:
case 140:
case 204:
case 1000:
case 1100:
return sec[3]+14; break;
default: break;
}
if (GB2_Center(sec) == NCEP) {
if (grid_def == 32768 || (grid_def == 32769)) return sec[3]+14;
}
return NULL;
}
int main(int argc, char **argv) {
#else
int wgrib2(int argc, const char **argv) {
#endif
//WNE FILE *in;
struct seq_file in_file;
unsigned char *msg, *sec[10]; /* sec[9] = last valid bitmap */
long int last_pos;
int file_arg, i, j, num_submsgs;
int n_arg;
unsigned int k, ndata;
int err_4_3_count;
float *data;
double ref;
// double *ddata, ref;
#ifdef USE_G2CLIB
float missing_c_val_1, missing_c_val_2;
g2int *bitmap, has_bitmap;
g2float *g2_data;
int ii;
#endif
struct ARGLIST arglist[N_ARGLIST];
int narglist;
const char *new_argv[N_ARGLIST];
void *local[N_ARGLIST];
int has_inv_option, last_submsg;
int err, new_GDS, center;
unsigned char dscale[2];
init_globals();
init_inv_out();
ndata = 0;
err_4_3_count = 0;
if (initial_call) { /* only done 1st time */
setup_user_gribtable();
// jas_init();
// gctpc initialiation
init(-1,-1,"gctpc_error,txt", "gctpc_param.txt");
initial_call = 0;
}
narglist = 0;
dscale[0] = dscale[1] = 0;
mode = 0;
if (fopen_file(&(rd_inventory_input), "-", "r")) fatal_error("opening stdin for rd_inventory","");
data = NULL;
// ddata = NULL;
#ifdef CALLABLE_WGRIB2
if (setjmp(fatal_err)) {
fprintf(stderr,"*** arg list to wgrib2:");
for (i=0; i < argc; i++) {
fprintf(stderr," %s", argv[i]);
}
fprintf(stderr,"\n\n");
if (ndata && data != NULL) free(data);
ndata=0;
return 1;
}
#endif
/* no arguments .. help screen */
if (argc == 1) {
mode = -1;
inv_out[0] = 0;
f_h(call_ARG0(inv_out,NULL));
// fprintf(inv_file, "%s\n", inv_out);
i = strlen(inv_out);
inv_out[i++] = '\n';
inv_out[i] = '\0';
fwrite_file(inv_out, 1, i, &inv_file);
err_bin(1); err_string(1);
return 8;
}
/* copy argv */
for (i = 0; i < argc; i++) {
new_argv[i] = argv[i];
}
/* scan for "inv" and input file */
has_inv_option = 0;
file_arg = 0;
for (i = 1; i < argc; i++) {
if (new_argv[i][0] != '-' || (strcmp(new_argv[i],"-") == 0) ) {
/* must be filename */
if (file_arg == 0) {
file_arg = i;
continue;
} else {
fatal_error_ss("too many grib files .. 1st=%s 2nd=%s", new_argv[file_arg], new_argv[i]);
}
}
/* must be an option */
for (j = 0; j < nfunctions; j++) {
if (strcmp(&(new_argv[i][1]),functions[j].name) == 0) {
if (functions[j].type == inv) has_inv_option = 1;
i += functions[j].nargs;
break;
}
}
if (j == nfunctions) {
fatal_error("undefined option: %s", new_argv[i]);
}
}
/* if no inv option, use default inventory .. put it at end */
if (has_inv_option == 0) new_argv[argc++] = "-s";
/* try opening the input file */
in_file.file_type = NOT_OPEN;
if (file_arg != 0) {
fopen_file(&in_file, new_argv[file_arg],"rb");
}
/* "compile" options */
#ifdef DEBUG
fprintf(stderr,"going to compile phase\n");
#endif
for (i = 1; i < argc; i++) {
/* if filename skip */
if (new_argv[i][0] != '-' || (strcmp(new_argv[i],"-") == 0)) continue;
/* must be an option */
for (j = 0; j < nfunctions; j++) {
if (strcmp(&(new_argv[i][1]),functions[j].name) == 0) break;
}
if (j == nfunctions) fatal_error("unknown option %s", new_argv[i]);
/* add to function argument list */
arglist[narglist].fn = j;
arglist[narglist].i_argc = i+1;
if (functions[j].type == inv) has_inv_option = 1;
i += functions[j].nargs;
if (i >= argc) fatal_error("missing arguments option=%s",functions[j].name);
narglist++;
}
if (has_inv_option == 0) fatal_error("missing arguments on last option","");
/* initialize options, mode = -1 */
#ifdef DEBUG
fprintf(stderr,"going to init options, narglist %d\n",narglist);
#endif
for (j = 0; j < narglist; j++) {
new_inv_out(); /* inv_out[0] = 0; */
n_arg = functions[arglist[j].fn].nargs;
err = 0;
#ifdef DEBUG
fprintf(stderr,"going to init option %s\n", functions[arglist[j].fn].name);
#endif
if (n_arg == 0) err = functions[arglist[j].fn].fn(-1,NULL,NULL,0, inv_out,local+j);
else if (n_arg == 1) err = functions[arglist[j].fn].fn(-1,NULL,NULL,0, inv_out,local+j,
new_argv[arglist[j].i_argc]);
else if (n_arg == 2) err = functions[arglist[j].fn].fn(-1,NULL,NULL,0, inv_out,local+j,
new_argv[arglist[j].i_argc],new_argv[arglist[j].i_argc+1]);
else if (n_arg == 3) err = functions[arglist[j].fn].fn(-1,NULL,NULL,0, inv_out,local+j,
new_argv[arglist[j].i_argc],new_argv[arglist[j].i_argc+1],new_argv[arglist[j].i_argc+2]);
else if (n_arg == 4) err = functions[arglist[j].fn].fn(-1,NULL,NULL,0, inv_out,local+j,
new_argv[arglist[j].i_argc],new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2], new_argv[arglist[j].i_argc+3]);
else if (n_arg == 5) err = functions[arglist[j].fn].fn(-1,NULL,NULL,0, inv_out,local+j,
new_argv[arglist[j].i_argc],new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2], new_argv[arglist[j].i_argc+3],
new_argv[arglist[j].i_argc+4]);
else if (n_arg == 6) err = functions[arglist[j].fn].fn(-1,NULL,NULL,0, inv_out,local+j,
new_argv[arglist[j].i_argc],new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2], new_argv[arglist[j].i_argc+3],
new_argv[arglist[j].i_argc+4], new_argv[arglist[j].i_argc+5]);
else if (n_arg == 7) err = functions[arglist[j].fn].fn(-1,NULL,NULL,0, inv_out,local+j,
new_argv[arglist[j].i_argc],new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2], new_argv[arglist[j].i_argc+3],
new_argv[arglist[j].i_argc+4], new_argv[arglist[j].i_argc+5],
new_argv[arglist[j].i_argc+6]);
else if (n_arg == 8) err = functions[arglist[j].fn].fn(-1,NULL,NULL,0, inv_out,local+j,
new_argv[arglist[j].i_argc],new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2], new_argv[arglist[j].i_argc+3],
new_argv[arglist[j].i_argc+4], new_argv[arglist[j].i_argc+5],
new_argv[arglist[j].i_argc+6], new_argv[arglist[j].i_argc+7]);
// if(inv_out[0] != 0) fprintf(inv_file, "%s", inv_out);
if(inv_out[0] != 0) fwrite_file(inv_out, 1, strnlen(inv_out,INV_BUFFER), &inv_file);
if (err) {
err_bin(1); err_string(1);
// cleanup
return 8;
}
}
/* error and EOF handlers have been initialized */
#ifdef DEBUG
fprintf(stderr,"initial error and EOF handlers\n");
#endif
if (has_inv_option == 0) fatal_error("missing arguments on last option","");
if (in_file.file_type == NOT_OPEN) {
if (file_arg == 0) fatal_error("no input file defined","");
else fatal_error("missing input file %s", new_argv[file_arg]);
}
if (latlon == 1 && output_order_wanted != wesn)
fatal_error("latitude-longitude information is only available with -order we:sn","");
if (input == inv_mode && (in_file.file_type != DISK && in_file.file_type != MEM))
fatal_error("wgrib2 cannot random access grib input file","");
#ifdef DEBUG
fprintf(stderr,"going to process data\n");
#endif
msg_no = 1;
inv_no = 0;
len = pos = 0;
submsg = 0;
msg = NULL;
last_pos = -1;
last_submsg = -1;
/* if dump mode .. position io stream */
if (input == dump_mode) {
while (msg_no < dump_msg) {
// // if (in_file.file_type == PIPE) msg = rd_grib2_msg_seq(sec, in, &pos, &len, &num_submsgs);
// if (in_file.file_type == PIPE) msg = rd_grib2_msg_seq_file(sec, &in_file, &pos, &len, &num_submsgs);
// else if (in_file.file_type == DISK) msg = rd_grib2_msg(sec, in, &pos, &len, &num_submsgs);
msg = rd_grib2_msg_seq_file(sec, &in_file, &pos, &len, &num_submsgs);
if (msg == NULL) fatal_error_i("record %d not found", dump_msg);
last_pos = pos;
pos += len;
msg_no++;
}
#ifdef DEBUG
printf("dump mode msg=%d\n", msg_no);
#endif
}
/*
* submsg = 0 .. beginning of unread record
* submsg = i .. start at ith submsg
* num_submsgs = number of submessages in grib message
*/
/* inventory loop */
for (;last_message == 0;) {
/* need position and submessage number of message */
if (input == inv_mode || input == dump_mode) {
if (input == inv_mode) {
if (rd_inventory(&msg_no,&submsg, &pos, &rd_inventory_input)) break;
if (fseek_file(&in_file, pos,SEEK_SET) != 0) fatal_error("fseek_file failed","");
}
else if (input == dump_mode) {
if (dump_msg == -1) break;
submsg = dump_submsg;
dump_msg = -1;
}
if (pos != last_pos) {
// // if (in_file.file_type == PIPE) msg = rd_grib2_msg_seq(sec, in, &pos, &len, &num_submsgs);
// if (in_file.file_type == PIPE) msg = rd_grib2_msg_seq_file(sec, &in_file, &pos, &len, &num_submsgs);
// else if (in_file.file_type == DISK) msg = rd_grib2_msg(sec, in, &pos, &len, &num_submsgs);
msg = rd_grib2_msg_seq_file(sec, &in_file, &pos, &len, &num_submsgs);
if (msg == NULL) {
fatal_error_i("grib message #%d not found", msg_no);
break;
}
last_pos = pos;
last_submsg = -1;
}
if (pos == last_pos && submsg == last_submsg + 1) {
/* read previous submessage */
if (parse_next_msg(sec) != 0) {
fprintf(stderr,"\n*** grib message #%d.%d not found ***\n\n", msg_no, submsg);
break;
}
}
else {
/* need to get desired submessage into sec */
if (parse_1st_msg(sec) != 0) {
fprintf(stderr,"\n*** grib message #%d.1 not found ***\n\n", msg_no);
break;
}
for (i = 2; i <= submsg; i++) {
if (parse_next_msg(sec) != 0) {
fprintf(stderr,"\n*** grib message #%d.%d not found ***\n\n", msg_no, i);
break;
}
}
}
last_submsg = submsg;
}
else if (input == all_mode) {
if (submsg == 0) {
// // if (in_file.file_type == PIPE) msg = rd_grib2_msg_seq(sec, in, &pos, &len, &num_submsgs);
// if (in_file.file_type == PIPE) msg = rd_grib2_msg_seq_file(sec, &in_file, &pos, &len, &num_submsgs);
// else if (in_file.file_type == DISK) msg = rd_grib2_msg(sec, in, &pos, &len, &num_submsgs);
msg = rd_grib2_msg_seq_file(sec, &in_file, &pos, &len, &num_submsgs);
if (msg == NULL) break;
submsg = 1;
}
else if (submsg > num_submsgs) {
pos += len;
msg_no++;
// // if (in_file.file_type == PIPE) msg = rd_grib2_msg_seq(sec, in, &pos, &len, &num_submsgs);
// if (in_file.file_type == PIPE) msg = rd_grib2_msg_seq_file(sec, &in_file, &pos, &len, &num_submsgs);
// else if (in_file.file_type == DISK) msg = rd_grib2_msg(sec, in, &pos, &len, &num_submsgs);
msg = rd_grib2_msg_seq_file(sec, &in_file, &pos, &len, &num_submsgs);
if (msg == NULL) break;
submsg = 1;
}
if (submsg == 1) {
if (parse_1st_msg(sec) != 0) {
fprintf(stderr,"illegal format: parsing 1st submessage\n");
}
}
else {
if (parse_next_msg(sec) != 0) {
fprintf(stderr,"illegal format: parsing submessages\n");
}
}
}
if (only_submsg > 0 && only_submsg != submsg) {
submsg++;
continue;
}
if (for_mode) {
if (msg_no < for_start || msg_no > for_end || ((msg_no - for_start) % for_step) != 0) {
if (msg_no > for_end && input != inv_mode) last_message = 1;
submsg++;
continue;
}
}
/* move inv_no++ before match_inv is made */
inv_no++;
if (match || match_fs) {
inv_out[0] = 0;
if (num_submsgs > 1) {
sprintf(inv_out,"%d.%d:", msg_no, submsg);
}
else {
sprintf(inv_out,"%d:", msg_no);
}
f_match_inv(call_ARG0(inv_out+strlen(inv_out), NULL));
if (is_match_fs(inv_out) != 0) {
submsg++;
inv_no--;
continue;
}
#ifdef USE_REGEX
if (is_match(inv_out) != 0) {
submsg++;
inv_no--;
continue;
}
#endif
}
match_flag = 0;
if (for_n_mode) {
if (inv_no < for_n_start || inv_no > for_n_end || ((inv_no - for_n_start) % for_n_step) != 0) {
if (inv_no > for_n_end) last_message = 1;
submsg++;
continue;
}
}
/* see if new GDS */
if ((i = (int) GB2_Sec3_size(sec)) != old_GDS_size) {
new_GDS = 1;
}
else {
new_GDS = 0;
for (j = 0; j < i; j++) {
if (old_gds[j] != sec[3][j]) new_GDS = 1;
}
}
if (new_GDS) {
GDS_change_no++;
if (i > GDS_max_size) {
if (GDS_max_size) free(old_gds);
GDS_max_size = i + 100; /* add 100 just to avoid excessive memory allocations */
if ((old_gds = (unsigned char *) malloc(GDS_max_size) ) == NULL) {
fatal_error("memory allocation problem old_gds in wgrib2.main","");
}
}
for (j = 0; j < i; j++) {
old_gds[j] = sec[3][j];
}
old_GDS_size = i;
/* update grid information */
get_nxny(sec, &nx, &ny, &npnts, &res, &scan); /* get nx, ny, and scan mode of grid */
get_nxny_(sec, &nx_, &ny_, &npnts, &res, &scan); /* get nx, ny, and scan mode of grid */
output_order = (nx_ < 1 || ny_ < 1) ? raw : output_order_wanted;
if (latlon) {
i = 1;
#ifdef USE_PROJ4
if (use_proj4 && i != 0) { /* use Proj4 to get lat lon values */
i = proj4_get_latlon(sec, &lon, &lat);
// if (i == 0) fprintf(stderr,"proj4_get_lat used\n");
}
#endif
if (use_gctpc && i != 0) { /* use gctpc to get lat lon values */
i = gctpc_get_latlon(sec, &lon, &lat);
// if (i == 0) fprintf(stderr,"gctpc_get_lat used\n");
}
if (i != 0) get_latlon(sec, &lon, &lat); /* get lat lon of grid points using built-in code */
}
}
/* Decode NDFD WxText */
if (WxText) mk_WxKeys(sec);
// any fixes to raw grib message before decode need to be placed here
if (fix_ncep_2_flag) fix_ncep_2(sec);
if (fix_ncep_3_flag) fix_ncep_3(sec);
if (fix_ncep_4_flag) fix_ncep_4(sec);
#ifdef CHECK
j = code_table_5_0(sec); // type of compression
/* yes this can be simplified but want to split it up in case other decoders have problems */
if (j == 0 && sec[5][19] == 0 && int2(sec[5] + 17) != 0 && ieee2flt(sec[5]+11) != 0.0)
fprintf(stderr,"Warning: g2lib/g2clib/grib-api simple encode/decode may differ from WMO standard, use -g2clib 0 for WMO standard\n");
if ((j == 2 || j == 3) && int2(sec[5]+17) != 0 && int4(sec[5] + 31) == 0 && ieee2flt(sec[5]+11) != 0.0)
fprintf(stderr,"Warning: g2lib/g2clib complex encode/decode may differ from WMO standard, use -g2clib 0 for WMO standard\n");
if (j == 40 && sec[5][19] == 0 && int2(sec[5] + 17) != 0 && ieee2flt(sec[5]+11) != 0.0)
fprintf(stderr,"Warning: g2lib/g2clib jpeg encode/deocde may differ from WMO standard, use -g2clib 0 for WMO standard\n");
if (j == 41 && sec[5][19] == 0 && int2(sec[5] + 17) != 0 && ieee2flt(sec[5]+11) != 0.0)
fprintf(stderr,"Warning: g2lib/g2clib/grib-api png encode/decode may differ from WMO standard, use -g2clib 0 for WMO standard\n");
/* check the size of Section 7 */
/* code to check the other sizes needs to be placed in decode routines */
j = code_table_5_0(sec); // type of compression
if (j == 0) { /* simple */
/* to avoid overflow on 32 bit machines */
/* old: k = (GB2_Sec5_nval(sec) * sec[5][19] + 7) / 8 + 5; */
k = 5 + (GB2_Sec5_nval(sec)/8) * sec[5][19] + (GB2_Sec5_nval(sec)%8) * (sec[5][19]/8)
+ ( (GB2_Sec5_nval(sec)%8) * (sec[5][19]%8) + 7) / 8;
if (k != GB2_Sec7_size(sec)) {
fprintf(stderr,"Detected a size mismatch, Section 7, wanted %d found %d\n", k, GB2_Sec7_size(sec));
if (decode) fatal_error("Section 7 size, mismatch, simple packing","");
}
}
else if (j == 4) { /* IEEE */
k = GB2_Sec5_nval(sec) * 4 + 5;
if (k != GB2_Sec7_size(sec)) {
fprintf(stderr,"Detected a size mismatch, Section 7, wanted %d found %d\n", k, GB2_Sec7_size(sec));
if (decode) fatal_error("Section 7 size, mismatch, IEEE packing","");
}
}
/* code table 4.3 can change units, warn if undefined */
if (err_4_3_count < 2) {
if (code_table_4_3(sec) == 255) {
fprintf(stderr,"** WARNING input Code Table 4.3 = 255 (undefined) **\n");
err_4_3_count++;
}
}
#endif
if (decode) {
#ifdef CHECK
if (code_table_6_0(sec) == 0) { // has bitmap
k = GB2_Sec3_npts(sec) - GB2_Sec5_nval(sec);
if (k != missing_points(sec[6]+6, GB2_Sec3_npts(sec)))
fatal_error_uu("inconsistent number of bitmap points sec3-sec5: %u sec6: %u",
k, missing_points(sec[6]+6, GB2_Sec3_npts(sec)));
}
else if (code_table_6_0(sec) == 255) { // no bitmap
if (GB2_Sec3_npts(sec) != GB2_Sec5_nval(sec))
fatal_error_uu("inconsistent number of data points sec3: %u sec5: %u",
GB2_Sec3_npts(sec), GB2_Sec5_nval(sec));
}
#endif
/* allocate data */
if (GB2_Sec3_npts(sec) != ndata) {
if (ndata) free(data);
ndata = GB2_Sec3_npts(sec);
if (ndata) {
data = (float *) malloc(sizeof(float) * (size_t) ndata);
if (data == NULL) {
ndata = 0;
fatal_error("main: memory allocation failed data","");
}
}
else { data = NULL; }
}
j = code_table_5_0(sec); // type of compression
/* USE G2CLIB */
#ifdef USE_G2CLIB
if (use_g2clib == 2) {
err = g2_getfld(msg,submsg,1,1,&grib_data);
if (err != 0) fatal_error_ii("Fatal g2clib decode err=%d msg=%d", err, msg_no);
free_gribfield = 1;
has_bitmap = grib_data->ibmap;
g2_data = &(grib_data->fld[0]);
if (has_bitmap == 0 || has_bitmap == 254) {
bitmap = grib_data->bmap;
for (k = 0; k < ndata; k++) {
data[k] = (bitmap[k] == 0) ? UNDEFINED : g2_data[k];
}
}
else {
for (k = 0; k < ndata; k++) {
data[k] = g2_data[k];
}
}
/* complex packing uses special values for undefined */
ii = sub_missing_values(sec, &missing_c_val_1, &missing_c_val_2);
if (ii == 1) {
for (k = 0; k < ndata; k++) {
if (data[k] == missing_c_val_1) data[k] = UNDEFINED;
}
}
else if (ii == 2) {
for (k = 0; k < ndata; k++) {
if (data[k] == missing_c_val_1) data[k] = UNDEFINED;
if (data[k] == missing_c_val_2) data[k] = UNDEFINED;
}
}
}
#endif
/* USE INTERNAL DECODER */
if (use_g2clib != 2) {
center = GB2_Center(sec);
if (use_g2clib == 1) { // introduce g2clib constant field error
/* g2clib ignores decimal scaling for constant fields make internal decoders look like g2clib*/
if ( (j == 0 && sec[5][19] == 0) || ((j == 2 || j == 3) && int4(sec[5] + 31) == 0) ||
(j == 40 && sec[5][19] == 0) || (j == 41 && sec[5][19] == 0) ||
(center == NCEP && j == 40000 && sec[5][19] == 0) ||
(center == NCEP && j == 40010 && sec[5][19] == 0) ) {
dscale[0] = sec[5][17];
dscale[1] = sec[5][18];
sec[5][17] = sec[5][18] = 0;
}
}
err = unpk_grib(sec, data);
if (err != 0) fatal_error_i("Fatal decode packing type %d",err);
if (use_g2clib == 1) { // fix up data
/* restore decimal scaling */
if ( (j == 0 && sec[5][19] == 0) || ((j == 2 || j == 3) && int4(sec[5] + 31) == 0) ||
(j == 40 && sec[5][19] == 0) || (j == 41 && sec[5][19]) ||
(center == NCEP && j == 40000 && sec[5][19] == 0) ||
(center == NCEP && j == 40010 && sec[5][19] == 0) ) {
sec[5][17] = dscale[0];
sec[5][18] = dscale[1];
}
}
}
/* convert to standard output order we:sn */
if (output_order_wanted == wesn) to_we_sn_scan(data, scan, npnts, nx, ny, save_translation);
else if (output_order_wanted == wens) to_we_ns_scan(data, scan, npnts, nx, ny, save_translation);
}
else {
if (ndata) free(data);
ndata = 0;
data = NULL;
}
/* get scaling parameters */
use_scale = input_scale = scaling(sec, &ref, &dec_scale, &bin_scale, &i) == 0;
/* make sure msg_no:pos is put in inv_out so that -last will work */
new_inv_out(); // inv_out[0] = 0;
if (num_submsgs > 1) {
sprintf(inv_out, "%d.%d%s%ld", msg_no, submsg, ":", pos);
}
else {
sprintf(inv_out, "%d%s%ld", msg_no, ":", pos);
}
// fprintf(inv_file, "%s", inv_out);
fwrite_file(inv_out, 1, strnlen(inv_out,INV_BUFFER), &inv_file);
for (j = 0; j < narglist; j++) {
/* skip execution if match_flag == 1 */
/* an output option acts as endif for match_flag */
if (match_flag == 1) {
if (functions[arglist[j].fn].type == output) match_flag = 0;
continue;
}
// if (functions[arglist[j].fn].type == inv) fprintf(inv_file, "%s", item_deliminator);
if (functions[arglist[j].fn].type == inv) fwrite_file(item_deliminator, 1, strlen(item_deliminator), &inv_file);
if (functions[arglist[j].fn].type != setup) {
new_inv_out(); // inv_out[0] = 0;
n_arg = functions[arglist[j].fn].nargs;
if (n_arg == 0)
functions[arglist[j].fn].fn(mode, sec, data, ndata, inv_out, local+j);
else if (n_arg == 1)
functions[arglist[j].fn].fn(mode, sec, data, ndata, inv_out, local+j,
new_argv[arglist[j].i_argc]);
else if (n_arg == 2)
functions[arglist[j].fn].fn(mode, sec, data, ndata, inv_out, local+j,
new_argv[arglist[j].i_argc], new_argv[arglist[j].i_argc+1]);
else if (n_arg == 3)
functions[arglist[j].fn].fn(mode, sec, data, ndata, inv_out, local+j,
new_argv[arglist[j].i_argc], new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2]);
else if (n_arg == 4)
functions[arglist[j].fn].fn(mode, sec, data, ndata, inv_out, local+j,
new_argv[arglist[j].i_argc], new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2], new_argv[arglist[j].i_argc+3]);
else if (n_arg == 5)
functions[arglist[j].fn].fn(mode, sec, data, ndata, inv_out, local+j,
new_argv[arglist[j].i_argc], new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2], new_argv[arglist[j].i_argc+3],
new_argv[arglist[j].i_argc+4]);
else if (n_arg == 6)
functions[arglist[j].fn].fn(mode, sec, data, ndata, inv_out, local+j,
new_argv[arglist[j].i_argc], new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2], new_argv[arglist[j].i_argc+3],
new_argv[arglist[j].i_argc+4], new_argv[arglist[j].i_argc+5]);
else if (n_arg == 7)
functions[arglist[j].fn].fn(mode, sec, data, ndata, inv_out, local+j,
new_argv[arglist[j].i_argc], new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2], new_argv[arglist[j].i_argc+3],
new_argv[arglist[j].i_argc+4], new_argv[arglist[j].i_argc+5],
new_argv[arglist[j].i_argc+6]);
else if (n_arg == 8)
functions[arglist[j].fn].fn(mode, sec, data, ndata, inv_out, local+j,
new_argv[arglist[j].i_argc], new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2], new_argv[arglist[j].i_argc+3],
new_argv[arglist[j].i_argc+4], new_argv[arglist[j].i_argc+5],
new_argv[arglist[j].i_argc+6], new_argv[arglist[j].i_argc+7]);
// if(inv_out[0] != 0) fprintf(inv_file, "%s", inv_out);
if(inv_out[0] != 0) fwrite_file(inv_out, 1, strnlen(inv_out,INV_BUFFER), &inv_file);
}
}
#ifdef CHECK
if (!decode) {
if (code_table_6_0(sec) == 0) { // has bitmap
k = GB2_Sec3_npts(sec) - GB2_Sec5_nval(sec);
if (k != missing_points(sec[6]+6, GB2_Sec3_npts(sec)))
fatal_error_uu("inconsistent number of bitmap points sec3-sec5: %u sec6: %u",
k, missing_points(sec[6]+6, GB2_Sec3_npts(sec)));
}
else if (code_table_6_0(sec) == 255) { // no bitmap
if (GB2_Sec3_npts(sec) != GB2_Sec5_nval(sec))
fatal_error_ii("inconsistent number of data points sec3: %d sec5: %d",
(int) GB2_Sec3_npts(sec), (int) GB2_Sec5_nval(sec));
}
}
#endif
submsg++;
#ifdef USE_G2CLIB
if (free_gribfield) { g2_free(grib_data); free_gribfield = 0;}
#endif
// fprintf(inv_file, "%s",end_inv);
fwrite_file(end_inv, 1, strlen(end_inv), &inv_file);
if (flush_mode) fflush_file(&inv_file);
if (dump_msg > 0) break;
}
/* finalize all functions, call with mode = -2 */
err = 0;
for (j = 0; j < narglist; j++) {
// if (functions[arglist[j].fn].type != setup) {
n_arg = functions[arglist[j].fn].nargs;
new_inv_out(); // inv_out[0] = 0;
if (n_arg == 0)
err |= functions[arglist[j].fn].fn(-2, NULL, NULL, 0, inv_out, local+j);
else if (n_arg == 1)
err |= functions[arglist[j].fn].fn(-2, NULL, NULL, 0, inv_out, local+j,
new_argv[arglist[j].i_argc]);
else if (n_arg == 2)
err |= functions[arglist[j].fn].fn(-2, NULL, NULL, 0, inv_out, local+j,
new_argv[arglist[j].i_argc], new_argv[arglist[j].i_argc+1]);
else if (n_arg == 3)
err |= functions[arglist[j].fn].fn(-2, NULL, NULL, 0, inv_out, local+j,
new_argv[arglist[j].i_argc], new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2]);
else if (n_arg == 4)
err |= functions[arglist[j].fn].fn(-2, NULL, NULL, 0, inv_out, local+j,
new_argv[arglist[j].i_argc], new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2], new_argv[arglist[j].i_argc+3]);
else if (n_arg == 5)
err |= functions[arglist[j].fn].fn(-2, NULL, NULL, 0, inv_out, local+j,
new_argv[arglist[j].i_argc], new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2], new_argv[arglist[j].i_argc+3],
new_argv[arglist[j].i_argc+4]);
else if (n_arg == 6)
err |= functions[arglist[j].fn].fn(-2, NULL, NULL, 0, inv_out, local+j,
new_argv[arglist[j].i_argc], new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2], new_argv[arglist[j].i_argc+3],
new_argv[arglist[j].i_argc+4], new_argv[arglist[j].i_argc+5]);
else if (n_arg == 7)
err |= functions[arglist[j].fn].fn(-2, NULL, NULL, 0, inv_out, local+j,
new_argv[arglist[j].i_argc], new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2], new_argv[arglist[j].i_argc+3],
new_argv[arglist[j].i_argc+4], new_argv[arglist[j].i_argc+5],
new_argv[arglist[j].i_argc+6]);
else if (n_arg == 8)
err |= functions[arglist[j].fn].fn(-2, NULL, NULL, 0, inv_out, local+j,
new_argv[arglist[j].i_argc], new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2], new_argv[arglist[j].i_argc+3],
new_argv[arglist[j].i_argc+4], new_argv[arglist[j].i_argc+5],
new_argv[arglist[j].i_argc+6], new_argv[arglist[j].i_argc+7]);
// if (inv_out[0]) fprintf(stderr, "%s\n", inv_out);
if (inv_out[0]) fprintf(stderr, "%s%s", inv_out, end_inv);
// }
}
err_bin(0); err_string(0);
fclose_file(&in_file);
if (ndata) {
ndata = 0;
free(data);
}
// return 0;
return err;
}
void set_mode(int new_mode) {
mode = new_mode;
}
int f_mysql_speed(ARG7) {
char temp_pathname[STRING_SIZE];
char sql[1500];
char server[100];
char user[100];
char password[100];
char database[100];
char table[100];
MYSQL_RES *res;
MYSQL_ROW row;
char name[100], desc[100], unit[100], level_buf[100];
int year, month, day, hour, minute, second;
char vt[20],rt[20];
// unsigned char *p;
int i,j;
double longitude,latitude,value;
char new_level[50];
char new_name[50];
char conv[50];
char precision[10];
char last;
char param[50];
int ctr;
int level;
unsigned int load_local_infile;
struct local_struct {
MYSQL *conn;
FILE *temp_fileptr;
unsigned int npts;
char *rows[MAX_NXNY];
char has_value[MAX_NXNY];
char *params;
unsigned int isset;
unsigned int wlon;
unsigned int remove_unlikely;
char runtime[20], validtime[20];
int last_GDS_change_no;
};
struct local_struct *save;
strcpy(server,arg1);
strcpy(user,arg2);
strcpy(password,arg3);
strcpy(database,arg4);
strcpy(table,arg5);
/* initialization phase */
if (mode == -1) {
decode = latlon = 1;
*local = save = (struct local_struct *) malloc( sizeof(struct local_struct));
if (save == NULL) fatal_error("mysql_speed memory allocation ","");
for (i = 0; i < MAX_NXNY; i++) {
save->rows[i] = NULL;
save->has_value[i] = 0;
}
save->last_GDS_change_no = 0;
save->conn = mysql_init(NULL);
save->temp_fileptr= NULL;
save->params = (char *) malloc(1500*sizeof(char));
sprintf(save->params,"%s","");
save->isset = 0;
save->runtime[0] = 0;
save->validtime[0] = 0;
load_local_infile = 1; // 1 = LOAD LOCAL INFILE, 0 = do not LOAD LOCAL INFILE
if (sscanf(arg6,"%d", &save->wlon) != 1) {
fatal_error("Argument 6, use western longitudes, has to be 0 or 1, error parsing %s", arg6);
}
if (sscanf(arg7,"%d", &save->remove_unlikely) != 1) {
fatal_error("Argument 7, remove unlikely values, has to be 0 or 1, error parsing %s", arg7);
}
/* Set options for database */
mysql_options(save->conn,MYSQL_OPT_LOCAL_INFILE, (char *) &load_local_infile);
/* Connect to database */
if (!mysql_real_connect(save->conn, server, user, password, database, 0, NULL, 0)) {
fatal_error("f_mysql_speed: could not connect to %s", mysql_error(save->conn));
}
return 0;
}
/* cleanup phase */
if (mode == -2) {
save = (struct local_struct *) *local;
strcpy(temp_pathname, "/tmp/wgrib2_mysqlXXXXXX");
if ( -1 == (load_local_infile = mkstemp(temp_pathname)) ) {
fatal_error("f_mysql_speed: error making temporary filename","");
}
if ( !(save->temp_fileptr = fdopen(load_local_infile, "w")) ) {
fatal_error("f_mysql_speed: error making temporary filename","");
}
fprintf(stdout, "Columns to insert: rt,vt,lat,lon%s\n", save->params);
fprintf(save->temp_fileptr,"rt,vt,lat,lon%s\n",save->params);
for (j = 0; j < save->npts; j++) {
if (save->has_value[j] == 1) fprintf(save->temp_fileptr,"%s\n",save->rows[j]);
}
fflush(save->temp_fileptr);
sprintf(sql,"LOAD DATA LOCAL INFILE '%s' INTO TABLE %s FIELDS TERMINATED BY ',' ENCLOSED BY '\"' LINES TERMINATED BY '\\n' IGNORE 1 LINES (rt,vt,lat,lon%s)",temp_pathname,table,save->params);
printf("Inserting into database, table=%s \n",table);
printf("%s\n",sql);
/* send SQL query */
if (mysql_query(save->conn, sql)) {
fatal_error("f_mysql_speed: connection error %s", mysql_error(save->conn));
}
fclose(save->temp_fileptr);
remove(temp_pathname);
mysql_close(save->conn);
return 0;
}
/* processing phase */
save = (struct local_struct *) *local;
// if (new_GDS && save->isset == 0 ) {
if ((save->last_GDS_change_no != GDS_change_no) && save->isset == 0 ) {
//save->npts = GB2_Sec3_npts(sec);
save->npts = ndata;
if (ndata > MAX_NXNY) fatal_error_i("f_mysql_speed: MAX_NXNY exceeded %d", MAX_NXNY);
for (i = 0; i < ndata; i++) {
if (save->rows[i] == NULL) {
save->rows[i] = (char *) malloc(1500*sizeof(char));
if (save->rows[i] == NULL) fatal_error("f_mysql_speed: memory allocation problem","");
}
save->rows[i][0] = '\0';
}
save->isset = 1;
// } else if (new_GDS) {
} else if (save->last_GDS_change_no != GDS_change_no) {
fatal_error("f_mysql_speed, grid definition has to be the same for all fields","");
}
save->last_GDS_change_no = GDS_change_no;
/*Collect runtime and validtime into vt and rt*/
reftime(sec, &year, &month, &day, &hour, &minute, &second);
// p = sec[1];
// sprintf(rt, "%4.4d-%2.2d-%2.2d %2.2d:%2.2d:%2.2d", (p[12]<<8)+p[13], p[14],p[15],p[16],p[17],p[18]);
sprintf(rt, "%4.4d-%2.2d-%2.2d %2.2d:%2.2d:%2.2d", year,month,day,hour,minute,second);
if (verftime(sec, &year, &month, &day, &hour, &minute, &second) == 0) {
sprintf(vt,"%4.4d-%2.2d-%2.2d %2.2d:%2.2d:%2.2d", year,month,day,hour,minute,second);
}
/*Check that runtimes are equal and validtimes are equal*/
if (save->runtime[0] != 0 && strcmp(save->runtime,rt) != 0) {
fprintf(stderr, "Error, runtime has to be equal for all fields\n");
}
strcpy(save->runtime,rt);
if (save->validtime[0] != 0 && strcmp(save->validtime,vt) != 0) {
fprintf(stderr, "Error, validtime has to be equal for all fields\n");
}
strcpy(save->validtime,vt);
/*Get levels, parameter name, description and unit*/
// f_lev(mode, sec, data, ndata, level_buf, local);
f_lev(call_ARG0(level_buf, NULL));
//if (ndata != save->npts && save->npts>0) fprintf(stderr,"ERROR: fields do not contain equally many gridpoints, %d , %d \n",save->npts,ndata);
if (strcmp(level_buf, "reserved") == 0) return(0);
getName(sec, mode, NULL, name, desc, unit);
fprintf(stderr,"Start processing of %s at %s, runtime=%s, validtime=%s \n", name, level_buf,rt,vt);
ctr = GB2_Center(sec);
sprintf(sql,"select * from wgrib2_parameter_mapping where center_id=%d and wgrib_param_name='%s' and wgrib_param_levelname='%s'", ctr, name, level_buf);
fprintf(stderr,"SQL: %s\n", sql);
if (mysql_query(save->conn, sql)) {
fatal_error("f_mysql_speed: mysql error %s", mysql_error(save->conn));
}
sprintf(new_level,"%s", "None");
sprintf(conv,"%s", "None");
sprintf(precision, "%s", "2");
printf("\nCenter_id: %d \n", ctr);
res = mysql_use_result(save->conn);
while ( (row = mysql_fetch_row(res) ) != NULL) {
printf("our_name: %s \n", row[3]);
sprintf(new_name,"%s", row[3]);
sprintf(new_level,"%s", row[4]);
sprintf(conv,"%s", row[5]);
sprintf(precision, "%s", row[6]);
}
printf("val_precision: %s \n", precision);
printf("conversion: %s \n", conv);
printf("our_levelname: %s \n", new_level);
if (strcmp(new_level,"None") != 0) {
fprintf(stderr,"Sets level %s as level %s\n", level_buf, new_level);
sprintf(param, "%s_%s", new_name, new_level);
} else {
if (strcmp(level_buf,"surface") == 0) {
sprintf(param, "%s_0", name);
} else {
sscanf(level_buf,"%d %s", &level, &last);
sprintf(param, "%s_%d", name, level);
}
}
mysql_free_result(res);
strcat(save->params,",");
strcat(save->params,param);
fprintf(stderr,"Gridpoints in data: %d\n", ndata);
fprintf(stderr, "Gridpoints for insert: %d\n", save->npts);
fprintf(stderr, "Remove unlikely: %d\n", save->remove_unlikely);
fprintf(stderr, "Western longitudes: %d\n", save->wlon);
for (j = 0; j < ndata; j++) {
longitude = lon[j];
if (longitude > 180 && save->wlon==1) longitude-=360;
latitude = lat[j];
value = data[j];
if (save->remove_unlikely == 1 && value>0 && value<10e-8) value=0;
if (save->remove_unlikely == 1 && (strcmp(name,"APCP")==0 || strcmp(name,"ASNOW")==0 || strcmp(name,"ACPCP")==0 ) && value==1) value=0;
value = convunit(value, conv);
if (strlen(save->rows[j]) < 2) {
if (!UNDEFINED_VAL(data[j])) {
save->has_value[j] = 1;
sprintf(save->rows[j],"\"%s\",\"%s\",%g,%g,%lg",rt,vt,latitude,longitude,value);
} else {
sprintf(save->rows[j],"\"%s\",\"%s\",%g,%g,NULL",rt,vt,latitude,longitude);
}
} else {
if (!UNDEFINED_VAL(data[j])) {
save->has_value[j] = 1;
sprintf(sql,",%lg",value);
} else {
sprintf(sql,",NULL");
}
strcat(save->rows[j],sql);
}
}
return 0;
}
/*
* HEADER:100:fix_CFSv2_fcst:misc:3:fixes CFSv2 monthly fcst X=daily or 00/06/12/18 Y=pert no. Z=number ens fcsts v1.0
*/
int f_fix_CFSv2_fcst(ARG3) {
int subcenter, n_time_ranges, fcst_time_0, fcst_time_1, i;
int dtime, unit, id;
int ref_year, ref_month, ref_day, ref_hour, ref_minute, ref_second;
int fcst0_year, fcst0_month, fcst0_day, fcst0_hour, fcst0_minute, fcst0_second;
int fcst1_year, fcst1_month, fcst1_day, fcst1_hour, fcst1_minute, fcst1_second;
static unsigned char new_sec4[61];
struct local_struct {
int fixed; // number of fields processed
int dt; // dt either 6 or 24
int hour0; // initial hour
int pert; // perturbation no, 0 = ctl
int num_ensemble_members;
};
struct local_struct *save;
if (mode == -1) {
*local = save = (struct local_struct *)malloc( sizeof(struct local_struct));
save->fixed = 0;
save->dt = 24;
if (strcmp(arg1,"daily") == 0) {
save->dt = 6;
save->hour0 = 0;
}
else if (strcmp(arg1,"00") == 0) {
save->hour0 = 0;
}
else if (strcmp(arg1,"06") == 0) {
save->hour0 = 6;
}
else if (strcmp(arg1,"12") == 0) {
save->hour0 = 12;
}
else if (strcmp(arg1,"18") == 0) {
save->hour0 = 1;
}
else fatal_error("fix_CFSv2_fcst: bad arg %s, wanted daily,00,06,12 or 18",arg1);
save->pert = atoi(arg2);
save->num_ensemble_members = atoi(arg2);
}
save = *local;
if (mode == -2) sprintf(inv_out,"fix_CFSRv2_fcst %d fields fixed" , save->fixed);
if (mode < 0) return 0;
// only process NCEP CFSv2 monthly forecasts
// must be NCEP generated
if (GB2_Center(sec) != NCEP) return 0;
// subcenter should be 0, 1 (reanalysis), 3 (EMC) or 4 (NCO)
subcenter = GB2_Subcenter(sec);
if (subcenter != 0 && subcenter != 1 && subcenter != 3 && subcenter != 4) return 0;
// must have process id = 82 or 98
id = analysis_or_forecast_generating_process_identifier(sec);
if (id != 82 && id != 98) return 0;
// product defn table must be 8 (fcst average)
if ( code_table_4_0(sec) != 8) return 0;
// n_time_ranges should be 1
if ( (n_time_ranges = sec[4][41]) != 1) {
fprintf(stderr,"unexected CFSv2 type forecast field .. pdt=8, n=%d\n", n_time_ranges);
return 0;
}
// forecast time units should be months
if (sec[4][17] != 3 || sec[4][48] != 3) return 0;
fcst_time_0 = int4(sec[4]+18);
fcst_time_1 = int4(sec[4]+49) + fcst_time_0;
// get reference time
reftime(sec, &ref_year, &ref_month, &ref_day, &ref_hour, &ref_minute, &ref_second);
if (ref_minute != 0 || ref_second != 0) {
fprintf(stderr,"unexected CFSv2 minute/second value != 0\n");
return 0;
}
if (fcst_time_0 != 0) { // start at day=1 hour=save->hour0 of proper month
fcst0_year = ref_year;
fcst0_month = ref_month + fcst_time_0;
if (fcst0_month > 12) {
i = (fcst0_month - 1) / 12;
fcst0_year += i;
fcst0_month -= (i*12);
}
fcst0_day = 1;
fcst0_hour = save->hour0;
fcst0_minute = 0;
fcst0_second= 0;
}
else { // start at current month
fcst0_year = ref_year;
fcst0_month = ref_month;
fcst0_day = ref_day;
fcst0_hour = ref_hour;
fcst0_minute = ref_minute;
fcst0_second= ref_second;
if (save->dt == 24) {
i = save->hour0 - fcst0_hour;
if (i < 0) i += 24;
add_time(&fcst0_year, &fcst0_month, &fcst0_day, &fcst0_hour, &fcst0_minute, &fcst0_second, i, HOUR);
}
}
sub_time(fcst0_year, fcst0_month, fcst0_day, fcst0_hour, fcst0_minute, fcst0_second,
ref_year, ref_month, ref_day, ref_hour, ref_minute, ref_second,
&dtime, &unit);
// set forecast time.
if (dtime == 0) unit = MONTH;
sec[4][17] = unit;
int_char(dtime, sec[4]+18);
if (fcst_time_1 == 0) fatal_error("fix-CFSv2_fcst: unexpected end_ft","");
if (save->dt == 6) { // ends at 18Z of last day of month
fcst1_year = ref_year;
fcst1_month = ref_month + fcst_time_1-1;
if (fcst1_month > 12) {
i = (fcst1_month - 1) / 12;
fcst1_year += i;
fcst1_month -= (i*12);
}
fcst1_day = num_days_in_month(fcst1_year, fcst1_month);
fcst1_hour = 18;
fcst1_minute = 0;
fcst1_second= 0;
// time increment = 6 hours
sec[4][53] = (unsigned char) HOUR;
int_char(6, sec[4]+54);
}
else {
fcst1_year = ref_year;
fcst1_month = ref_month + fcst_time_1-1;
if (fcst1_month > 12) {
i = (fcst1_month - 1) / 12;
fcst1_year += i;
fcst1_month -= (i*12);
}
fcst1_day = num_days_in_month(fcst1_year, fcst1_month);
fcst1_hour = save->hour0;
fcst1_minute = 0;
fcst1_second= 0;
// time increment = 24 hours
sec[4][53] = (unsigned char) HOUR;
int_char(24, sec[4]+54);
}
// find stat processing time
sub_time(fcst1_year, fcst1_month, fcst1_day, fcst1_hour, fcst1_minute, fcst1_second,
fcst0_year, fcst0_month, fcst0_day, fcst0_hour, fcst0_minute, fcst0_second,
&dtime, &unit);
// change length of processing
sec[4][48] = (unsigned char) unit;
int_char(dtime, sec[4]+49);
// fprintf(stderr,"length of processing %d unit (%d)\n", dtime, unit);
// save end-of-processing time
save_time(fcst1_year, fcst1_month, fcst1_day, fcst1_hour, fcst1_minute, fcst1_second, sec[4]+34);
// make sure that some basic info is correct
sec[4][47] = 2; // fcst_time++
// now to add ensemble information
for (i = 0; i < 34; i++) new_sec4[i] = sec[4][i];
for (i = 34; i <= 57; i++) new_sec4[i+3] = sec[4][i];
uint_char(61, new_sec4); // length of new sec[4[
new_sec4[7] = 0; // pdt = 11
new_sec4[8] = 11;
// add perturbation info
if (save->pert == 0) { // pert == 0 means control forecast
new_sec4[34] = 0;
new_sec4[35] = 0;
}
else {
new_sec4[34] = 3;
new_sec4[35] = save->pert;
}
new_sec4[35] = save->num_ensemble_members;
sec[4] = new_sec4;
save->fixed = save->fixed + 1;
return 0;
}
/*
* HEADER:100:misc:inv:0:variable name qualifiers like chemical, ensemble, probability, etc
*/
int f_misc(ARG0) {
const char *string;
int need_space = 0;
int pdt, val;
static int error_count = 0;
if (mode < 0) return 0;
pdt = GB2_ProdDefTemplateNo(sec);
inv_out += strlen(inv_out);
f_ens(call_ARG0(inv_out,NULL) );
if (strlen(inv_out)) {
if (need_space) strcat(inv_out,":");
need_space = 1;
}
inv_out += strlen(inv_out);
f_prob(call_ARG0(inv_out,NULL));
if (strlen(inv_out)) {
if (need_space) strcat(inv_out,":");
need_space = 1;
}
inv_out += strlen(inv_out);
f_spatial_proc(call_ARG0(inv_out,NULL));
if (strlen(inv_out)) {
if (need_space) strcat(inv_out,":");
need_space = 1;
}
inv_out += strlen(inv_out);
f_wave_partition(call_ARG0(inv_out,NULL) );
if (strlen(inv_out)) {
if (need_space) strcat(inv_out,":");
need_space = 1;
}
inv_out += strlen(inv_out);
val = code_table_4_3(sec);
if (val == 5) {
if (need_space) strcat(inv_out,":");
strcat(inv_out,"probability forecast");
}
else if (val == 6 || val == 7) {
if (need_space) strcat(inv_out,":");
strcat(inv_out,"analysis/forecast error");
need_space = 1;
}
else if (val == 9) {
if (need_space) strcat(inv_out,":");
strcat(inv_out,"climatological");
need_space = 1;
}
else if (GB2_Center(sec) == 7 && val == 192) {
if (need_space) strcat(inv_out,":");
strcat(inv_out,"Confidence Indicator");
need_space = 1;
}
if (pdt == 7) {
if (need_space) strcat(inv_out,":");
strcat(inv_out,"analysis/forecast error");
need_space = 1;
}
else if (pdt == 6 || pdt == 10) {
if (need_space) strcat(inv_out,":");
f_percent(call_ARG0(inv_out,NULL) );
strcat(inv_out," level");
need_space = 1;
}
if ( (val = code_table_4_230(sec)) != -1) {
if (need_space) strcat(inv_out,":");
strcat(inv_out,"chemical=");
val = code_table_4_230(sec);
if (val >= 0) {
if (GB2_MasterTable(sec) <= 4) {
if (error_count++ <= 10) {
if (GB2_Center(sec) == ECMWF)
fprintf(stderr, "Warning: possible incompatible chemistry table .. table turned off.\n");
else if (GB2_Center(sec) != NCEP)
fprintf(stderr,
"Warning: if file made with ECMWF API, possible incompatible chemistry table\n");
}
}
string = NULL;
switch(val) {
#include "CodeTable_4.230.dat"
}
if (GB2_MasterTable(sec) <= 4 && GB2_Center(sec) == ECMWF) {
string = NULL;
}
if (string != NULL) strcat(inv_out,string);
else {
inv_out += strlen(inv_out);
sprintf(inv_out,"chemical_%d",val);
}
}
need_space = 1;
}
if ( (val = code_table_4_233(sec)) != -1) {
if (need_space) strcat(inv_out,":");
strcat(inv_out,"aerosol=");
val = code_table_4_233(sec);
if (val >= 0) {
string = NULL;
switch(val) {
#include "CodeTable_4.233.dat"
}
if (string != NULL) strcat(inv_out,string);
else {
inv_out += strlen(inv_out);
sprintf(inv_out,"chemical_%d",val);
}
}
need_space = 1;
}
if (pdt == 44) {
if (need_space) strcat(inv_out,":");
inv_out += strlen(inv_out);
f_aerosol_size(call_ARG0(inv_out,NULL));
need_space = 1;
}
if (pdt == 48) {
if (need_space) strcat(inv_out,":");
inv_out += strlen(inv_out);
f_aerosol_size(call_ARG0(inv_out,NULL));
strcat(inv_out,":");
inv_out += strlen(inv_out);
f_aerosol_wavelength(call_ARG0(inv_out,NULL));
need_space = 1;
}
return 0;
}
int main(int argc, char **argv) {
FILE *in;
unsigned char *msg, *sec[9];
long int last_pos;
struct stat stat_buf; /* for type of grib input file */
int file_arg, i, j, num_submsgs;
int n_arg;
unsigned int k, ndata;
float *data;
double *ddata, ref;
#ifdef USE_G2CLIB
float missing_c_val_1, missing_c_val_2;
g2int *bitmap, has_bitmap;
g2float *g2_data;
int ii;
#endif
struct ARGLIST arglist[N_ARGLIST];
int narglist = 0;
const char *new_argv[N_ARGLIST];
void *local[N_ARGLIST];
int has_inv_option, last_submsg;
int err, new_GDS, gdt, pdt, center;
unsigned char dscale[2] = {0,0};
inv_file = stdout;
// jas_init();
// gctpc initialiation
init(-1,-1,"gctpc_1,txt", "gctpc_2.txt");
data = NULL;
ndata = 0;
/* no arguments .. help screen */
if (argc == 1) {
// f_help(-1,NULL,NULL,0,inv_out,local,"most");
mode = -1;
data = NULL;
ndata = 0;
*inv_out = 0;
f_h(call_ARG0(inv_out,NULL));
fprintf(inv_file, "%s\n", inv_out);
eof_bin(); eof_string();
exit(8);
}
setup_user_gribtable();
/* copy argv */
for (i = 0; i < argc; i++) {
new_argv[i] = argv[i];
}
/* scan for "inv" and input file */
has_inv_option = 0;
file_arg = 0;
for (i = 1; i < argc; i++) {
if (new_argv[i][0] != '-') {
/* must be filename */
file_arg = i;
continue;
}
/* must be an option */
for (j = 0; j < nfunctions; j++) {
if (strcmp(&(new_argv[i][1]),functions[j].name) == 0) {
if (functions[j].type == inv) has_inv_option = 1;
i += functions[j].nargs;
break;
}
}
}
/* if no inv option, use default inventory .. put it at end */
if (has_inv_option == 0) {
for (i = 0; i < argc; i++) {
new_argv[i] = new_argv[i];
}
new_argv[argc++] = "-s";
}
/* parse parameters */
file_arg = 0;
for (i = 1; i < argc; i++) {
if (new_argv[i][0] != '-' || (strcmp(new_argv[i],"-") == 0) ) {
/* must be filename */
if (file_arg == 0) {
file_arg = i;
continue;
} else {
fatal_error("too many grib files .. 2nd=%s", new_argv[i]);
}
}
/* must be an option */
for (j = 0; j < nfunctions; j++) {
if (strcmp(&(new_argv[i][1]),functions[j].name) == 0) {
#ifdef DEBUG
fprintf(stderr,"match .. -%s %d args\n", functions[j].name, functions[j].nargs);
#endif
/* add to function argument list */
arglist[narglist].fn = j;
arglist[narglist].i_argc = i+1;
if (functions[j].type == inv) has_inv_option = 1;
i += functions[j].nargs;
if (i >= argc) fatal_error("missing arguments option=%s",functions[j].name);
narglist++;
break;
}
}
if (j == nfunctions) {
fatal_error("unknown option %s", new_argv[i]);
}
}
if (has_inv_option == 0) {
fatal_error("missing arguments on last option","");
}
/* initialize options mode = -1 */
#ifdef DEBUG
fprintf(stderr,"init options narglist %d\n",narglist);
#endif
for (j = 0; j < narglist; j++) {
inv_out[0] = 0;
n_arg = functions[arglist[j].fn].nargs;
err = 0;
if (n_arg == 0) err = functions[arglist[j].fn].fn(-1,NULL,NULL,0, inv_out,local+j);
else if (n_arg == 1) err = functions[arglist[j].fn].fn(-1,NULL,NULL,0, inv_out,local+j,
new_argv[arglist[j].i_argc]);
else if (n_arg == 2) err = functions[arglist[j].fn].fn(-1,NULL,NULL,0, inv_out,local+j,
new_argv[arglist[j].i_argc],new_argv[arglist[j].i_argc+1]);
else if (n_arg == 3) err = functions[arglist[j].fn].fn(-1,NULL,NULL,0, inv_out,local+j,
new_argv[arglist[j].i_argc],new_argv[arglist[j].i_argc+1],new_argv[arglist[j].i_argc+2]);
else if (n_arg == 4) err = functions[arglist[j].fn].fn(-1,NULL,NULL,0, inv_out,local+j,
new_argv[arglist[j].i_argc],new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2], new_argv[arglist[j].i_argc+3]);
else if (n_arg == 5) err = functions[arglist[j].fn].fn(-1,NULL,NULL,0, inv_out,local+j,
new_argv[arglist[j].i_argc],new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2], new_argv[arglist[j].i_argc+3],
new_argv[arglist[j].i_argc+4]);
else if (n_arg == 6) err = functions[arglist[j].fn].fn(-1,NULL,NULL,0, inv_out,local+j,
new_argv[arglist[j].i_argc],new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2], new_argv[arglist[j].i_argc+3],
new_argv[arglist[j].i_argc+4], new_argv[arglist[j].i_argc+5]);
else if (n_arg == 7) err = functions[arglist[j].fn].fn(-1,NULL,NULL,0, inv_out,local+j,
new_argv[arglist[j].i_argc],new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2], new_argv[arglist[j].i_argc+3],
new_argv[arglist[j].i_argc+4], new_argv[arglist[j].i_argc+5],
new_argv[arglist[j].i_argc+6]);
else if (n_arg == 8) err = functions[arglist[j].fn].fn(-1,NULL,NULL,0, inv_out,local+j,
new_argv[arglist[j].i_argc],new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2], new_argv[arglist[j].i_argc+3],
new_argv[arglist[j].i_argc+4], new_argv[arglist[j].i_argc+5],
new_argv[arglist[j].i_argc+6], new_argv[arglist[j].i_argc+7]);
if(inv_out[0] != 0) fprintf(inv_file, "%s", inv_out);
if (err) {
err_bin(); err_string();
exit(8);
}
}
if (file_arg == 0 && argc > 1) fatal_error("no input file", "");
if (latlon == 1 && output_order_wanted != wesn)
fatal_error("latitude-longitude information is only available with -order we:sn","");
/* open input file */
seq_input = 0;
if (strcmp(new_argv[file_arg],"-") == 0) {
seq_input = 1;
in = stdin;
if (mode == 98) fprintf(stderr, "grib input is stdin\n");
}
else {
if (stat(new_argv[file_arg], &stat_buf) != -1) {
if (S_ISREG(stat_buf.st_mode)) {
if (mode == 98) fprintf(stderr, "grib input is a regular file\n");
}
else if (S_ISDIR(stat_buf.st_mode)) {
fatal_error("grib input is a directory: %s",new_argv[file_arg]);
}
else if (S_ISCHR(stat_buf.st_mode)) {
seq_input = 1;
if (mode == 98) fprintf(stderr, "grib input is a char device\n");
}
else if (S_ISBLK(stat_buf.st_mode)) {
seq_input = 1;
if (mode == 98) fprintf(stderr, "grib input is a block device\n");
}
else if (S_ISFIFO(stat_buf.st_mode)) {
seq_input = 1;
if (mode == 98) fprintf(stderr, "grib input is a fifo device\n");
}
else {
if (mode == 98) fprintf(stderr, "grib input has an unknown type\n");
}
}
if ((in = fopen(new_argv[file_arg],"rb")) == NULL) {
fatal_error("could not open file: %s", new_argv[file_arg]);
}
}
/* sequential input - can not do random access */
if (seq_input && input == inv_mode) fatal_error("wgrib2 cannot random access grib input file","");
ndata = 0;
data = NULL;
ddata = NULL;
msg_no = 1;
inv_no = 0;
len = pos = 0;
submsg = 0;
msg = NULL;
if ((old_gds = (unsigned char *) malloc(GDS_max_size * sizeof(char)) ) == NULL) {
fatal_error("memory allocation problem old_gds in wgrib2.main","");
}
last_pos = -1;
last_submsg = -1;
/* if dump mode .. position io stream */
if (input == dump_mode) {
while (msg_no < dump_msg) {
msg = seq_input ? rd_grib2_msg_seq(in, &pos, &len, &num_submsgs) :
rd_grib2_msg(in, &pos, &len, &num_submsgs);
if (msg == NULL) fatal_error_i("record %d not found", dump_msg);
last_pos = pos;
pos += len;
msg_no++;
}
#ifdef DEBUG
printf("dump mode msg=%d\n", msg_no);
#endif
}
/*
* submsg = 0 .. beginning of unread record
* submsg = i .. start at ith submsg
* num_submsgs = number of submessages in grib message
*/
/* inventory loop */
for (;last_message == 0;) {
/* need position and submessage number of message */
if (input == inv_mode || input == dump_mode) {
if (input == inv_mode) {
if (rd_inventory(&msg_no,&submsg, &pos)) break;
}
else if (input == dump_mode) {
if (dump_msg == -1) break;
submsg = dump_submsg;
dump_msg = -1;
}
if (pos != last_pos) {
msg = seq_input ? rd_grib2_msg_seq(in, &pos, &len, &num_submsgs) :
rd_grib2_msg(in, &pos, &len, &num_submsgs);
if (msg == NULL) {
fatal_error_i("grib message #%d not found", msg_no);
break;
}
last_pos = pos;
last_submsg = -1;
}
if (pos == last_pos && submsg == last_submsg + 1) {
/* read previous submessage */
if (parse_next_msg(sec) != 0) {
fprintf(stderr,"\n*** grib message #%d.%d not found ***\n\n", msg_no, submsg);
break;
}
}
else {
/* need to get desired submessage into sec */
if (parse_1st_msg(sec) != 0) {
fprintf(stderr,"\n*** grib message #%d.1 not found ***\n\n", msg_no);
break;
}
for (i = 2; i <= submsg; i++) {
if (parse_next_msg(sec) != 0) {
fprintf(stderr,"\n*** grib message #%d.%d not found ***\n\n", msg_no, i);
break;
}
}
}
last_submsg = submsg;
}
else if (input == all_mode) {
if (submsg == 0) {
msg = seq_input ? rd_grib2_msg_seq(in, &pos, &len, &num_submsgs) :
rd_grib2_msg(in, &pos, &len, &num_submsgs);
if (msg == NULL) break;
submsg = 1;
}
else if (submsg > num_submsgs) {
pos += len;
msg_no++;
msg = seq_input ? rd_grib2_msg_seq(in, &pos, &len, &num_submsgs) :
rd_grib2_msg(in, &pos, &len, &num_submsgs);
if (msg == NULL) break;
submsg = 1;
}
if (submsg == 1) {
if (parse_1st_msg(sec) != 0) {
fprintf(stderr,"illegal format: parsing 1st submessage\n");
}
}
else {
if (parse_next_msg(sec) != 0) {
fprintf(stderr,"illegal format: parsing submessages\n");
}
}
}
if (only_submsg > 0 && only_submsg != submsg) {
submsg++;
continue;
}
if (for_mode) {
if (msg_no < for_start || msg_no > for_end || ((msg_no - for_start) % for_step) != 0) {
if (msg_no > for_end && input != inv_mode) last_message = 1;
submsg++;
continue;
}
}
#ifdef USE_REGEX
/* move inv_no++ before match_inv is made */
inv_no++;
if (match) {
inv_out[0] = 0;
if (num_submsgs > 1) {
sprintf(inv_out,"%d.%d:", msg_no, submsg);
}
else {
sprintf(inv_out,"%d:", msg_no);
}
// f_match_inv(0, sec, NULL, 0, inv_out+strlen(inv_out), NULL);
f_match_inv(call_ARG0(inv_out+strlen(inv_out), NULL));
if (is_match(inv_out) != 0) {
submsg++;
inv_no--;
continue;
}
}
#endif
match_flag = 0;
if (for_n_mode) {
if (inv_no < for_n_start || inv_no > for_n_end || ((inv_no - for_n_start) % for_n_step) != 0) {
if (inv_no > for_n_end) last_message = 1;
submsg++;
continue;
}
}
/* see if new GDS */
if ((i = GB2_Sec3_size(sec)) != old_GDS_size) {
new_GDS = 1;
}
else {
new_GDS = 0;
for (j = 0; j < i; j++) {
if (old_gds[j] != sec[3][j]) new_GDS = 1;
}
}
if (new_GDS) {
GDS_change_no++;
if (i > GDS_max_size) {
free(old_gds);
GDS_max_size = i;
if ((old_gds = (unsigned char *) malloc(GDS_max_size) ) == NULL) {
fatal_error("memory allocation problem old_gds in wgrib2.main","");
}
}
for (j = 0; j < i; j++) {
old_gds[j] = sec[3][j];
}
old_GDS_size = i;
/* update grid information */
get_nxny(sec, &nx, &ny, &npnts, &res, &scan); /* get nx, ny, and scan mode of grid */
output_order = (nx == -1 || ny == -1) ? raw : output_order_wanted;
if (latlon) {
i = 1;
if (use_gctpc && output_order == wesn) { /* use gctpc to get lat lon values */
i = gctpc_get_latlon(sec, &lon, &lat);
}
if (i) get_latlon(sec, &lon, &lat); /* get lat lon of grid points */
}
}
/* Decode NDFD WxText */
if (WxText) mk_WxKeys(sec);
// any fixes to raw grib message before decode need to be placed here
if (fix_ncep_2_flag) fix_ncep_2(sec);
if (fix_ncep_3_flag) fix_ncep_3(sec);
if (fix_ncep_4_flag) fix_ncep_4(sec);
#ifdef CHECK
j = code_table_5_0(sec); // type of compression
/* yes this can be simplified but want to split it up in case other decoders have problems */
if (j == 0 && sec[5][19] == 0 && int2(sec[5] + 17) != 0 && ieee2flt(sec[5]+11) != 0.0)
fprintf(stderr,"Warning: g2lib/g2clib/grib-api simple decode may differ from WMO standard, use -g2clib 0 for WMO standard\n");
if ((j == 2 || j == 3) && int2(sec[5]+17) != 0 && int4(sec[5] + 31) == 0 && ieee2flt(sec[5]+11) != 0.0)
fprintf(stderr,"Warning: g2lib/g2clib complex decode may differ from WMO standard, use -g2clib 0 for WMO standard\n");
if (j == 40 && sec[5][19] == 0 && int2(sec[5] + 17) != 0 && ieee2flt(sec[5]+11) != 0.0)
fprintf(stderr,"Warning: g2lib/g2clib jpeg deocde may differ from WMO standard, use use -g2clib 0 for WMO standard\n");
if (j == 41 && sec[5][19] == 0 && int2(sec[5] + 17) != 0 && ieee2flt(sec[5]+11) != 0.0)
fprintf(stderr,"Warning: g2lib/g2clib/grib-api png decode may differ from WMO standard, use -g2clib 0 for WMO standard\n");
#endif
#ifdef CHECK
/* check the size of Section 7 */
/* code to check the other sizes needs to be placed in decode routines */
j = code_table_5_0(sec); // type of compression
if (j == 0) { /* simple */
/* to avoid overflow on 32 bit machines */
/* old: k = (GB2_Sec5_nval(sec) * sec[5][19] + 7) / 8 + 5; */
k = 5 + (GB2_Sec5_nval(sec)/8) * sec[5][19] + (GB2_Sec5_nval(sec)%8) * (sec[5][19]/8)
+ ( (GB2_Sec5_nval(sec)%8) * (sec[5][19]%8) + 7) / 8;
if (k != GB2_Sec7_size(sec)) {
fprintf(stderr,"Detected a size mismatch, Section 7, wanted %d found %d\n", k, GB2_Sec7_size(sec));
if (decode) fatal_error("Section 7 size, mismatch, simple packing","");
}
}
else if (j == 4) { /* IEEE */
k = GB2_Sec5_nval(sec) * 4 + 5;
if (k != GB2_Sec7_size(sec)) {
fprintf(stderr,"Detected a size mismatch, Section 7, wanted %d found %d\n", k, GB2_Sec7_size(sec));
if (decode) fatal_error("Section 7 size, mismatch, IEEE packing","");
}
}
#endif
if (decode) {
#ifdef CHECK
if (code_table_6_0(sec) == 0) { // has bitmap
k = GB2_Sec3_npts(sec) - GB2_Sec5_nval(sec);
if (k != missing_points(sec[6]+6, GB2_Sec3_npts(sec)))
fatal_error_ii("inconsistent number of bitmap points sec3-sec5: %d sec6: %d",
k, missing_points(sec[6]+6, GB2_Sec3_npts(sec)));
}
else if (code_table_6_0(sec) == 255) { // no bitmap
if (GB2_Sec3_npts(sec) != GB2_Sec5_nval(sec))
fatal_error_ii("inconsistent number of data points sec3: %d sec5: %d",
GB2_Sec3_npts(sec), GB2_Sec5_nval(sec));
}
#endif
/* allocate data */
if (GB2_Sec3_npts(sec) != ndata) {
if (ndata) free(data);
ndata = GB2_Sec3_npts(sec);
if (ndata) {
data = (float *) malloc(ndata * sizeof(float));
if (data == NULL) fatal_error("main: memory allocation failed data","");
}
else { data = NULL; }
}
j = code_table_5_0(sec); // type of compression
gdt = code_table_3_1(sec); // grid type
pdt = GB2_ProdDefTemplateNo(sec); // product defintion template
/* USE G2CLIB */
#ifdef USE_G2CLIB
if (use_g2clib == 2) {
err = g2_getfld(msg,submsg,1,1,&grib_data);
if (err != 0) fatal_error_ii("Fatal g2clib decode err=%d msg=%d", err, msg_no);
free_gribfield = 1;
has_bitmap = grib_data->ibmap;
g2_data = &(grib_data->fld[0]);
if (has_bitmap == 0 || has_bitmap == 254) {
bitmap = grib_data->bmap;
for (i = 0; i < ndata; i++) {
data[i] = (bitmap[i] == 0) ? UNDEFINED : g2_data[i];
}
}
else {
for (i = 0; i < ndata; i++) {
data[i] = g2_data[i];
}
}
/* complex packing uses special values for undefined */
ii = sub_missing_values(sec, &missing_c_val_1, &missing_c_val_2);
if (ii == 1) {
for (i = 0; i < ndata; i++) {
if (data[i] == missing_c_val_1) data[i] = UNDEFINED;
}
}
else if (ii == 2) {
for (i = 0; i < ndata; i++) {
if (data[i] == missing_c_val_1) data[i] = UNDEFINED;
if (data[i] == missing_c_val_2) data[i] = UNDEFINED;
}
}
}
#endif
/* USE INTERNAL DECODER */
if (use_g2clib != 2) {
center = GB2_Center(sec);
if (use_g2clib == 1) { // introduce g2clib constant field error
/* g2clib ignores decimal scaling for constant fields make internal decoders look like g2clib*/
if ( (j == 0 && sec[5][19] == 0) || ((j == 2 || j == 3) && int4(sec[5] + 31) == 0) ||
(j == 40 && sec[5][19] == 0) || (j == 41 && sec[5][19] == 0) ||
(center == NCEP && j == 40000 && sec[5][19] == 0) ||
(center == NCEP && j == 40010 && sec[5][19] == 0) ) {
dscale[0] = sec[5][17];
dscale[1] = sec[5][18];
sec[5][17] = sec[5][18] = 0;
}
}
err = unpk_grib(sec, data);
if (err != 0) fatal_error_i("Fatal decode packing type %d",err);
if (use_g2clib == 1) { // fix up data
/* restore decimal scaling */
if ( (j == 0 && sec[5][19] == 0) || ((j == 2 || j == 3) && int4(sec[5] + 31) == 0) ||
(j == 40 && sec[5][19] == 0) || (j == 41 && sec[5][19]) ||
(center == NCEP && j == 40000 && sec[5][19] == 0) ||
(center == NCEP && j == 40010 && sec[5][19] == 0) ) {
sec[5][17] = dscale[0];
sec[5][18] = dscale[1];
}
}
}
/* convert to standard output order we:sn */
if (output_order_wanted == wesn) to_we_sn_scan(data);
else if (output_order_wanted == wens) to_we_ns_scan(data);
}
else {
if (ndata) free(data);
ndata = 0;
data = NULL;
}
/* get scaling parameters */
use_scale = scaling(sec, &ref, &dec_scale, &bin_scale, &i) == 0;
if (num_submsgs > 1) {
fprintf(inv_file, "%d.%d%s%ld", msg_no, submsg, ":", pos);
}
else {
fprintf(inv_file, "%d%s%ld", msg_no, ":", pos);
}
for (j = 0; j < narglist; j++) {
/* skip execution if match_flag == 1 */
/* an output option acts as endif for match_flag */
if (match_flag == 1) {
if (functions[arglist[j].fn].type == output) match_flag = 0;
continue;
}
// if (functions[arglist[j].fn].type == inv) fprintf(inv_file, item_deliminator);
if (functions[arglist[j].fn].type == inv) fprintf(inv_file, "%s", item_deliminator);
if (functions[arglist[j].fn].type != setup) {
inv_out[0] = 0;
n_arg = functions[arglist[j].fn].nargs;
if (n_arg == 0)
functions[arglist[j].fn].fn(mode, sec, data, ndata, inv_out, local+j);
else if (n_arg == 1)
functions[arglist[j].fn].fn(mode, sec, data, ndata, inv_out, local+j,
new_argv[arglist[j].i_argc]);
else if (n_arg == 2)
functions[arglist[j].fn].fn(mode, sec, data, ndata, inv_out, local+j,
new_argv[arglist[j].i_argc], new_argv[arglist[j].i_argc+1]);
else if (n_arg == 3)
functions[arglist[j].fn].fn(mode, sec, data, ndata, inv_out, local+j,
new_argv[arglist[j].i_argc], new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2]);
else if (n_arg == 4)
functions[arglist[j].fn].fn(mode, sec, data, ndata, inv_out, local+j,
new_argv[arglist[j].i_argc], new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2], new_argv[arglist[j].i_argc+3]);
else if (n_arg == 5)
functions[arglist[j].fn].fn(mode, sec, data, ndata, inv_out, local+j,
new_argv[arglist[j].i_argc], new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2], new_argv[arglist[j].i_argc+3],
new_argv[arglist[j].i_argc+4]);
else if (n_arg == 6)
functions[arglist[j].fn].fn(mode, sec, data, ndata, inv_out, local+j,
new_argv[arglist[j].i_argc], new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2], new_argv[arglist[j].i_argc+3],
new_argv[arglist[j].i_argc+4], new_argv[arglist[j].i_argc+5]);
else if (n_arg == 7)
functions[arglist[j].fn].fn(mode, sec, data, ndata, inv_out, local+j,
new_argv[arglist[j].i_argc], new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2], new_argv[arglist[j].i_argc+3],
new_argv[arglist[j].i_argc+4], new_argv[arglist[j].i_argc+5],
new_argv[arglist[j].i_argc+6]);
else if (n_arg == 8)
functions[arglist[j].fn].fn(mode, sec, data, ndata, inv_out, local+j,
new_argv[arglist[j].i_argc], new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2], new_argv[arglist[j].i_argc+3],
new_argv[arglist[j].i_argc+4], new_argv[arglist[j].i_argc+5],
new_argv[arglist[j].i_argc+6], new_argv[arglist[j].i_argc+7]);
if(inv_out[0] != 0) fprintf(inv_file, "%s", inv_out);
}
}
#ifdef CHECK
if (!decode) {
if (code_table_6_0(sec) == 0) { // has bitmap
k = GB2_Sec3_npts(sec) - GB2_Sec5_nval(sec);
if (k != missing_points(sec[6]+6, GB2_Sec3_npts(sec)))
fatal_error_ii("inconsistent number of bitmap points sec3-sec5: %d sec6: %d",
k, missing_points(sec[6]+6, GB2_Sec3_npts(sec)));
}
else if (code_table_6_0(sec) == 255) { // no bitmap
if (GB2_Sec3_npts(sec) != GB2_Sec5_nval(sec))
fatal_error_ii("inconsistent number of data points sec3: %d sec5: %d",
GB2_Sec3_npts(sec), GB2_Sec5_nval(sec));
}
}
#endif
submsg++;
#ifdef USE_G2CLIB
if (free_gribfield) { g2_free(grib_data); free_gribfield = 0;}
#endif
// fprintf(inv_file, "\n");
fprintf(inv_file, "%s",end_inv);
if (flush_mode) fflush(inv_file);
if (dump_msg > 0) break;
}
/* finalize all functions, call with mode = -2 */
for (j = 0; j < narglist; j++) {
if (functions[arglist[j].fn].type != setup) {
n_arg = functions[arglist[j].fn].nargs;
inv_out[0] = 0;
if (n_arg == 0)
functions[arglist[j].fn].fn(-2, NULL, NULL, 0, inv_out, local+j);
else if (n_arg == 1)
functions[arglist[j].fn].fn(-2, NULL, NULL, 0, inv_out, local+j,
new_argv[arglist[j].i_argc]);
else if (n_arg == 2)
functions[arglist[j].fn].fn(-2, NULL, NULL, 0, inv_out, local+j,
new_argv[arglist[j].i_argc], new_argv[arglist[j].i_argc+1]);
else if (n_arg == 3)
functions[arglist[j].fn].fn(-2, NULL, NULL, 0, inv_out, local+j,
new_argv[arglist[j].i_argc], new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2]);
else if (n_arg == 4)
functions[arglist[j].fn].fn(-2, NULL, NULL, 0, inv_out, local+j,
new_argv[arglist[j].i_argc], new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2], new_argv[arglist[j].i_argc+3]);
else if (n_arg == 5)
functions[arglist[j].fn].fn(-2, NULL, NULL, 0, inv_out, local+j,
new_argv[arglist[j].i_argc], new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2], new_argv[arglist[j].i_argc+3],
new_argv[arglist[j].i_argc+4]);
else if (n_arg == 6)
functions[arglist[j].fn].fn(-2, NULL, NULL, 0, inv_out, local+j,
new_argv[arglist[j].i_argc], new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2], new_argv[arglist[j].i_argc+3],
new_argv[arglist[j].i_argc+4], new_argv[arglist[j].i_argc+5]);
else if (n_arg == 7)
functions[arglist[j].fn].fn(-2, NULL, NULL, 0, inv_out, local+j,
new_argv[arglist[j].i_argc], new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2], new_argv[arglist[j].i_argc+3],
new_argv[arglist[j].i_argc+4], new_argv[arglist[j].i_argc+5],
new_argv[arglist[j].i_argc+6]);
else if (n_arg == 8)
functions[arglist[j].fn].fn(-2, NULL, NULL, 0, inv_out, local+j,
new_argv[arglist[j].i_argc], new_argv[arglist[j].i_argc+1],
new_argv[arglist[j].i_argc+2], new_argv[arglist[j].i_argc+3],
new_argv[arglist[j].i_argc+4], new_argv[arglist[j].i_argc+5],
new_argv[arglist[j].i_argc+6], new_argv[arglist[j].i_argc+7]);
// if (inv_out[0]) fprintf(stderr, "%s\n", inv_out);
if (inv_out[0]) fprintf(stderr, "%s%s", inv_out, end_inv);
}
}
eof_bin(); eof_string();
exit(0);
}
int getName(unsigned char **sec, int mode, char *inv_out, char *name, char *desc, char *unit) {
int discipline, center, mastertab, localtab, parmcat, parmnum;
int pdt;
struct gribtable_s *p;
const char *p_unit;
p = NULL;
if (user_gribtable != NULL) p = search_gribtable(user_gribtable, sec);
#ifdef USE_TIGGE
if (tigge && p == NULL) p = search_gribtable(tigge_gribtable, sec); /* tigge is default table */
#endif
if (p == NULL) p = search_gribtable(gribtable, sec);
#ifdef USE_TIGGE
/* if undefined and a tigge file */
if (p == NULL && !tigge && (code_table_1_3(sec) == 4 || code_table_1_3(sec) == 5)) p = search_gribtable(tigge_gribtable, sec);
#endif
p_unit = "unit";
if (p) {
p_unit = p->unit;
pdt = code_table_4_0(sec);
if (pdt == 5 || pdt == 9) p_unit = "prob";
}
if (p) {
if (name) strcpy(name, p->name);
if (desc) strcpy(desc, p->desc);
if (unit) strcpy(unit, p_unit);
if (inv_out) {
sprintf(inv_out, "%s", p->name);
inv_out += strlen(inv_out);
if (mode) sprintf(inv_out," %s [%s]", p->desc, p_unit);
}
}
else {
discipline = GB2_Discipline(sec);
center = GB2_Center(sec);
mastertab = GB2_MasterTable(sec);
localtab = GB2_LocalTable(sec);
parmcat = GB2_ParmCat(sec);
parmnum = GB2_ParmNum(sec);
if (name) sprintf(name,"var%d_%d_%d",discipline,parmcat,parmnum);
if (desc) strcpy(desc,"desc");
if (unit) strcpy(unit,p_unit);
if (inv_out) {
if (parmnum < 192 && parmcat < 192) {
sprintf(inv_out,"var discipline=%d master_table=%d parmcat=%d parm=%d",
discipline, mastertab, parmcat, parmnum);
}
else {
sprintf(inv_out,"var discipline=%d center=%d local_table=%d parmcat=%d parm=%d",
discipline, center, localtab, parmcat, parmnum);
}
}
}
return 0;
}
int f_ens(ARG0) {
int type, typefcst, pdt, center;
const char *string;
if (mode >= 0) {
typefcst = code_table_4_7(sec);
type = code_table_4_6(sec);
center = GB2_Center(sec);
pdt = GB2_ProdDefTemplateNo(sec);
if (type >= 0) {
switch(type) {
case 0: sprintf(inv_out,"ENS=hi-res ctl"); break;
case 1: sprintf(inv_out,"ENS=low-res ctl"); break;
case 2:
sprintf(inv_out,"ENS=-%d", perturbation_number(sec)); break;
case 3:
sprintf(inv_out,"ENS=+%d", perturbation_number(sec)); break;
case 4:
sprintf(inv_out,"MM-ENS=%d", perturbation_number(sec)); break;
default:
sprintf(inv_out,"ENS=? table4.6=%d pert=%d",type,perturbation_number(sec)); break;
}
inv_out += strlen(inv_out);
if (typefcst >= 0) {
*inv_out++=' ';
*inv_out=0;
}
}
if (typefcst >= 0) {
string = "unknown derived fcst";
switch(typefcst) {
case 0: string = "ens mean"; break;
case 1: string = "wt ens mean"; break;
case 2: if ((pdt == 2) || (pdt == 12)) {
string = "ens std dev"; break;
}
string = "cluster std dev"; break;
case 3: if ((pdt == 2) || (pdt == 12)) {
string = "normalized ens std dev"; break;
}
string = "normalized cluster std dev"; break;
case 4: string = "ens spread"; break;
case 5: string = "ens large anom index"; break;
case 6: if ((pdt == 2) || (pdt == 12)) {
string = "unwt ens mean"; break;
}
string = "unwt cluster mean"; break;
case 7: string = "25%-75% range"; break;
case 8: string = "min all members"; break;
case 9: string = "max all members"; break;
case 192: if (center == NCEP) string = "unwt mode all members"; break;
case 193: if (center == NCEP) string = "10% all members"; break;
case 194: if (center == NCEP) string = "50% all members"; break;
case 195: if (center == NCEP) string = "90% all members"; break;
case 196: if (center == NCEP) string = "stat. weight for each members"; break;
case 197: if (center == NCEP) string = "percentile from climate distribution"; break;
}
sprintf(inv_out,"%s", string);
inv_out += strlen(inv_out);
}
}
return 0;
}
/*
* HEADER:400:JMA:inv:0:inventory for JMA locally defined PDT
*/
int f_JMA(ARG0) {
int pdt, i, nr;
int center;
if (mode < 0) return 0;
center = GB2_Center(sec);
if (center != JMA1 && center != JMA2) return 0;
pdt = GB2_ProdDefTemplateNo(sec);
if (pdt < 50000) return 0;
if (pdt == 51022) {
sprintf(inv_out,"site=%c%c%c%c", sec[4][24], sec[4][25], sec[4][26], sec[4][27]);
inv_out += strlen(inv_out);
if (mode >= 1) {
sprintf(inv_out,":site_id#=%u", uint2(sec[4]+28));
inv_out += strlen(inv_out);
sprintf(inv_out,":site_lon=%lf:site_lat=%lf:site_elev=%.1lfm",int4(sec[4]+18)*1e-6,
int4(sec[4]+14)*1e-6,int2(sec[4]+22)*0.1) ;
inv_out += strlen(inv_out);
sprintf(inv_out,":mag_dec=%.2fdeg", int2(sec[4]+30)*0.001);
inv_out += strlen(inv_out);
sprintf(inv_out,":freq=%ukHz", uint4(sec[4]+32));
inv_out += strlen(inv_out);
sprintf(inv_out,":pol=%d:opn_mode=%d:reflec=%d", (int) sec[4][36], (int) sec[4][37], (int) sec[4][38]);
inv_out += strlen(inv_out);
sprintf(inv_out,":qc=%d:cluster_filter=%d", (int) sec[4][39], (int) sec[4][40]);
inv_out += strlen(inv_out);
sprintf(inv_out,":angle_elev_constant=%.2f", int2(sec[4]+41)*0.01);
inv_out += strlen(inv_out);
}
nr = GB2_Sec4_size(sec);
nr = (nr % 4) ? -1 : (nr - 60) / 4;
if (mode >= 1) {
sprintf(inv_out,":Ntheta=%d", nr);
inv_out += strlen(inv_out);
i = uint_n(sec[4]+55, 3);
if (i != 0xffffff) {
sprintf(inv_out,":bin_size=%dm", i);
inv_out += strlen(inv_out);
}
if ((sec[4][58] != 255) || (sec[4][59] != 255)) {
sprintf(inv_out,":d_theta=%.1lfdeg", int2(sec[4]+58)*0.1);
inv_out += strlen(inv_out);
}
}
if (mode < 2) {
sprintf(inv_out,":ele(1)=%.2fdeg", int2(sec[4]+60)*0.01);
inv_out += strlen(inv_out);
sprintf(inv_out,":ele(%d)=%.2fdeg", nr/2,int2(sec[4]+60+((nr-1)/2)*4)*0.01);
inv_out += strlen(inv_out);
}
if (mode >= 2) {
sprintf(inv_out,"\n:elevation angle(1..Ntheta)=");
inv_out += strlen(inv_out);
for (i = 0; i < nr; i++) {
sprintf(inv_out,"%.2f ", int2(sec[4]+60+i*4)*0.01);
inv_out += strlen(inv_out);
}
sprintf(inv_out,"\n:pulse repetion freq(1..Ntheta)=");
inv_out += strlen(inv_out);
for (i = 0; i < nr; i++) {
sprintf(inv_out,"%.2f ", int2(sec[4]+62+i*4)*0.01);
inv_out += strlen(inv_out);
}
}
}
return 0;
}
