Disaggregation calculator

The Disaggregation Calculator allows calculating contributions (in terms of probabilities) to a specific level of seismic hazard (represented by a ground motion value corresponding to a probability of exceedence in given time span), in terms of:

 - latitude and longitude of the surface projection of the closest point of
  the rupture area

  - magnitude

  - epsilon

  - tectonic region type (active shallow crust, stable tectonic, subduction
  interface, subduction intraslab, volcanic)

The implementation of the disaggregation calculator requires the definition of a new CALCULATION_MODE = DISAGGREGATION

Given a region or a set of site where to compute disaggregation results, the following new parameters must be defined:

LATITUDE_BIN_LIMITS = list of latitude values representing bin limits used for disaggregation, values must be in increasing order
LONGITUDE_BIN_LIMITS = list of longitude values representing bin limits for disaggregation, values must be in increasing order
MAGNITUDE_BIN_LIMITS = list of magnitude values representing bin limits for disaggregation, values must be in increasing order
EPSILON_BIN_LIMITS = list of epsilon values representing bin limits for disaggregation, values must be in increasing order
DISTANCE_BIN_LIMITS = list of distance values representing bin limits for disaggregation, values must be in increasing order
POES_TO_DISAGGREGATE = list of probabilities of exceedence to disaggregate

The input for the disaggregation calculator is represented by a source model logic tree and a gmpe logic tree (represented by the keys SOURCE_MODEL_LOGIC_TREE_FILE and GMPE_LOGIC_TREE_FILE)

The parameters needed to compute the earthquake rupture forecast from a source model and to set the GMPE parameters must be also defined:

MINIMUM_MAGNITUDE
INVESTIGATION_TIME
MAXIMUM_DISTANCE
WIDTH_OF_MFD_BIN

COMPONENT
INTENSITY_MEASURE_TYPE
PERIOD
DAMPING
INTENSITY_MEASURE_LEVELS
GMPE_TRUNCATION_TYPE
TRUNCATION_LEVEL
STANDARD_DEVIATION_TYPE
REFERENCE_VS30_VALUE
REFERENCE_DEPTH_TO_2PT5KM_PER_SEC_PARAM
SADIGH_SITE_TYPE

INCLUDE_AREA_SOURCES
TREAT_AREA_SOURCE_AS
AREA_SOURCE_DISCRETIZATION
AREA_SOURCE_MAGNITUDE_SCALING_RELATIONSHIP

INCLUDE_GRID_SOURCES
TREAT_GRID_SOURCE_AS
GRID_SOURCE_MAGNITUDE_SCALING_RELATIONSHIP

INCLUDE_FAULT_SOURCE
FAULT_RUPTURE_OFFSET
FAULT_SURFACE_DISCRETIZATION
FAULT_MAGNITUDE_SCALING_RELATIONSHIP
FAULT_MAGNITUDE_SCALING_SIGMA
RUPTURE_ASPECT_RATIO
RUPTURE_FLOATING_TYPE

INCLUDE_SUBDUCTION_FAULT_SOURCE
SUBDUCTION_FAULT_RUPTURE_OFFSET
SUBDUCTION_FAULT_SURFACE_DISCRETIZATION
SUBDUCTION_FAULT_MAGNITUDE_SCALING_RELATIONSHIP
SUBDUCTION_FAULT_MAGNITUDE_SCALING_SIGMA
SUBDUCTION_RUPTURE_ASPECT_RATIO
SUBDUCTION_RUPTURE_FLOATING_TYPE

A prototype of the disaggregation calculator has been implemented and can be found here: https://github.com/monellid/shaCalcs/blob/master/src/DisaggregationCalculator.java
A test class checking the scientific validity of the calculator is here:
https://github.com/monellid/shaCalcs/blob/master/src/DisaggregationCalculatorTest.java

For a given a site of interest, the disaggregation calculator computes probability mass functions (PMF) for:
- magnitude
- distance
- tectonic region type
- magnitude- tectonic region type
- latitute-longitude
- latitude-longitude-magnitude
- latitude-longitude-epsilon
- latitude-longitude-magnitude-epsilon
- latitude-longitude-tectonic region type

The above listed results must be calculated and saved for each site of interest and for each logic tree sample. The option to save the full disaggregation matrix must be also given.

Aggregated results must be also produced. From the set of results produced by the logic tree sampling, the mean disaggregation results must be computed and also the disaggregation results associated to quantile levels of the ground motion values.