for field_name in digitised_test_fields:
t, var = problem1_test.mSuite.resultsList[run_index].\
getFieldHistoryAtCell(field_name, obs_cell_index)
plt.semilogx(t, var / scale[field_name], '-', label='Waiwera', zorder=4)
t, var = AUTOUGH2_result.getFieldHistoryAtCell(field_name, obs_cell_index)
plt.semilogx(t, var / scale[field_name], 's', label='AUTOUGH2', zorder=3)
for sim in digitised_simulators:
result = digitised_result[sim]
t, var = result.getFieldHistoryAtCell(field_name, obs_cell_index)
plt.semilogx(t, var / scale[field_name], symbol[sim], label=sim)
t, var = analytical_result.getFieldHistoryAtCell(field_name,
obs_cell_index)
plt.semilogx(t,
var / scale[field_name],
'k--',
label='analytical',
zorder=1)
plt.xlabel('time (s)')
plt.ylabel(field_name + ' (' + unit[field_name] + ')')
plt.legend(loc='best')
plt.title(' '.join((model_name, field_name.lower(), 'results at', obspt)))
img_filename_base = '_'.join((model_name, obspt, field_name))
img_filename_base = img_filename_base.replace(' ', '_')
img_filename = os.path.join(problem1_test.mSuite.runs[run_index].basePath,
problem1_test.mSuite.outputPathBase,
img_filename_base)
getFieldHistoryAtCell(field_name, obs_cell_index)
plt.plot(t / yr,
var / scale[field_name],
'-',
label='Waiwera',
zorder=3)
t, var = AUTOUGH2_result.getFieldHistoryAtCell(field_name,
obs_cell_index)
plt.plot(t[::3] / yr,
var[::3] / scale[field_name],
's',
label='AUTOUGH2',
zorder=2)
t, var = ref_result.getFieldHistoryAtCell(field_name, obs_cell_index)
plt.plot(t / yr,
var / scale[field_name],
symbol[ref_sim],
label=ref_sim,
zorder=1)
plt.xlabel('time (years)')
plt.ylabel(field_name + ' (' + unit[field_name] + ')')
plt.legend(loc='best')
plt.title(' '.join(
(model_name, field_name.lower(), 'results at depth', obspt, 'm')))
img_filename_base = '_'.join((model_name, obspt, field_name))
img_filename_base = img_filename_base.replace(' ', '_')
img_filename = os.path.join(
problem4_test.mSuite.runs[run_index].basePath,