Run this script get_fault_numrels.py with 2 input files as below. all_fault_list.txt contains name of 484 faults
Some faults listed in input 1 are in subduction zones and are not listed in input 2. which are omitted in the output file.
These omitted faults are : ['HikRaukmax', 'HikHBaymin', 'Lachlan1-2', 'HikHBaymax', 'HikRaukmin', 'Wairaka05', 'HikWgtnmax', 'HawkeBay6-12', 'Waimarama1-2', 'Waimarama3-4', 'FiordSZ09', 'HawkeBay5-11', 'FiordMar1-2', 'FiordSZ03', 'HikWgtnmin']
This is the list of fault and realisations we used for running Cybershakes : CS_fault_numrels_jason.txt
We noted the output produced by Python code above is not identical to this file. 21 faults have 1 more/less realisations - for the consistency, we are sticking to this file.
The following script can be used to calculate how many simulation jobs are needed for each fault given a model based on rupture magnitude.
The file names to be used for input and output can be found at the top and modified as needed.This can be changed to take command line input if required.
Works with both Python 2 and Python 3, but Python 3 is preferred.
This file provides the same output for the matlab file found here, given the input files provided on the page.
Dependencies: matplotlib
Inputs: A whitespace delimited file with source names and corresponding magnitudes, and a source list file with one source per line.
Outputs: A whitespace delimited file with source names and corresponding number of simulation counts and an image displaying the number of simulations for a given magnitude.
Run with
python3 run_C18p6_NumSim.py
Given the input files found on the Step 2 page
CS_list_v18p6.txtFault_LeonardMag_Length_Width_Rake.txt
We get the following output files, which are identical to the output when they are run with the Matlab script
fault_selection_c18p6_for_SRF.txt