Changeset 686 for XIOS/trunk/src/io/netCdfInterface.cpp

Timestamp:

09/15/15 15:15:57 (9 years ago)

Author:

rlacroix

Message:

Use the NetCDF wrapper for inputs for better error checking.

File:

• XIOS/trunk/src/io/netCdfInterface.cpp (modified) (29 diffs)

XIOS/trunk/src/io/netCdfInterface.cpp

@@ -22 +22 @@
 int CNetCdfInterface::create(const StdString& fileName, int cMode, int& ncId)
 {
-  int status = nc_create((fileName.c_str()), cMode, &ncId);
-  if (NC_NOERR != status)
-  {
-    StdString errormsg(nc_strerror(status));
-    StdStringStream sstr;
-    sstr << "Error in calling function : nc_create((fileName.c_str()), cMode, &ncId) " << std::endl
-         << errormsg << std::endl
-         << "Unable to create file, given its name : " << fileName
+  int status = nc_create(fileName.c_str(), cMode, &ncId);
+  if (NC_NOERR != status)
+  {
+    StdString errormsg(nc_strerror(status));
+    StdStringStream sstr;
+    sstr << "Error when calling function: nc_create(fileName.c_str(), cMode, &ncId) " << std::endl
+         << errormsg << std::endl
+         << "Unable to create file, given its name: " << fileName
          << "and its creation mode " << creationMode2String(cMode) << std::endl;
     StdString e = sstr.str();
@@ -49 +49 @@
 int CNetCdfInterface::createPar(const StdString& fileName, int cMode, MPI_Comm comm, MPI_Info info, int& ncId)
 {
-  int status = xios::nc_create_par((fileName.c_str()), cMode, comm, info, &ncId);
-  if (NC_NOERR != status)
-  {
-    StdString errormsg(nc_strerror(status));
-    StdStringStream sstr;
-    sstr << "Error in calling function : nc_create_par((fileName.c_str()), cMode, comm, info, &ncId) " << std::endl
-         << errormsg << std::endl
-         << "Unable to create file on parallel filesys, given its name : " << std::endl
+  int status = xios::nc_create_par(fileName.c_str(), cMode, comm, info, &ncId);
+  if (NC_NOERR != status)
+  {
+    StdString errormsg(nc_strerror(status));
+    StdStringStream sstr;
+    sstr << "Error when calling function: nc_create_par(fileName.c_str(), cMode, comm, info, &ncId) " << std::endl
+         << errormsg << std::endl
+         << "Unable to create file on parallel file system, given its name: " << std::endl
          << "and its creation mode " << creationMode2String(cMode) << std::endl;
     StdString e = sstr.str();
@@ -74 +74 @@
 int CNetCdfInterface::open(const StdString& fileName, int oMode, int& ncId)
 {
-  int status = nc_open((fileName.c_str()), oMode, &ncId);
-  if (NC_NOERR != status)
-  {
-    StdString errormsg(nc_strerror(status));
-    StdStringStream sstr;
-    sstr << "Error in calling function : nc_open((fileName.c_str()), oMode, &ncId) "<< std::endl
-         << errormsg << std::endl
-         << "Unable to open file, given its name : " << fileName
+  int status = nc_open(fileName.c_str(), oMode, &ncId);
+  if (NC_NOERR != status)
+  {
+    StdString errormsg(nc_strerror(status));
+    StdStringStream sstr;
+    sstr << "Error when calling function: nc_open(fileName.c_str(), oMode, &ncId) "<< std::endl
+         << errormsg << std::endl
+         << "Unable to open file, given its name: " << fileName
          << "and its open mode " << openMode2String(oMode) << std::endl;
     StdString e = sstr.str();
@@ -102 +102 @@
 int CNetCdfInterface::openPar(const StdString& fileName, int oMode, MPI_Comm comm, MPI_Info info, int& ncId)
 {
-  int status = xios::nc_open_par((fileName.c_str()), oMode, comm, info, &ncId);
-  if (NC_NOERR != status)
-  {
-    StdString errormsg(nc_strerror(status));
-    StdStringStream sstr;
-    sstr << "Error in calling function " << "nc_open_par((fileName.c_str()), oMode, comm, info, &ncId) " << std::endl
-         << errormsg << std::endl
-         << "Unable to open file on parallel filesys, given its name : " << fileName
+  int status = xios::nc_open_par(fileName.c_str(), oMode, comm, info, &ncId);
+  if (NC_NOERR != status)
+  {
+    StdString errormsg(nc_strerror(status));
+    StdStringStream sstr;
+    sstr << "Error when calling function nc_open_par(fileName.c_str(), oMode, comm, info, &ncId) " << std::endl
+         << errormsg << std::endl
+         << "Unable to open file on parallel file system, given its name: " << fileName
          << "and its open mode " << openMode2String(oMode) << std::endl;
     StdString e = sstr.str();
@@ -130 +130 @@
     StdString errormsg(nc_strerror(status));
     StdStringStream sstr;
-    sstr << "Error in calling function " << "nc_close(ncId)" << std::endl
-         << errormsg << std::endl
-         << "Unable to close file, given its id : " << ncId << std::endl;
+    sstr << "Error when calling function nc_close(ncId)" << std::endl
+         << errormsg << std::endl
+         << "Unable to close file, given its id: " << ncId << std::endl;
     StdString e = sstr.str();
     throw CNetCdfException(e);
@@ -152 +152 @@
     StdString errormsg(nc_strerror(status));
     StdStringStream sstr;
-    sstr << "Error in calling function " << "nc_redef(ncId)" << std::endl
+    sstr << "Error when calling function nc_redef(ncId)" << std::endl
       << errormsg << std::endl
-      << "Unable to put this file into define mode given its id : " << ncId << std::endl;
+      << "Unable to put this file into define mode given its id: " << ncId << std::endl;
     StdString e = sstr.str();
     throw CNetCdfException(e);
@@ -175 +175 @@
     StdStringStream sstr;
 
-    sstr << "Error in calling function " << "nc_enddef(ncId)" << std::endl
-         << errormsg << std::endl
-         << "Unable to end define mode of this file, given its id : " << ncId << std::endl;
+    sstr << "Error when calling function nc_enddef(ncId)" << std::endl
+         << errormsg << std::endl
+         << "Unable to end define mode of this file, given its id: " << ncId << std::endl;
     StdString e = sstr.str();
     throw CNetCdfException(e);
@@ -194 +194 @@
 int CNetCdfInterface::inqNcId(int ncid, const StdString& grpName, int& grpId)
 {
-  int status = nc_inq_ncid(ncid, (grpName.c_str()), &grpId);
-  if (NC_NOERR != status)
-  {
-    StdString errormsg(nc_strerror(status));
-    StdStringStream sstr;
-
-    sstr << "Error in calling function " << "nc_inq_ncid(ncid, (grpName.c_str()), &grpId)" << std::endl
-         << errormsg << std::endl
-         << "Unable to get id of a group (File), given its name : " << grpName << std::endl;
+  int status = nc_inq_ncid(ncid, grpName.c_str(), &grpId);
+  if (NC_NOERR != status)
+  {
+    StdString errormsg(nc_strerror(status));
+    StdStringStream sstr;
+
+    sstr << "Error when calling function nc_inq_ncid(ncid, grpName.c_str(), &grpId)" << std::endl
+         << errormsg << std::endl
+         << "Unable to get id of a group (File), given its name: " << grpName << std::endl;
     StdString e = sstr.str();
     throw CNetCdfException(e);
@@ -218 +218 @@
 \return Status code
 */
-int CNetCdfInterface::inqVarId(int ncid,const StdString& varName, int& varId)
-{
-  int status = nc_inq_varid(ncid, (varName.c_str()), &varId);
-  if (NC_NOERR != status)
-  {
-    StdString errormsg(nc_strerror(status));
-    StdStringStream sstr;
-
-    sstr << "Error in calling function : nc_inq_varid(ncid, (varName.c_str()), &varId)" << std::endl
+int CNetCdfInterface::inqVarId(int ncid, const StdString& varName, int& varId)
+{
+  int status = nc_inq_varid(ncid, varName.c_str(), &varId);
+  if (NC_NOERR != status)
+  {
+    StdString errormsg(nc_strerror(status));
+    StdStringStream sstr;
+
+    sstr << "Error when calling function: nc_inq_varid(ncid, varName.c_str(), &varId)" << std::endl
          << (errormsg) << std::endl
-         << "Unable to get id of variable with name : " << varName << std::endl;
+         << "Unable to get id of variable with name: " << varName << std::endl;
     StdString e = sstr.str();
     throw CNetCdfException(e);
@@ -243 +243 @@
 \return Status code
 */
-int CNetCdfInterface::inqDimId(int ncid,const StdString& dimName, int& dimId)
-{
-  int status = nc_inq_dimid(ncid, (dimName.c_str()), &dimId);
-  if (NC_NOERR != status)
-  {
-    StdString errormsg(nc_strerror(status));
-    StdStringStream sstr;
-
-    sstr << "Error in calling function " << "nc_inq_dimid(ncid, (dimName.c_str()), &dimId)" << std::endl
-         << errormsg << std::endl
-         << "Unable to get id of dimension, given its name : " << dimName << std::endl;
-    StdString e = sstr.str();
-    throw CNetCdfException(e);
-  }
-
+int CNetCdfInterface::inqDimId(int ncid, const StdString& dimName, int& dimId)
+{
+  int status = nc_inq_dimid(ncid, dimName.c_str(), &dimId);
+  if (NC_NOERR != status)
+  {
+    StdString errormsg(nc_strerror(status));
+    StdStringStream sstr;
+
+    sstr << "Error when calling function nc_inq_dimid(ncid, dimName.c_str(), &dimId)" << std::endl
+         << errormsg << std::endl
+         << "Unable to get id of dimension, given its name: " << dimName << std::endl;
+    StdString e = sstr.str();
+    throw CNetCdfException(e);
+  }
+
+  return status;
+}
+
+/*!
+This function queries the name of a variable given its id.
+\param [in] ncid Groupd id (or File Id)
+\param [in] varId Id of desired variable
+\param [out] varName name of desired variable
+\return Status code
+*/
+int CNetCdfInterface::inqVarName(int ncid, int varId, StdString& varName)
+{
+  char varNameBuff[NC_MAX_NAME + 1];
+  int status = nc_inq_varname(ncid, varId, varNameBuff);
+  if (NC_NOERR != status)
+  {
+    StdString errormsg(nc_strerror(status));
+    StdStringStream sstr;
+
+    sstr << "Error when calling function nc_inq_varname(ncid, varId, varNameBuff)" << std::endl
+         << errormsg << std::endl
+         << "Unable to get variable name given its id: " << varId << std::endl;
+    StdString e = sstr.str();
+    throw CNetCdfException(e);
+  }
+  varName = varNameBuff;
   return status;
 }
@@ -275 +301 @@
     StdStringStream sstr;
 
-    sstr << "Error in calling function " << "nc_inq_dimid" << std::endl
+    sstr << "Error when calling function nc_inq_dimid" << std::endl
       << errormsg << std::endl
       << "Unable to get id of unlimited dimension " << std::endl;
@@ -289 +315 @@
 \param [in] ncid Groupd id (or File Id)
 \param [in] dimId Id of desired dimension
-\param [in/out] dimName Name of desired dimension
+\param [out] dimName Name of desired dimension
 \return Status code
 */
 int CNetCdfInterface::inqDimName(int ncid, int dimId, StdString& dimName)
 {
-  char fullNameIn[NC_MAX_NAME +1];
+  char fullNameIn[NC_MAX_NAME + 1];
   int status = nc_inq_dimname(ncid, dimId, fullNameIn);
   if (NC_NOERR != status)
@@ -301 +327 @@
     StdStringStream sstr;
 
-    sstr << "Error in calling function " << "nc_inq_dimname(ncid, dimId, fullNameIn)" << std::endl
-         << errormsg << std::endl
-         << "Unable to get dimension name from its id : " << dimId << std::endl;
+    sstr << "Error when calling function nc_inq_dimname(ncid, dimId, fullNameIn)" << std::endl
+         << errormsg << std::endl
+         << "Unable to get dimension name given its id: " << dimId << std::endl;
     StdString e = sstr.str();
     throw CNetCdfException(e);
@@ -326 +352 @@
     StdStringStream sstr;
 
-    sstr << "Error in calling function " << "nc_inq_dimlen(ncid, dimId, &dimLen)" << std::endl
-         << errormsg << std::endl
-         << "Unable to get dimension length from its id: " << dimId << std::endl;
+    sstr << "Error when calling function nc_inq_dimlen(ncid, dimId, &dimLen)" << std::endl
+         << errormsg << std::endl
+         << "Unable to get dimension length given its id: " << dimId << std::endl;
     StdString e = sstr.str();
     throw CNetCdfException(e);
@@ -351 +377 @@
     StdStringStream sstr;
 
-    sstr << "Error in calling function " << "nc_inq_varndims(ncid, varId, &nDims)" << std::endl
-         << errormsg << std::endl
-         << "Unable to get the number of dimension of variable with Id : " << varId << std::endl;
+    sstr << "Error when calling function nc_inq_varndims(ncid, varId, &nDims)" << std::endl
+         << errormsg << std::endl
+         << "Unable to get the number of dimension of variable with Id: " << varId << std::endl;
     StdString e = sstr.str();
     throw CNetCdfException(e);
@@ -376 +402 @@
     StdStringStream sstr;
 
-    sstr << "Error in calling function " << "nc_inq_vardimid(ncid, varId, dimIds)" << std::endl
+    sstr << "Error when calling function nc_inq_vardimid(ncid, varId, dimIds)" << std::endl
          << errormsg << std::endl
          << "Unable to get list of dimension id of the variable with id " << varId << std::endl;
@@ -402 +428 @@
     StdStringStream sstr;
 
-    sstr << "Error in calling function " << "nc_inq_dimids(ncid, &nDims, dimIds, includeParents)" << std::endl;
-    sstr << errormsg << std::endl;
-    sstr << "Unable to retrieve number of dimension in the group with id : " << ncid << std::endl;
+    sstr << "Error when calling function nc_inq_dimids(ncid, &nDims, dimIds, includeParents)" << std::endl;
+    sstr << errormsg << std::endl;
+    sstr << "Unable to retrieve number of dimension in the group with id: " << ncid << std::endl;
     sstr << "With number of Parents " << includeParents << std::endl;
+    StdString e = sstr.str();
+    throw CNetCdfException(e);
+  }
+
+  return status;
+}
+
+/*!
+This function queries the full name of a group given its id.
+\param [in] ncid Groupd id (or File Id)
+\param [in/out] grpFullName the full name of the group
+\return Status code
+*/
+int CNetCdfInterface::inqGrpFullName(int ncid, StdString& grpFullName)
+{
+  StdSize strlen = 0;
+  std::vector<char> buff;
+  int status = nc_inq_grpname_full(ncid, &strlen, NULL);
+  if (NC_NOERR == status)
+  {
+    buff.resize(strlen + 1);
+    status = nc_inq_grpname_full(ncid, NULL, &buff[0]);
+  }
+
+  if (NC_NOERR != status)
+  {
+    StdString errormsg(nc_strerror(status));
+    StdStringStream sstr;
+
+    sstr << "Error when calling function nc_inq_grpname_full(ncid, &strlen, &buff[0])" << std::endl
+         << errormsg << std::endl
+         << "Unable to get the full group name given its id: " << ncid << std::endl;
+    StdString e = sstr.str();
+    throw CNetCdfException(e);
+  }
+
+  grpFullName.assign(buff.begin(), buff.end());
+
+  return status;
+}
+
+/*!
+This function queries the list of group ids given a location id.
+\param [in] ncid Groupd id (or File Id)
+\param [in/out] numgrps number of groups
+\param [in/out] ncids list of group ids
+\return Status code
+*/
+int CNetCdfInterface::inqGrpIds(int ncid, int& numgrps, int* ncids)
+{
+  int status = nc_inq_grps(ncid, &numgrps, ncids);
+  if (NC_NOERR != status)
+  {
+    StdString errormsg(nc_strerror(status));
+    StdStringStream sstr;
+
+    sstr << "Error when calling function nc_inq_grps(ncid, &numgrps, ncids)" << std::endl;
+    sstr << errormsg << std::endl;
+    sstr << "Unable to retrieve the list of groups for location id: " << ncid << std::endl;
+    StdString e = sstr.str();
+    throw CNetCdfException(e);
+  }
+
+  return status;
+}
+
+/*!
+This function queries the list of variable ids given a location id.
+\param [in] ncid Groupd id (or File Id)
+\param [in/out] nvars number of variables
+\param [in/out] varids list of variable ids
+\return Status code
+*/
+int CNetCdfInterface::inqVarIds(int ncid, int& nvars, int* varids)
+{
+  int status = nc_inq_varids(ncid, &nvars, varids);
+  if (NC_NOERR != status)
+  {
+    StdString errormsg(nc_strerror(status));
+    StdStringStream sstr;
+
+    sstr << "Error when calling function nc_inq_varids(ncid, &nvars, varids)" << std::endl;
+    sstr << errormsg << std::endl;
+    sstr << "Unable to retrieve the list of variables for location id: " << ncid << std::endl;
+    StdString e = sstr.str();
+    throw CNetCdfException(e);
+  }
+
+  return status;
+}
+
+/*!
+This function queries the type and the size of an attribute given its name and the id of the variable to which it is attached.
+\param [in] ncid Groupd id (or File Id)
+\param [in] varid the id of the variable to which the attribute is attached
+\param [in] name the name of the attribute
+\param [out] type the type of the attribute
+\param [out] len the size of the attribute
+\return Status code
+*/
+int CNetCdfInterface::inqAtt(int ncid, int varid, const StdString& name, nc_type& type, size_t& len)
+{
+  int status = nc_inq_att(ncid, varid, name.c_str(), &type, &len);
+  if (NC_NOERR != status)
+  {
+    StdString errormsg(nc_strerror(status));
+    StdStringStream sstr;
+
+    sstr << "Error when calling function nc_inq_att(ncid, varid, name.c_str(), &type, &len)" << std::endl;
+    sstr << errormsg << std::endl;
+    sstr << "Unable to query the attribute information given its name: " << name << " and its variable id:" << varid << std::endl;
+    StdString e = sstr.str();
+    throw CNetCdfException(e);
+  }
+
+  return status;
+}
+
+/*!
+This function queries the number of global attributes given a location id.
+\param [in] ncid Groupd id (or File Id)
+\param [out] ngatts the number of global attributes
+\return Status code
+*/
+int CNetCdfInterface::inqNAtts(int ncid, int& ngatts)
+{
+  int status = nc_inq_natts(ncid, &ngatts);
+  if (NC_NOERR != status)
+  {
+    StdString errormsg(nc_strerror(status));
+    StdStringStream sstr;
+
+    sstr << "Error when calling function nc_inq_natts(ncid, &ngatts)" << std::endl;
+    sstr << errormsg << std::endl;
+    sstr << "Unable to query the number of global attributes given the location id:" << ncid << std::endl;
+    StdString e = sstr.str();
+    throw CNetCdfException(e);
+  }
+
+  return status;
+}
+
+/*!
+This function queries the number of global attributes given a location id and a variable id.
+\param [in] ncid Groupd id (or File Id)
+\param [in] varid the id of the variable
+\param [out] natts the number of global attributes
+\return Status code
+*/
+int CNetCdfInterface::inqVarNAtts(int ncid, int varid, int& natts)
+{
+  int status = nc_inq_varnatts(ncid, varid, &natts);
+  if (NC_NOERR != status)
+  {
+    StdString errormsg(nc_strerror(status));
+    StdStringStream sstr;
+
+    sstr << "Error when calling function nc_inq_varnatts(ncid, varid, &natts)" << std::endl;
+    sstr << errormsg << std::endl;
+    sstr << "Unable to query the number of attributes given the location id:" << ncid << " and the variable id:" << varid << std::endl;
+    StdString e = sstr.str();
+    throw CNetCdfException(e);
+  }
+
+  return status;
+}
+
+
+//! Query the name of an attribute given a location id, a variable id and the attribute number
+int CNetCdfInterface::inqAttName(int ncid, int varid, int attnum, StdString& name)
+{
+  char attName[NC_MAX_NAME + 1];
+  int status = nc_inq_attname(ncid, varid, attnum, attName);
+  if (NC_NOERR != status)
+  {
+    StdString errormsg(nc_strerror(status));
+    StdStringStream sstr;
+
+    sstr << "Error when calling function nc_inq_attname(ncid, varid, attnum, attName)" << std::endl;
+    sstr << errormsg << std::endl;
+    sstr << "Unable to query the name of attribute " << attnum << " given the location id:" << ncid << " and the variable id:" << varid << std::endl;
     StdString e = sstr.str();
     throw CNetCdfException(e);
@@ -422 +629 @@
 int CNetCdfInterface::defGrp(int parentNcid, const StdString& grpName, int& grpId)
 {
-  int status = nc_def_grp(parentNcid, (grpName.c_str()), &grpId);
-  if (NC_NOERR != status)
-  {
-    StdString errormsg(nc_strerror(status));
-    StdStringStream sstr;
-
-    sstr << "Error in calling function " << "nc_def_grp(parentNcid, (grpName.c_str()), &grpId)" << std::endl;
-    sstr << errormsg << std::endl;
-    sstr << "Unable to create group Id, given its name : " << grpName << std::endl;
+  int status = nc_def_grp(parentNcid, grpName.c_str(), &grpId);
+  if (NC_NOERR != status)
+  {
+    StdString errormsg(nc_strerror(status));
+    StdStringStream sstr;
+
+    sstr << "Error when calling function nc_def_grp(parentNcid, grpName.c_str(), &grpId)" << std::endl;
+    sstr << errormsg << std::endl;
+    sstr << "Unable to create group Id, given its name: " << grpName << std::endl;
     StdString e = sstr.str();
     throw CNetCdfException(e);
@@ -447 +654 @@
 int CNetCdfInterface::defDim(int ncid, const StdString& dimName, StdSize dimLen, int& dimId)
 {
-  int status = nc_def_dim(ncid, (dimName.c_str()), dimLen, &dimId);
-  if (NC_NOERR != status)
-  {
-    StdString errormsg(nc_strerror(status));
-    StdStringStream sstr;
-
-    sstr << "Error in calling function " << "nc_def_dim(ncid, (dimName.c_str()), dimLen, &dimId)" << std::endl;
-    sstr << errormsg << std::endl;
-    sstr << "Unable to create dimension with name : " << dimName
+  int status = nc_def_dim(ncid, dimName.c_str(), dimLen, &dimId);
+  if (NC_NOERR != status)
+  {
+    StdString errormsg(nc_strerror(status));
+    StdStringStream sstr;
+
+    sstr << "Error when calling function nc_def_dim(ncid, dimName.c_str(), dimLen, &dimId)" << std::endl;
+    sstr << errormsg << std::endl;
+    sstr << "Unable to create dimension with name: " << dimName
          << " and with length " << dimLen << std::endl;
     StdString e = sstr.str();
@@ -475 +682 @@
 \return Status code
 */
-int CNetCdfInterface::defVar(int ncid,const StdString& varName, nc_type xtype,
-                            int nDims, const int dimIds[], int& varId)
-{
-  int status = nc_def_var(ncid, (varName.c_str()), xtype, nDims, dimIds, &varId);
-  if (NC_NOERR != status)
-  {
-    StdString errormsg(nc_strerror(status));
-    StdStringStream sstr;
-
-    sstr << "Error in calling function " << " nc_def_var(ncid, (varName.c_str()), xtype, nDims, dimIds, &varId)" << std::endl;
-    sstr << errormsg << std::endl;
-    sstr << "Unable to add a new variable with name : " << varName
+int CNetCdfInterface::defVar(int ncid, const StdString& varName, nc_type xtype,
+                             int nDims, const int dimIds[], int& varId)
+{
+  int status = nc_def_var(ncid, varName.c_str(), xtype, nDims, dimIds, &varId);
+  if (NC_NOERR != status)
+  {
+    StdString errormsg(nc_strerror(status));
+    StdStringStream sstr;
+
+    sstr << "Error when calling function  nc_def_var(ncid, varName.c_str(), xtype, nDims, dimIds, &varId)" << std::endl;
+    sstr << errormsg << std::endl;
+    sstr << "Unable to add a new variable with name: " << varName
          << " with type " << xtype
          << " and number of dimension " << nDims << std::endl;
@@ -501 +708 @@
 \param [in] ncid Id groupd(or File Id)
 \param [in] varId Id of the variable
-\param [in] storage Type of storage (It can be : NC_CONTIGUOUS, NC_CHUNKED)
+\param [in] storage Type of storage (It can be: NC_CONTIGUOUS, NC_CHUNKED)
 \param [in/out] chunkSize array list of chunk sizes
 \return Status code
@@ -513 +720 @@
     StdStringStream sstr;
 
-    sstr << "Error in calling function " << "nc_def_var_chunking(ncid, varId, storage, chunkSize)" << std::endl;
-    sstr << errormsg << std::endl;
-    sstr << "Unable to set chunk size of the variable with id : " << varId
+    sstr << "Error when calling function nc_def_var_chunking(ncid, varId, storage, chunkSize)" << std::endl;
+    sstr << errormsg << std::endl;
+    sstr << "Unable to set chunk size of the variable with id: " << varId
       << " and storage type " << storage << std::endl;
     StdString e = sstr.str();
@@ -539 +746 @@
     StdStringStream sstr;
 
-    sstr << "Error in calling function " << "nc_def_var_deflate(ncid, varId, false, (compressionLevel > 0), compressionLevel)" << std::endl;
+    sstr << "Error when calling function nc_def_var_deflate(ncid, varId, false, (compressionLevel > 0), compressionLevel)" << std::endl;
     sstr << errormsg << std::endl;
     sstr << "Unable to set the compression level of the variable with id: " << varId
@@ -559 +766 @@
 {
   int old_fill_mode;
-  int status = nc_set_fill(ncid, fill ? NC_FILL : NC_NOFILL, &old_fill_mode);
-  if (NC_NOERR != status)
-  {
-    StdString errormsg(nc_strerror(status));
-    StdStringStream sstr;
-
-    sstr << "Error when calling function nc_set_fill(ncid, fill ? NC_FILL : NC_NOFILL, &old_fill_mode)" << std::endl;
-    sstr << errormsg << std::endl;
-    sstr << "Unable to set the fill mode to : " << (fill ? "NC_FILL" : "NC_NOFILL") << std::endl;
+  int status = nc_set_fill(ncid, fill ? NC_FILL: NC_NOFILL, &old_fill_mode);
+  if (NC_NOERR != status)
+  {
+    StdString errormsg(nc_strerror(status));
+    StdStringStream sstr;
+
+    sstr << "Error when calling function nc_set_fill(ncid, fill ? NC_FILL: NC_NOFILL, &old_fill_mode)" << std::endl;
+    sstr << errormsg << std::endl;
+    sstr << "Unable to set the fill mode to: " << (fill ? "NC_FILL": "NC_NOFILL") << std::endl;
     StdString e = sstr.str();
     throw CNetCdfException(e);
@@ -592 +799 @@
     StdStringStream sstr;
 
-    sstr << "Error in calling function " << "nc_def_var_fill(ncid, varId, noFill, fillValue)" << std::endl;
-    sstr << errormsg << std::endl;
-    sstr << "Unable to set fill parameters of the variable with id : " << varId
+    sstr << "Error when calling function nc_def_var_fill(ncid, varId, noFill, fillValue)" << std::endl;
+    sstr << errormsg << std::endl;
+    sstr << "Unable to set fill parameters of the variable with id: " << varId
       << " and fill option " << noFill << std::endl;
     StdString e = sstr.str();
@@ -620 +827 @@
     StdStringStream sstr;
 
-    sstr << "Error in calling function " << "nc_var_par_access(ncid, varId, access)" << std::endl;
-    sstr << errormsg << std::endl;
-    sstr << "Unable to change read/write option of the variable with id : " << varId << std::endl;
+    sstr << "Error when calling function nc_var_par_access(ncid, varId, access)" << std::endl;
+    sstr << errormsg << std::endl;
+    sstr << "Unable to change read/write option of the variable with id: " << varId << std::endl;
     StdString e = sstr.str();
     throw CNetCdfException(e);
@@ -643 +850 @@
     StdStringStream sstr;
 
-    sstr << "Error in calling function " << "nc_sync(ncid)" << std::endl;
-    sstr << errormsg << std::endl;
-    sstr << "Unable to make a synchronization of a netCDF file with id : " << ncid << std::endl;
-    StdString e = sstr.str();
-    throw CNetCdfException(e);
-  }
-
-  return status;
-}
-
-/*!
-This function makes a request to netcdf with its id, to add or change a variable attribute or gloabl attribute,
-given its name, type, number of values provided for attribute
-\param [in] ncid Id of groupd(or File Id)
-\param [in] varId Id of the variable
-\param [in] attrName Name of the attribute
-\param [in] xtypes One of the set of predefined netCDF data types
-\param [in] numVal Number of values
-\param [in] op Array of values provided for attribute
-\return Status code
-*/
-int CNetCdfInterface::putAtt(int ncid, int varId, const StdString& attrName, nc_type xtype,
-                            StdSize numVal, const void* op)
-{
-  int status = nc_put_att(ncid, varId, (attrName.c_str()), xtype, numVal, op);
-  if (NC_NOERR != status)
-  {
-    StdString errormsg(nc_strerror(status));
-    StdStringStream sstr;
-
-    sstr << "Error in calling function " << "nc_put_att(ncid, varId, (attrName.c_str()), xtype, numVal, op)" << std::endl;
-    sstr << errormsg << std::endl;
-    sstr << "Unable to set attribute " << attrName << " for a variable with id : " << varId
-         << " with number of attribute  " << numVal
-         << " with type " << xtype << std::endl;
-    StdString e = sstr.str();
-    throw CNetCdfException(e);
-  }
-
-  return status;
-}
-
- // Some specilization of putAttributeType
+    sstr << "Error when calling function nc_sync(ncid)" << std::endl;
+    sstr << errormsg << std::endl;
+    sstr << "Unable to make a synchronization of a netCDF file with id: " << ncid << std::endl;
+    StdString e = sstr.str();
+    throw CNetCdfException(e);
+  }
+
+  return status;
+}
+
+// Some specializations of getAttributeType
+template<>
+int CNetCdfInterface::ncGetAttType(int ncid, int varid, const char* attrName, double* data)
+{
+  return nc_get_att_double(ncid, varid, attrName, data);
+}
+
+template<>
+int CNetCdfInterface::ncGetAttType(int ncid, int varid, const char* attrName, float* data)
+{
+  return nc_get_att_float(ncid, varid, attrName, data);
+}
+
+template<>
+int CNetCdfInterface::ncGetAttType(int ncid, int varid, const char* attrName, int* data)
+{
+  return nc_get_att_int(ncid, varid, attrName, data);
+}
+
+template<>
+int CNetCdfInterface::ncGetAttType(int ncid, int varid, const char* attrName, long* data)
+{
+  return nc_get_att_long(ncid, varid, attrName, data);
+}
+
+template<>
+int CNetCdfInterface::ncGetAttType(int ncid, int varid, const char* attrName, short* data)
+{
+  return nc_get_att_short(ncid, varid, attrName, data);
+}
+
+template<>
+int CNetCdfInterface::ncGetAttType(int ncid, int varid, const char* attrName, char* data)
+{
+  return nc_get_att_text(ncid, varid, attrName, data);
+}
+
+// Some specializations of putAttributeType
 template<>
 int CNetCdfInterface::ncPutAttType(int ncid, int varid, const char* attrName,
-                                 StdSize numVal, const double* op)
-{
-  return (nc_put_att_double(ncid, varid, attrName, NC_DOUBLE, numVal, op));
+                                   StdSize numVal, const double* data)
+{
+  return nc_put_att_double(ncid, varid, attrName, NC_DOUBLE, numVal, data);
 }
 
 template<>
 int CNetCdfInterface::ncPutAttType(int ncid, int varid, const char* attrName,
-                                   StdSize numVal, const float* op)
-{
-  return (nc_put_att_float(ncid, varid, attrName, NC_FLOAT, numVal, op));
+                                   StdSize numVal, const float* data)
+{
+  return nc_put_att_float(ncid, varid, attrName, NC_FLOAT, numVal, data);
 }
 
 template<>
 int CNetCdfInterface::ncPutAttType(int ncid, int varid, const char* attrName,
-                                   StdSize numVal, const int* op)
-{
-  return (nc_put_att_int(ncid, varid, attrName, NC_INT, numVal, op));
+                                   StdSize numVal, const int* data)
+{
+  return nc_put_att_int(ncid, varid, attrName, NC_INT, numVal, data);
 }
 
 template<>
 int CNetCdfInterface::ncPutAttType(int ncid, int varid, const char* attrName,
-                                   StdSize numVal, const long* op)
-{
-  return (nc_put_att_long(ncid, varid, attrName, NC_LONG, numVal, op));
+                                   StdSize numVal, const long* data)
+{
+  return nc_put_att_long(ncid, varid, attrName, NC_LONG, numVal, data);
 }
 
 template<>
 int CNetCdfInterface::ncPutAttType(int ncid, int varid, const char* attrName,
-                                   StdSize numVal, const short* op)
-{
-  return (nc_put_att_short(ncid, varid, attrName, NC_SHORT, numVal, op));
-}
-
-
-// Some specilization of putVariableType
-template<>
-int CNetCdfInterface::ncPutVaraType(int ncid, int varid, const StdSize* start, const StdSize* count, const double* op)
-{
-  return (nc_put_vara_double(ncid, varid, start, count, op));
-}
-
-template<>
-int CNetCdfInterface::ncPutVaraType(int ncid, int varid, const StdSize* start, const StdSize* count, const float* op)
-{
-  return (nc_put_vara_float(ncid, varid, start, count, op));
-}
-
-template<>
-int CNetCdfInterface::ncPutVaraType(int ncid, int varid, const StdSize* start, const StdSize* count, const int* op)
-{
-  return (nc_put_vara_int(ncid, varid, start, count, op));
+                                   StdSize numVal, const short* data)
+{
+  return nc_put_att_short(ncid, varid, attrName, NC_SHORT, numVal, data);
+}
+
+template<>
+int CNetCdfInterface::ncPutAttType(int ncid, int varid, const char* attrName,
+                                   StdSize numVal, const char* data)
+{
+  return nc_put_att_text(ncid, varid, attrName, numVal, data);
+}
+
+// Some specializations of getVariableType
+template<>
+int CNetCdfInterface::ncGetVaraType(int ncid, int varid, const StdSize* start, const StdSize* count, double* data)
+{
+  return nc_get_vara_double(ncid, varid, start, count, data);
+}
+
+template<>
+int CNetCdfInterface::ncGetVaraType(int ncid, int varid, const StdSize* start, const StdSize* count, float* data)
+{
+  return nc_get_vara_float(ncid, varid, start, count, data);
+}
+
+template<>
+int CNetCdfInterface::ncGetVaraType(int ncid, int varid, const StdSize* start, const StdSize* count, int* data)
+{
+  return nc_get_vara_int(ncid, varid, start, count, data);
+}
+
+// Some specializations of putVariableType
+template<>
+int CNetCdfInterface::ncPutVaraType(int ncid, int varid, const StdSize* start, const StdSize* count, const double* data)
+{
+  return nc_put_vara_double(ncid, varid, start, count, data);
+}
+
+template<>
+int CNetCdfInterface::ncPutVaraType(int ncid, int varid, const StdSize* start, const StdSize* count, const float* data)
+{
+  return nc_put_vara_float(ncid, varid, start, count, data);
+}
+
+template<>
+int CNetCdfInterface::ncPutVaraType(int ncid, int varid, const StdSize* start, const StdSize* count, const int* data)
+{
+  return nc_put_vara_int(ncid, varid, start, count, data);
 }
 
@@ -760 +997 @@
   {
   case NC_NOWRITE:
-    modeMes = StdString("NC_NOWRITE : Opening netCDF file with read-only access with buffering and caching access");
+    modeMes = StdString("NC_NOWRITE: Opening netCDF file with read-only access with buffering and caching access");
     break;
   case NC_SHARE:
-    modeMes = StdString("NC_SHARE : Several processes can read the file concurrently");
+    modeMes = StdString("NC_SHARE: Several processes can read the file concurrently");
     break;
   case NC_WRITE:
-    modeMes = StdString("NC_WRITE : NetCDF file is readable and writable");
+    modeMes = StdString("NC_WRITE: NetCDF file is readable and writable");
     break;
   default:
@@ -781 +1018 @@
   {
   case NC_NOCLOBBER:
-    modeMes = StdString("NC_NOCLOBBER : Not overwrite an exisiting netCDF file ");
+    modeMes = StdString("NC_NOCLOBBER: Not overwrite an exisiting netCDF file ");
     break;
   case NC_SHARE:
-    modeMes = StdString("NC_SHARE : Several processes can read from and write into the file concurrently");
+    modeMes = StdString("NC_SHARE: Several processes can read from and write into the file concurrently");
     break;
   case NC_64BIT_OFFSET:
-    modeMes = StdString("NC_64BIT_OFFSET : NetCDF file is 64-bit offset");
+    modeMes = StdString("NC_64BIT_OFFSET: NetCDF file is 64-bit offset");
     break;
   case NC_NETCDF4:
-    modeMes = StdString("NC_NETCDF4 : NetCDF file is HDF5/NetCDF-4");
+    modeMes = StdString("NC_NETCDF4: NetCDF file is HDF5/NetCDF-4");
     break;
   case NC_CLASSIC_MODEL:
-    modeMes = StdString("NC_CLASSIC_MODEL : NetCDF file is classical model");
+    modeMes = StdString("NC_CLASSIC_MODEL: NetCDF file is classical model");
     break;
   default: