Method operator()

IO/NetCDF/vtkNetCDFUGRIDReader.cxx:831–869

{
  template <typename OutArray>
  void operator()(OutArray* output, int NcId, int var, std::size_t time, std::size_t size,
    bool isTemporal, bool replaceFill, int& result)
  {
    using T = vtk::GetAPIType<OutArray>;

    output->SetNumberOfComponents(1);
    output->SetNumberOfTuples(size);

    if (isTemporal)
    {
      const std::array<std::size_t, 2> start{ time, 0 };
      const std::array<std::size_t, 2> count{ 1, size };
      result = nc_get_vara(NcId, var, start.data(), count.data(), output->GetPointer(0));
    }
    else
    {
      const std::array<std::size_t, 1> start{ 0 };
      const std::array<std::size_t, 1> count{ size };
      result = nc_get_vara(NcId, var, start.data(), count.data(), output->GetPointer(0));
    }

    if (result != NC_NOERR)
    {
      return;
    }

    if (replaceFill && (std::is_same<T, float>::value || std::is_same<T, double>::value))
    {
      T fillValue{};
      if (nc_get_att(NcId, var, "_FillValue", &fillValue) != NC_NOERR)
      {
        vtkDebugWithObjectMacro(output, "No fill value defined");
        return;
      }

      auto range(vtk::DataArrayValueRange(output));
      std::replace(range.begin(), range.end(), fillValue, static_cast<T>(vtkMath::Nan()));
    }
  }
};

//--------------------------------------------------------------------------------------------------