Gives details into what plots are available within cybershake_investigation and their corresponding scripts with inputs

Plot Examples

plot_avg_event_residual_spectra

Produce a single plot involving all events for a cybershake run.
This plot will show the residual spectra for each event across all sites as shown above.
Note: Only shows the lines where 50% of the line is above or below the 16th or 84th percentile to reduce clutter


python cybershake_investigation/scripts/plot_avg_event_residual_spectra.py
    <path to the root directory of simulation data> (Can be either a simulation or archive structure)
    <path to the empirical directory for comparison>
    <path to the output directory for plot>
	[--sites <The sites to use during computing, defaults to all sites>]
    [--component <Default rotd50>]
 	[--faults <The faults to plot>] 
 	[--from_archive <Specify this flag if the root directory is from an archive directory>]
    [--n_procs <To use multiprocessing for computing, default 1>]



plot_avg_site_residual_spectra

produce a single plot involving all events for a cybershake run.
This plot will show the residual spectra for each site across all events as shown above.
Note: Only shows the lines where 50% of the line is above or below the 16th or 84th percentile to reduce some clutter


python cybershake_investigation/scripts/plot_avg_site_residual_spectra.py
    <path to the root directory of simulation data> (Can be either a simulation or archive structure)
    <path to the empirical directory for comparison>
    <path to the output directory for plot>
	[--sites <The sites to plot>]
    [--component <Default rotd50>]
 	[--faults <The faults to use during computing, defaults to all faults>] 
 	[--from_archive <Specify this flag if the root directory is from an archive directory>]
    [--n_procs <To use multiprocessing for computing, default 1>]

plot_avg_site_residual

Creates Full Spatial plots that show the log ratio between the empirical results and the IM csvs produced by simulations averaged across all events for all sites across NZ.

python cybershake_investigation/scripts/plot_avg_site_residual.py
    <path to the root directory of simulation data> (Can be either a simulation or archive structure)
    <path to the empirical directory for comparison>
	<path to the qcore/data directory for map plot data>
	<path to the station ll file>
    <path to the output directory for the plots>
    [--component <Default rotd50>]
 	[--ims <The ims to plot, defaults to all ims>] 
    [--faults <The faults to limit when extracting data, defaults to all faults>] 
 	[--from_archive <Specify this flag if the root directory is from an archive directory>]
    [--high_res <Specify this flag if you want high res plots, will take significantly longer for a large amount of plots>]
    [--n_procs <To use multiprocessing for computing, default 1>]

plot_spatial_residual

          


Produces ratio plots based on cybershake results.
Produces a plot per Fault, Realisation, IM specified in the desired output path

python cybershake_investigation/scripts/plot_spatial_residual.py
    <path to the root directory of simulation data> (Can be either a simulation or archive structure)
    <path to the empirical directory for comparison>
	<path to the qcore/data directory for map plot data>
	<path to the station ll file>
	<path to the nhm file>
    <path to the output directory for the plots>
    [--component <Default rotd50>]
 	[--ims <The ims to plot, defaults to all ims>] 
    [--faults <The faults to plot, defaults to all faults>] 
 	[--rels <Only produce these realisation plots. Can be set to 'all', default of None>]
 	[--no_combined <Specify this flag to not produce an average plot over all realisations>]   
    [--from_archive <Specify this flag if the root directory is from an archive directory>]
    [--high_res <Specify this flag if you want high res plots, will take significantly longer for a large amount of plots>]
    [--n_procs <To use multiprocessing for computing, default 1>]

plot_spatial_im


Produces IM Spatial plots based on cybershake results.
Produces a plot per Fault, Realisation, IM specified in the desired output path for empirical or simulation data.

python cybershake_investigation/scripts/plot_spatial_im.py
    <path to the root directory of simulation data> (Can be either a simulation or archive structure)
    <path to the empirical directory for comparison>
	<path to the qcore/data directory for map plot data>
	<path to the station ll file>
	<path to the nhm file>
    <path to the output directory for the plots>
    [--component <Default rotd50>]
 	[--ims <The ims to plot, defaults to all ims>] 
    [--faults <The faults to plot, defaults to all faults>] 
 	[--rels <Only produce these realisation plots. Can be set to 'all', default of None>]
 	[--no_combined <Specify this flag to not produce an average plot over all realisations>]  
	[--plot_emp <Specify this flag to produce an empirical plot along with all simulation plots that get plotted>]   
    [--from_archive <Specify this flag if the root directory is from an archive directory>]
    [--high_res <Specify this flag if you want high res plots, will take significantly longer for a large amount of plots>]
    [--n_procs <To use multiprocessing for computing, default 1>]

plot_hist_residual


Produce histogram plots based on residual values between cybershake and empirical results.
Produces a plot per Realisation, IM specified in the desired output path.

python cybershake_investigation/scripts/plot_hist_residual.py
    <path to the root directory of simulation data> (Can be either a simulation or archive structure)
    <path to the empirical directory for comparison>
    <path to the output directory for the plots>
    [--component <Default rotd50>]
 	[--ims <The ims to plot, defaults to all ims>] 
    [--faults <The faults to plot, defaults to all faults>] 
	[--rels <Only produce these realisation plots. Can be set to 'all', default of None>]
 	[--no_combined <Specify this flag to not produce an average plot over all realisations>] 
    [--from_archive <Specify this flag if the root directory is from an archive directory>]
    [--n_procs <To use multiprocessing for plotting, default 1>]

plot_single_site_residual

Produce a single site residual histogram plot based on cybershake results.
Produces a plot per Site, IM specified in the desired output path

python cybershake_investigation/scripts/plot_single_site_residual.py
    <path to the root directory of simulation data> (Can be either a simulation or archive structure)
    <path to the empirical directory for comparison>
    <path to the output directory for the plots>
	<The sites to plot>
    [--component <Default rotd50>]
 	[--ims <The ims to plot, defaults to all ims>] 
    [--faults <The faults to plot, defaults to all faults>] 
    [--from_archive <Specify this flag if the root directory is from an archive directory>]
    [--n_procs <To use multiprocessing for plotting, default 1>]

plot_single_site_spectrum

Produces a single site spatial histogram plot based on cybershake results.
Produces a plot per Site, Fault specified in the desired output path

python cybershake_investigation/scripts/plot_hist_residual.py
    <path to the root directory of simulation data> (Can be either a simulation or archive structure)
    <path to the empirical directory for comparison>
    <path to the output directory for the plots>
 	<The sites to plot> 
    [--component <Default rotd50>]
    [--faults <The faults to plot, defaults to all faults>] 
    [--from_archive <Specify this flag if the root directory is from an archive directory>]
    [--n_procs <To use multiprocessing for plotting, default 1>]


  • No labels