@@ -90,7 +90,6 @@ def detect_plates(step, vrms_surface, fids, time):
9090 pom2 , np .less , order = trench_span , mode = 'wrap' )[0 ]
9191 trench = ph_coord [argless_dv ]
9292 velocity_trench = vph2 [argless_dv , indsurf ]
93- dv_trench = dvph2 [argless_dv ]
9493
9594 # finding ridges
9695 pom2 = np .copy (dvph2 )
@@ -111,9 +110,7 @@ def detect_plates(step, vrms_surface, fids, time):
111110 if argdel :
112111 print ('deleting from ridge' , trench , ridge [argdel ])
113112 ridge = np .delete (ridge , np .array (argdel ))
114- arggreat_dv = np .delete (arggreat_dv , np .array (argdel ))
115113
116- dv_ridge = dvph2 [arggreat_dv ]
117114 if 'age' in conf .plates .plot :
118115 agefld = step .fields ['age' ].values [0 , :, :, 0 ]
119116 age_surface = np .ma .masked_where (agefld [:, indsurf ] < 0.00001 ,
@@ -134,7 +131,7 @@ def detect_plates(step, vrms_surface, fids, time):
134131 agetrench [itrench ]
135132 ))
136133
137- return trench , ridge , agetrench , dv_trench , dv_ridge
134+ return trench , ridge , agetrench
138135
139136
140137def plot_plate_limits (axis , ridges , trenches ):
@@ -364,7 +361,7 @@ def main_plates(sdat):
364361
365362 time = step .time * vrms_surface * \
366363 conf .scaling .ttransit / conf .scaling .yearins / 1.e6
367- trenches , ridges , agetrenches , _ , _ = \
364+ trenches , ridges , agetrenches = \
368365 detect_plates (step , vrms_surface , fids , time )
369366 _write_trench_diagnostics (step , time , trenches , agetrenches , fids )
370367
0 commit comments