This is now included in the PR you just approved. Two main things to be fixed on plotting:
color palette should differ between subplots (gives better contrast between subplots)
there's nothing in the legend to indicate that the blue strips are fault locations. Would be good to include that in the legend. I'm not super familiar with the nitty gritty of matplotlib, so not sure how to just add something to the legend.
Bonus: considering changing this to something like seaborn, which has nicer high-level plotting features
@christopherDT thanks for finding this...
https://stackoverflow.com/questions/3681872/shading-an-area-between-two-points-in-a-matplotlib-plot
Final solution would just be one combined plot:
https://gist.github.com/christopherDT/129883b6dda40216aabd8a8b74beb818