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

ECC-445: Add count of global when ECCODES_DEBUG set

This commit is contained in:
Shahram Najm 2018-06-21 15:32:07 +01:00
parent 4220f3fffe
commit d4a4cd1406
1 changed files with 14 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
src/grib_iterator_class_gaussian_reduced.c
View File

@ -104,6 +104,14 @@ static int next(grib_iterator* i, double *lat, double *lon, double *val)
typedef void (*get_reduced_row_proc)(long pl, double lon_first, double lon_last, long* npoints, long* ilon_first, long* ilon_last); typedef void (*get_reduced_row_proc)(long pl, double lon_first, double lon_last, long* npoints, long* ilon_first, long* ilon_last);
static size_t count_global_points(long* pl, size_t plsize)
{
long i, count=0;
for (i=0;i<plsize;i++) {
count += pl[i];
}
return count;
}
static size_t count_subarea_points(grib_handle* h, get_reduced_row_proc get_reduced_row, static size_t count_subarea_points(grib_handle* h, get_reduced_row_proc get_reduced_row,
long* pl, size_t plsize, double lon_first, double lon_last) long* pl, size_t plsize, double lon_first, double lon_last)
{ {
@ -140,7 +148,7 @@ static int iterate_reduced_gaussian_subarea(grib_iterator* iter, grib_handle* h,
if (h->context->debug) { if (h->context->debug) {
const size_t np = count_subarea_points(h, get_reduced_row, pl, plsize, lon_first, lon_last); const size_t np = count_subarea_points(h, get_reduced_row, pl, plsize, lon_first, lon_last);
printf("ECCODES DEBUG iterate_reduced_gaussian_subarea: num points=%ld\n", np); printf("ECCODES DEBUG grib_iterator_class_gaussian_reduced: sub-area num points=%ld\n", np);
} }
/*find starting latitude */ /*find starting latitude */
@ -274,6 +282,11 @@ static int init(grib_iterator* iter, grib_handle* h, grib_arguments* args)
} else { } else {
/*global*/ /*global*/
iter->e=0; iter->e=0;
if (h->context->debug) {
const size_t np = count_global_points(pl, plsize);
printf("ECCODES DEBUG grib_iterator_class_gaussian_reduced: global num points=%ld\n", np);
}
for (j=0;j<plsize;j++) { for (j=0;j<plsize;j++) {
row_count=pl[j]; row_count=pl[j];
for (i=0;i<row_count;i++) { for (i=0;i<row_count;i++) {