Skip to content

Tensor interpolation for non-stationary geostat fields #641

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
wkitlasten wants to merge 37 commits into
base: develop
Choose a base branch
from
Open

Tensor interpolation for non-stationary geostat fields #641

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
37 commits
Select commit Hold shift + click to select a range
069b8ad
added an option for full tensor interpolation in apply_ppu_hyperpars_…
weskitlasten Jul 2, 2025
b3f275c
functioning tensor interpolation. play with conceptual_point params s…
weskitlasten Jul 4, 2025
fb6513f
working version, minor issues with array shape, working on making eve…
weskitlasten Jul 8, 2025
210824c
working version, smooth boundary, increase variance at boundar
weskitlasten Jul 9, 2025
404e14f
minor tweaks, good version
weskitlasten Jul 9, 2025
e4b9eda
external iid files, small tweaks
weskitlasten Jul 10, 2025
17ee758
create 3d tensors
weskitlasten Jul 17, 2025
ec39ff2
3d without zones
weskitlasten Jul 18, 2025
a244824
3d with zones
weskitlasten Jul 18, 2025
bdf26fb
3d tensor interpolation, some functions added to plot_utils
weskitlasten Jul 22, 2025
436dd62
removing some redundancies
weskitlasten Jul 22, 2025
a58cebe
input data all non-transformed
weskitlasten Aug 5, 2025
2e03035
moved variogram def out of for loop, other efficiencies
weskitlasten Aug 5, 2025
6771a99
fixed transform bug and cleaned up some transform bits, much better r…
weskitlasten Aug 5, 2025
7fcbf69
moving example out of tensor_utils
weskitlasten Aug 5, 2025
aea823d
pypestutils implementation with tensors
weskitlasten Aug 13, 2025
0f8c1b4
pypestutils implementation todo:clean up file handling and figures
weskitlasten Aug 14, 2025
3c86600
proper back transformation, circle example, todo fieldgen auto filena…
weskitlasten Aug 15, 2025
7ca5bb9
massive refactoring to use modelgrid object, dealing with layers better
weskitlasten Aug 19, 2025
663204a
plot 1 layer at a time
weskitlasten Aug 19, 2025
a8e8ab3
useage tweaks
weskitlasten Aug 20, 2025
df44e57
model files for HBHP, dis only
weskitlasten Aug 30, 2025
3d9e834
model_name as parameter
weskitlasten Aug 30, 2025
f06bcbd
additional modflow files required for dis
weskitlasten Aug 30, 2025
1283c51
forgot shared_props folder from hbhp model structure
weskitlasten Aug 30, 2025
bd62b4b
results vs field bug fixed, tweak some plots
weskitlasten Sep 2, 2025
30ca942
cleaned up some figures
weskitlasten Sep 2, 2025
5fee6db
figure arrow tweaks
weskitlasten Sep 2, 2025
55a50ac
no need for points at the tail
weskitlasten Sep 2, 2025
e13f979
function to output raster file
weskitlasten Sep 3, 2025
99869b0
save tif to figures dir
weskitlasten Sep 3, 2025
c9a6b91
pull from upstream
weskitlasten Nov 23, 2025
acca030
new python implementation of fieldgen, good on tangential but wairau …
weskitlasten Nov 27, 2025
7d0bc8d
working python version, moving average options, multiple realizations
weskitlasten Nov 27, 2025
509c510
fixed bug in iid generation, python similar to fieldgen for all teste…
weskitlasten Nov 28, 2025
d4ee72d
cleaned up a bit
weskitlasten Nov 29, 2025
bc68240
model_ws change
weskitlasten Dec 4, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
270 changes: 270 additions & 0 deletions examples/HBHP/hpm_mf6/Bridgpa_HU2_sva.csv
View file
Open in desktop

Large diffs are not rendered by default.

151,303 changes: 151,303 additions & 0 deletions examples/HBHP/hpm_mf6/grid.dat
View file
Open in desktop

Large diffs are not rendered by default.

42 changes: 42 additions & 0 deletions examples/HBHP/hpm_mf6/hpm_mf6.dis
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,42 @@
# File generated by Flopy version 3.4.2 on 04/27/2024 at 08:12:37.
BEGIN options
LENGTH_UNITS meters
XORIGIN 1.89146000E+06
YORIGIN 5.59440000E+06
END options

BEGIN dimensions
NLAY 9
NROW 302
NCOL 501
END dimensions

