...
Post-processing related plots
note: this step requires observation generated, please follow the wiki page referring How to create observation files
1.Create (or copy from template) a config file: post-processing_config.cfgCode Block title post-processing_config.cfg [files] outputDirectory = /home/nesi00213/RealTime/YYYYMMDD_EVENT_NAME_mXpX_YYYYMMDDHHmm/GM/Validation/ srfFile = /home/nesi00213/RealTime/YYYYMMDD_EVENT_NAME_mXpX_YYYYMMDDHHmm/Src/Model/YYYYMMDD_EVENT_NAME_mXpX/t2/Ian02_s103245.srf stationFile = /home/nesi00213/RealTime/YYYYMMDD_EVENT_NAME_mXpX_YYYYMMDDHHmm/Stat/YYYYMMDD_EVENT_NAME_mXpX/YYYYMMDD_EVENT_NAME_mXpX.ll velSimDir = /home/nesi00213/RealTime/YYYYMMDD_EVENT_NAME_mXpX_YYYYMMDDHHmm/GM/Sim/Data/BB/Vel/ dbLocation = /home/nesi00213/RealTime/YYYYMMDD_EVENT_NAME_mXpX_YYYYMMDDHHmm/GM/Validation/database.db siteClassFile = /nesi/projects/nesi00213/Pre-processing/RegionalSeismicityTectonics/Sept42010GmMetadata.ll vs30File = /home/nesi00213/RealTime/YYYYMMDD_EVENT_NAME_mXpX_YYYYMMDDHHmm/Stat/YYYYMMDD_EVENT_NAME_mXpX/YYYYMMDD_EVENT_NAME_mXpX.vs30 velObsDir = /home/nesi00213/RealTime/RealTimeYYYYMMDD_EVENT_NAME_mXpX_YYYYMMDDHHmm/GM/Obs/Data/velBB runName = RealTimeYYYYMMDD_EVENT_NAME_mXpX_YYYYMMDDHHmm [plotGMPE] DistMin = 1.0 DistMax = 100 nval = 51 vs30Default = 250.0, rake = ERROR:PLEASE_SET Rtvz = 0 Mw = ERROR:PLEASE_SET Ztor = ERROR:PLEASE_SET model = 1 dip = ERROR:PLEASE_SET [plotSeismogram] normalize = True magnify = 3.0 [calcIM] plotGMPE = True recalculateGMPE = False recalculateIMs = False calculateRRups = True calcIM = 3 [plotTimeSeries] gapStations = 5 nLimit = 5 tmax = 200 stationsToPlot = [] subdir = waveforms gapObsSim = 2 [plotIM] ImName = PGA,pSA,PGV,Ds575,Ds595,AI,CAV pSAStations = [] pSAPlotCommonAxisLimits = False unknown = D pSAPeriod = 0.1,0.2,0.5,1.0,3.0,10.0 siteClassMarkers = False
2.Edit all the file directory to the correct folder. – it will error if you don't change all of the paths above.
3a) If you are running a RealTime simulation:python /home/nesi00213/post-processing/examples/run_default_plots.py path_to_config_file
Code Block python /home/nesi00213/post-processing/examples/run_documentationdefault_plots.py /home/nesi00213/RealTime/20161213_Kaikoura_m7p8_201705301501/GM/Validation/post-processing_config.cfg
3b) If you are running a future event simulation:
python /home/nesi00213/post-processing/examples/run_simulation_plots.py path_to_config_file
Code Block python /home/nesi00213/post-processing/examples/run_simulation_plots.py /home/nesi00213/RealTime/20161213_Kaikoura_m7p8_201705301501/GM/Validation/post-processing_config.cfg
* important: must edit the stationsToPlot if simulation used non-uniform grid. Otherwise it will plot all the virtual stations
* for future events, which mean there is not observation data
Spatial Distribution
note:this must be done after post-processing
1.cd GM/Validation/Data. There should exist two .ll files PGV and pSA
2.plot the files using plot_stations.pyCode Block plot_stations.py stat_ratios_PGV_20161213_Amberley_m7p8_201705301501.ll IMs # this will write the images to IMs relative directory plot_stations.py stat_ratios_pSA_20161213_Amberley_m7p8_201705301501.ll pSA # this will write the images to pSA relative directory
Note: new formate format of config will have IM instead of PGV, 6 file instead of 1. pick the first one and ignore the rest.
...