Skip to content

Commit

Permalink
pipt tutorial (#67)
Browse files Browse the repository at this point in the history
* Small update to handle ensemble keywords

* Added tutorial for pipt

* Update README.md

* Update README.md
  • Loading branch information
rolfjl authored Aug 5, 2024
1 parent 7c3cd13 commit ace6ec2
Show file tree
Hide file tree
Showing 16 changed files with 3,351 additions and 15 deletions.
3 changes: 2 additions & 1 deletion README.md
Original file line number Diff line number Diff line change
Expand Up @@ -67,7 +67,8 @@ python3 -m pip install -e .
PET needs to be set up with a configuration file. See the example folder for inspiration.

## Tutorials
A POPT tutorial is found [here](https://github.com/Python-Ensemble-Toolbox/PET/blob/main/docs/tutorials/popt/tutorial_popt.ipynb)
- A POPT tutorial is found [here](https://github.com/Python-Ensemble-Toolbox/PET/blob/main/docs/tutorials/popt/tutorial_popt.ipynb)
- A PIPT tutorial is found [here](https://github.com/Python-Ensemble-Toolbox/PET/blob/main/docs/tutorials/pipt/tutorial_pipt.ipynb)

## Suggested readings:

Expand Down
39 changes: 39 additions & 0 deletions docs/tutorials/pipt/3D_ESMDA.toml
Original file line number Diff line number Diff line change
@@ -0,0 +1,39 @@
[ensemble]
ne = 50.0
state = "permx"
prior_permx = [["vario", "sph"], ["mean", "priormean.npz"], ["var", 1.0], ["range", 10.0], ["aniso", 1.0],
["angle", 0.0], ["grid", [10.0, 10.0, 2.0]]]

[dataassim]
daalg = ["esmda", "esmda"]
analysis = "approx"
energy = 98.0
obsvarsave = "yes"
restartsave = "no"
analysisdebug = ["pred_data", "state", "data_misfit", "prev_data_misfit"]
restart = "no"
obsname = "days"

truedataindex = [400, 800, 1200, 1600, 2000, 2400, 2800, 3200, 3600, 4000]
truedata = "true_data.csv"
assimindex = [0,1,2,3,4,5,6,7,8,9]
datatype = ["WOPR PRO1", "WOPR PRO2", "WOPR PRO3", "WWPR PRO1", "WWPR PRO2",
"WWPR PRO3", "WWIR INJ1", "WWIR INJ2", "WWIR INJ3"]
staticvar = "permx"

datavar = "var.csv"

mda = [ ["tot_assim_steps", 3], ['inflation_param', [2, 4, 4]] ]

[fwdsim]
reporttype = "days"
reportpoint = [400, 800, 1200, 1600, 2000, 2400, 2800, 3200, 3600, 4000]
replace = "yes"
saveforecast = "yes"
sim_limit = 300.0
rerun = 1
runfile = "runfile"
datatype = ["WOPR PRO1", "WOPR PRO2", "WOPR PRO3", "WWPR PRO1", "WWPR PRO2",
"WWPR PRO3", "WWIR INJ1", "WWIR INJ2", "WWIR INJ3"]
parallel = 4
startdate = "1/1/2022"
Loading

0 comments on commit ace6ec2

Please sign in to comment.