subscript 1.5.0__tar.gz → 1.6.1__tar.gz

{subscript-1.5.0 → subscript-1.6.1}/.github/workflows/codecov.yml

@@ -27,6 +27,7 @@ jobs:
       run: |
         sudo apt-get install software-properties-common
         sudo apt-add-repository ppa:opm/ppa
+        sudo apt-add-repository ppa:opm/testing
         sudo apt-get update
         sudo apt-get install mpi-default-bin
         sudo apt-get install libopm-simulators-bin
@@ -47,5 +48,6 @@ jobs:
         pytest tests --disable-warnings --cov=subscript --cov-report=xml
 
     - name: Upload coverage to Codecov
-      run: |
-        bash <(curl -s https://codecov.io/bash)
+      uses: codecov/codecov-action@v4
+      with:
+        token: ${{ secrets.CODECOV_TOKEN }}

{subscript-1.5.0 → subscript-1.6.1}/.github/workflows/subscript.yml

@@ -39,6 +39,7 @@ jobs:
         run: |
           sudo apt-get install software-properties-common
           sudo apt-add-repository ppa:opm/ppa
+          sudo apt-add-repository ppa:opm/testing
           sudo apt-get update
           sudo apt-get install mpi-default-bin
           sudo apt-get install libopm-simulators-bin

{subscript-1.5.0 → subscript-1.6.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: subscript
-Version: 1.5.0
+Version: 1.6.1
 Summary: Equinor's collection of subsurface reservoir modelling scripts
 Author-email: Equinor <rnyb@equinor.com>
 License:                     GNU GENERAL PUBLIC LICENSE
@@ -700,7 +700,7 @@ Description-Content-Type: text/markdown
 License-File: LICENSE
 Requires-Dist: resdata
 Requires-Dist: res2df
-Requires-Dist: ert>=2.38.0b7
+Requires-Dist: ert
 Requires-Dist: fmu-tools
 Requires-Dist: matplotlib
 Requires-Dist: numpy<2
@@ -730,7 +730,7 @@ Requires-Dist: types-python-dateutil; extra == "tests"
 Provides-Extra: docs
 Requires-Dist: autoapi; extra == "docs"
 Requires-Dist: sphinx; extra == "docs"
-Requires-Dist: sphinx-argparse; extra == "docs"
+Requires-Dist: sphinx-argparse<0.5.0; extra == "docs"
 Requires-Dist: sphinx-autodoc-typehints; extra == "docs"
 Requires-Dist: sphinx-copybutton; extra == "docs"
 Requires-Dist: sphinx-togglebutton; extra == "docs"

subscript-1.6.1/docs/scripts/grav_subs_maps.rst

@@ -0,0 +1,10 @@
+
+GRAV_SUBS_MAPS
+==============
+
+.. argparse::
+   :module: subscript.grav_subs_maps.grav_subs_maps
+   :func: get_parser
+   :prog: grav_subs_maps
+	
+   

{subscript-1.5.0 → subscript-1.6.1}/docs/scripts/summaryplot.rst

@@ -3,7 +3,7 @@ SUMMARYPLOT
 ===========
 
 ``summaryplot`` is a command line utility to generate plots from Eclipse
-simulations, based on `libecl <http://github.com/equinor/libecl>`_
+simulations, based on `resdata <http://github.com/equinor/resdata>`_
 for processing Eclipse output files and
 `matplotlib <http://matplotlib.sourceforge.net>`_ for plotting.
 

{subscript-1.5.0 → subscript-1.6.1}/pyproject.toml

@@ -38,7 +38,7 @@ dynamic = ["version"]
 dependencies = [
     "resdata",
     "res2df",
-    "ert>=2.38.0b7",
+    "ert",
     "fmu-tools",
     "matplotlib",
     "numpy<2",
@@ -72,7 +72,7 @@ tests = [
 docs = [
     "autoapi",
     "sphinx",
-    "sphinx-argparse",
+    "sphinx-argparse<0.5.0",
     "sphinx-autodoc-typehints",
     "sphinx-copybutton",
     "sphinx-togglebutton",
@@ -96,6 +96,7 @@ eclcompress = "subscript.eclcompress.eclcompress:main"
 ecldiff2roff = "subscript.ecldiff2roff.ecldiff2roff:main"
 fmu_copy_revision = "subscript.fmu_copy_revision.fmu_copy_revision:main"
 fmuobs = "subscript.fmuobs.fmuobs:main"
+grav_subs_maps = "subscript.grav_subs_maps.grav_subs_maps:main"
 interp_relperm = "subscript.interp_relperm.interp_relperm:main"
 merge_rft_ertobs = "subscript.merge_rft_ertobs.merge_rft_ertobs:main"
 merge_unrst_files = "subscript.merge_unrst_files.merge_unrst_files:main"

{subscript-1.5.0 → subscript-1.6.1}/src/subscript/check_swatinit/check_swatinit.py

@@ -326,7 +326,7 @@ def qc_flag(qc_frame: pd.DataFrame) -> pd.DataFrame:
 
     contact = "OWC" if "OWC" in qc_frame else "GWC"
 
-    # Eclipse and libecl does not calculate cell centres to the same decimals.
+    # Eclipse and resdata does not calculate cell centres to the same decimals.
     # Add some tolerance when testing towards fluid contacts.
     contacttolerance = 1e-4
 

subscript-1.6.1/src/subscript/config_jobs/GRAV_SUBS_MAPS

@@ -0,0 +1,9 @@
+EXECUTABLE  grav_subs_maps
+
+DEFAULT  <ROOT_PATH> "./"
+
+ARGLIST  "--configfile" <GRAVMAPS_CONFIG> "--root-path" <ROOT_PATH> "--outputdir" <OUTPUT_DIR> <UNRST_FILE> 
+
+MIN_ARG  2
+MAX_ARG  4
+ARG_TYPE 0 STRING

{subscript-1.5.0 → subscript-1.6.1}/src/subscript/eclcompress/eclcompress.py

@@ -312,7 +312,7 @@ def find_keyword_sets(filelines: List[str]) -> List[Tuple[int, int]]:
         /
 
     we are not able to detect anything but the first record (line) without
-    having a full Eclipse parser (OPM). This means we we only compress the
+    having a full Eclipse parser (OPM). This means we only compress the
     first line. These type of keywords are not important to compress, and we
     could just as well avoid compressing them altogether.
 
@@ -336,7 +336,7 @@ def find_keyword_sets(filelines: List[str]) -> List[Tuple[int, int]]:
             continue
         # Remove embracing quotes if in a multi-keyword
         keyword = line.split(" ")[0].strip("'")
-        if keyword in ALLOWLIST_KEYWORDS:
+        if (keyword in ALLOWLIST_KEYWORDS) or keyword.startswith("FIP"):
             kwstart = lineidx
             if "/" in line:
                 keywordsets.append((kwstart, lineidx))

subscript-1.6.1/src/subscript/grav_subs_maps/grav_subs_maps.py

@@ -0,0 +1,324 @@
+import argparse
+import logging
+import os
+import sys
+from datetime import date
+from pathlib import Path
+from typing import Any, Dict, List, Optional
+
+import xtgeo
+import yaml
+from pydantic import BaseModel, Field, FilePath, field_validator
+from resdata.gravimetry import ResdataGrav, ResdataSubsidence
+from resdata.grid import Grid
+from resdata.resfile import ResdataFile
+from typing_extensions import Annotated
+
+import subscript
+
+logger = subscript.getLogger(__name__)
+
+# Constant for subsidence modelling, not influencing results
+# since subsidence is calculated from porevolume change
+# therefore defaulted
+DUMMY_YOUNGS = 0.5
+
+PREFIX_GRAVSURF = "all--delta_gravity_"
+PREFIX_SUBSSURF = "all--subsidence"
+
+DESCRIPTION = """
+Modelling maps of gravity change and subsidence from flow
+simulation output (EGRID, INIT and UNRST files).
+
+The script reads flow simulation results and a yaml configuration file specifying input
+and calculation parameters. Output is surfaces in irap binary format.
+"""
+
+EPILOGUE = """
+.. code-block:: yaml
+
+  # Example config file for grav_subs_maps
+
+  input:
+      diffdates:
+        - [2020-07-01, 2018-01-01] # Difference date to model. Must exist in UNRST file.
+      seabed_map: seabed.gri  # Path to file with seabed, irap binary format.
+                              # Also used as map template
+
+    calculations:
+      poisson_ratio: 0.45  # For subsidence calulcations, used in Geertsma model
+      coarsening: 8        # Coarsening factor for maps to speed up calculations
+      phases: ['gas', 'oil','water', 'total']  # One map for each phase specified
+
+"""
+
+CATEGORY = "modelling.reservoir"
+
+EXAMPLES = """
+.. code-block:: console
+
+ FORWARD_MODEL GRAV_SUBS_MAPS(<UNRST_FILE>=<ECLBASE>.UNRST, <GRAV_CONFIG>=grav_subs_maps.yml, <ROOT_PATH>=<CONFIG_PATH>, <OUTPUTDIR>=share/results/maps)
+
+
+where ``ECLBASE`` is already defined in your ERT config, pointing to the flowsimulator
+basename relative to ``RUNPATH``, grav_subs_maps.yml is a YAML file defining
+the inputs and modelling parameters and ``OUTPUTDIR`` is the path to the output folder.
+``ROOT_PATH`` is optinal and defaulted to "./". This is the rooth path assumed for
+relative paths in the yml config.
+
+The directory to export maps to must exist.
+"""  # noqa
+
+
+class GravInput(BaseModel):
+    diffdates: List[List[date]]
+    seabed_map: FilePath
+
+
+class GravCalc(BaseModel):
+    poisson_ratio: Annotated[float, Field(strict=True, ge=0, le=0.5)]
+    coarsening: Optional[Annotated[int, Field(strict=True, ge=1)]] = None
+    phases: List[str]
+
+    @field_validator("phases")
+    @classmethod
+    def check_phases(cls, phases: List[str]) -> List[str]:
+        allowed_phases = ["oil", "gas", "water", "total"]
+        for item in phases:
+            assert item in allowed_phases, f"allowed phases are {str(allowed_phases)}"
+        return phases
+
+
+class GravMapsConfig(BaseModel):
+    input: GravInput
+    calculations: GravCalc
+
+
+def get_parser() -> argparse.ArgumentParser:
+    """Function to create the argument parser that is going to be served to the user.
+
+    Returns:
+        argparse.ArgumentParser: The argument parser to be served
+
+    """
+    parser = argparse.ArgumentParser(
+        prog="grav_subs_maps.py",
+        description=DESCRIPTION,
+        epilog=EPILOGUE,
+        formatter_class=argparse.RawTextHelpFormatter,
+    )
+
+    parser.add_argument("UNRSTfile", type=str, help="Path to flowsimulator UNRST file")
+    parser.add_argument(
+        "-c",
+        "-C",
+        "--configfile",
+        type=str,
+        help="Name of YAML config file",
+        required=True,
+    )
+    parser.add_argument(
+        "-r",
+        "--root-path",
+        type=str,
+        default="./",
+        help=("Root path assumed for relative paths" " in config file."),
+    )
+    parser.add_argument(
+        "-o",
+        "--outputdir",
+        type=str,
+        help="Path to directory for output maps. Directory must exist.",
+        default="./",
+    )
+    parser.add_argument(
+        "--version",
+        action="version",
+        version="%(prog)s (subscript version " + subscript.__version__ + ")",
+    )
+    return parser
+
+
+def main() -> None:
+    """Invocated from the command line, parsing command line arguments"""
+    parser = get_parser()
+    args = parser.parse_args()
+
+    logger.setLevel(logging.INFO)
+
+    # parse the config file
+    if not Path(args.configfile).exists():
+        sys.exit("No such file:" + args.configfile)
+    config = yaml.safe_load(Path(args.configfile).read_text(encoding="utf8"))
+
+    # cfg = GravMapsConfig.model_validate(config).model_dump()
+
+    if not Path(args.outputdir).exists():
+        sys.exit("Output folder does not exist:" + args.outputdir)
+    if not Path(args.UNRSTfile).exists():
+        sys.exit("UNRST file does not exist:" + args.UNRSTfile)
+
+    main_gravmaps(args.UNRSTfile, config, Path(args.root_path), Path(args.outputdir))
+
+
+def prepend_root_path_to_relative_files(
+    cfg: Dict[str, Any], root_path: Path
+) -> Dict[str, Any]:
+    """Prepend root_path to relative files found paths in a configuration
+    dictionary.
+
+    Note: This function is before prior to validation of the configuration!
+
+    Will look for filename in the key "input["seabed_map"]"
+
+    Args:
+        cfg: grav_subs_maps configuration dictionary
+        root_path: An relative or absolute path to be prepended
+
+    Returns:
+        Modified configuration for interp_relperm
+    """
+    if (
+        "input" in cfg
+        and "seabed_map" in cfg["input"]
+        and not os.path.isabs(cfg["input"]["seabed_map"])
+    ):
+        cfg["input"]["seabed_map"] = str(root_path / Path(cfg["input"]["seabed_map"]))
+    return cfg
+
+
+def main_gravmaps(
+    unrst_file: str,
+    config: Dict[str, Any],
+    root_path: Optional[Path],
+    output_folder: Path,
+) -> None:
+    """
+    Process a configuration, model gravity and subsidence surfaces and write to disk.
+
+    Args:
+        resdata: Path to flow simulation UNRST file
+        config: Configuration for modelling
+    """
+
+    if root_path is not None:
+        config = prepend_root_path_to_relative_files(config, root_path)
+
+    cfg = GravMapsConfig.model_validate(config).model_dump()
+
+    # Read inputs and calculation parameters
+    input_diffdates = cfg["input"]["diffdates"]
+    map_template = cfg["input"]["seabed_map"]
+    coarsening = cfg["calculations"]["coarsening"]
+    phases = cfg["calculations"]["phases"]
+    poisson_ratio = cfg["calculations"]["poisson_ratio"]
+
+    # Read seabed map and coarsen
+    seabed = xtgeo.surface_from_file(map_template)
+    seabed.coarsen(coarsening)
+
+    if isinstance(unrst_file, str):
+        restart_file = unrst_file[:-6] + ".UNRST"
+        egrid_file = unrst_file[:-6] + ".EGRID"
+        init_file = unrst_file[:-6] + ".INIT"
+        grid = Grid(egrid_file)
+        init = ResdataFile(init_file)
+        rest = ResdataFile(restart_file)
+
+    restart_index = {}
+
+    # From restart datetime format to YYYYMMDD as key
+    for i, restart_date in enumerate(rest.dates):
+        restart_index[restart_date.strftime("%Y%m%d")] = i
+
+    diffdates = []
+    # Convert dates from datetime format to strings
+    logger.info("Will do modelling for diffdates: ")
+    for diffdate in input_diffdates:
+        diff = [diffdate[0].strftime("%Y%m%d"), diffdate[1].strftime("%Y%m%d")]
+        diffdates.append(diff)
+        logger.info(f"{diffdate[0]}_{diffdate[1]}")
+
+    grav = ResdataGrav(grid, init)
+    subsidence = ResdataSubsidence(grid, init)
+
+    added_dates = []
+
+    for diffdate in diffdates:
+        for singledate in diffdate:  # base and monitor
+            rsb = rest.restartView(0)
+            if singledate not in added_dates:
+                if singledate in restart_index:
+                    rsb = rest.restartView(restart_index[singledate])
+                    grav.add_survey_RFIP(singledate, rsb)
+                    subsidence.add_survey_PRESSURE(singledate, rsb)
+                    added_dates.append(singledate)
+                else:
+                    logger.error(
+                        f"Date {singledate} specified but not found in UNRST file."
+                    )
+                    sys.exit(1)
+    phase_code = {"oil": 1, "gas": 2, "water": 4, "total": 7}
+
+    # Gravity
+    for diffdate in diffdates:
+        for phase in phases:
+            logger.info(
+                f"Calculating delta gravity map from {phase} "
+                f"for {diffdate[0]}_{diffdate[1]}"
+            )
+            dgsim = seabed.copy()
+            df_dgsim = dgsim.get_dataframe()
+            dgsim_series = []
+            for index, row in df_dgsim.iterrows():
+                dgsim_series.append(
+                    grav.eval(
+                        diffdate[1],
+                        diffdate[0],
+                        (row["X_UTME"], row["Y_UTMN"], row["VALUES"]),
+                        phase_mask=phase_code[phase],
+                    )
+                )
+            dgsim.values = dgsim_series
+            filename = (
+                PREFIX_GRAVSURF
+                + phase
+                + "--"
+                + diffdate[0]
+                + "_"
+                + diffdate[1]
+                + ".gri"
+            )
+            dgsim.to_file(os.path.join(output_folder, filename))
+
+    # Subsidence
+    for diffdate in diffdates:
+        logger.info(f"Calculating subsidence map for {diffdate[0]}_{diffdate[1]}")
+        dzsim = seabed.copy()
+        df_dzsim = dzsim.get_dataframe()
+        dzsim_series = []
+        for index, row in df_dzsim.iterrows():
+            dzsim_series.append(
+                subsidence.eval_geertsma_rporv(
+                    diffdate[1],
+                    diffdate[0],
+                    (row["X_UTME"], row["Y_UTMN"], row["VALUES"]),
+                    DUMMY_YOUNGS,
+                    poisson_ratio,
+                    row["VALUES"],
+                )
+            )
+
+        dzsim.values = [i * 100 for i in dzsim_series]  # From m to cms
+
+        filename = PREFIX_SUBSSURF + "--" + diffdate[0] + "_" + diffdate[1] + ".gri"
+        dzsim.to_file(os.path.join(output_folder, filename))
+
+    logger.info(
+        "Done; All gravity and subsidence maps written to folder: %s",
+        str(output_folder),
+    )
+
+
+if __name__ == "__main__":
+    main()

{subscript-1.5.0 → subscript-1.6.1}/src/subscript/interp_relperm/interp_relperm.py

@@ -167,7 +167,10 @@ def make_wateroilgas(dframe: pd.DataFrame, delta_s: float) -> pyscal.WaterOilGas
             The data must be restricted to only one SATNUM.
     """
     dframe = dframe.copy()  # Copy since we will modify it.
-    wog = pyscal.WaterOilGas(swl=dframe["SW"].min(), h=delta_s)
+    wog = pyscal.WaterOilGas(
+        swl=dframe["SW"].min() if "SW" in dframe.columns else 1.0 - dframe["SG"].max(),
+        h=delta_s,
+    )
     if "PCOW" not in dframe:
         dframe = dframe.assign(PCOW=0)
     if "PCOG" not in dframe:
@@ -191,6 +194,7 @@ def make_wateroilgas(dframe: pd.DataFrame, delta_s: float) -> pyscal.WaterOilGas
             .drop_duplicates()
             .reset_index()
         )
+        wog.wateroil.add_fromtable(wo_dframe)
         go_dframe = (
             dframe[["SG", "KRG", "KROG", "PCOG"]]
             .set_index("SG")
@@ -202,23 +206,33 @@ def make_wateroilgas(dframe: pd.DataFrame, delta_s: float) -> pyscal.WaterOilGas
             .drop_duplicates()
             .reset_index()
         )
+        wog.gasoil.add_fromtable(go_dframe)
     else:
-        wo_dframe = dframe[["SW", "KRW", "KROW", "PCOW"]].dropna().reset_index()
-        go_dframe = dframe[["SG", "KRG", "KROG", "PCOG"]].dropna().reset_index()
-
-    wog.wateroil.add_fromtable(wo_dframe)
-    wog.gasoil.add_fromtable(go_dframe)
+        if {"SW", "KRW", "KROW", "PCOW"}.issubset(dframe.columns):
+            wog.wateroil.add_fromtable(
+                dframe[["SW", "KRW", "KROW", "PCOW"]].dropna().reset_index()
+            )
+        else:
+            wog.wateroil = None
+        if {"SG", "KRG", "KROG", "PCOG"}.issubset(dframe.columns):
+            wog.gasoil.add_fromtable(
+                dframe[["SG", "KRG", "KROG", "PCOG"]].dropna().reset_index()
+            )
+        else:
+            wog.gasoil = None
 
     # socr can for floating point reasons become estimated to be larger than
     # sorw, which means we are in an oil paleo zone setting. This is not
     # supported by interp_relperm. Reset the property to ensure interpolation
     # is not affected:
-    wog.wateroil.socr = wog.wateroil.sorw
+    if wog.wateroil:
+        wog.wateroil.socr = wog.wateroil.sorw
 
     # If sgro > 0, it is a gas condensate object, which cannot be
     # mixed with non-gas condensate (during interpolation). Avoid pitfalls
     # in the estimated sgro by always setting it to zero:
-    wog.gasoil.sgro = 0.0
+    if wog.gasoil:
+        wog.gasoil.sgro = 0.0
     return wog
 
 
@@ -248,7 +262,9 @@ def make_interpolant(
     high = make_wateroilgas(high_df.loc[satnum], delta_s)
     rec = pyscal.SCALrecommendation(low, base, high, "SATNUM " + str(satnum), h=delta_s)
 
-    return rec.interpolate(interp_param["param_w"], interp_param["param_g"], h=delta_s)
+    return rec.interpolate(
+        interp_param.get("param_w", 0.0), interp_param.get("param_g", 0), h=delta_s
+    )
 
 
 def get_parser() -> argparse.ArgumentParser:

{subscript-1.5.0 → subscript-1.6.1}/src/subscript/legacy/runeclipse

@@ -8,15 +8,26 @@ use Cwd;
 use IO::File;
 
 # ------------------------------------------------------
+my $hostname           = `hostname`;
+my $_balance     = `date +%s`;
 
+$_balance = $_balance % 2;
+my $lisens = "7321\@eclipse1-lic-no.statoil.no:7321\@eclipse2-lic-no.statoil.no";
+if ($_balance == 1) {
+    $lisens = "7321\@eclipse2-lic-no.statoil.no:7321\@eclipse1-lic-no.statoil.no";
+}
+
+if (index($hostname, "stjohn-") == 0) {
+    $lisens = "7321\@stjohn-grid01.stjohn.statoil.no:$lisens";
+} elsif (index($hostname, "rio-") == 0) {
+    $lisens = "7321\@rio-grid01.rio.statoil.no:$lisens";
+}
 
 my $ecldir       = "/prog/res/ecl/grid";
 my $eclolddir    = "/prog/ecl/grid";
 my $ecltestdir   = "/prog/res/ecl/testgrid";
 my $mpiintelroot = "/prog/res/ecl/grid/tools/linux_x86_64/intel/mpi/2018.1.163";
 
-my $lisens = "7321\@eclipse-lic-no.statoil.no";
-
 my $extconfig   = "macros/CONFIG.ECL";
 my $extmpirun64 = "tools/linux_x86_64/mpich_x86_64/bin/mpirun";
 
@@ -73,7 +84,7 @@ my $workdir = cwd;
 my $scriptname = "runeclipse";
 
 my $usage =
-  "usage $scriptname  [-p program] [-q que_name] [-v version] [-x]  [-i] root";
+    "usage $scriptname  [-p program] [-q que_name] [-v version] [-x]  [-i] root";
 
 #
 # Get command line options:
@@ -400,7 +411,7 @@ sub build_batch {
     if ($testversion) {
         print BATCHFIL "echo \"*** RUNNING TEST VERSION $version ****\"\n";
     }
-    print BATCHFIL "export LM_LICENSE_FILE=$lisens\n";
+    print BATCHFIL "export SLBSLS_LICENSE_FILE=$lisens\n";
     print BATCHFIL "export F_UFMTENDIAN=big\n";
     print BATCHFIL "export ARCH=\`uname -m\`\n";
     print BATCHFIL "cd $workdir\n";
@@ -437,7 +448,7 @@ sub build_parallel_batch {
       or die "$scriptname: Cannot create batch file $batchfilename\n";
 
     print BATCHFIL "#!/bin/sh \n";
-    print BATCHFIL "export LM_LICENSE_FILE=$lisens\n";
+    print BATCHFIL "export SLBSLS_LICENSE_FILE=$lisens\n";
 
     #    print BATCHFIL
     #        "export FLEXLM_NO_CKOUT_INSTALL_LIC=1\n" ;

{subscript-1.5.0 → subscript-1.6.1}/src/subscript/restartthinner/restartthinner.py

@@ -30,9 +30,9 @@ written to the same filename (keeping the original is optional)
 def find_resdata_app(toolname: str) -> str:
     """Locate path of apps in resdata.
 
-    These have varying suffixes due through the history of libecl Makefiles.
+    These have varying suffixes due through the history of resdata Makefiles.
 
-    Depending on libecl-version, it has the .x or the .c.x suffix
+    Depending on resdata-version, it has the .x or the .c.x suffix
     We prefer .x.
 
     Returns:
@@ -75,7 +75,7 @@ def rd_repacker(rstfilename: str, slicerstindices: list, quiet: bool) -> None:
     modify the original filename.
     """
     out = " >/dev/null" if quiet else ""
-    # Error early if libecl tools are not available
+    # Error early if resdata tools are not available
     try:
         find_resdata_app("rd_unpack")
         find_resdata_app("rd_pack")

{subscript-1.5.0 → subscript-1.6.1}/src/subscript/version.py

@@ -12,5 +12,5 @@ __version__: str
 __version_tuple__: VERSION_TUPLE
 version_tuple: VERSION_TUPLE
 
-__version__ = version = '1.5.0'
-__version_tuple__ = version_tuple = (1, 5, 0)
+__version__ = version = '1.6.1'
+__version_tuple__ = version_tuple = (1, 6, 1)

{subscript-1.5.0 → subscript-1.6.1}/src/subscript/welltest_dpds/welltest_dpds.py

@@ -7,7 +7,6 @@ from pathlib import Path
 import numpy as np
 import pandas as pd
 from resdata.summary import Summary
-from scipy.interpolate import interp1d
 
 from subscript import __version__
 
@@ -419,19 +418,7 @@ def genobs_vec(filen, vec, time):
     dframe = pd.read_csv(filen, sep="\t")
     obs_time = dframe["dTime"][1:None].dropna().to_numpy(dtype=float)
 
-    gen_data = np.zeros(len(obs_time))
-
-    interp = interp1d(time, vec)
-
-    for idx, timepoint in enumerate(obs_time):
-        if timepoint < time[0]:
-            gen_data[idx] = vec[0]
-        elif timepoint > time[-1]:
-            gen_data[idx] = vec[-1]
-        else:
-            gen_data[idx] = interp(timepoint)
-
-    return gen_data
+    return np.interp(obs_time, time, vec)
 
 
 def main():

{subscript-1.5.0 → subscript-1.6.1}/src/subscript.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: subscript
-Version: 1.5.0
+Version: 1.6.1
 Summary: Equinor's collection of subsurface reservoir modelling scripts
 Author-email: Equinor <rnyb@equinor.com>
 License:                     GNU GENERAL PUBLIC LICENSE
@@ -700,7 +700,7 @@ Description-Content-Type: text/markdown
 License-File: LICENSE
 Requires-Dist: resdata
 Requires-Dist: res2df
-Requires-Dist: ert>=2.38.0b7
+Requires-Dist: ert
 Requires-Dist: fmu-tools
 Requires-Dist: matplotlib
 Requires-Dist: numpy<2
@@ -730,7 +730,7 @@ Requires-Dist: types-python-dateutil; extra == "tests"
 Provides-Extra: docs
 Requires-Dist: autoapi; extra == "docs"
 Requires-Dist: sphinx; extra == "docs"
-Requires-Dist: sphinx-argparse; extra == "docs"
+Requires-Dist: sphinx-argparse<0.5.0; extra == "docs"
 Requires-Dist: sphinx-autodoc-typehints; extra == "docs"
 Requires-Dist: sphinx-copybutton; extra == "docs"
 Requires-Dist: sphinx-togglebutton; extra == "docs"