Page History

...

1. Source generation (Viktor/Hoby)
   1. Srf Generation

             - Manually extract location, dip,strike, rake from Mw tensor (Ristau) into setSrfParams.py

             - Change TYPE=2 (point source to finite fault)

             - Run createSRF.py

        b. Plot Srf

             -  change  Change  SRF name and output name in plotSRF.py

             - Run  Run plotSRF.py

        c. Plot Srf on Map (For verification purpose)

              /nesi/projects/nesi00213/gm_sim_workflow/devel/plot_srf.py file.srf

   
   

2. Velocity model generation
   1. Ethan
3. Observed ground motions
   1. Ahsan
      

      NOTE: Installation is required only once, unless there is a need to update. geoNet
      package is self contained and should work on your pcs and laptops too.

      0) INSTALL geoNet.
      From 
      https://github.com/ucgmsim/Pre-processing/tree/master/geoNet
      clone geoNet and follow the installation instructions.

      1) GET THE REQUIRED SCRIPTS.
      From geoNet/examples copy

      getData.py
      event_statsll.py
      processData.py
      plot_accvel.py
      plot_psa.py

      to your working directory

      2) GET THE DATA FOR THE EARTHQUAKE.

      In getData.py set BASE_URL to the location to the geoNet URL e.g
      BASE_URL="ftp://ftp.geonet.org.nz/strong/processed/Proc/2017/02_Feb/2017-02-01_102129/Vol1/data/"

      save and run.

      3) GET THE LIST OF STATIONS FOR THE EARTHQUAKE
      Run event_statsll.py which creates event_stats.ll and places it
      in your working directory.

      4)PROCESS THE GEONET VOL1 DATA TO GET NORTH, EAST AND VERTICAL COMPONENTS
      Run processData.py which filters and rotates the SMS data placing them
      in Vol1/data/accBB etc.
      Note: It is better to redirect the std output to a file like so
      python processData.py > processData_output

      5) CREATE PLOTS FOR OBSEREVED DATA
      Run plot_accvel.py and plot_psa.py

4. Ground motion simulation computation (Sung)
   1. Upload Rupture model to /nesi/projects/nesi00213/RupModel/Realtime/MonYEAR where MonYear is like Dec2016, Jan2017. Under this directory, you need Srf and Stoch.
   2. Upload Velocity model to /nesi/projects/nesi00213/VelocityModels/. Currently we have Canterbury or South_Island. Unless it is a Canterbury specific model, place it under South_Island. We will be doing some clean up later.
      Your velocity model should include rho3dfile.d, vp3dfile.p, vs3dfile.s, params_vel.py and outputs from gen_cords (gridfile, gridout, model_bounds, model_coords, model_params).
   3. Upload station files to /nesi/projects/nesi00213/StationInfo. This should include geonet_stations_xxx.ll, .vs30 and .vs30ref.
   4. Go to /nesi/projects/nesi00213/RunFolder and execute ./install.sh. and follow the instruction.
   5. During install.sh, it will advise you expected wallclock time. Enter a sensible wall clock time.
   6. Follow the remaining instruction.
      
      
5. simulation post-processing
   1. Viktor/Hoby
   2. General Plotting Instructions:
      run on hypocentre, not fitzroy
      make sure the scripts are in PATH (hypocentre):
      export PATH=$PATH:/nesi/projects/nesi00213/gm_sim_workflow/devel
   3. Plotting Station based data such as IMs, Vs30, Observed PGA, pSA:
      plot_stations.py datafile.ll
      Run on hypocentre, datafile has 6 line header as described at the top of plot_stations.py
   4. Plotting Timeslice based data (includes PGV, MMI, timeslice animation) on hypocentre:
      plot_transfer.py auto /nesi/projects/nesi00213/RunFolder/baes/RT2017Jan17_VMSI_20170117_drill_200m-h0p200_EMODv3p0p4_170116 .
      will copy the remote fitzroy folder to the current directory (.) and start plotting using defaults
      run plot_ts_sum.py to re-run PGV and MMI plots
      run plot_ts.py to re-run timeslice image creation 

...