Hi, I am not sure which versions of gmpe-smtk and OQ engine you are using, but for me it worked (slightly modified below - cannot attach due to company policies).
# Import the modules
%matplotlib inline
import numpy as np
import smtk.trellis.trellis_plots as trpl
# Set up the configuration
gmpe_list = ["AkkarBommer2010", "AkkarCagnan2010", "AkkarEtAlRjb2014", "BooreAtkinson2008", "ChiouYoungs2008", "ZhaoEtAl2006Asc"]
imts = ["PGA", "SA(0.2)", "SA(1.0)", "SA(2.0)"]
params = {"ztor": 5.0, # Top of rupture depth
"hypo_depth": 10.0, # Hypocentral depth
"distance_type": "rjb",
"vs30": 800.0, # Vs30 for all sites
"vs30measured": True, # Vs30 value is measured
"z1pt0": 100.0, # Depth (m) to the 1.0 km/s Vs interface
"dip": 90.0, # Vertical Fault
"rake": 0.0 # Strike-slip fault
}
magnitudes = 6.5
distances = {"repi": np.arange(0.0, 151.0, 1.0)}
distances["rhypo"] = np.sqrt(distances["repi"] ** 2.0 + params["hypo_depth"] ** 2)
distances["rjb"] = distances["repi"]
distances["rrup"] = np.sqrt(distances["rjb"] ** 2.0 + params["ztor"] ** 2)
distances["rx"] = distances["rjb"]
tp=trpl.DistanceIMTTrellis(magnitudes, distances, gmpe_list, imts, params,
distance_type="rjb", plot_type="loglog", dpi=400, figure_size=(20,20),
filename="distance_imt_trellis_simple.pdf",
filetype="pdf")
tp.plot()
Produced both a figure in jupyter and the pdf.
Peter