Method ParseHeader

IO/NetCDF/vtkNetCDFUGRIDReader.cxx:289–529

--------------------------------------------------------------------------------------------------

//--------------------------------------------------------------------------------------------------
bool vtkNetCDFUGRIDReader::ParseHeader()
{
  int varCount{};
  if (!this->CheckError(nc_inq_nvars(this->NcId, &varCount)))
  {
    return false;
  }

  if (varCount == 0)
  {
    vtkErrorMacro("No variable defined in file.");
    return false;
  }

  std::vector<int> vars{};
  vars.resize(static_cast<std::size_t>(varCount));
  if (!this->CheckError(nc_inq_varids(this->NcId, &varCount, vars.data())))
  {
    return false;
  }

  std::vector<int> meshIds{};
  std::vector<int> faceIds{};
  std::vector<int> nodeIds{};

  for (int var : vars)
  {
    int attCount{};
    if (!this->CheckError(nc_inq_varnatts(this->NcId, var, &attCount)))
    {
      return false;
    }

    for (int i = 0; i < attCount; ++i)
    {
      const auto name(this->GetAttributeName(var, i));

      nc_type type{ NC_NAT };
      if (!this->CheckError(nc_inq_atttype(this->NcId, var, name.c_str(), &type)))
      {
        return false;
      }

      if (type != NC_CHAR)
      {
        continue;
      }

      // if the cf_role attribute is "mesh_topology" then the var is a mesh
      if (name == "cf_role")
      {
        const auto value = this->GetAttributeString(var, name);
        if (value == "mesh_topology")
        {
          meshIds.emplace_back(var);
        }
      }
      else if (name == "location")