BEGIN griddata
delr
OPEN/CLOSE 'hpm_mf6.dis_delr.txt' FACTOR 1.0
delc
OPEN/CLOSE 'hpm_mf6.dis_delc.txt' FACTOR 1.0
top
OPEN/CLOSE '..\shared_props\top.txt' FACTOR 1.0
botm LAYERED
OPEN/CLOSE '..\shared_props\botm_0.txt' FACTOR 1.0
OPEN/CLOSE '..\shared_props\botm_1.txt' FACTOR 1.0
OPEN/CLOSE '..\shared_props\botm_2.txt' FACTOR 1.0
OPEN/CLOSE '..\shared_props\botm_3.txt' FACTOR 1.0
OPEN/CLOSE '..\shared_props\botm_4.txt' FACTOR 1.0
OPEN/CLOSE '..\shared_props\botm_5.txt' FACTOR 1.0
OPEN/CLOSE '..\shared_props\botm_6.txt' FACTOR 1.0
OPEN/CLOSE '..\shared_props\botm_7.txt' FACTOR 1.0
OPEN/CLOSE '..\shared_props\botm_8.txt' FACTOR 1.0
idomain LAYERED
OPEN/CLOSE '..\shared_props\idomain_0.txt' FACTOR 1
OPEN/CLOSE '..\shared_props\idomain_1.txt' FACTOR 1
OPEN/CLOSE '..\shared_props\idomain_2.txt' FACTOR 1
OPEN/CLOSE '..\shared_props\idomain_3.txt' FACTOR 1
OPEN/CLOSE '..\shared_props\idomain_4.txt' FACTOR 1
OPEN/CLOSE '..\shared_props\idomain_5.txt' FACTOR 1
OPEN/CLOSE '..\shared_props\idomain_6.txt' FACTOR 1
OPEN/CLOSE '..\shared_props\idomain_7.txt' FACTOR 1
OPEN/CLOSE '..\shared_props\idomain_8.txt' FACTOR 1
END griddata

1 change: 1 addition & 0 deletions examples/HBHP/hpm_mf6/hpm_mf6.dis_delc.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000
1 change: 1 addition & 0 deletions examples/HBHP/hpm_mf6/hpm_mf6.dis_delr.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000 100.00000000
28 changes: 28 additions & 0 deletions examples/HBHP/hpm_mf6/hpm_mf6.ims
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,28 @@
# File generated by Flopy version 3.4.2 on 04/27/2024 at 08:12:37.
BEGIN options
PRINT_OPTION all
COMPLEXITY complex
END options

BEGIN nonlinear
OUTER_DVCLOSE 0.00100000
OUTER_MAXIMUM 500
UNDER_RELAXATION dbd
UNDER_RELAXATION_GAMMA 0.20000000
UNDER_RELAXATION_THETA 0.70000000
UNDER_RELAXATION_KAPPA 0.10000000
UNDER_RELAXATION_MOMENTUM 0.00100000
BACKTRACKING_NUMBER 20
BACKTRACKING_TOLERANCE 1.10000000
BACKTRACKING_REDUCTION_FACTOR 0.20000000
BACKTRACKING_RESIDUAL_LIMIT 50.00000000
END nonlinear

BEGIN linear
INNER_MAXIMUM 100
INNER_DVCLOSE 1.00000000E-04
inner_rclose 0.01000000
PRECONDITIONER_LEVELS 8
PRECONDITIONER_DROP_TOLERANCE 1.00000000E-04
END linear

21 changes: 21 additions & 0 deletions examples/HBHP/hpm_mf6/hpm_mf6.nam
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,21 @@
# File generated by Flopy version 3.4.2 on 04/27/2024 at 08:12:37.
BEGIN options
SAVE_FLOWS
NEWTON UNDER_RELAXATION
END options

BEGIN packages
DIS6 hpm_mf6.dis dis
IC6 hpm_mf6.ic ic
SFR6 hpm_mf6.sfr sfr_0
WEL6 hpm_mf6.wel wel_0
NPF6 hpm_mf6.npf npf
OC6 hpm_mf6.oc oc
RCH6 hpm_mf6.rcha rcha_0
GHB6 hpm_mf6_coast.ghb ghb_0
GHB6 hpm_mf6_bound.ghb ghb_1
DRN6 hpm_mf6.drn drn_0
MVR6 hpm_mf6.mvr mvr
OBS6 hpm_mf6.obs obs_4
END packages

14 changes: 14 additions & 0 deletions examples/HBHP/hpm_mf6/hpm_mf6.tdis
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,14 @@
# File generated by Flopy version 3.4.2 on 04/27/2024 at 08:12:37.
BEGIN options
TIME_UNITS days
START_DATE_TIME 01-01-1940
END options

BEGIN dimensions
NPER 1
END dimensions

BEGIN perioddata
1.00000000 1 1.00000000
END perioddata

19 changes: 19 additions & 0 deletions examples/HBHP/hpm_mf6/mfsim.nam
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,19 @@
# File generated by Flopy version 3.4.2 on 04/27/2024 at 08:12:37.
BEGIN options
END options

BEGIN timing
TDIS6 hpm_mf6.tdis
END timing

BEGIN models
gwf6 hpm_mf6.nam hpm_mf6
END models

BEGIN exchanges
END exchanges

BEGIN solutiongroup 1
ims6 hpm_mf6.ims hpm_mf6
END solutiongroup 1

Loading
Loading