Non-contiguous active FIPNUM

src/subscript/prtvol2csv/prtvol2csv.py

@@ -379,10 +379,16 @@ def prtvol2df(
     # Remove extra empty 'regions' (from the reservoir volume table in .PRT)
     # Concatenate dataframes horizontally. Both are/must be indexed by value
     # of fipname (FIPNUM):
+
+    # Get maximum FIPNUM
+    sim_max_fipnum = simvolumes_df.index.max()
+    res_max_fipnum = resvolumes_df.loc[(resvolumes_df != 0).any(axis=1)].index.max()
+    max_fipnum = max(sim_max_fipnum, res_max_fipnum)
     volumes = (
-        pd.concat([simvolumes_df, resvolumes_df[: len(simvolumes_df)]], axis=1)
+        pd.concat([simvolumes_df, resvolumes_df[:max_fipnum]], axis=1)
         .apply(pd.to_numeric)
         .fillna(value=0.0)
+        .sort_index()
     )
 
     # Rename the index to "FIPNUM", as required by webviz-subsurface, if requested

tests/test_prtvol2csv.py

@@ -554,6 +554,19 @@ def test_fipxxx(tmp_path, mocker):
     expected["FIPNAME"] = "FIPNUM"
     pd.testing.assert_frame_equal(dframe, expected)
 
+def test_inactive_fipnum(tmp_path, mocker):
+    """Test the case with non-contiguous active FIPNUM"""
+
+    prtfile = TEST_PRT_DATADIR / "DROGON_INACTIVE_FIPNUM.PRT"
+
+    os.chdir(tmp_path)
+    with pytest.warns(FutureWarning, match="Output directories"):
+        mocker.patch("sys.argv", ["prtvol2csv", "--debug", str(prtfile)])
+        prtvol2csv.main()
+    dframe = pd.read_csv("share/results/volumes/simulator_volume_fipnum.csv")
+
+    assert dframe.loc[dframe["FIPNUM"].idxmax(), "HCPV_TOTAL"] != 0
+
 
 def test_find_prtfile(tmp_path):
     """Test location service for PRT files"""