Versions Compared

Old Version 45

changes.mady.by.user Jason Motha

Saved on

compared with

New Version 46

changes.mady.by.user Jason Motha

Saved on

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

There are 4 types of SRF files that can be created, all SRF files are created directly or indirectly via createSRF.py.

Type 1point source. appears as 1 plane with 1 subfault.
Type 2finite fault. a single plane created from point source attributes, it has an are determined by a magnitude scaling relation.
Type 3finite fault. a single plane created with specified parameters.
Type 4finite fault. a multiple plane fault. planes can be internally generated together and cut up (like plane extensions that continue at a different strike angle) or:
separate cases where multiple SRFs are combined. where separate cases are used, each can have a delay time to synchronise rupture initiation time between cases.

Common Attributes (always available)

attributedescription
type1: point source, 2: finite fault from point source info, 3: single plane, 4: multiple plane
dtslip time step (s)
rakerake angles (type 1-3: single value, type 4: 2D list for cases, fault planes)
magmagnitude (type 1-3: single value, type 4: list for each case)
hlonhypocentre longitude for first fault (degrees)
hlathypocentre latitude for first fault (degrees)
hdepthhypocentre depth for first fault (km)
corners4 corners (longitude, latitude) for each plane. derived parameter (shape of array is nplane, 4 corners, 2 lonlat values).

SRF Plane Attributes (lists with a length matching number of fault planes)

naming taken from qcore/srf.py

attributedescription
centretop centre longitude, latitude
nstrikenumber of subfaults along strike
ndipnumber of subfaults along dip
lengthlength (along strike) of fault plane (km)
widthwidth (along dip) of fault plane (km)
strikestrike angle (degrees)
dipdip angle (degrees)
shyphypocentre location along strike from top centre (km), -999.9 for plane extentions
dhyphypocentre location along dip from top edge (km), -999.9 for plane extentions
dtopdepth of top edge (km)
dbottomdepth to bottom of plane, derived parameter (km)

Type 1 (point source) Attributes

attributedescription
vsvs at centroid depth (km/s)
rhodensity

Type 2 (finite fault from point source data) Attributes

attributedescription
mwsrmagnitude scaling relation name

Type 2+ (finite faults) Attributes

attributedescription
vmvelocity model used (basename)
shyp0distance along strike (km) to hypocentre for each case from top left corner (strike = 0, dip = 0)
dhyp0distance along dip (km) to hypocentre for each case from top left corner (strike = 0, dip = 0)

Optional Attributes (need to check if they exist)

attributedescription
tect_typetectonic type, taken from the NHM file in nhm2srf.py. known options:
ACTIVE_SHALLOW
VOLCANIC
SUBDUCTION_INTERFACE

LF/HF/BB binary format

These files store timeseries data. All formats follow a style derived from the LF seis format produced by EMOD3D:

...

HF and BB have gaps between the first 2 sections to allow future additions to section 1 without breaking backwards compatibility.

LFHFBB

i4 number of stations

TOTAL 4 BYTES

i4 number of stations
i4 number of timesteps
i4 seed used in simulation
i4 site amplification used (bool)
i4 path duration method

  • 0: GP2010 formulation
  • 1: WUS modification trial/error
  • 2: ENA modification trial/error
  • 11: WUS formulation of BT2014
    over predicts for multiple rays
  • 12: ENA modification trial/error
    over predicts for multiple rays

i4 number of ray methods
known options for first-fourth type below:

  • 1: direct
  • 2: moho

i4 first ray method used
i4 second ray method used
i4 third ray method used
i4 fourth ray method used
i4 nbu parameter
i4 ift parameter
i4 nlskip parameter
i4 icflag parameter (bool)
i4 individually run stations (bool)
i4 site specific VMs used (bool)
f4 duration of timeseries (s)
f4 timestep of timeseries (s)
f4 start time of timeseries (s)
f4 stress drop average (bars)
f4 kappa parameter
f4 q frequency exponent
f4 max sim frequency (Hz)
f4 flo parameter (Hz)
f4 fhi parameter (Hz)
f4 rupture velocity factor (rupture : Vs)
f4 rvfac shallow fault multiplier
f4 rvfac deep fault multiplier
f4 czero coefficient, -1: used binary default
f4 calpha coefficient, -1: used binary default
f4 seismic moment, -1: used rupture model
f4 rupture velocity, -1: used rupture model
f4 depth to moho, -1: used binary default (999.9)
f4 vp_sig parameter
sig parameters are likely uncertainties
f4 vsh_sig parameter
f4 rho_sig parameter
f4 qs_sig parameter
f4 fourier amplitude uncertainty (1)
f4 fourier amplitude uncertainty (2)
f4 rupture velocity uncertainty
s64 stoch file used (basename)
s64 velocity model used (basename)

TOTAL 288 BYTES

i4 number of stations
i4 number of timesteps
f4 duration of timeseries (s)
f4 timestep of timeseries (s)
f4 start time of timeseries (s)
s256 LF directory path used
s256 LF VM directory path used
s256 HF file path used
possibly add vsite file path used here?

TOTAL 788 BYTES

START OFFSET 4 BYTES

i4 index of station in input file
i4 x gridpoint of station
i4 y gridpoint of station
i4 z gridpoint of station
i4 simulation number of timesteps
f4 simulation timestep (s)
f4 simulation grid spacing (km)
f4 simulation grid rotation (degrees)
f4 latitude of station (degrees)
f4 longitude of station (degrees)
s8 name of station

TOTAL 48 BYTES * NUM_STATIONS

START OFFSET 512 BYTES

f4 longitude of station (degrees)
f4 latitude of station (degrees)
s8 name of station
f4 epicentre distance to station (km)
f4 vs30 at station (m/s)

TOTAL 24 BYTES * NUM_STATIONS

START OFFSET 1280 BYTES

f4 longitude of station (degrees)
f4 latitude of station (degrees)
s8 name of station
i4 x gridpoint of station
i4 y gridpoint of station
i4 z gridpoint of station
f4 epicentre distance to station (km)
f4 HF

vs30

vs30ref (m/s)
f4 LF

vs30

vs30ref (m/s)
f4 BB vs30 (m/s)

TOTAL 44 BYETS * NUM_STATIONS

START OFFSET 0 FROM ABOVE

f4 velocity (cm/s) timeseries in array dimensions:
timestep, station, component (x, y, z, ..., 9)

TOTAL 4 BYTES * PRODUCT_OF_DIMENTIONS

START OFFSET 0 FROM ABOVE

f4 acceleration (cm/s^2) timeseries in array dimensions:
station, timestep, component (x, y, z)

TOTAL 4 BYTES * PRODUCT_OF_DIMENTIONS

START OFFSET 0 FROM ABOVE

f4 acceleration (g) timeseries in array dimensions:
station, timestep, component (x, y, z)

TOTAL 4 BYTES * PRODUCT_OF_DIMENTIONS

XYTS.e3d binary format

This file is produced by EMOD3D and contains a timeseries of ground motions on the XY plane. Unlike the LF seis files, this contains data at all grid points and may have a decimated resolution specified when running EMOD3D through the e3d.par file with the parameters dxts and dyts.

...