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:40] wphase [RUNA scripts] |
wphase:documentation [2022/01/10 07:39] (current) wphase [Extracting data from a gzipped tar archive] |
||
|---|---|---|---|
| Line 27: | Line 27: | ||
| - gcc and gfortran | - gcc and gfortran | ||
| - python2.7 (or later) | - python2.7 (or later) | ||
| - | - You have to install numpy, matplotlib, basemap and netCDF4 to run some python scripts which make figures. | + | - You have to install numpy, matplotlib, cartopy and netCDF4 to run some python scripts which make figures. |
| - rdseed | - rdseed | ||
| Line 90: | Line 90: | ||
| ==== Extracting data from a gzipped tar archive ==== | ==== Extracting data from a gzipped tar archive ==== | ||
| - | As SEED volumes are no longer supported, waveforms and instrument response can also be retrived from a gzipped tar archive (hereafter referred to as "tgz") including SAC and SAC_PZ (poles and zeros) files. Instead of using the "extract.csh" script to extract data from SEED, we can then extract data from a tgz archive using: | + | As SEED volumes are no longer supported, waveforms and instrument response can also be extracted from a gzipped tar archive (hereafter referred to as "tgz") including both SAC and SAC_PZ (poles and zeros) files. |
| + | Instead of using the "extract.csh" script to extract data from SEED, we can then extract data from a tgz archive using: | ||
| <code>${WPHASE_HOME}/bin/tgz2DATA_org.py path_to_tgz_archive.tgz</code> | <code>${WPHASE_HOME}/bin/tgz2DATA_org.py path_to_tgz_archive.tgz</code> | ||
| - | where "path_to_tgz_archive.tgz" is the path to the tgz archive including SAC and SAC_PZ files. | + | where "path_to_tgz_archive.tgz" is the path to the tgz archive including SAC and SAC_PZ files. This script will extract SAC and SAC_PZ files a the directory named "DATA_org". |
| 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 150: | Line 156: | ||
| ----- | ----- | ||
| ==== RUNA scripts ==== | ==== RUNA scripts ==== | ||
| - | Data screening and inversion can be done using the RUNA3_lite.csh script: | + | Data screening and inversion can all be done at once using the RUNA3_lite.csh script: |
| <code> | <code> | ||
| ${WPHASE_HOME}/bin/RUNA3_lite.csh | ${WPHASE_HOME}/bin/RUNA3_lite.csh | ||
| Line 187: | 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 203: | 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 391: | 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.1634920841.txt.gz · Last modified: 2021/10/22 18:40 by wphase
Except where otherwise noted, content on this wiki is licensed under the following license: CC Attribution-Share Alike 4.0 International
