core
/
eccodes
mirror of https://github.com/ecmwf/eccodes.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

ECC-713 and ECC-880: grib_to_netcdf and stepUnits in minutes

This commit is contained in:
Shahram Najm 2019-02-01 18:12:55 +00:00
parent 1802a2acb9
commit 8c9c0eefd8
1 changed files with 21 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

25
tools/grib_to_netcdf.c
View File

@ -597,7 +597,6 @@ static err handle_to_request(request *r, grib_handle* g)
while(grib_keys_iterator_next(ks))
{
strcpy(name, grib_keys_iterator_get_name(ks));
if((e = grib_keys_iterator_get_string(ks, value, &len)) != GRIB_SUCCESS)
@ -607,6 +606,11 @@ static err handle_to_request(request *r, grib_handle* g)
len = sizeof(value);
}
strcpy(name, "stepUnits");
if((e = grib_get_string(g, name, value, &len)) == GRIB_SUCCESS) {
set_value(r, name, "%s", value);
}
/*
Assert(grib_get_long(g, "validityDate", &l ) == 0);
set_value(r, "validityDate", "%ld", l);
@ -1937,11 +1941,12 @@ static void validation_time(request *r)
{
long date = 0;
long time = 0;
long step = 0;
double step = 0;
long fcmonthdays = 0;
long fcmonth = 0;
double v;
long julian = 0;
const char* step_units = NULL;
long nstep = count_values(r, "step");
long ndate = count_values(r, "date");
@ -1965,7 +1970,7 @@ static void validation_time(request *r)
}
if(nstep)
step = atol(get_value(r, "step", 0));
step = atof(get_value(r, "step", 0));
if(ndate)
{
@ -2017,8 +2022,14 @@ static void validation_time(request *r)
}
julian = grib_date_to_julian(date);
step_units = get_value(r, "stepUnits", 0);
if (step_units){
if(strcmp("m", step_units)==0) {
step /= 60;
}
}
v = julian * 24.0 + fcmonthdays * 24.0 + time / 100.0 + step * 1.0;
grib_context_log(ctx, GRIB_LOG_DEBUG, "grib_to_netcdf: date=%ld, julian=%ld, fcmonthdays=%ld, time=%ld, step=%ld, validation=%ld", date, julian, fcmonthdays, time, step, v);
grib_context_log(ctx, GRIB_LOG_DEBUG, "grib_to_netcdf: date=%ld, julian=%ld, fcmonthdays=%ld, time=%ld, step=%g, validation=%ld", date, julian, fcmonthdays, time, step, v);
set_value(r, "_validation", "%lf", v);
set_value(r, "_juliandate", "%ld", julian);
@ -2932,6 +2943,12 @@ static int define_netcdf_dimensions(hypercube *h, fieldset *fs, int ncid, datase
grib_context_log(ctx, GRIB_LOG_ERROR, "Please split input data into different files. Exiting!\n");
exit(1);
}
if(count_values(data_r,"stepUnits") > 1)
{
grib_context_log(ctx, GRIB_LOG_ERROR, "Cannot handle fields for different stepUnits.\n");
grib_context_log(ctx, GRIB_LOG_ERROR, "Please split input data into different files. Exiting!\n");
exit(1);
}
if(strcmp(axis, "levelist") == 0)
{