Extra command line argument write-lodsmry for creating LODSMRY file

- default is false
 - file LODSRMY optimized for fast loading of selected vectors

ebos/eclgenericwriter.cc

@@ -163,7 +163,8 @@ EclGenericWriter(const Schedule& schedule,
                  const GridView& gridView,
                  const Dune::CartesianIndexMapper<Grid>& cartMapper,
                  const Dune::CartesianIndexMapper<EquilGrid>* equilCartMapper,
-                 bool enableAsyncOutput)
+                 bool enableAsyncOutput,
+                 bool enableEsmry )
     : collectToIORank_(grid,
                        equilGrid,
                        gridView,
@@ -182,7 +183,7 @@ EclGenericWriter(const Schedule& schedule,
         eclIO_.reset(new EclipseIO(eclState_,
                                    UgGridHelpers::createEclipseGrid(*equilGrid, eclState_.getInputGrid()),
                                    schedule_,
-                                   summaryConfig_));
+                                   summaryConfig_, "", enableEsmry));
 
         const auto& wbp_calculators = eclIO_->summary().wbp_calculators( schedule.size() - 1 );
         wbp_index_list_ = wbp_calculators.index_list();

ebos/eclgenericwriter.hh

@@ -69,7 +69,8 @@ public:
                      const GridView& gridView,
                      const Dune::CartesianIndexMapper<Grid>& cartMapper,
                      const Dune::CartesianIndexMapper<EquilGrid>* equilCartMapper,
-                     bool enableAsyncOutput);
+                     bool enableAsyncOutput,
+                     bool enableEsmry);
 
     const EclipseIO& eclIO() const;
 

ebos/eclproblem.hh

@@ -437,6 +437,11 @@ template<class TypeTag>
 struct EnableAsyncEclOutput<TypeTag, TTag::EclBaseProblem> {
     static constexpr bool value = true;
 };
+// Write ESMRY file for fast loading of summary data
+template<class TypeTag>
+struct EnableEsmry<TypeTag, TTag::EclBaseProblem> {
+    static constexpr bool value = false;
+};
 
 // By default, use single precision for the ECL formated results
 template<class TypeTag>
@@ -2544,7 +2549,7 @@ private:
             if (FluidSystem::phaseIsActive(FluidSystem::waterPhaseIdx))
                 dofFluidState.setSaturation(FluidSystem::waterPhaseIdx,
                                             waterSaturationData[dofIdx]);
-                                            
+
             if (FluidSystem::phaseIsActive(FluidSystem::gasPhaseIdx)){
                 if (!FluidSystem::phaseIsActive(FluidSystem::oilPhaseIdx)){
                     dofFluidState.setSaturation(FluidSystem::gasPhaseIdx,

ebos/eclwriter.hh

@@ -53,6 +53,10 @@ template<class TypeTag, class MyTypeTag>
 struct EclOutputDoublePrecision {
     using type = UndefinedProperty;
 };
+template<class TypeTag, class MyTypeTag>
+struct EnableEsmry {
+    using type = UndefinedProperty;
+};
 
 } // namespace Opm::Properties
 
@@ -108,6 +112,8 @@ public:
 
         EWOMS_REGISTER_PARAM(TypeTag, bool, EnableAsyncEclOutput,
                              "Write the ECL-formated results in a non-blocking way (i.e., using a separate thread).");
+        EWOMS_REGISTER_PARAM(TypeTag, bool, EnableEsmry,
+                             "Write ESMRY file for fast loading of summary data.");
     }
 
     // The Simulator object should preferably have been const - the
@@ -122,7 +128,7 @@ public:
                    simulator.vanguard().gridView(),
                    simulator.vanguard().cartesianIndexMapper(),
                    simulator.vanguard().grid().comm().rank() == 0 ? &simulator.vanguard().equilCartesianIndexMapper() : nullptr,
-                   EWOMS_GET_PARAM(TypeTag, bool, EnableAsyncEclOutput))
+                   EWOMS_GET_PARAM(TypeTag, bool, EnableAsyncEclOutput), EWOMS_GET_PARAM(TypeTag, bool, EnableEsmry))
         , simulator_(simulator)
     {
         this->eclOutputModule_ = std::make_unique<EclOutputBlackOilModule<TypeTag>>(simulator, this->wbp_index_list_, this->collectToIORank_);