[cig-commits] [commit] master: Bug fix (cdf74ff)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/vq

On branch  : master
Link       : https://github.com/geodynamics/vq/compare/f26b7fa26e56a9309aa9658f7441dfd63922dee6...0aca7bf57ea151a65788e0106b736313ec48b926

>---------------------------------------------------------------

commit cdf74ffd6a9cf50896017662eeca2433dc7a0638
Author: Eric Heien <emheien at ucdavis.edu>
Date:   Fri Jan 9 14:50:52 2015 -0800

    Bug fix


>---------------------------------------------------------------

cdf74ffd6a9cf50896017662eeca2433dc7a0638
 quakelib/src/QuakeLibIO.cpp | 9 +--------
 1 file changed, 1 insertion(+), 8 deletions(-)

diff --git a/quakelib/src/QuakeLibIO.cpp b/quakelib/src/QuakeLibIO.cpp
index b8f507a..3715afe 100644
--- a/quakelib/src/QuakeLibIO.cpp
+++ b/quakelib/src/QuakeLibIO.cpp
@@ -1632,21 +1632,17 @@ int quakelib::ModelWorld::read_files_eqsim(const std::string &geom_file_name, co
 
             new_element.set_lame_lambda(friction_data.get_lame_lambda());
             new_element.set_lame_mu(friction_data.get_lame_mu());
-            //new_element.set_static_strength(friction_data.get_static_strength(it->first));
-            //new_element.set_dynamic_strength(friction_data.get_dynamic_strength(it->first));
+            new_element.set_max_slip(0);    // Set a temporary maximum slip of 0 (this will be changed below)
 
             // Insert partially finished element into the eqsim_world
             eqsim_world.insert(new_element);
 
             // Compute area of the current element, add it to the total for this section
             fault_areas[sit->second.sid()] += eqsim_world.create_sim_element(new_element.id()).area();
-
-
         }
 
         // Go through the created elements and assign maximum slip based on fault section area
         for (eit=eqsim_world.begin_element(); eit!=eqsim_world.end_element(); ++eit) {
-
             // From Table 2A in Wells Coppersmith 1994
             double moment_magnitude = 4.07+0.98*log10(conv.sqm2sqkm(fault_areas[eit->section_id()]));
 
@@ -1655,10 +1651,7 @@ int quakelib::ModelWorld::read_files_eqsim(const std::string &geom_file_name, co
 
             // Set the max slip for the current element
             eit->set_max_slip(max_slip);
-
         }
-
-
     }
 
     insert(eqsim_world);