Comments (7)
@basaks @alexgorb
Fixed and just sent you guys the code.
Note: you seed to do 'make all' inside 'tau' first before 'ellip'
from hiperseis.
When I tried to compile the travel time code, I get the following error:
tau:$ pwd
/home/sudipta/Downloads/rses.anu.edu.au/pub/ak135/tau
tau:$ tree
.
├── ak135.tvel
├── emtvelin.f
├── iasp91.tvel
├── libsun.f
├── libtau.f
├── limits.inc
├── Makefile
├── remodl.f
├── setbrn.f
├── setbrn.o
├── travel.help
├── ttimes.f
├── ttimes.o
└── ttlim.inc
0 directories, 14 files
tau:$ make ttimes
f77 -c -o libtau.o libtau.f
libtau.f:1405:9:
if(iargc(i).lt.n) go to 1
1
Error: Too many arguments in call to ‘iargc’ at (1)
libtau.f:602:72: Error: Too many arguments in call to ‘abort’ at (1)
libtau.f:255:72: Error: Too many arguments in call to ‘abort’ at (1)
<builtin>: recipe for target 'libtau.o' failed
make: *** [libtau.o] Error 1
@alexgorb @rh-downunder Can help compiling this code?
from hiperseis.
@alexgorb Can you confirm that we just need to implement ellipticity correction for P and S waves only?
from hiperseis.
from hiperseis.
Thanks for that @alexgorb. A typical output is the following for a sc3 captured event origin/arrival combination:
Source latitude: -4.1000332
Source depth (km): 10
Azimuth from source: 202.9688263
Enter delta: 10.72813511
delta # code time time(el) dT/dD dT/dh d2T/dD2
PbPb is not available
PgPg is not available
SbSb is not available
SgSg is not available
10.73 1 Pn 153.66 153.82 13.6907 -1.20E-01 -2.43E-02
2 pPn 156.07 156.25 13.6930 1.20E-01 -2.47E-02
3 sPn 157.48 157.64 13.6924 2.61E-01 -2.46E-02
4 PnPn 161.35 161.87 13.7406 -1.20E-01 -1.04E-01
5 PbPb 188.38 188.90 17.0386 -7.86E-02 -9.86E-02
6 PgPg 205.54 206.06 19.1396 -2.40E-03 -1.11E-01
7 Sn 273.80 274.07 24.5187 -1.86E-01 -9.18E-03
8 sSn 277.53 277.93 24.5257 1.86E-01 -9.37E-03
9 SnSn 285.63 286.57 24.6497 -1.85E-01 -4.02E-02
10 SbSb 317.72 318.66 28.7667 -1.28E-01 -5.87E-02
11 SgSg 344.55 345.49 32.0837 -4.03E-03 -6.64E-02
12 PcP 515.45 516.12 1.0155 -1.72E-01 3.32E-03
13 ScP 727.95 728.92 1.3115 -2.89E-01 2.59E-03
14 PcS 729.12 730.04 1.3118 -1.72E-01 2.59E-03
15 ScS 943.02 944.23 1.8733 -2.89E-01 1.80E-03
16 PKiKP 994.39 995.44 0.2395 -1.72E-01 1.37E-02
17 pPKiKP 997.83 998.89 0.2395 1.72E-01 1.37E-02
18 sPKiKP 999.00 1000.06 0.2395 2.89E-01 1.37E-02
19 SKiKP 1205.34 1206.70 0.2529 -2.89E-01 1.30E-02
20 PKKPdf 1910.36 1914.25 -0.2380 -1.72E-01 -1.38E-02
21 SKKPdf 2121.31 2124.98 -0.2261 -2.89E-01 -1.45E-02
22 PKKSdf 2122.48 2127.92 -0.2261 -1.72E-01 -1.45E-02
23 SKKSdf 2333.41 2338.53 -0.2153 -2.89E-01 -1.52E-02
24 P'P'df 2421.67 2421.24 -0.3095 -1.72E-01 -1.06E-02
25 P'P'ab 2608.82 2608.39 -4.4433 -1.68E-01 4.43E-01
26 S'S'df 3269.05 3266.74 -0.2428 -2.89E-01 -1.35E-02
from hiperseis.
We have an adaptation for ellipticity correction here: https://github.com/GeoscienceAustralia/ellip-corr
Results match that of the previous post. Examples:
In [1]: import ellipcorr
In [2]: ellipcorr.ellipticity_corr('sP', 65, 124, 45, 39)
Out[2]: -0.38976147770881653
In [3]: ellipcorr.ellipticity_corr('P', 65, 124, 45, 39)
Out[3]: -0.3774765431880951
In [4]: ellipcorr.ellipticity_corr('pPKiKP', 65, 124, 45, 39)
Out[4]: -0.7800000309944153
from hiperseis.
We now have tests for ellipticity correction
showing that the ellipticity correction is computed correctly for each of the phases in the list above.
There is also basic python packaging implemented in that repo.
from hiperseis.
Related Issues (20)
- uncover learn had a new problem HOT 10
- can not find the code about rf 1D Earth structure inversion
- Build of ellipcorr fails HOT 2
- receiver function handpick
- Some bug during calculating receiver function using ITER method HOT 7
- try to compile the rf-inversion code from i-earth HOT 2
- run iLoc on the events generated for temporary surveys HOT 5
- travel time tomographic inversion steps HOT 1
- improve the picking code by implementing parallel execution HOT 3
- improve relevant passive-seismic code quality by adding documentation in docstrings
- convert temporary stations metadata to a format compatible with cluster code HOT 1
- run the picking code on antelope events and check against analyst picks, for both p-picks and s-picks HOT 1
- run s-phase receiver function analysis on mtIsa HOT 1
- incorporate rotation from NE to RT before picking for s-phase arrivals HOT 5
- inversion of seismogram to below the sediment layer in basins to improve receiver function results HOT 1
- Improvements/enhancements for the picking of P and S phases HOT 4
- picking of Australian and other relevant permanent stations cwb waveform data for ISC events HOT 4
- receiver function analysis for mtIsa data with new set of "parameters" HOT 2
- Calculating shifts in the GPS times of temporary station picks using nearby permanent station picks HOT 4
- fix bugs in picking code HOT 1
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from hiperseis.