Page History

...

• Source_map

  plot_srf_map.py src_model.srf

  Code Block
  language bash
  $plot_srf_map.py Ian02_s103245.srf

• Source_square

  plot_srf_square.py src_model.srf

  Code Block
  language bash
  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
  *uses non_uniform_grid and prune out the virtual stations manually if you are making document for future event.

  
  plot_vs30.py root_simulation_dir

  Code Block
  language bash
  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

  Code Block
  language bash
  plot_obs.py .

• PNG_timeslices

  plot_ts.py processes

  Code Block
  language bash
  plot_ts.py 30

  note: if files missing, try looking at the realtime sim wikipage referencing "plot_transfer.py auto"

  Code Block
  plot_transfer.py auto folder_path_on_fitz local_path user_on_fitzroy $plot_transer.py auto /nesi/projects/nesi00213/RunFolder/baes/AlpineFault_Rupture_VMSI_v1p65_200m-h0p200_EMODv3p0p4_170531_s2n/ . ykh22

  
  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 time slices.

  Code Block
  title 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 the movie file.
  * you will need at least 10 time slices for the document

  Code Block
  title 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

  The below command needs to be run from the same folder as params_base.py is in. It is normally located in $root_dir/GM/Sim/Data

  plot_ts_sum.py

  Code Block
  language bash
  plot_ts_sum.py

...