Orthogonal convergent velocity

Q&A
#1

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

#2

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

#3

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

#4

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)