forked from NoisyLeon/NoisePy
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathdebug_xcorr_azi.py
44 lines (32 loc) · 2.57 KB
/
debug_xcorr_azi.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
import eikonaltomo
import raytomo
import numpy as np
#-----------------------
#initialization
#-----------------------
# dset = eikonaltomo.EikonalTomoDataSet('/work1/leon/ALASKA_work/hdf5_files/debug2_xcorr_60.h5')
# dset.set_input_parameters(minlon=188, maxlon=238, minlat=52, maxlat=72, pers=np.array([28., 30., 32.]))
# dset.eikonal_stack_mp(anisotropic=True, spacing_ani=2., N_bin=20, workingdir='./debug2_xcorr_2deg_60', enhanced=True)
# dset.compute_azi_aniso()
# dset = eikonaltomo.EikonalTomoDataSet('/work1/leon/ALASKA_work/hdf5_files/debug2_xcorr_90.h5')
# dset = eikonaltomo.EikonalTomoDataSet('/work1/leon/ALASKA_work/hdf5_files/comp_xcorr_2deg.h5')
# dset.set_input_parameters(minlon=188, maxlon=238, minlat=52, maxlat=72, pers=np.array([28., 30., 32.]))
# dset.eikonal_stack_mp(anisotropic=True, spacing_ani=2., N_bin=20, azi_amp_tresh=0.1, workingdir='./comp_xcorr_2deg', enhanced=False)
# dset.compute_azi_aniso()
# dset = eikonaltomo.EikonalTomoDataSet('/work1/leon/ALASKA_work/hdf5_files/comp_xcorr_2deg_0.05.h5')
# dset.set_input_parameters(minlon=188, maxlon=238, minlat=52, maxlat=72, pers=np.array([30.]))
# dset.set_input_parameters(minlon=188, maxlon=238, minlat=52, maxlat=72, pers=np.array([28., 30., 32.]))
# dset.eikonal_stack_mp(anisotropic=True, spacing_ani=2., N_bin=20, azi_amp_tresh=0.05, workingdir='./comp_xcorr_2deg_0.05', enhanced=False)
# dset.compute_azi_aniso()
# dset = eikonaltomo.EikonalTomoDataSet('/work1/leon/ALASKA_work/hdf5_files/comp_xcorr_2deg_0.05_untest.h5')
# dset.compute_azi_aniso_enhanced()
dset = eikonaltomo.EikonalTomoDataSet('/work1/leon/ALASKA_work/hdf5_files/azi_xcorr_all_2deg_0.05.h5')
# dset.set_input_parameters(minlon=188, maxlon=238, minlat=52, maxlat=72, pers=np.array([20., 22., 24., 26., 28., 30., 32., 34., 36., 38., 40.]))
# dset.eikonal_stack_mp(anisotropic=True, spacing_ani=2., N_bin=20, azi_amp_tresh=0.05, workingdir='./azi_xcorr_all_2deg_0.05')
# dset.compute_azi_aniso()
# dset.plot_fast_axis(runid=0, period=30., factor=5, ampref = 1.)
# dset.diff_fast_axis(inh5fname='/work1/leon/ALASKA_work/hdf5_files/azi_quake_2deg_0.05.h5', vmin=0., vmax=90., period=30., runid=0)
# dset.combine_amp(inh5fname='/work1/leon/ALASKA_work/hdf5_files/azi_quake_2deg_0.05.h5', period=30., runid=0)
# dset.combine_fast_axis(inh5fname='/work1/leon/ALASKA_work/hdf5_files/azi_quake_2deg_0.05.h5', period=30., runid=0)
# dset.plot_fast_axis(0, 30., datatype='vel_iso', vmin=3.6, vmax=3.9, ampref=1., normv=1.5, scaled=True, factor=5)
# dset.diff_amp(inh5fname='/work1/leon/ALASKA_work/hdf5_files/azi_quake_2deg_0.05.h5', period=30., runid=0)