PC2ReadDomains.C

Go to the documentation of this file.

#ifndef __PLASCOM2_WRITE_DOMAINS_H__
#define __PLASCOM2_WRITE_DOMAINS_H__
#include "PlasCom2.H"
#include "PC2IO.H"

namespace plascom2 {

  int application::ReadDomains()
  {
    const std::string &inputFileName("ugh");
    std::ostringstream messageStream;
    

    int numDomains = appDomains.size();
    for(int iDomain = 0;iDomain < numDomains;iDomain++){
      std::string &domainName(domainInfo.domainNames[iDomain]);
      DomainBaseType &outDomain(appDomains[iDomain]);
      //      double domainTime = outDomain.Time();
      DomainBaseType::gridvector  &domainGrids(outDomain.Grids());
      DomainBaseType::statevector &gridStates(outDomain.States());
      DomainBaseType::statevector &gridParams(outDomain.Params());
      const std::vector<std::string> &geometryNames(outDomain.GeometryNames());
      const std::vector<std::string> &gridNames(outDomain.GridNames());
      int numGrids = domainGrids.size();
      for(int iGrid = 0;iGrid < numGrids;iGrid++){
        GridType  &outGrid(*domainGrids[iGrid]);
        StateType &outState(*gridStates[iGrid]);
        StateType &outParam(*gridParams[iGrid]);
        const std::string &geometryName(geometryNames[iGrid]);
        std::string gridName(gridNames[iGrid]);
        std::string::size_type x = gridName.find_last_of(":");
        if(x != std::string::npos){
          gridName = gridName.substr(x+1);
        }
#ifdef ENABLE_HDF5
        if(plascom2::io::hdf5::InputSingle(inputFileName,domainName,geometryName,gridName,
                                            outGrid,outState,outParam,appConfig,messageStream)){
          ErrOut("plascom2::io::hdf5::InputSingle failed.\n");
          return(1);
        }
#else
        ErrOut("WARNING!> HDF5 not enabled.  Cannot perform IO.\n");
        return(1);
#endif
      }
    }
    return(0);
  };
}

#endif