bes Updated for version 3.20.10
HDF5CFInt16.cc
Go to the documentation of this file.
1// This file is part of hdf5_handler: an HDF5 file handler for the OPeNDAP
2// data server.
3
4// Copyright (c) 2011-2016 The HDF Group, Inc. and OPeNDAP, Inc.
5//
6// This is free software; you can redistribute it and/or modify it under the
7// terms of the GNU Lesser General Public License as published by the Free
8// Software Foundation; either version 2.1 of the License, or (at your
9// option) any later version.
10//
11// This software is distributed in the hope that it will be useful, but
12// WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
13// or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
14// License for more details.
15//
16// You should have received a copy of the GNU Lesser General Public
17// License along with this library; if not, write to the Free Software
18// Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
19//
20// You can contact OPeNDAP, Inc. at PO Box 112, Saunderstown, RI. 02874-0112.
21// You can contact The HDF Group, Inc. at 1800 South Oak Street,
22// Suite 203, Champaign, IL 61820
23
32
33#include "config_hdf5.h"
34
35#include <libdap/InternalErr.h>
36#include "HDF5CFInt16.h"
37#include <BESDebug.h>
38#include "h5common.h"
39
40using namespace std;
41using namespace libdap;
42
43HDF5CFInt16::HDF5CFInt16(const string &n, const string &d) : Int16(n, d)
44{
45
46}
47
48HDF5CFInt16::HDF5CFInt16(const string &n, const string &d,const string &d_f) : Int16(n, d),filename(d_f)
49{
50}
51
52HDF5CFInt16::~HDF5CFInt16()
53{
54}
55BaseType *HDF5CFInt16::ptr_duplicate()
56{
57 return new HDF5CFInt16(*this);
58}
59
60bool HDF5CFInt16::read()
61{
62 BESDEBUG("h5","Coming to HDF5CFInt16 read "<<endl);
63
64 if (read_p())
65 return true;
66
67 hid_t file_id = H5Fopen(filename.c_str(),H5F_ACC_RDONLY,H5P_DEFAULT);
68 if(file_id < 0) {
69 throw InternalErr(__FILE__,__LINE__, "Fail to obtain the HDF5 file ID .");
70 }
71
72 hid_t dset_id = -1;
73 dset_id = H5Dopen2(file_id,dataset().c_str(),H5P_DEFAULT);
74 if(dset_id < 0) {
75 H5Fclose(file_id);
76 throw InternalErr(__FILE__,__LINE__, "Fail to obtain the dataset .");
77 }
78
79 hid_t dtypeid = H5Dget_type(dset_id);
80 if(dtypeid < 0) {
81 H5Dclose(dset_id);
82 H5Fclose(file_id);
83 throw InternalErr(__FILE__,__LINE__, "Fail to obtain the datatype .");
84 }
85
86 hid_t memtype = H5Tget_native_type(dtypeid, H5T_DIR_ASCEND);
87 if (memtype < 0){
88 H5Tclose(dtypeid);
89 H5Dclose(dset_id);
90 H5Fclose(file_id);
91 throw InternalErr(__FILE__, __LINE__, "Cannot obtain the native datatype.");
92 }
93
94 try {
95
96 dods_int16 buf;
97 if(1 == H5Tget_size(memtype) && H5T_SGN_2 == H5Tget_sign(memtype)) {
98 signed char buf2;
99 get_data(dset_id,(void*)&buf2);
100 buf =(short)buf2;
101 }
102 else
103 get_data(dset_id, (void *) &buf);
104
105 set_read_p(true);
106 set_value(buf);
107
108 if(H5Tclose(memtype) < 0) {
109 throw InternalErr(__FILE__, __LINE__, "Unable to close the memory datatype.");
110 }
111 if(H5Tclose(dtypeid) < 0) {
112 throw InternalErr(__FILE__, __LINE__, "Unable to close the datatype id.");
113 }
114 // Release the handles.
115 if (H5Dclose(dset_id) < 0) {
116 throw InternalErr(__FILE__, __LINE__, "Unable to close the dset.");
117 }
118 H5Fclose(file_id);
119 }
120 catch(...) {
121 H5Tclose(memtype);
122 H5Tclose(dtypeid);
123 H5Dclose(dset_id);
124 H5Fclose(file_id);
125 throw;
126 }
127
128 return true;
129
130}
131
This class provides a way to map HDF5 int16 to DAP int16 for the CF option.
void get_data(hid_t dset, void *buf)
Definition: h5common.cc:50