|
|
|
@ -103,11 +103,13 @@ if PLOT:
|
|
|
|
y = stats.cauchy(loc=row['loc'], scale=row['scale']).pdf(x)
|
|
|
|
y = stats.cauchy(loc=row['loc'], scale=row['scale']).pdf(x)
|
|
|
|
ax.plot(x, y * row['scale'], c=c[j])
|
|
|
|
ax.plot(x, y * row['scale'], c=c[j])
|
|
|
|
if j % 2 == 0:
|
|
|
|
if j % 2 == 0:
|
|
|
|
ax.annotate(f' {i}', (x[y.argmax()], 1),
|
|
|
|
ax.annotate(f' {i}', (x[y.argmax()], 0.32),
|
|
|
|
ha='center',
|
|
|
|
ha='center',
|
|
|
|
va='bottom',
|
|
|
|
va='bottom',
|
|
|
|
|
|
|
|
clip_on=False,
|
|
|
|
rotation=90)
|
|
|
|
rotation=90)
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ax.set_ylim(top=0.35)
|
|
|
|
ax.set_xlabel('SLR (m)', labelpad=10)
|
|
|
|
ax.set_xlabel('SLR (m)', labelpad=10)
|
|
|
|
ax.set_ylabel('Probability density (-)', labelpad=10)
|
|
|
|
ax.set_ylabel('Probability density (-)', labelpad=10)
|
|
|
|
|
|
|
|
|
|
|
|
|
