#include #include #include #include #include #include #include #include using namespace Opm::EclIO; using EclEntry = EclFile::EclEntry; template void write(EclOutput& outFile, EclFile& file1, const std::string& name, int index) { auto vect = file1.get(index); outFile.write(name, vect); } template void write(EclOutput& outFile, ERst& file1, const std::string& name, int index, int reportStepNumber) { auto vect = file1.getRestartData(index, reportStepNumber); outFile.write(name, vect); } template void write(EclOutput& outFile, ERst& file1, const std::string& name, int index) { auto vect = file1.get(index); outFile.write(name, vect); } template void writeArray(std::string name, eclArrType arrType, T& file1, int index, EclOutput& outFile) { if (arrType == INTE) { write(outFile, file1, name, index); } else if (arrType == REAL) { write(outFile, file1, name,index); } else if (arrType == DOUB) { write(outFile, file1, name, index); } else if (arrType == LOGI) { write(outFile, file1, name, index); } else if (arrType == CHAR) { write(outFile, file1, name, index); } else if (arrType == C0NN) { write(outFile, file1, name, index); } else if (arrType == MESS) { outFile.message(name); } else { std::cout << "unknown array type " << std::endl; exit(1); } } template void writeArray(std::string name, eclArrType arrType, T& file1, int index, int reportStepNumber, EclOutput& outFile) { if (arrType == INTE) { write(outFile, file1, name, index, reportStepNumber); } else if (arrType == REAL) { write(outFile, file1, name, index, reportStepNumber); } else if (arrType == DOUB) { write(outFile, file1, name, index, reportStepNumber); } else if (arrType == LOGI) { write(outFile, file1, name, index, reportStepNumber); } else if (arrType == CHAR) { write(outFile, file1, name, index, reportStepNumber); } else if (arrType == MESS) { outFile.message(name); } else { std::cout << "unknown array type " << std::endl; exit(1); } } void writeC0nnArray(std::string name, int elementSize, EclFile& file1, int index, EclOutput& outFile) { auto vect = file1.get(index); outFile.write(name, vect, elementSize); } void writeArrayList(std::vector& arrayList, std::vector& elementSizeList, EclFile file1, EclOutput& outFile) { for (size_t index = 0; index < arrayList.size(); index++) { std::string name = std::get<0>(arrayList[index]); eclArrType arrType = std::get<1>(arrayList[index]); if (arrType == Opm::EclIO::C0NN){ writeC0nnArray(name, elementSizeList[index], file1, index, outFile); } else writeArray(name, arrType, file1, index, outFile); } } void writeArrayList(std::vector& arrayList, ERst file1, int reportStepNumber, EclOutput& outFile) { for (size_t index = 0; index < arrayList.size(); index++) { std::string name = std::get<0>(arrayList[index]); eclArrType arrType = std::get<1>(arrayList[index]); writeArray(name, arrType, file1, index , reportStepNumber, outFile); } } static void printHelp() { std::cout << "\nconvertECL needs one argument which is the input file to be converted. If this is a binary file the output file will be formatted. If the input file is formatted the output will be binary. \n" << "\nIn addition, the program takes these options (which must be given before the arguments):\n\n" << "-h Print help and exit.\n" << "-l List report step numbers in the selected restart file.\n" << "-i Enforce IX standard on output file.\n" << "-r Extract and convert a spesific report time step number from a unified restart file. \n\n"; } int main(int argc, char **argv) { int c = 0; int reportStepNumber = -1; bool specificReportStepNumber = false; bool listProperties = false; bool enforce_ix_output = false; while ((c = getopt(argc, argv, "hr:li")) != -1) { switch (c) { case 'h': printHelp(); return 0; case 'l': listProperties=true; break; case 'i': enforce_ix_output=true; break; case 'r': specificReportStepNumber=true; reportStepNumber = atoi(optarg); break; default: return EXIT_FAILURE; } } int argOffset = optind; // start reading auto start = std::chrono::system_clock::now(); std::string filename = argv[argOffset]; EclFile file1(filename); bool formattedOutput = file1.formattedInput() ? false : true; int p = filename.find_last_of("."); int l = filename.length(); std::string rootN = filename.substr(0,p); std::string extension = filename.substr(p,l-p); std::string resFile; if (listProperties) { if (extension==".UNRST") { ERst rst1(filename); rst1.loadData("INTEHEAD"); std::vector reportStepList=rst1.listOfReportStepNumbers(); for (auto seqn : reportStepList) { std::vector inteh = rst1.getRestartData("INTEHEAD", seqn, 0); std::cout << "Report step number: " << std::setfill(' ') << std::setw(4) << seqn << " Date: " << inteh[66] << "/" << std::setfill('0') << std::setw(2) << inteh[65] << "/" << std::setfill('0') << std::setw(2) << inteh[64] << std::endl; } std::cout << std::endl; } else { std::cout << "\n!ERROR, option -l only only available for unified restart files (*.UNRST) " << std::endl; exit(1); } return 0; } std::map to_formatted = {{".EGRID", ".FEGRID"}, {".INIT", ".FINIT"}, {".SMSPEC", ".FSMSPEC"}, {".UNSMRY", ".FUNSMRY"}, {".UNRST", ".FUNRST"}, {".RFT", ".FRFT"}, {".ESMRY", ".FESMRY"}}; std::map to_binary = {{".FEGRID", ".EGRID"}, {".FINIT", ".INIT"}, {".FSMSPEC", ".SMSPEC"}, {".FUNSMRY", ".UNSMRY"}, {".FUNRST", ".UNRST"}, {".FRFT", ".RFT"}, {".FESMRY", ".ESMRY"}}; if (formattedOutput) { auto search = to_formatted.find(extension); if (search != to_formatted.end()){ resFile = rootN + search->second; } else if (extension.substr(1,1)=="X"){ resFile = rootN + ".F" + extension.substr(2); } else if (extension.substr(1,1)=="S"){ resFile = rootN + ".A" + extension.substr(2); } else { std::cout << "\n!ERROR, unknown file type for input file '" << rootN + extension << "'\n" << std::endl; exit(1); } } else { auto search = to_binary.find(extension); if (search != to_binary.end()){ resFile = rootN + search->second; } else if (extension.substr(1,1)=="F"){ resFile = rootN + ".X" + extension.substr(2); } else if (extension.substr(1,1)=="A"){ resFile = rootN + ".S" + extension.substr(2); } else { std::cout << "\n!ERROR, unknown file type for input file '" << rootN + extension << "'\n" << std::endl; exit(1); } } std::cout << "\033[1;31m" << "\nconverting " << argv[argOffset] << " -> " << resFile << "\033[0m\n" << std::endl; EclOutput outFile(resFile, formattedOutput); if ((file1.is_ix()) || (enforce_ix_output)) { std::cout << "setting IX flag on output file \n"; outFile.set_ix(); } if (specificReportStepNumber) { if (extension!=".UNRST") { std::cout << "\n!ERROR, option -r only can only be used with unified restart files (*.UNRST) " << std::endl; exit(1); } ERst rst1(filename); if (!rst1.hasReportStepNumber(reportStepNumber)) { std::cout << "\n!ERROR, selected unified restart file doesn't have report step number " << reportStepNumber << "\n" << std::endl; exit(1); } rst1.loadReportStepNumber(reportStepNumber); auto arrayList = rst1.listOfRstArrays(reportStepNumber); writeArrayList(arrayList, rst1, reportStepNumber, outFile); } else { file1.loadData(); auto arrayList = file1.getList(); std::vector elementSizeList = file1.getElementSizeList(); writeArrayList(arrayList, elementSizeList, file1, outFile); } auto end = std::chrono::system_clock::now(); std::chrono::duration elapsed_seconds = end-start; std::cout << "runtime : " << argv[argOffset] << ": " << elapsed_seconds.count() << " seconds\n" << std::endl; return 0; }