When tessellating subduction zones, GPlately now only looks for subduction polarity (not subduction feature type). This means features of type OrogenicBelt also get processed, when only SubductionZone should (by default).
I had mistakenly assumed only a subduction zone would have a subductionPolarity. But I see the GPGIM allows an OrogenicBelt to have one too - it says “Similar to a subduction zone, one plate has often overridden the other, so we use a similar model to that of gpml:SubductionZone”. So I guess there’s potentially an overriding plate there too (but maybe not in all cases).
Solution is to add a parameter to enable sampling all subducting boundaries (including orogenic belts) and have it default to False (to sample only SubductionZones).
When tessellating subduction zones, GPlately now only looks for subduction polarity (not subduction feature type). This means features of type
OrogenicBeltalso get processed, when onlySubductionZoneshould (by default).I had mistakenly assumed only a subduction zone would have a subductionPolarity. But I see the GPGIM allows an OrogenicBelt to have one too - it says “Similar to a subduction zone, one plate has often overridden the other, so we use a similar model to that of gpml:SubductionZone”. So I guess there’s potentially an overriding plate there too (but maybe not in all cases).
Solution is to add a parameter to enable sampling all subducting boundaries (including orogenic belts) and have it default to False (to sample only SubductionZones).