bes  Updated for version 3.20.6
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 "InternalErr.h"
36 #include "HDF5CFInt16.h"
37 #include <BESDebug.h>
38 #include "h5common.h"
39 
40 using namespace std;
41 using namespace libdap;
42 
43 HDF5CFInt16::HDF5CFInt16(const string &n, const string &d) : Int16(n, d)
44 {
45 
46 }
47 
48 HDF5CFInt16::HDF5CFInt16(const string &n, const string &d,const string &d_f) : Int16(n, d),filename(d_f)
49 {
50 }
51 
52 HDF5CFInt16::~HDF5CFInt16()
53 {
54 }
55 BaseType *HDF5CFInt16::ptr_duplicate()
56 {
57  return new HDF5CFInt16(*this);
58 }
59 
60 bool 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 
HDF5CFInt16.h
This class provides a way to map HDF5 int16 to DAP int16 for the CF option.
get_data
void get_data(hid_t dset, void *buf)
Definition: h5common.cc:50
HDF5CFInt16
Definition: HDF5CFInt16.h:40
libdap
Definition: BESDapFunctionResponseCache.h:35
h5common.h