PC2WriteDomains.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 {

  std::string IORootName(){
    return(std::string("PlasCom2")); 
  };

  int application::WriteDomains()
  {
    std::ostringstream messageStream;
    std::string fileNameRoot;
    if(appConfig.IsSet("PlasCom2:OutputName"))
      fileNameRoot = appConfig.GetValue("PlasCom2:OutputName");
    if(fileNameRoot.empty()){
      fileNameRoot = std::string("PlasCom2");
    }
    std::string outputFileSuffix(std::string("_")+
                                 pcpp::io::StepStamp(iStep));
    std::string outputFileNameRoot(fileNameRoot+outputFileSuffix);
    
    //std::string("_")+pcpp::io::StepStamp(iStep));
    
    int numDomains = appDomains.size();
    
    for(int iDomain = 0;iDomain < numDomains;iDomain++){
      std::string &domainName(domainInfo.domainNames[iDomain]);
      DomainBaseType &outDomain(appDomains[iDomain]);

      std::ostringstream wsdMessageStream;
      //      std::string outputFileName2(outDomain.Name()+outputFileSuffix+".h5");
      std::string outputFileName2(outputFileNameRoot+".h5");
      plascom2::io::hdf5::WriteSingleDomain(outputFileName2,outDomain,appConfig,wsdMessageStream);
    }
    return(0);
    //      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 numLocalGrids = outDomain.NumberOfLocalGrids();
//       for(int iLocalGrid = 0;iLocalGrid < numLocalGrids;iLocalGrid++){
//         int iGrid = outDomain.LocalGridID(iLocalGrid);

//         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
//         std::string outputFileName(outputFileNameRoot+".h5");

//         if(plascom2::io::hdf5::OutputSingle(outputFileName,domainName,geometryName,gridName,
//                                             outGrid,outState,outParam,appConfig,messageStream,iStep)){
//           ErrOut("plascom2::io::hdf5::OutputSingle failed.\n");
//           return(1);
//         }
// #else
//         ErrOut("WARNING!> HDF5 not enabled.  Cannot perform IO.\n");
//         return(1);
// #endif
//       }
//     }
//     return(0);
  };
}

#endif