fix: 🩹 strata_trapper: initialize output arrays with NaN

ImperialCollegeLondon · Aug 9, 2024 · 0f13d2c · 0f13d2c
1 parent c6dc14c
commit 0f13d2c
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/strata_trapper.m b/strata_trapper.m
@@ -13,9 +13,9 @@
 
 saturations = linspace(params.sw_resid,1,options.sat_num_points);
 
-cap_pres_upscaled = zeros(grid.cells.num,length(saturations));
-krw = zeros(grid.cells.num,3,length(saturations));
-krg = zeros(grid.cells.num,3,length(saturations));
+cap_pres_upscaled = nan(grid.cells.num,length(saturations));
+krw = nan(grid.cells.num,3,length(saturations));
+krg = nan(grid.cells.num,3,length(saturations));
 
 cells_num = min(length(mask),grid.cells.num);
 mask = mask(1:cells_num);