AkB2010 extended and adjusted for use in the new Swiss Hazard model [2014]

AkB2010 extended and adjusted for use in the new Swiss Hazard model [2014]
I'm not sure about how you implemented this method. 'C1_rrup' should depend on rrup, but depending on the if case it is sometimes a scalar and sometimes an array