The purpose of this bug to describe the calculation of the occurrence rate of earthquakes on a fault using the Anderson & Luco Type 1 method.
There are three types of methods to calculate the occurrence rate made by Anderson & Luco (1983). They share the same set of input parameters, and differ only in the equation used for the calculation.
The parameters used are:
1) Rate of slip (mm/yr) on the fault (positive float)
2) Aseismic coefficient - a positive float in the range 0.0 to 1.0 that describes the proportion of the fault slip that is not "used" by earthquakes.
3) Rake (degrees, i.e. positive float in the range -180.0 to 180.0) of the fault
4) Fault config, which contains the Shear modulus, magnitude scaling relation model and displacement to length ratio (see tectonic region)
5) Fault geometry ( https://bugs.launchpad.net/openquake/+bug/975104)
6) Minimum magnitude (compulsory - input by the user)
7) Magnitude interval (bin width)
8) b-value
9) Maximum magnitude: this can be input directly by the user. However, if it is not specifed by the user then the maximum magnitude is calculated by the magnitude scaling relation(WC1994 get_median_mag or get_mag), which takes as input the magnitude scaling relation model (implemented from the nhlib) selection, the fault geometry area and the rake.
A possible implementation is defined for calculating the occurence rate:
https://github.com/g-weatherill/gw_misc_code/blob/master/fault_calculator/new_fault_class.py#L257 see also:
https://github.com/g-weatherill/gw_misc_code/blob/master/fault_calculator/new_fault_class_2.py
