ReadNCDF.hpp

Go to the documentation of this file.

//-------------------------------------------------------------------------
// Filename      : ReadNCDF.hpp
//
// Purpose       : ExodusII reader
//
// Special Notes : Lots of code taken from verde implementation
//
// Creator       : Tim Tautges & Corey Ernst
//
// Date          : 3/02
//
// Owner         : Tim Tautges & Corey Ernst
//-------------------------------------------------------------------------

#ifndef READNCDF_HPP
#define READNCDF_HPP

#ifndef IS_BUILDING_MB
//#error "ReadNCDF.hpp isn't supposed to be included into an application"
#endif

#include <vector>
#include <string>

#include "moab/Forward.hpp"
#include "moab/ReaderIface.hpp"
#include "moab/ExoIIInterface.hpp"
#include "moab/Range.hpp"

namespace moab {

class ReadUtilIface;

struct ReadBlockData
{
  int blockId;
  int startExoId; 
  EntityHandle startMBId; 
  int numElements;
  bool reading_in;
  ExoIIElementType elemType;
};

class ReadNCDF : public ReaderIface
{
   
public:
  
  static ReaderIface* factory( Interface* );
  
  void tokenize( const std::string& str,
                 std::vector<std::string>& tokens,
                 const char* delimiters );

  ErrorCode load_file( const char* file_name,
                       const EntityHandle* file_set,
                       const FileOptions& opts,
                       const SubsetList* subset_list = 0,
                       const Tag* file_id_tag = 0 );

  ErrorCode read_tag_values( const char* file_name,
                             const char* tag_name,
                             const FileOptions& opts,
                             std::vector<int>& tag_values_out,
                             const SubsetList* subset_list = 0 );
  
   ReadNCDF(Interface* impl = NULL);

  virtual ~ReadNCDF();

  //update the coords for deformed mesh according to FileOptions
  ErrorCode update(const char *exodus_file_name, const FileOptions& opts,
                     const int num_blocks, const int *blocks_to_load,
                     const EntityHandle file_set );

private:

  ReadUtilIface* readMeshIface;

  bool dimension_exists(const char *attrib_name);
  
  void reset();

  ErrorCode read_exodus_header();
  
  ErrorCode read_nodes(const Tag* file_id_tag);
  
  ErrorCode read_block_headers(const int *blocks_to_load,
                                  const int num_blocks);
  
  ErrorCode read_elements(const Tag* file_id_tag);
  
  ErrorCode read_global_ids();

  ErrorCode read_nodesets();
  
  ErrorCode read_sidesets();

  int exodus_file();

  int number_dimensions();

  ErrorCode get_type(char *exo_element_type,
                        EntityType &elem_type);
 
  ErrorCode get_type(EntityType &elem_type,
                        std::string &exo_element_type);

  /* 
  int get_int_tag(const MB_MeshSet *this_ms,
                  const TagHandle tag_id);
 */

  //qa record stuff 
  ErrorCode read_qa_records(EntityHandle file_set);
  ErrorCode read_qa_information( std::vector<std::string> &qa_record_list);

  ErrorCode read_qa_string(char *string,
                              int record_number,
                              int record_position); 

  ErrorCode create_ss_elements( int *element_ids, int *side_list,
                                   int num_sides, int num_dist_factors,
                                   std::vector<EntityHandle> &entities_to_add,
                                   std::vector<EntityHandle> &reverse_entities,
                                   std::vector<double> &dist_factor_vector,
                                   int ss_seq_id);

  ErrorCode find_side_element_type( const int element_id, ExoIIElementType &type, 
                                       ReadBlockData &block_data, int &df_index, int side_id );
  
 /* ErrorCode assign_block_ids_to_ssets(EntityHandle ss_handle,
                                         MB_MeshSet *ss_mesh_set);
                                         */

  ErrorCode create_sideset_element( const std::vector<EntityHandle>&, EntityType, EntityHandle&, int&);

  int get_number_nodes( EntityHandle handle );



  //------------member variables ------------//

  Interface* mdbImpl;
  
  int ncFile;        // netcdf/exodus file

  int CPU_WORD_SIZE;
  int IO_WORD_SIZE;

  int vertexOffset;

  std::string exodusFile;

  int numberNodes_loading;

  int numberElements_loading;

  int numberElementBlocks_loading; 

  int numberNodeSets_loading; 

  int numberSideSets_loading; 

  int numberDimensions_loading;

  EntityHandle mCurrentMeshHandle;

  //keeps track of the exodus ids of the elements and nodes just loaded
  std::vector<char> nodesInLoadedBlocks;
  //note- vector<bool> has limited capabilities

  //vector of blocks that are loading 
  std::vector< ReadBlockData > blocksLoading;

  Tag mMaterialSetTag;
  Tag mDirichletSetTag;
  Tag mNeumannSetTag;
  Tag mHasMidNodesTag;
  Tag mDistFactorTag;
  Tag mGlobalIdTag;
  Tag mQaRecordTag;

  int max_line_length, max_str_length;

  Range initRange;
};

// inline functions
inline int ReadNCDF::number_dimensions() 
{
   return numberDimensions_loading;
}

} // namespace moab

#endif