Merge pull request #193 from johnmaxwilson/master

Fixed major error in mesher interpolation between trace points
geodynamics · Dec 7, 2017 · fe9d0b8 · fe9d0b8
2 parents 717b4b3 + 1ffcd85
commit fe9d0b8
Showing 1 changed file with 8 additions and 1 deletion.
diff --git a/quakelib/src/QuakeLibIO.cpp b/quakelib/src/QuakeLibIO.cpp
@@ -647,7 +647,9 @@ void quakelib::ModelWorld::create_section(std::vector<unsigned int> &unused_trac
 
         double          inner_t;
         // And the element and inner_t corresponding to it
-        spline.get_element(cur_t, next_elem_ind, inner_t);
+        spline.get_element(cur_t, cur_elem_ind, inner_t);
+
+        next_elem_ind = cur_elem_ind + 1;
 
         // Mark the element corresponding to this point as having been used
         unused_trace_pts.erase(next_elem_ind);
@@ -659,10 +661,15 @@ void quakelib::ModelWorld::create_section(std::vector<unsigned int> &unused_trac
 
 
         elem_slip_rate = conv.cm_per_yr2m_per_sec(inner_t *trace.at(next_elem_ind).slip_rate()+(1.0-inner_t)*trace.at(cur_elem_ind).slip_rate());
+
         elem_aseismic = inner_t *trace.at(next_elem_ind).aseismic()+(1.0-inner_t)*trace.at(cur_elem_ind).aseismic();
+
         elem_dip = conv.deg2rad(inner_t *trace.at(next_elem_ind).dip()+(1.0-inner_t)*trace.at(cur_elem_ind).dip());
+
         elem_rake = conv.deg2rad(inner_t *trace.at(next_elem_ind).rake()+(1.0-inner_t)*trace.at(cur_elem_ind).rake());
+
         elem_lame_mu = inner_t *trace.at(next_elem_ind).lame_mu()+(1.0-inner_t)*trace.at(cur_elem_ind).lame_mu();
+
         elem_lame_lambda = inner_t *trace.at(next_elem_ind).lame_lambda()+(1.0-inner_t)*trace.at(cur_elem_ind).lame_lambda();
 
         // Ensure square elements, vertical number for best fit.