computational-seismology / pytomo3d Goto Github PK

It might be good to have the option of padding zeroes here. One example were this would be beneficial would be when filtering specfem synthetics (that start from time=0) for near-field stations.

The coordinates of event and receiver are in opposite order, so the back azimuth is actually azimuth.
in pytomo3d/adjoint/adjsrc.py line 203 change
_, baz, _ = gps2DistAzimuth(elat, elon, slat, slon)
to
_, baz, _ = gps2DistAzimuth(slat, slon, elat, elon)
should fix the issue

Hi,

Using the default value for water_level (60) in instrument response removal can cause problems as discussed in obspy#1594.

Corresponding lines are in processed.py L246-L250.

Adding water_level=100 parameter to the remove_response calls solves the issue for my tests. However, water_level value may need to be decided more carefully.

For example, if stationxml provided, then sanity check, like check Z, N, E, azimuth and dip should be done.
If NE is not exact NE, then rotate should be done on the base of 12.

I think the new adjoint source has a default location id "" (null)
https://github.com/wjlei1990/pytomo3d/blob/master/pytomo3d/adjoint/sum_adjoint.py#L129

the check will not let the code proceed because the adj_base.id has default loc_id while adj.id doesn't
https://github.com/wjlei1990/pytomo3d/blob/master/pytomo3d/adjoint/sum_adjoint.py#L29-L31

Hi,

I think there is a typo in source_weights.py#L109. param[flag] should be param[plot]. This causes to call the plot function even if the plot is set to False.

Here, mentioned by Ridvan.
https://github.com/wjlei1990/pytomo3d/blob/master/pytomo3d/window/window.py#L159-L163

Issue 1:
Error during the processing of station 'TA.K56A' on rank 8:
Traceback (At max 3 levels - most recent call last):
File "/autofs/nccs-svm1_home1/youyir/src/pypaw/src/pypaw/procbase.py", line 200, in smart_run
self._core(path, param)
File "/autofs/nccs-svm1_home1/youyir/src/pypaw/src/pypaw/adjoint.py", line 236, in _core
output_filename)
File "/autofs/nccs-svm1_home1/youyir/src/pyasdf/pyasdf/asdf_data_set.py", line 1423, in process_two_files
getattr(other_ds.waveforms, station))
File "/autofs/nccs-svm1_home1/youyir/src/pypaw/src/pypaw/adjoint.py", line 109, in adjoint_wrapper
**postproc_param)
File "/autofs/nccs-svm1_home1/youyir/src/pytomo3d/pytomo3d/adjoint/process_adjsrc.py", line 361, in process_adjoint
add_missing_components(adj_stream, component_list=["Z", "R", "T"])
File "/autofs/nccs-svm1_home1/youyir/src/pytomo3d/pytomo3d/adjoint/process_adjsrc.py", line 249, in add_missing_components
missinglist.remove(tr.stats.channel[-1])

ValueError: list.remove(x): x not in list

Issue 2:
Error extracting station('TA.G59A') info:list index out of range
Error extracting station('AK.DHY') info:list index out of range

When using --user flag for pip installation, the packages (and the dependencies) are installed inside the ~/.local directory and not the conda environment.

(ASDF6) chinook00 % python
Python 2.7.11 |Continuum Analytics, Inc.| (default, Dec  6 2015, 18:08:32) 
[GCC 4.4.7 20120313 (Red Hat 4.4.7-1)] on linux2
Type "help", "copyright", "credits" or "license" for more information.
Anaconda is brought to you by Continuum Analytics.
Please check out: http://continuum.io/thanks and https://anaconda.org
>>> import obspy
>>> obspy
<module 'obspy' from '/u1/uaf/vsilwal/conda_environments/ASDF6/lib/python2.7/site-packages/obspy/__init__.pyc'>
>>> obspy.__version__
u'1.0.3'
(ASDF6) chinook00 % cd pyflex
(ASDF6) chinook00 % pip install -v -e . --user
(ASDF6) chinook00 % python
Python 2.7.11 |Continuum Analytics, Inc.| (default, Dec  6 2015, 18:08:32) 
[GCC 4.4.7 20120313 (Red Hat 4.4.7-1)] on linux2
Type "help", "copyright", "credits" or "license" for more information.
Anaconda is brought to you by Continuum Analytics.
Please check out: http://continuum.io/thanks and https://anaconda.org
ImportError: No module named __version__
>>> import obspy
>>> obspy
<module 'obspy' from '/u1/uaf/vsilwal/.local/lib/python2.7/site-packages/obspy/__init__.pyc'>
>>> obspy.__version__
u'1.1.0'

Should I not be using --user flag?
I am using the system module conda and pip comes along with it.

(ASDF6) chinook00 % module list
Currently Loaded Modulefiles:
  1) lang/Anaconda2/4.0.0                       5) compiler/ifort/2016.3.210-GCC-5.4.0-2.26   9) tools/XZ/5.2.2-pic-intel-2016b
  2) compiler/GCCcore/5.4.0                     6) toolchain/pic-iompi/2016b                 10) lib/zlib/1.2.8-pic-intel-2016b
  3) tools/binutils/2.26-GCCcore-5.4.0          7) numlib/imkl/11.3.3.210-pic-iompi-2016b    11) tools/Szip/2.1-pic-intel-2016b
  4) compiler/icc/2016.3.210-GCC-5.4.0-2.26     8) toolchain/pic-intel/2016b                 12) data/HDF5/1.8.15-pic-intel-2016b

It seems that:

https://github.com/wjlei1990/pytomo3d/blob/master/pytomo3d/adjoint/adjsrc.py#L137

is not consistent with

https://github.com/wjlei1990/pytomo3d/blob/master/pytomo3d/adjoint/adjsrc.py#L82

windows are different objects in these two functions.