32#include "config_hdf5.h"
37#include <libdap/InternalErr.h>
40#include "h5apicompatible.h"
45static int visit_obj_cb(hid_t group_id,
const char *name,
const H5O_info_t *oinfo,
void *_op_data);
46static herr_t attr_info(hid_t loc_id,
const char *name,
const H5A_info_t *ainfo,
void *_op_data);
48BaseType *HDF5GMCFMissLLArray::ptr_duplicate()
53bool HDF5GMCFMissLLArray::read()
56 BESDEBUG(
"h5",
"Coming to HDF5GMCFMissLLArray read "<<endl);
58 if (
nullptr == HDF5RequestHandler::get_lrdata_mem_cache())
59 read_data_NOT_from_mem_cache(
false,
nullptr);
63 vector<string> cur_lrd_non_cache_dir_list;
64 HDF5RequestHandler::get_lrd_non_cache_dir_list(cur_lrd_non_cache_dir_list);
69 if ((cur_lrd_non_cache_dir_list.empty())
70 || (
"" == check_str_sect_in_list(cur_lrd_non_cache_dir_list, filename,
'/'))) {
72 vector<string> cur_cache_dlist;
73 HDF5RequestHandler::get_lrd_cache_dir_list(cur_cache_dlist);
74 string cache_dir = check_str_sect_in_list(cur_cache_dlist, filename,
'/');
75 if (cache_dir !=
"") {
77 cache_key = cache_dir + varname;
80 cache_key = filename + varname;
84 vector<size_t> dim_sizes;
85 Dim_iter i_dim = dim_begin();
86 Dim_iter i_enddim = dim_end();
87 while (i_dim != i_enddim) {
88 dim_sizes.push_back(dimension_size(i_dim));
92 size_t total_elems = 1;
93 for (
unsigned int i = 0; i < dim_sizes.size(); i++)
94 total_elems = total_elems * dim_sizes[i];
96 handle_data_with_mem_cache(dtype, total_elems, cache_flag, cache_key,
false);
99 read_data_NOT_from_mem_cache(
false,
nullptr);
105void HDF5GMCFMissLLArray::obtain_aqu_obpg_l3_ll(
int* offset,
int* step,
int nelms,
bool add_cache,
void* buf)
108 BESDEBUG(
"h5",
"Coming to obtain_aqu_obpg_l3_ll read "<<endl);
114 throw InternalErr(__FILE__, __LINE__,
"The number of dimension for Aquarius Level 3 map data must be 1");
116 bool check_pass_fileid_key = HDF5RequestHandler::get_pass_fileid();
117 if (
false == check_pass_fileid_key) {
118 if ((fileid = H5Fopen(filename.c_str(), H5F_ACC_RDONLY, H5P_DEFAULT)) < 0) {
120 eherr <<
"HDF5 File " << filename <<
" cannot be opened. " << endl;
121 throw InternalErr(__FILE__, __LINE__, eherr.str());
126 if ((rootid = H5Gopen(fileid,
"/", H5P_DEFAULT)) < 0) {
127 HDF5CFUtil::close_fileid(fileid, check_pass_fileid_key);
129 eherr <<
"HDF5 dataset " << varname <<
" cannot be opened. " << endl;
130 throw InternalErr(__FILE__, __LINE__, eherr.str());
133 float LL_first_point = 0.0;
135 int LL_total_num = 0;
137 if (CV_LAT_MISS == cvartype) {
138 string Lat_SWP_name = (Aqu_L3 == product_type) ?
"SW Point Latitude" :
"sw_point_latitude";
139 string Lat_step_name = (Aqu_L3 == product_type) ?
"Latitude Step" :
"latitude_step";
140 string Num_lines_name = (Aqu_L3 == product_type) ?
"Number of Lines" :
"number_of_lines";
142 float Lat_step = 0.0;
144 vector<char> dummy_str;
146 obtain_ll_attr_value(fileid, rootid, Lat_SWP_name, Lat_SWP, dummy_str);
147 obtain_ll_attr_value(fileid, rootid, Lat_step_name, Lat_step, dummy_str);
148 obtain_ll_attr_value(fileid, rootid, Num_lines_name, Num_lines, dummy_str);
149 if (Num_lines <= 0) {
151 HDF5CFUtil::close_fileid(fileid, check_pass_fileid_key);
152 throw InternalErr(__FILE__, __LINE__,
"The number of line must be >0");
156 LL_first_point = Lat_SWP + (Num_lines - 1) * Lat_step;
157 LL_step = (float)(Lat_step * (-1.0));
158 LL_total_num = Num_lines;
161 if (CV_LON_MISS == cvartype) {
163 string Lon_SWP_name = (Aqu_L3 == product_type) ?
"SW Point Longitude" :
"sw_point_longitude";
164 string Lon_step_name = (Aqu_L3 == product_type) ?
"Longitude Step" :
"longitude_step";
165 string Num_columns_name = (Aqu_L3 == product_type) ?
"Number of Columns" :
"number_of_columns";
167 float Lon_step = 0.0;
170 vector<char> dummy_str_value;
172 obtain_ll_attr_value(fileid, rootid, Lon_SWP_name, Lon_SWP, dummy_str_value);
173 obtain_ll_attr_value(fileid, rootid, Lon_step_name, Lon_step, dummy_str_value);
174 obtain_ll_attr_value(fileid, rootid, Num_columns_name, Num_cols, dummy_str_value);
177 HDF5CFUtil::close_fileid(fileid, check_pass_fileid_key);
178 throw InternalErr(__FILE__, __LINE__,
"The number of line must be >0");
182 LL_first_point = Lon_SWP;
184 LL_total_num = Num_cols;
190 if (nelms > LL_total_num) {
192 HDF5CFUtil::close_fileid(fileid, check_pass_fileid_key);
193 throw InternalErr(__FILE__, __LINE__,
194 "The number of elements exceeds the total number of Latitude or Longitude");
197 for (
int i = 0; i < nelms; ++i)
198 val[i] = LL_first_point + (offset[0] + i * step[0]) * LL_step;
200 if (
true == add_cache) {
201 vector<float> total_val;
202 total_val.resize(LL_total_num);
203 for (
int total_i = 0; total_i < LL_total_num; total_i++)
204 total_val[total_i] = LL_first_point + total_i * LL_step;
205 memcpy(buf, total_val.data(), 4 * LL_total_num);
208 set_value((dods_float32 *) val.data(), nelms);
210 HDF5CFUtil::close_fileid(fileid, check_pass_fileid_key);
214void HDF5GMCFMissLLArray::obtain_gpm_l3_ll(
int* offset,
int* step,
int nelms,
bool add_cache,
void*buf)
218 throw InternalErr(__FILE__, __LINE__,
"The number of dimension for GPM Level 3 map data must be 1");
220 bool check_pass_fileid_key = HDF5RequestHandler::get_pass_fileid();
222 if (
false == check_pass_fileid_key) {
223 if ((fileid = H5Fopen(filename.c_str(), H5F_ACC_RDONLY, H5P_DEFAULT)) < 0) {
225 eherr <<
"HDF5 File " << filename <<
" cannot be opened. " << endl;
226 throw InternalErr(__FILE__, __LINE__, eherr.str());
230 vector<char> grid_info_value;
232 float lon_start = 0.;
240 if(GPMM_L3 == product_type || GPMS_L3 == product_type) {
241 hid_t grid_grp_id = -1;
242 string grid_grp_name;
244 if ((name() ==
"nlat") || (name() ==
"nlon")) {
246 string temp_grid_grp_name(GPM_GRID_GROUP_NAME1, strlen(GPM_GRID_GROUP_NAME1));
247 temp_grid_grp_name =
"/" + temp_grid_grp_name;
248 if (H5Lexists(fileid, temp_grid_grp_name.c_str(), H5P_DEFAULT) > 0)
249 grid_grp_name = temp_grid_grp_name;
251 string temp_grid_grp_name2(GPM_GRID_GROUP_NAME2, strlen(GPM_GRID_GROUP_NAME2));
252 temp_grid_grp_name2 =
"/" + temp_grid_grp_name2;
253 if (H5Lexists(fileid, temp_grid_grp_name2.c_str(), H5P_DEFAULT) > 0)
254 grid_grp_name = temp_grid_grp_name2;
256 throw InternalErr(__FILE__, __LINE__,
"Unknown GPM grid group name ");
262 string temp_grids_group_name(GPM_GRID_MULTI_GROUP_NAME, strlen(GPM_GRID_MULTI_GROUP_NAME));
263 if (name() ==
"lnH" || name() ==
"ltH")
264 grid_grp_name = temp_grids_group_name +
"/G2";
265 else if (name() ==
"lnL" || name() ==
"ltL") grid_grp_name = temp_grids_group_name +
"/G1";
272 if(name() ==
"lnH" || name() ==
"ltH" ||
273 name() ==
"lnL" || name() ==
"ltL") {
274 string temp_grids_group_name(GPM_GRID_MULTI_GROUP_NAME,strlen(GPM_GRID_MULTI_GROUP_NAME));
277 size_t grids_group_pos = varname.find(temp_grids_group_name);
278 if(string::npos == grids_group_pos) {
279 throw InternalErr (__FILE__, __LINE__,
280 "Cannot find group Grids.");
283 string grids_cgroup_path = varname.substr(grids_group_pos+1);
284 size_t grids_cgroup_pos = varname.find_first_of(
"/");
285 if(string::npos == grids_cgroup_pos) {
286 throw InternalErr (__FILE__, __LINE__,
287 "Cannot find child group of group Grids.");
290 string temp_sub_grp_name = grids_cgroup_path.substr(0,grids_cgroup_pos);
291 if(name() ==
"lnH" || name() ==
"ltH")
292 sub_grp1_name = temp_sub_grp_name;
293 else if(name() ==
"lnL" || name() ==
"ltL")
294 sub_grp2_name = temp_sub_grp_name;
296 grid_grp_name = temp_grids_group_name +
"/" + temp_sub_grp_name;
301 if ((grid_grp_id = H5Gopen(fileid, grid_grp_name.c_str(), H5P_DEFAULT)) < 0) {
302 HDF5CFUtil::close_fileid(fileid, check_pass_fileid_key);
304 eherr <<
"HDF5 dataset " << varname <<
" cannot be opened. " << endl;
305 throw InternalErr(__FILE__, __LINE__, eherr.str());
309 string grid_info_name(GPM_ATTR2_NAME, strlen(GPM_ATTR2_NAME));
310 if (name() ==
"lnL" || name() ==
"ltL")
311 grid_info_name =
"G1_" + grid_info_name;
312 else if (name() ==
"lnH" || name() ==
"ltH") grid_info_name =
"G2_" + grid_info_name;
314 float dummy_value = 0.0;
316 obtain_ll_attr_value(fileid, grid_grp_id, grid_info_name, dummy_value, grid_info_value);
317 HDF5CFUtil::parser_gpm_l3_gridheader(grid_info_value, latsize, lonsize, lat_start, lon_start, lat_res, lon_res,
320 H5Gclose(grid_grp_id);
323 H5Gclose(grid_grp_id);
332 vector<char> grid_info_value1;
333 vector<char> grid_info_value2;
334 obtain_gpm_l3_new_grid_info(fileid,grid_info_value1,grid_info_value2);
335 obtain_lat_lon_info(grid_info_value1,grid_info_value2,latsize,lonsize,lat_start,lon_start,lat_res,lon_res);
338 HDF5CFUtil::close_fileid(fileid, check_pass_fileid_key);
341 send_gpm_l3_ll_to_dap(latsize,lonsize,lat_start,lon_start,lat_res,lon_res,offset,step,nelms,add_cache, buf);
348 float lon_start = 0.;
355 HDF5CFUtil::parser_gpm_l3_gridheader(grid_info_value, latsize, lonsize, lat_start, lon_start, lat_res, lon_res,
358 if (0 == latsize || 0 == lonsize) {
359 HDF5CFUtil::close_fileid(fileid, check_pass_fileid_key);
360 throw InternalErr(__FILE__, __LINE__,
"Either latitude or longitude size is 0. ");
366 if (CV_LAT_MISS == cvartype) {
368 if (nelms > latsize) {
369 H5Gclose(grid_grp_id);
370 HDF5CFUtil::close_fileid(fileid, check_pass_fileid_key);
371 throw InternalErr(__FILE__, __LINE__,
"The number of elements exceeds the total number of Latitude ");
374 for (
int i = 0; i < nelms; ++i)
375 val[i] = lat_start + offset[0] * lat_res + lat_res / 2 + i * lat_res * step[0];
377 if (add_cache ==
true) {
378 vector<float> total_val;
379 total_val.resize(latsize);
380 for (
int total_i = 0; total_i < latsize; total_i++)
381 total_val[total_i] = lat_start + lat_res / 2 + total_i * lat_res;
382 memcpy(buf, total_val.data(), 4 * latsize);
385 else if (CV_LON_MISS == cvartype) {
387 if (nelms > lonsize) {
388 H5Gclose(grid_grp_id);
389 HDF5CFUtil::close_fileid(fileid, check_pass_fileid_key);
390 throw InternalErr(__FILE__, __LINE__,
"The number of elements exceeds the total number of Longitude");
394 for (
int i = 0; i < nelms; ++i)
395 val[i] = lon_start + offset[0] * lon_res + lon_res / 2 + i * lon_res * step[0];
397 if (add_cache ==
true) {
398 vector<float> total_val;
399 total_val.resize(lonsize);
400 for (
int total_i = 0; total_i < lonsize; total_i++)
401 total_val[total_i] = lon_start + lon_res / 2 + total_i * lon_res;
402 memcpy(buf, total_val.data(), 4 * lonsize);
407 set_value((dods_float32 *) val.data(), nelms);
409 H5Gclose(grid_grp_id);
410 HDF5CFUtil::close_fileid(fileid, check_pass_fileid_key);
418 if (nelms > LL_total_num) {
421 throw InternalErr (__FILE__, __LINE__,
422 "The number of elements exceeds the total number of Latitude or Longitude");
425 for (
int i = 0; i < nelms; ++i)
426 val[i] = LL_first_point + (offset[0] + i*step[0])*LL_step;
428 set_value ((dods_float32 *) val.data(), nelms);
439void HDF5GMCFMissLLArray::obtain_ll_attr_value(hid_t , hid_t s_root_id,
const string & s_attr_name,
440 T& attr_value, vector<char> & str_attr_value)
443 BESDEBUG(
"h5",
"Coming to obtain_ll_attr_value"<<endl);
444 hid_t s_attr_id = -1;
445 if ((s_attr_id = H5Aopen_by_name(s_root_id,
".", s_attr_name.c_str(),
446 H5P_DEFAULT, H5P_DEFAULT)) < 0) {
447 string msg =
"Cannot open the HDF5 attribute ";
450 throw InternalErr(__FILE__, __LINE__, msg);
453 hid_t attr_type = -1;
454 if ((attr_type = H5Aget_type(s_attr_id)) < 0) {
455 string msg =
"cannot get the attribute datatype for the attribute ";
459 throw InternalErr(__FILE__, __LINE__, msg);
462 hid_t attr_space = -1;
463 if ((attr_space = H5Aget_space(s_attr_id)) < 0) {
464 string msg =
"cannot get the hdf5 dataspace id for the attribute ";
469 throw InternalErr(__FILE__, __LINE__, msg);
472 hssize_t num_elm = H5Sget_simple_extent_npoints(attr_space);
475 string msg =
"cannot get the number for the attribute ";
479 H5Sclose(attr_space);
481 throw InternalErr(__FILE__, __LINE__, msg);
485 string msg =
"The number of attribute must be 1 for Aquarius level 3 data ";
489 H5Sclose(attr_space);
491 throw InternalErr(__FILE__, __LINE__, msg);
494 size_t atype_size = H5Tget_size(attr_type);
495 if (atype_size <= 0) {
496 string msg =
"cannot obtain the datatype size of the attribute ";
500 H5Sclose(attr_space);
502 throw InternalErr(__FILE__, __LINE__, msg);
505 if (H5T_STRING == H5Tget_class(attr_type)) {
506 if (H5Tis_variable_str(attr_type)) {
509 H5Sclose(attr_space);
511 throw InternalErr(__FILE__, __LINE__,
512 "Currently we assume the attributes we use to retrieve lat and lon are NOT variable length string.");
515 str_attr_value.resize(atype_size);
516 if (H5Aread(s_attr_id, attr_type, str_attr_value.data()) < 0) {
517 string msg =
"cannot retrieve the value of the attribute ";
521 H5Sclose(attr_space);
523 throw InternalErr(__FILE__, __LINE__, msg);
529 else if (H5Aread(s_attr_id, attr_type, &attr_value) < 0) {
530 string msg =
"cannot retrieve the value of the attribute ";
534 H5Sclose(attr_space);
536 throw InternalErr(__FILE__, __LINE__, msg);
541 H5Sclose(attr_space);
545void HDF5GMCFMissLLArray::obtain_gpm_l3_new_grid_info(hid_t file,
546 vector<char>& grid_info_value1,
547 vector<char>& grid_info_value2){
555 attr_na.name =
nullptr;
556 attr_na.value =
nullptr;
558 herr_t ret_o= H5OVISIT(file, H5_INDEX_NAME, H5_ITER_INC, visit_obj_cb, (
void*)&attr_na);
561 throw InternalErr(__FILE__, __LINE__,
"H5OVISIT failed. ");
571 vector<char> grid_info_value_1(attr_na.value,attr_na.value+strlen(attr_na.value));
572 vector<char> grid_info_value_2;
575 BESDEBUG(
"h5",
"Found the GPM level 3 Grid_info attribute."<<endl);
576 grid_info_value1.resize(strlen(attr_na.value));
577 memcpy(grid_info_value1.data(),attr_na.value,strlen(attr_na.value));
579 string tv(grid_info_value1.begin(),grid_info_value1.end());
580 cerr<<
"grid_info_value1 is "<<tv <<endl;
581 printf(
"attr_name 1st is %s\n",attr_na.name);
582 printf(
"attr_value 1st is %s\n",attr_na.value);
588 herr_t ret_o2= H5OVISIT(file, H5_INDEX_NAME, H5_ITER_INC, visit_obj_cb, (
void*)&attr_na);
591 throw InternalErr(__FILE__, __LINE__,
"H5OVISIT failed again. ");
605 grid_info_value2.resize(strlen(attr_na.value));
606 memcpy(grid_info_value2.data(),attr_na.value,strlen(attr_na.value));
608 string tv(grid_info_value2.begin(),grid_info_value2.end());
617void HDF5GMCFMissLLArray::obtain_lat_lon_info(
const vector<char>& grid_info_value1,
618 const vector<char>& grid_info_value2,
619 int& latsize,
int& lonsize,
620 float& lat_start,
float& lon_start,
621 float& lat_res,
float& lon_res){
623 float lat1_start = 0;
624 float lon1_start = 0.;
630 float lat2_start = 0;
631 float lon2_start = 0.;
637 HDF5CFUtil::parser_gpm_l3_gridheader(grid_info_value1, lat1size, lon1size, lat1_start, lon1_start,
638 lat1_res, lon1_res,
false);
640 HDF5CFUtil::parser_gpm_l3_gridheader(grid_info_value2, lat2size, lon2size, lat2_start, lon2_start,
641 lat2_res, lon2_res,
false);
643 bool pick_gv1 =
true;
646 if (name() ==
"lnL" || name() ==
"ltL") {
647 if(lat1_res <lat2_res)
650 else if (name() ==
"lnH" || name() ==
"ltH") {
651 if(lat1_res >lat2_res)
655 if(
true == pick_gv1) {
658 lat_start = lat1_start;
659 lon_start = lon1_start;
666 lat_start = lat2_start;
667 lon_start = lon2_start;
675attr_info(hid_t loc_id,
const char *name,
const H5A_info_t *ainfo,
void *_op_data)
684 attr_info_t *op_data = (attr_info_t *)_op_data;
687 if(strstr(name,GPM_ATTR2_NAME)!=
nullptr) {
690 attr = H5Aopen(loc_id, name, H5P_DEFAULT);
693 atype = H5Aget_type(attr);
698 if(H5T_STRING == H5Tget_class(atype)){
704 if(strncmp(name,op_data->name,strlen(name))!=0) {
705 hid_t aspace = H5Aget_space(attr);
711 hsize_t num_elms = H5Tget_size(atype)*H5Sget_simple_extent_npoints(aspace);
713 char *attr_value = op_data->value;
714 attr_value = malloc(num_elms+1);
715 H5Aread(attr,atype,attr_value);
716 printf(
"attr_value is %s\n",attr_value);
718 char *cur_attr_value = (
char*)malloc(num_elms+1);
719 if(H5Aread(attr,atype,(
void*)cur_attr_value)<0) {
723 free(cur_attr_value);
728 if(strncmp(cur_attr_value,op_data->value,strlen(op_data->value))!=0) {
730 op_data->name =
nullptr;
731 op_data->name = (
char*)malloc(strlen(name)+1);
732 strncpy(op_data->name,name,strlen(name));
734 free(op_data->value);
735 op_data->value =
nullptr;
736 op_data->value=(
char*)malloc(num_elms+1);
737 strncpy(op_data->value,cur_attr_value,strlen(cur_attr_value));
740 free(cur_attr_value);
745 hid_t aspace = H5Aget_space(attr);
752 hsize_t num_elms = H5Tget_size(atype)*H5Sget_simple_extent_npoints(aspace);
753 op_data->name = (
char*)malloc(strlen(name)+1);
754 strncpy(op_data->name,name,strlen(name));
757 char *attr_value = op_data->value;
758 attr_value = malloc(num_elms+1);
759 H5Aread(attr,atype,attr_value);
760 printf(
"attr_value is %s\n",attr_value);
762 op_data->value = (
char*)malloc(num_elms+1);
763 if(H5Aread(attr,atype,(
void*)op_data->value)<0) {
767 free(op_data->value);
786visit_obj_cb(hid_t group_id,
const char *name,
const H5O_info_t *oinfo,
799 attr_info_t *op_data = (attr_info_t *)_op_data;
802 if(oinfo->type == H5O_TYPE_GROUP) {
805 grp = H5Gopen2(group_id,name,H5P_DEFAULT);
808 ret = H5Aiterate2(grp, H5_INDEX_NAME, H5_ITER_INC,
nullptr, attr_info, op_data);
811 printf(
"object: attr name is %s\n",op_data->name);
812 printf(
"object: attr value is %s\n",op_data->value);
828void HDF5GMCFMissLLArray::send_gpm_l3_ll_to_dap(
const vector<char>& grid_info_value,
int* offset,
int* step,
829 int nelms,
bool add_cache,
void*buf) {
832 float lon_start = 0.;
840void HDF5GMCFMissLLArray::send_gpm_l3_ll_to_dap(
const int latsize,
const int lonsize,
const float lat_start,
const float lon_start,
841 const float lat_res,
const float lon_res,
const int* offset,
const int* step,
842 const int nelms,
const bool add_cache,
void*buf) {
845 if (0 == latsize || 0 == lonsize) {
846 throw InternalErr(__FILE__, __LINE__,
"Either latitude or longitude size is 0. ");
852 if (CV_LAT_MISS == cvartype) {
854 if (nelms > latsize) {
855 throw InternalErr(__FILE__, __LINE__,
"The number of elements exceeds the total number of Latitude ");
858 for (
int i = 0; i < nelms; ++i)
859 val[i] = lat_start + offset[0] * lat_res + lat_res / 2 + i * lat_res * step[0];
861 if (add_cache ==
true) {
862 vector<float> total_val;
863 total_val.resize(latsize);
864 for (
int total_i = 0; total_i < latsize; total_i++)
865 total_val[total_i] = lat_start + lat_res / 2 + total_i * lat_res;
866 memcpy(buf, total_val.data(), 4 * latsize);
869 else if (CV_LON_MISS == cvartype) {
875 if (nelms > lonsize) {
876 throw InternalErr(__FILE__, __LINE__,
"The number of elements exceeds the total number of Longitude");
879 for (
int i = 0; i < nelms; ++i)
880 val[i] = lon_start + offset[0] * lon_res + lon_res / 2 + i * lon_res * step[0];
882 if (add_cache ==
true) {
883 vector<float> total_val;
884 total_val.resize(lonsize);
885 for (
int total_i = 0; total_i < lonsize; total_i++)
886 total_val[total_i] = lon_start + lon_res / 2 + total_i * lon_res;
887 memcpy(buf, total_val.data(), 4 * lonsize);
892 set_value((dods_float32 *) val.data(), nelms);
896void HDF5GMCFMissLLArray::read_data_NOT_from_mem_cache(
bool add_cache,
void*buf)
899 BESDEBUG(
"h5",
"Coming to HDF5GMCFMissLLArray: read_data_NOT_from_mem_cache "<<endl);
911 int nelms = format_constraint(offset.data(), step.data(), count.data());
913 if (GPMM_L3 == product_type || GPMS_L3 == product_type || GPM_L3_New == product_type)
914 obtain_gpm_l3_ll(offset.data(), step.data(), nelms, add_cache, buf);
915 else if (Aqu_L3 == product_type || OBPG_L3 == product_type)
916 obtain_aqu_obpg_l3_ll(offset.data(), step.data(), nelms, add_cache, buf);
This class specifies the retrieval of the missing lat/lon values for general HDF5 products.
include the entry functions to execute the handlers