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:49] wphase [Dependencies] |
wphase:documentation [2022/01/10 07:30] 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. 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 188: | 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 204: | 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 392: | 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. |