plot_subduction_teeth

[thomas-schouten]

Hi,

I am using GPlately to plot my tectonic reconstruction. I downloaded GPlately today, so I have the latest version of GPlately and pyGPlates. I have a topological closed plate polygon with subduction zones loaded as a pgyplates feature collection. All subduction zone segments have a polarity. I am able to plot trenches, ridges, faults without problem. However, I am unable to plot the subduction teeth, even though I receive no error message.

Sample code: fig, ax = plt.subplots(subplot_kw={'projection': cartopy.crs.PlateCarree()}, figsize=(12, 8)) ax.set_extent([50, 105, -10, 40], crs=cartopy.crs.PlateCarree()) plate_model = gplately.PlateReconstruction(rotation_model, plates) plot_plates = gplately.plot.PlotTopologies(plate_model, time) plot_plates.time = i plot_plates.plot_ridges_and_transforms(ax, color="r") plot_plates.plot_trenches(ax, color="b") plot_plates.plot_faults(ax, color="k") plot_plates.plot_subduction_teeth(ax, color="k")

What could be the issue?

Many thanks, Thomas Schouten

[brmather]

Hi Thomas,

Good on you for being an early adopter of gplately. So that we can diagnose your issue, can you check the following:

• Does your topology features (in your case I assume this is plates) contain “subductionPolarity”?
• Does plot_subduction_teeth work for other plate models on your machine?
• Have you tried adjusting the scale parameters?
• What is your OS and Python version?

Thanks!

[michaelchin]

The size of triangles is in meters. But the lines are in degrees. The mismatch in measure units caused no triangle has been plotted.

[michaelchin]

fixed in 1d60b76

[brmather]

@michaelchin we need to update the notebooks in case the subduction teeth have changed. Perhaps @lauren-ilano @laurilano can re-run the notebooks to make sure everything works...?