Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

  • 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.cfg

    Code Block
    titlepost-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_documentation_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.py

    Code Block
    plot_stations.py stat_ratios_PGV_20161213_Amberley_m7p8_201705301501.ll [path_to_output_directory_optional]IMs # this will write the images to IMs relative directory
    plot_stations.py stat_ratios_pSA_20161213_Amberley_m7p8_201705301501.ll [path_to_output_directory_optional] pSA # this will write the images to pSA relative directory

    Note: new formate of config will have IM instead of PGV, 6 file instead of 1. pick the first one and ignore the rest.

...