Skip to content
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.

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

Issues with restarting #202

Open
ArMOdr opened this issue Aug 26, 2022 · 1 comment
Open

Issues with restarting #202

ArMOdr opened this issue Aug 26, 2022 · 1 comment

Comments

@ArMOdr
Copy link

ArMOdr commented Aug 26, 2022

Hi
I am trying to re-start my PEST run. I have these lines to my pst file.

++ ies_restart_observation_ensemble(03.0.obs.csv)
++ ies_restart_parameter_ensemble(03.0.par.csv)

Both 03.0.obs.csv and 03.0.par.csv look fine.

When I try to restart, I get the below error. Any idea what the issue could be?

Cheers

--- initializing ---
...using REDSVD for truncated svd solve
...maxsing: 1000000
...eigthresh: 1e-06
...initializing localizer
...using automatic adaptive localization
with autoadaloc_sigma_dist 1
...using lambda multipliers: 1 , 5 , 10 ,
...using lambda scaling factors: 0.5 , 0.8 , 0.95 ,
...acceptable phi factor: 1.05
...lambda increase factor: 10
...lambda decrease factor: 0.75
...max run fail: 1

--- sanity_check errors - uh oh ---
...ies_par_en is empty but ies_restart_par_en is not - how can this work?
...ies_par_en is empty but ies_restart_obs_en is not - how can this work?


IterEnsembleSmoother error: sanity_check() found some problems - please review rec file
@jtwhite79
Copy link
Collaborator

The restart for ies is confusing and a lot of it has to do with allowing par and obs (and noise) realizations that have different names. The short answer is you need to pass ies_parameter_ensemble(03.0.par.csv), ies_observation_ensemble(03.obs+noise.csv) and ies_restart_observation_ensemble(03.0.obs.csv). You dont need to pass ies_restart_parameter_ensemble in most cases (its only needed if you want to track parameter changes from the initial/prior ensemble, but you are restarting with the prior ensemble, so its irrelevant).

Also, just a heads up: if your control file is named "03.pst", then when you restart, ies will stomp on "03.0.obs.csv", "03.obs+noise.csv" and "03.0.par.csv"

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants