Skip to content

Commit 0ba479c

Browse files
author
bodlukas
committed
Changed Figures 1 and 3 and docs
1 parent 54c69c0 commit 0ba479c

6 files changed

+5661
-5648
lines changed

.gitignore

+5-1
Original file line numberDiff line numberDiff line change
@@ -1 +1,5 @@
1-
.ipynb_checkpoints
1+
.ipynb_checkpoints
2+
3+
__pycache__/
4+
*.py[cod]
5+
*$py.class

data/intro_virtualregion.csv

+5,625-5,625
Large diffs are not rendered by default.

figures.ipynb

+23-17
Large diffs are not rendered by default.

introductory_example_estimation.ipynb

+1-1
Large diffs are not rendered by default.

main_estimate_postprocess.py

+3-1
Original file line numberDiff line numberDiff line change
@@ -11,7 +11,9 @@
1111
1212
For a thorough explanation of the estimation and
1313
post-processing computations check out the following notebook:
14-
...
14+
15+
introductory_example_estimation.ipynb
16+
1517
This can be easily imported to colab and you don't have
1618
to set up a special environment locally.
1719

utils_plotting.py

+4-3
Original file line numberDiff line numberDiff line change
@@ -127,7 +127,7 @@ def fig1_e_to_h(fig, grid_spec, df, ref_idx, kE, kEA, kEAS):
127127
lev = np.arange(0,1.2,0.2)
128128

129129
x_ref, y_ref = (df.x[ref_idx], df.y[ref_idx])
130-
x_event, y_event = (0, 0)
130+
x_event, y_event = (-5, 10)
131131
for i, val in enumerate(vals):
132132
ax = fig.add_subplot(grid_spec[0,i])
133133
if i == 0:
@@ -146,7 +146,7 @@ def fig1_e_to_h(fig, grid_spec, df, ref_idx, kE, kEA, kEAS):
146146
ax.set_xlabel('Easting [km]')
147147
if i==0: ax.set_ylabel('Northing [km]')
148148
if i==1:
149-
ax.text(0, 3.5, 'epicenter', fontsize=8, ha='center')
149+
ax.text(-5, 6.5, 'epicenter', fontsize=8, ha='center')
150150
ax.annotate("reference site",
151151
xy=(x_ref, y_ref-0.5), xycoords='data',
152152
xytext=(x_ref, -11.3), textcoords='data',
@@ -256,7 +256,8 @@ def fig3(axss, dfs, titles):
256256
ax.set_xlabel(r'Lengthscale $\ell_{\mathrm{E}}$ [km]')
257257
ax.set_xlim([4, 46])
258258
ax.set_ylim([0.26, 1.24])
259-
ax.set_title(titles[i], fontsize=9)
259+
ax.set_title(titles[0][i], fontsize=9)
260+
ax.text(0.5, 0.98, titles[1][i], fontsize=9, transform=ax.transAxes, ha='center', va='top')
260261
if i==1:
261262
ax.text(0.68, 0.3, r'sample $\mathbf{\psi}_{\mathrm{E},r}$',
262263
transform=ax.transAxes,

0 commit comments

Comments
 (0)