Orthogonal convergent velocity

Hi Gplates community,

I’m working on a research project that require to use pyGplate to calculate the orthogonal convergent velocity for a subducting plate toward a trench line. I find “subduction_convergence.py” might be the right one to use, but I do not know how, like how to define the subducting plate id and the trench line id in it.

Could anyone help me out?

Thanks,
Jeremy

Do you have the rotation and topology files? You need them to use subduction_convergence.py. The plate ids should be assigned to subducting plates and trench in the topology file. You can find a topology example(DynamicPolygons) in the sample data GPlates - Browse /gplates/2.2 at SourceForge.net.

@john.cannon

Hi!

Yes, I have the rotation and topology files from the original install package. I will dp the rotation file from Torsvik et al., 2019 as well.

Thank you for your kindly Reply!
Jeremy

Here are some example code. The results will be put in folder “result”.

rotation_files = ['my_rotation.rot']
topology_files = ['my_topology_1.gpml', ['my_topology_2.gpml']

kwargs = {    
    'output_distance_to_nearest_edge_of_trench':True,
    'output_distance_to_start_edge_of_trench':True,
    'output_convergence_velocity_components':True,
    'output_trench_absolute_velocity_components':True,
    'output_subducting_absolute_velocity':True,
    'output_subducting_absolute_velocity_components':True}

return_code = subduction_convergence_over_time(
        'result/subdata',
        'csv',
        rotation_files,
        topology_files,
        math.radians(0.2),
        0,
        100,
        1,
        1,
        0,
        output_gpml_filename = None,
        **kwargs)