This page provides rough outlines on the current progress on Cybershake in NZ.
Note that his page is for internal sharing purposes only, it is likely inaccurate and out-of-date, and therefore it is advised to conduct researchers directly if you want more reliable information
| Version | Num sources | SRF generation | VM domain | VM generation | Simulation specs | Estimated core-hours | Actual core hours | Notes |
|---|---|---|---|---|---|---|---|---|
| v17.8 | 15 (dominant in Canterbury) | vX.XX rupgenerator 3 hypo and 2 slip dist per source | automated based on PGV>5cm/s; 15kmRup, 5km land cutoff Default depth and duration scaling | 0.4km regular grid, Vs_min=500 | Transition freq = 0.25Hz | ~3k Fitzroy | First implementation; Focus on running workflow and comparison with empirically-derived hazard curves | |
| v17.9 (planning) | ~100 (dominant in Canterbury, but most of Sth Island) | hypo every X km along strike, Y slip dist per source | automated for PGV>5cm/s ... Improved depth scaling | as for v17.8 | as for v17.8 | ~30k Fitzroy | Focus on extending number of sources and srf uncertianties | |
| v17.10 | ~100 (Canterbury) | ? | as for v17.9 | 0.2km regular grid, Vs_min=500 | Transition frequency = 0.5 Hz | ~170k Fitzroy | Focus on improved freq transition in South Island | |
| v17.11 | ~540 (all in NZ) | ? | as for v17.9 | as for v17.8 | as for v17.8 | ~500k Fitzroy(?) | Focus on producing NZ-wide output at lower freq transition |
Verification process:
Fault rupture (SRF) files:
- Source plot on the map
- Length-dependent hypocenter location (done)
- Different slip distribution realizations (done)
Velocity model (VM) files:
- Velocity model plot on the map (done)
- Rotation and reduction in size to minimize the ocean coverage (done)
- Generate a code to check VM files (done)
- Check if given VM parameter (folder path) exists
- .p, .s, .d and params_vel.py files must exist
- warns if model_params/coords/bounds etc don't exist using params_vel sufx in filename
- params_vel matching hh/xlen vs nx for x, y, z
- file size = nx * ny * nz * 4 bytes, checked for .p, .s, .d
- if numpy available: checks first xz slice for >0 and not NaN in .p, .s, .d
Simulation process:
- Verification code for param files for SRFs and VMs (e.g., transition frequency, grid size) (in-progress)
- Config file for batch submission of the simulations (in-progress)
- Wall clock and number of cores estimations (in-progress)
- Automated check to see if the runs for LF, HF, and BB are done successfully (in-progress)
Post-processing:
- Automated IM plot on the map for all simulations (in-progress)
- Automated standard post-processing for all simulations (in-progress)
- Hazard curve calculations from simulated ground motions (in-progress)
- Empirical hazard calculations (in-progress)
- Comparison and combination of the simulated and empirical hazard analyses (in-progress)
Prioritization:
To do now:
- Automated IM plot on the map for all simulations (in-progress)
- Automated standard post-processing for all simulations (in-progress)
- Hazard curve calculations from simulated ground motions (in-progress)
- Empirical hazard calculations (in-progress)
To do before v17.9:
- Config file for batch submission of the simulations (in-progress)
- Verification code for param files for SRFs and VMs (e.g., transition frequency, grid size) (in-progress)
- Wall clock and number of cores estimations (in-progress)
- Automated check to see if the runs for LF, HF, and BB are done successfully (in-progress)
- Comparison and combination of the simulated and empirical hazard analyses (in-progress)
To do, but no timeline for completion:
- ...