This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision Last revision Both sides next revision | ||
wphase:documentation [2021/10/22 18:40] wphase [RUNA scripts] |
wphase:documentation [2022/01/10 07:30] 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. 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</code> | ||
+ | Although this is the default format, extraction can also be stated explicitly by running | ||
<code>${WPHASE_HOME}/bin/tgz2DATA_org.py path_to_tgz_archive.tgz wilber</code> | <code>${WPHASE_HOME}/bin/tgz2DATA_org.py path_to_tgz_archive.tgz wilber</code> | ||
Line 150: | Line 153: | ||
----- | ----- | ||
==== 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 190: | ||
==== 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 206: | ||
<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 394: | ||
${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. |