To receive notifications about scheduled maintenance, please subscribe to the mailing-list gitlab-operations@sympa.ethz.ch. You can subscribe to the mailing-list at https://sympa.ethz.ch

Commit 605ef7ca authored by Feliks Kiszkurno's avatar Feliks Kiszkurno
Browse files

added reporting of chi2 to invert_data.py

parent 87ebb58a
......@@ -32,12 +32,20 @@ def invert_data(profile_name, *, lambda_param=20, z_weight_param=0.2):
model_inverted = ert_manager.invert(lam=lambda_param, paraDX=0.25, paraMaxCellSize=2, zWeight=z_weight_param, # paraDepth=2 * max_depth,
quality=34, zPower=0.4)
log_file_name = settings.settings['log_file_name']
log_file = open(os.path.join(settings.settings['results_folder'], log_file_name), 'a')
log_file.write('\n')
log_file.write('Inverted profile: {}'.format(profile_name))
log_file.write('\n')
log_file.write('Last chi2 = {}'.format(ert_manager.inv.chi2()))
log_file.write('\n')
result_full = ert_manager.inv.model
result_array = result_full.array()
limits = pg.utils.interperc(ert_manager.inv.model, trimval=25.0)
limits = [50, 2000]
limits = [4, 50]
result_array[result_array <= limits[0]] = limits[0]
result_array[result_array >= limits[1]] = limits[1]
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment