wphase:documentation
Differences
This shows you the differences between two versions of the page.
| Both sides previous revision Previous revision Next revision | Previous revision | ||
|
wphase:documentation [2021/10/22 18:49] wphase [Dependencies] |
wphase:documentation [2022/01/10 07:39] (current) wphase [Extracting data from a gzipped tar archive] |
||
|---|---|---|---|
| Line 97: | Line 97: | ||
| Two tgz archive formats are supported: WILBER (default) and NIED. | Two tgz archive formats are supported: WILBER (default) and NIED. | ||
| - | Using the WILBER format, SAC files should be named with a ".SAC" or ".sac" extension and SAC_PZ file names should start with "SAC_PZ". Although this is the default format, extraction from a WILBER archive can be stated explicitly by running | + | Using the WILBER format, the tgz archive should include SAC files named with a ".SAC" or ".sac" extension and SAC_PZ file names starting with "SAC_PZ". If data is ordered using the [[https://ds.iris.edu/wilber3|IRIS Wilber 3 system]], this format can be obtained directly by selecting "Output format: SAC binary" (little or big-endian) and "Bundle As: gzipped tar archive" in the final Request Data page (see below). |
| - | <code>${WPHASE_HOME}/bin/tgz2DATA_org.py path_to_tgz_archive.tgz wilber</code> | + | |
| + | {{:wilber3_screenshot.png?300|{{:wphase:wilber3_screenshot.png?300|}} | ||
| + | |||
| + | Once the tgz archive is prepared, extraction from a WILBER formatted archive is done by running: | ||
| + | <code>${WPHASE_HOME}/bin/tgz2DATA_org.py path_to_tgz_archive.tgz wilber</code> | ||
| + | As this is the default format, the "wilber" option can be removed from the command line. | ||
| Using the NIED tgz archive format, SAC_PZ files should end with ".zp" and SAC file names can have any format (but should not end with ".zp"). Extraction from a NIED archive can be done using | Using the NIED tgz archive format, SAC_PZ files should end with ".zp" and SAC file names can have any format (but should not end with ".zp"). Extraction from a NIED archive can be done using | ||
| Line 188: | Line 193: | ||
| ==== Plot routines ==== | ==== Plot routines ==== | ||
| - | All the plotting scripts are coded using python and the module pylab which have to be installed before using them. The module basemap is also needed for plotting maps but it is optional (even if we recommend to install it). | + | All the plotting scripts are coded using python and the module pylab which have to be installed before using them. The module cartopy is also needed for plotting maps but it is optional (even if we recommend to install it). |
| There are 3 script which can be used to plot the W phase inversion results: | There are 3 script which can be used to plot the W phase inversion results: | ||
| * The first script is <code>{WPHASE_HOME}/bin/make_grids.py</code> which be used to display grid-searches results. Use <code>{WPHASE_HOME}/bin/make_grids.py -h</code> to have more details on available options and arguments. The 2 other scripts plot observed and synthetic seismograms after W phase inversion. | * The first script is <code>{WPHASE_HOME}/bin/make_grids.py</code> which be used to display grid-searches results. Use <code>{WPHASE_HOME}/bin/make_grids.py -h</code> to have more details on available options and arguments. The 2 other scripts plot observed and synthetic seismograms after W phase inversion. | ||
| - | * In order to plot complete seismograms individually and place station on a map: <code>${WPHASE_HOME}/bin/traces.py</code> which draw complete seismograms and show station location on a map (if basemap is available). Please use <code>${WPHASE_HOME}/bin/traces.py -h</code> for more details on available options. | + | * In order to plot complete seismograms individually and place station on a map: <code>${WPHASE_HOME}/bin/traces.py</code> which draw complete seismograms and show station location on a map (if cartopy is available). Please use <code>${WPHASE_HOME}/bin/traces.py -h</code> for more details on available options. |
| * To plot concatenated waveforms:<code>{WPHASE_HOME}/bin/make_cwp.py</code> For more details on available options and arguments see <code>{WPHASE_HOME}/bin/make_cwp.py -h</code> to see the options. | * To plot concatenated waveforms:<code>{WPHASE_HOME}/bin/make_cwp.py</code> For more details on available options and arguments see <code>{WPHASE_HOME}/bin/make_cwp.py -h</code> to see the options. | ||
| Line 204: | Line 209: | ||
| <http://www.ngdc.noaa.gov/mgg/global/global.html> | <http://www.ngdc.noaa.gov/mgg/global/global.html> | ||
| - | If Basemap is installed, you can optionally draw ETOPO01 topography and | + | If Cartopy is installed, you can optionally draw ETOPO01 topography and |
| bathymetry. To do so, the path to the ETOPO01 netCDF file must be | bathymetry. To do so, the path to the ETOPO01 netCDF file must be | ||
| assigned to the environment variable $ETOPOFILE which can be | assigned to the environment variable $ETOPOFILE which can be | ||
| Line 392: | Line 397: | ||
| ${WPHASE_HOME}/bin/make_grids.py –p -b | ${WPHASE_HOME}/bin/make_grids.py –p -b | ||
| </code> | </code> | ||
| - | to create the file grid_search_xy.pdf containing Fig. 1 which corresponds to the above example. The option –b is optional: it allows to use the basemap module in order to plot coastlines and topography. | + | to create the file grid_search_xy.pdf containing Fig. 1 which corresponds to the above example. The option –b is optional: it allows to use the cartopy module in order to plot coastlines and topography. |
wphase/documentation.1634921358.txt.gz · Last modified: 2021/10/22 18:49 by wphase
Except where otherwise noted, content on this wiki is licensed under the following license: CC Attribution-Share Alike 4.0 International
