/* * (C) Copyright 2005- ECMWF. * * This software is licensed under the terms of the Apache Licence Version 2.0 * which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. * * In applying this licence, ECMWF does not waive the privileges and immunities granted to it by * virtue of its status as an intergovernmental organisation nor does it submit to any jurisdiction. */ /************************************** * Enrico Fucile **************************************/ #include "grib_api_internal.h" /* This is used by make_class.pl START_CLASS_DEF CLASS = accessor SUPER = grib_accessor_class_double IMPLEMENTS = unpack_double;pack_double;is_missing IMPLEMENTS = pack_long IMPLEMENTS = init MEMBERS=const char* value MEMBERS=const char* multiplier MEMBERS=const char* divisor MEMBERS=const char* truncating END_CLASS_DEF */ /* START_CLASS_IMP */ /* Don't edit anything between START_CLASS_IMP and END_CLASS_IMP Instead edit values between START_CLASS_DEF and END_CLASS_DEF or edit "accessor.class" and rerun ./make_class.pl */ static int is_missing(grib_accessor*); static int pack_double(grib_accessor*, const double* val, size_t* len); static int pack_long(grib_accessor*, const long* val, size_t* len); static int unpack_double(grib_accessor*, double* val, size_t* len); static void init(grib_accessor*, const long, grib_arguments*); typedef struct grib_accessor_scale { grib_accessor att; /* Members defined in gen */ /* Members defined in double */ /* Members defined in scale */ const char* value; const char* multiplier; const char* divisor; const char* truncating; } grib_accessor_scale; extern grib_accessor_class* grib_accessor_class_double; static grib_accessor_class _grib_accessor_class_scale = { &grib_accessor_class_double, /* super */ "scale", /* name */ sizeof(grib_accessor_scale), /* size */ 0, /* inited */ 0, /* init_class */ &init, /* init */ 0, /* post_init */ 0, /* destroy */ 0, /* dump */ 0, /* next_offset */ 0, /* get length of string */ 0, /* get number of values */ 0, /* get number of bytes */ 0, /* get offset to bytes */ 0, /* get native type */ 0, /* get sub_section */ 0, /* pack_missing */ &is_missing, /* is_missing */ &pack_long, /* pack_long */ 0, /* unpack_long */ &pack_double, /* pack_double */ 0, /* pack_float */ &unpack_double, /* unpack_double */ 0, /* unpack_float */ 0, /* pack_string */ 0, /* unpack_string */ 0, /* pack_string_array */ 0, /* unpack_string_array */ 0, /* pack_bytes */ 0, /* unpack_bytes */ 0, /* pack_expression */ 0, /* notify_change */ 0, /* update_size */ 0, /* preferred_size */ 0, /* resize */ 0, /* nearest_smaller_value */ 0, /* next accessor */ 0, /* compare vs. another accessor */ 0, /* unpack only ith value (double) */ 0, /* unpack only ith value (float) */ 0, /* unpack a given set of elements (double) */ 0, /* unpack a given set of elements (float) */ 0, /* unpack a subarray */ 0, /* clear */ 0, /* clone accessor */ }; grib_accessor_class* grib_accessor_class_scale = &_grib_accessor_class_scale; /* END_CLASS_IMP */ static void init(grib_accessor* a, const long l, grib_arguments* c) { grib_accessor_scale* self = (grib_accessor_scale*)a; int n = 0; self->value = grib_arguments_get_name(grib_handle_of_accessor(a), c, n++); self->multiplier = grib_arguments_get_name(grib_handle_of_accessor(a), c, n++); self->divisor = grib_arguments_get_name(grib_handle_of_accessor(a), c, n++); self->truncating = grib_arguments_get_name(grib_handle_of_accessor(a), c, n++); } static int unpack_double(grib_accessor* a, double* val, size_t* len) { grib_accessor_scale* self = (grib_accessor_scale*)a; int ret = 0; long value = 0; long multiplier = 0; long divisor = 0; if (*len < 1) { ret = GRIB_ARRAY_TOO_SMALL; grib_context_log(a->context, GRIB_LOG_ERROR, "Accessor %s cannot gather value for %s and/or %s", a->name, self->multiplier, self->divisor); return ret; } if ((ret = grib_get_long_internal(grib_handle_of_accessor(a), self->divisor, &divisor)) != GRIB_SUCCESS) return ret; if ((ret = grib_get_long_internal(grib_handle_of_accessor(a), self->multiplier, &multiplier)) != GRIB_SUCCESS) return ret; if ((ret = grib_get_long_internal(grib_handle_of_accessor(a), self->value, &value)) != GRIB_SUCCESS) return ret; if (value == GRIB_MISSING_LONG) *val = GRIB_MISSING_DOUBLE; else *val = ((double)(value * multiplier)) / divisor; if (ret == GRIB_SUCCESS) *len = 1; return ret; } static int pack_long(grib_accessor* a, const long* val, size_t* len) { const double dval = (double)*val; return pack_double(a, &dval, len); } static int pack_double(grib_accessor* a, const double* val, size_t* len) { grib_accessor_scale* self = (grib_accessor_scale*)a; int ret = 0; long value = 0; long divisor = 0; long multiplier = 0; long truncating = 0; double x; ret = grib_get_long_internal(grib_handle_of_accessor(a), self->divisor, &divisor); if (ret != GRIB_SUCCESS) return ret; ret = grib_get_long_internal(grib_handle_of_accessor(a), self->multiplier, &multiplier); if (ret != GRIB_SUCCESS) return ret; if (self->truncating) { ret = grib_get_long_internal(grib_handle_of_accessor(a), self->truncating, &truncating); if (ret != GRIB_SUCCESS) return ret; } if (multiplier == 0) { grib_context_log(a->context, GRIB_LOG_ERROR, "Accessor %s: cannot divide by a zero multiplier %s", a->name, self->multiplier); return GRIB_ENCODING_ERROR; } x = *val * (double)divisor / (double)multiplier; if (*val == GRIB_MISSING_DOUBLE) value = GRIB_MISSING_LONG; else if (truncating) { value = (long)x; } else { value = x > 0 ? (long)(x + 0.5) : (long)(x - 0.5); } ret = grib_set_long_internal(grib_handle_of_accessor(a), self->value, value); if (ret) grib_context_log(a->context, GRIB_LOG_ERROR, "Accessor %s: cannot pack value for %s (%s)\n", a->name, self->value, grib_get_error_message(ret)); if (ret == GRIB_SUCCESS) *len = 1; return ret; } static int is_missing(grib_accessor* a) { grib_accessor_scale* self = (grib_accessor_scale*)a; grib_accessor* av = grib_find_accessor(grib_handle_of_accessor(a), self->value); if (!av) return GRIB_NOT_FOUND; return grib_is_missing_internal(av); // int ret=0; // long value=0; // if((ret = grib_get_long_internal(grib_handle_of_accessor(a),self->value, &value))!= GRIB_SUCCESS){ // grib_context_log(a->context, GRIB_LOG_ERROR, // "Accessor %s cannot gather value for %s error %d \n", a->name, // self->value, ret); // return 0; // } // return (value == GRIB_MISSING_LONG); }