Added handling of missing SES in seismicity forecast
The seismicity forecasts provided as input may not contain events for all realisations for which the forecast was run. This is the case when the software used to produce the seismicity forecast is set to output earthquakes above a certain threshold that results in some stochastic event sets being empty and, thus, non-existent in the forecast file. However, the user might want to take into account these "empty" or "missing" stochastic event sets when calculating the mean damage and losses associated with the whole forecast, as these stochastic event sets can be considered as not having produced any damage or losses.
The present MR addresses this issue by allowing the user to input two parameters in the config.yml
file that control the IDs of the stochastic event sets to be taken into account:
-
continuous_ses_numbering
: True/False (i.e. use IDs defined with consecutive numbering as specified inses_range
). -
ses_range
: ID_min, ID_max.
E.g., if continuous_ses_numbering=True
and ses_range=1, 10
, the software will consider 10 stochastic event sets with IDs [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]
.
Closes https://git.gfz-potsdam.de/cnievas/real-time-loss-tools/-/issues/10.