You are viewing an old version of this page. View the current version.

Compare with Current View Page History

« Previous Version 25 Next »

Introduction

Creating MetaData

  • Go you the simulation folder on Fitzroy

  • Run the genCSV.py

    python /nesi/projects/nesi00213/workflow/genCSV.py -location Christchurch -datetime 2010-09-04-04:30:00 -exetime 201705301501 -mag 6.8

    -datetime is the time when the earthquake happened, and must have the formate : %Y-%M-%D-hh-mm-ss
    -exetime is when the simulation is ran, and must fit the formate: %Y%M%%Dhhmm, e.g. 201705301501

    Note: In default , the meta data should be located under CSV folder. But if you are creating metadata for other user's simulation you can pass an extra argument to genCSV for it to save at specific place.

    python /nesi/projects/nesi00213/workflow/genCSV.py -location Christchurch -datetime 2010-09-04-04:30:00 -exetime 201705301501 -mag 6.8 -save ~/metadata/20100904

Images needed

  • Source_map

    plot_srf_map.py src_model.srf

    $plot_srf_map.py Ian02_s103245.srf

  • Source_square

    plot_srf_square.py src_model.srf

    plot_srf_square.py Ian02_s103245.srf

  • Vs30StationMap

    note: this step requires observation generated, please follow the wiki page referring How to create observation files
    plot_vs30.py sim_folder

    plot_vs30.py /home/nesi00213/RealTime/20161213_Kaikoura_m7p8_201705301501

  • ObsTimeSeriesOnMap

    note: this step requires observation generated, please follow the wiki page referring How to create observation files
    plot_obs.py sim_folder

    plot_obs.py .

  • PNG_timeslices

    plot_ts.py processes

    plot_ts.py 30

    note: if files missing, try looking at the realtime sim wikipage referencing "plot_transfer.py auto"
    important: the original plot_ts.py will combine all ts.png into a mov file and delete all the png files. MUST modify the code so that it will not delete the files.

    Comment out the rm line
    # temporary files can be quite large
rmtree(gmt_temp)
#rmtree(png_dir)

    *below is a example if you just want to plot specific ts without producing PGV and MMI.

    Plot specific TS
    pool = mp.Pool(processes)
#pool.map(render_slice, xrange(xyts.t0, xyts.nt - xyts.t0))
pool.map(render_slice, [120,330,550,700,800,910,1000,1200,1500,1600])
print('FINISHED TIMESLICE SEGMENT IN %.2fs' % (time() - ts0))
print('AVERAGE TIMESLICE TIME %.2fs' % \ 
        ((time() - ts0) / (xyts.nt - xyts.t0)))
# images -> animation
#make_movie('%s/ts%%04d.png' % (png_dir), \
#        os.path.join(base_dir, 'animation.mov'), fps = 20)
# temporary files can be quite large
rmtree(gmt_temp)
#rmtree(png_dir)

  • PGV and MMI

    plot_ts_sum.py

    plot_ts_sum.py
  • Spacial Distribution

 

  • 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

    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/20161213_Kaikoura_m7p8/t2/Ian02_s103245.srf
stationFile = /home/nesi00213/RealTime/YYYYMMDD_EVENT_NAME_mXpX_YYYYMMDDHHmm/Stat/20161213_Kaikoura_m7p8/20161213_Kaikoura_m7p8.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/20161213_Kaikoura_m7p8/20161213_Kaikoura_m7p8.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 = 120
Rtvz = ERROR:PLEASE_SET
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.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.
    3. python /home/nesi00213/post-processing/examples/run_default_plots.py path_to_config_file

    run_default_plots.py /home/nesi00213/RealTime/20161213_Kaikoura_m7p8_201705301501/GM/Validation/post-processing_config.cfg

Running publish

  • cd GM/Sim

  • create of copy the config file: publish.cfg

    publish.cfg
    [event]
title:13 Nov 2016 Mw7.8 Kaikoura
location:/home/nesi00213/RealTime/20161213_Amberley_m7p8_201705301501
csv:/home/nesi00213/RealTime/20161213_Amberley_m7p8_201705301501/GM/Sim/MetaData
output: 20161213_Amberley_m7p8_201705301501

[image]
SrfMapView:Src/Figures/source_map.png
srfSlipRiseTimeRake:Src/Figures/source_square.png
#
Vs30StationMap:VM/Figures/Vs30StationMap.png
ObsTimeSeriesOnMap:GM/Obs/Figures/ObsTimeSeriesOnMap.png
#
timesliceLocation:GM/Sim/Figures/PNG_timeslices/
timesliceFiles:ts0120.png,ts0330.png,ts0550.png,ts0700.png,ts0800.png,ts0910.png,ts1000.png,ts1200.png,ts1500.png

validationLocation:GM/Validation/Figures/IM/IM_vs_source_site_distance/
validationFiles:PGA_comp_geom_with_Rrup_with_GMPE_with_Station_ClassesVel.png,pSA(0p2)_comp_geom_with_Rrup_with_GMPE_20161213_Amberley_m7p8_201705301501.png,pSA(0p5)_comp_geom_with_Rrup_with_GMPE_20161213_Amberley_m7p8_201705301501.png,pSA(1)_comp_geom_with_Rrup_with_GMPE_20161213_Amberley_m7p8_201705301501.png,pSA(3)_comp_geom_with_Rrup_with_GMPE_20161213_Amberley_m7p8_201705301501.png,pSA(10)_comp_geom_with_Rrup_with_GMPE_20161213_Amberley_m7p8_201705301501.png

pSAWithPeriod:GM/Validation/Figures/IM/spectra_ratio/pSAWithPeriod_compgeom_20161213_Amberley_m7p8_201705301501.png
#missing
spatialLocation:GM/Validation/Figures/ratios
spatialFiles:PGV/p000.png,pSA/p000.png,pSA/p001.png,pSA/p002.png,pSA/p003.png,pSA/p004.png
#
simPgv:GM/Sim/Figures/PNG_tssum/PGV.png
simMmi:GM/Sim/Figures/PNG_tssum/MMI.png
#
responseSpectraLocation:GM/Validation/Figures/IM/PSA
responseSpectraFiles:pSA_comp_geom_vs_Period_Vel_AMBC.png,pSA_comp_geom_vs_Period_Vel_CMHS.png,pSA_comp_geom_vs_Period_Vel_DALS.png,pSA_comp_geom_vs_Period_Vel_FTPS.png,pSA_comp_geom_vs_Period_Vel_GLWS.png,pSA_comp_geom_vs_Period_Vel_INGS.png,pSA_comp_geom_vs_Period_Vel_KOWC.png,pSA_comp_geom_vs_Period_Vel_MCAS.png,pSA_comp_geom_vs_Period_Vel_NCBS.png,pSA_comp_geom_vs_Period_Vel_WVFS.png
#
waveformLocation:GM/Validation/Figures/waveforms
waveformFiles:velts_ADCSAKSSAMBCAPPSARKS_2016Nov13_Ian02_s103245_VMCant_Amberly_200m-h0p200_EMODv3p0p4_161221.png,velts_SMTCSOCSSOMSSPRSSTAS_2016Nov13_Ian02_s103245_VMCant_Amberly_200m-h0p200_EMODv3p0p4_161221.png
  • Edit all path to the correct directory

  • Run the publish script

    python /home/nesi00213/gm_publish/publish.py publish.cfg

 

  • No labels