NcDataReader.h

Go to the documentation of this file.

// This file os part of FVM
// Copyright (c) 2012 FVM Authors
// See LICENSE file for terms.

#ifndef NCDATAREADER_H
#define NCDATAREADER_H

#include <string>
#include <vector>

#include "Mesh.h"


using namespace std;

class NcFile;
class NcDim;
class NcVar;


typedef   shared_ptr< StorageSite >  StorageSitePtr; 
typedef   shared_ptr< CRConnectivity > CRConnectivityPtr;


class NcDataReader {

public :
    typedef   shared_ptr< Array<int> >     ArrayIntPtr;

    NcDataReader( const string& fname );
    
    MeshList    getMeshList();

    static void  destroyMeshList( MeshList meshList );
    ~NcDataReader();

  void  createMappers( const MeshList& );

private :

    void  init();
    void  setNcFile();
    void  getDims();
    void  getVars();

    void  get_var_values();
    void  allocate_vars();
    void  get_bndry_vals();
    void  get_interface_vals();
    void  get_coord_vals();
    void  get_connectivity_vals();
    void  get_mapper_vals();

    MeshList  meshList();
    void  storage_sites ( int id, const MeshList&  meshList );
    void  boundary_faces( int id, const MeshList&  meshList );
    void  interfaces    ( int id, const MeshList&  meshList );
    void  coords        ( int id, const MeshList&  meshList );
    void  face_cells    ( int id, const MeshList&  meshList );
    void  face_nodes    ( int id, const MeshList&  meshList );

    string _fname;


    //netcdf variables    
    NcFile   *_ncFile;
    int  _nmesh;
    int  _nBoun;
    int  _charSize;
    int  _nNeighMesh;
    int  _nnodes;

    int   _nfaceRow;
    int   _nfaceCellsCol;
    int   _nfaceNodesCol;
    int   _nInterface;



    NcVar *_dimension;
    NcVar *_meshID;
    NcVar *_facesCount;
    NcVar *_cellsCount;
    NcVar *_ghostCellsCount;
    NcVar *_nodesCount;
    NcVar *_mapCount;
    NcVar *_interiorFacesGroup;

    NcVar*  _boundaryGroup;
    NcVar*  _boundarySize;
    NcVar*  _boundaryOffset;
    NcVar*  _boundaryID;
    NcVar*  _boundaryType;


    NcVar*  _interfaceGroup;
    NcVar*  _interfaceSize;
    NcVar*  _interfaceOffset;
    NcVar*  _interfaceID;
  
    NcVar*  _x;
    NcVar*  _y;
    NcVar*  _z;

    NcVar* _faceCellsRowCount;
    NcVar* _faceCellsColCount;
    NcVar* _faceNodesRowCount;
    NcVar* _faceNodesColCount;

    NcVar* _faceCellsRow;
    NcVar* _faceCellsCol;
    NcVar* _faceNodesRow;
    NcVar* _faceNodesCol;

    NcVar* _gatherIndices;
    NcVar* _scatterIndices;

    int*  _dimensionVals;
    int*  _meshIDVals;

    int *_facesCountVals;
    int *_cellsCountVals;
    int *_ghostCellsCountVals;
    int *_nodesCountVals;
    int *_mapCountVals;
    int *_interiorFacesGroupVals;


    int *_boundaryGroupVals;
    int *_boundarySizeVals;
    int *_boundaryOffsetVals;
    int * _boundaryIDVals;
    vector<char*>  _boundaryTypeVals;

    int  *_interfaceGroupVals;
    int  *_interfaceSizeVals;
    int  *_interfaceOffsetVals;
    int  *_interfaceIDVals;

    double  *_xVals;
    double  *_yVals;
    double  *_zVals;

    int  *_faceCellsRowCountVals;
    int  *_faceCellsColCountVals;
    int  *_faceNodesRowCountVals;
    int  *_faceNodesColCountVals;
    
    int  *_faceCellsRowVals;
    int  *_faceCellsColVals;
    int  *_faceNodesRowVals;
    int  *_faceNodesColVals;
    
    int  *_gatherIndicesVals;
    int  *_scatterIndicesVals;



};

#endif