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

Compare with Current View Page History

« Previous Version 16 Next »

Introduction

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

    plot_vs30.py sim_folder

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

  • ObsTimeSeriesOnMap

    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

    1.Create (or copy from template) a config file: post-processing_config.cfg

    post-processing_config.cfg
    [files]
outputDirectory = /home/nesi00213/RealTime/20161213_Kaikoura_m7p8_201705301501/GM/Validation/
srfFile = /home/nesi00213/RealTime/20161213_Kaikoura_m7p8_201705301501/Src/Model/20161213_Kaikoura_m7p8/t2/Ian02_s103245.srf
stationFile = /home/nesi00213/RealTime/20161213_Kaikoura_m7p8_201705301501/Stat/20161213_Kaikoura_m7p8/20161213_Kaikoura_m7p8.ll
velSimDir = /home/nesi00213/RealTime/20161213_Kaikoura_m7p8_201705301501/GM/Sim/Data/BB/Vel/
dbLocation = /home/nesi00213/RealTime/20161213_Kaikoura_m7p8_201705301501/GM/Validation/database2.db
siteClassFile = /nesi/projects/nesi00213/RunFolder/leer/Pre-processing/RegionalSeismicityTectonics/Sept42010GmMetadata.ll
velObsDir = /home/nesi00213/RealTime/20161213_Kaikoura_m7p8_201705301501/GM/Obs/Data/velBB
runName = 20161213_Amberley_m7p8_201705301501

[plotGMPE]
DistMin = 1.0
DistMax = 200
nval = 51
V30measured = False
SiteClassForPrediction = D
rake = 0
Rtvz = 0
Mw = 7.8
Vs30 = 1500,800,450,250,200
Ztor = 0
model = 1
dip = 82

[plotSeismogram]
normalize = True # used to be 1.0
magnify = 3.0

[flag]
plotGMPE = 1
plotIM = 3
plotSeismogram = 3
plotTimeSeries = 3
calcIM = 3

[plotTimeSeries]
gapStations = 5
nLimit = 5
tmax = 300
stationsToPlot= HAVS,DAVS,WRCS,PWES,BMTS,MOLS,NWFS,WANS,WTES,MOTS,NELS,WSTS,NBSS,WDFS,NEWS,SOMS,PAPS,KEKS,TEPS,LIRS,WNKS,FAIS,PGMS,SEVS,KHLS,PIPS,BTWS,UHSS,WNAS,WTYS,SEAS,HSES,NNZ,HSSS,TRMS,LRSS,MKBS,SEDS,TRTS,KIKS,WEL,WDAS,WAZ,MISS,UHCS,VUWS,LNBS,MATS,KIRS,POTS,TAIS,MKVS,WEMS,BOWS,BWRS,POKS,PHHS,PVCS,MAVS,PTOS,THOB
subdir = waveforms
gapObsSim = 2

[calcIM]
IMsToCalculate = AI,CAV,Ds575,Ds595,PGA,PGV,pSA
calculateRRups = True

[plotIM]
ImName = PGA,PGV,pSA
saveData = True
component = Geom
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.
    3.edit the stationsToPlot - to the list of real station with both obs and sim data.
    4.cd /home/nesi00213/post-processing/examples
    5. python 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