View Source

SUMMARY

Old HF workflow
- losing 10% core hours for formatting text
- losing 250% storage space for storing txt
New BB workflow
- checks LF/HF compatibility, previous case of different `nt` no longer undetected
- kupe 80 processes = 55 seconds, hypocentre 28 processes = 15 seconds (siteamp not fully vectorised)
- quickly finished, not fully tested yet
- No longer losing 600% storage space for storing both Vel and Acc as text

How to run

HF

mpirun -n 80 python2 hf_sim.py source.stoch stations.ll HF.bin \ # HF.bin is the output file
-m /home/nesi00213/VelocityModel/Mod-1D/Cant1D_v1-midQ.1d \ # override vm
--dt 0.02
-i # individually run stations

BB

mpirun -n 80 python2 bb_sim.py LF/bevan2012_v3_s103252/OutBin /path/to/vm/ HF.bin BB.bin # BB.bin is the output file

BB gets hf_vs_ref from HF file (stored as part of HF), lf_vs_ref from vs3d.file.s

Python interface

HF (and BB)

from qcore.timeseries import HFSeis
# load file
hf = HFSeis('HF.bin')
# retrieve timeseries (x,y,z)
hf.acc('CACS')
# only x
hf.acc('CACS', comp = <0, or hf.X or hf.COMP['090']>)
# store station as txt
hf.acc2txt('CACS')
# store all stations as txt, same as if ran standard txt version
hf.all2txt(prefix = 'Acc/hf_')

LF

# just like for vel
lf.vel('CACS')
# you can access transparently converted accelerations
lf.acc('CACS')