It is recommended to run this on NeSI to avoid excessive waiting time in the queue.
The supplied GeoNet_CMT_solutions.csv has the event info, where the bold-faced triplet is the (strike1, dip1, rake1) we will be using to generate the source model from GCMT
PublicID,Date,Latitude,Longitude,strike1,dip1,rake1,strike2,dip2,rake2,ML,Mw,Mo,CD,NS,DC,Mxx,Mxy,Mxz,Myy,Myz,Mzz,VR,Tva,Tpl,Taz,Nva,Npl,Naz,Pva,Ppl,Paz,Method
....
2023p707798,20230919211400,-43.6988,171.0682,63,56,156,172,47,46,6.1,5.6,3.20e+24,18,12,85,5193.00,13350.00,13600.00,-27870.00,-929.10,22670.00,62,30676.89,59,11,2464.09,31,205,-33140.98,6,111,2
This triplet should be consistent with GeoNet Moment Tensor solution published by John Ristau
Best Double-Couple
Mo = 3.26e+24 dyne.cm Mw = 5.6
NP1: Strike = 167 Dip = 70 Rake = 36
NP2: Strike = 63 Dip = 56 Rake = 156
Go to /nesi/nobackup/nesi00213/RunFolder/{your username}. Let's create a directory RangitataGorge, and go to that directory.
The latest GeoNet_CMT_solution.csv can be downloaded from here: https://raw.githubusercontent.com/GeoNet/data/main/moment-tensor/GeoNet_CMT_solutions.csv
wget https://raw.githubusercontent.com/GeoNet/data/main/moment-tensor/GeoNet_CMT_solutions.csv
You will also need
cp /nesi/project/nesi00213/VelocityModel/Mod-1D/Cant1D_v3-midQ_OneRay.1d . cp /nesi/project/nesi00213/StationInfo/non_uniform_whole_nz_with_real_stations-hh400_v20p3_land.* .
Let's create a file list.txt that just contains the event PublicID, which is 2023p707798
$ cat list.txt 2023p707798
We will generate Type 2 (point source finite fault) source model,
python $gmsim/Pre-processing/srf_generation/source_parameter_generation/generate_realisations_from_gcmt.py list.txt GeoNet_CMT_solutions.csv 2 --vel_mod_1d Cant1D_v3-midQ_OneRay.1d --cybershake_root .
This should be very quick.
Let's install the workflow.
python $gmsim/workflow/workflow/automation/install_scripts/install_cybershake.py . list.txt 22.2.2.1 --keep_dup_station --stat_file_path non_uniform_whole_nz_with_real_stations-hh400_v20p3_land.ll
If a nicely interpolated distribution of IM plot map is desired, it is best to use a non-uniform grid instead of geoNet stations list. The non-uniform grid contains most real stations, but double check it contains the stations you are interested in for validation purposes. Also bear in mind that HF, BB and IM calculation can be significantly slower as non-uniform grid contains a lot of virtual stations,
We need to execute run_cybershake.py twice.
For the first run (=pre-processing step), let's create a task_config_prepro.yaml. For details on the meaning of ALL, MEDIAN , have a look at this page: Changes to automated workflow from the inclusion of median events
The complete pipeline of workflow can be found here: Automated workflow pipelines
SRF_GEN: ALL INSTALL_REALISATION: ALL VM_PARAMS: MEDIAN VM_GEN: MEDIAN INSTALL_FAULT: MEDIAN VM_PERT: NONE # ALL # REL_ONLY NO_VM_PERT: ALL
python $gmsim/workflow/workflow/automation/execution_scripts/run_cybershake.py . ./task_config_prepro.yaml
When VM gen is completed, check the VM domain with the PNG image produced in Data/VMs directory.
When all the input models are successfully produced, run it again with task_config.yaml configured as below.
EMOD3D: MEDIAN HF: MEDIAN BB: MEDIAN IM_calc: MEDIAN clean_up: MEDIAN merge_ts: MEDIAN plot_ts: MEDIAN IM_plot: MEDIAN rrup: NONE LF2BB: NONE HF2BB: NONE plot_srf: NONE advanced_IM: NONE
python $gmsim/workflow/workflow/automation/execution_scripts/run_cybershake.py . ./task_config.yaml
This took about 5 hours on NeSI, including plot_ts (2D animation) and IM plots under verification directory.
└── verification
├── 2023p707798.m4v
└── IM_PLOT
├── geom
│ ├── im_order.txt
│ ├── non_uniform_im
│ │ ├── plot_items_0.png
│ │ ├── plot_items_10.png
│ │ ├── plot_items_11.png
│ │ ├── plot_items_12.png
│ │ ├── plot_items_13.png
│ │ ├── plot_items_14.png
│ │ ├── plot_items_15.png
...
Congratulations. You have completed your first ground motion simulations.
Have a look at plots of PGV, PGA, pSA etc (done by IM plot) and 2D animation (done by plot_ts) under verification folder.
Output data from each step can be also found under folders representing each step of simulation (eg. BB, IM_calc etc.)
GMSimViz animation
If you wish to produce a GMSimViz animation, copy SRF, XYTS file to RCC's visualisation VM (10.195.0.37) using rccuser account and run the following command (written by Jonney)
cd /data/ykh22/Atlas_animation/tmp mkdir RangitataGorge (and place SRF and XTRS files) /data/ykh22/Atlas_animation/scripts/run_viz.sh RangitataGorge/2023p707798.srf RangitataGorge/2023p707798_xyts.e3d "Rangitata Gorge Mw5.6"
By default, the video will show SeisTech logo on the top-left corner, and QuakeCoRE logo on the top-right corner.
The logo files are at /data/GMSimViz/gmsimviz/seistech-logo.png and quakecore-logo-logo.png.
Make symbolic links to display the logo or remove one not to show it.
ln -s quakecore-logo.png logo-right.png lrwxrwxrwx 1 rccuser rccuser 18 Sep 22 09:18 logo-right.png -> quakecore-logo.png
The rendering can take 4~5 hours (but worth the effort!)
The final video can be viewable from here: https://www.dropbox.com/scl/fi/rfmi9fl3khdd71vegz14y/2023p707798_viktor_v3.m4v?rlkey=tyt97jf5ny49ndkg8uj0mrxi7&st=n3nto3y2&dl=0