m2.add_perturb_layer(0, 20., 2, 6.06, False)
m2.add_perturb_layer(0, 20., 3, 6.24, False)
m2.add_perturb_layer(0, 20., 4, 0.72, False)
m2.add_perturb_layer(0, 20., 5, 2.73, False)
m2.init_tilt()
m2.dipArr[-1] = 27; m2.dipArr[-2] = 27
m2.strikeArr[-1] = 110; m2.strikeArr[-2] = 110
m2.rot_dip_strike()
m2.decompose()
rsolver1 = ref.ref_solver(m1)
rsolver1.init_default_2()
for baz in np.arange(12)*30.:
rsolver1.solve_anirec(baz=baz)
rsolver1.plot_baz_rf(comp='T', showfig=False)
rsolver2 = ref.ref_solver(m2)
rsolver2.init_default_2()
for baz in np.arange(12)*30.:
rsolver2.solve_anirec(baz=baz)
rsolver2.plot_baz_rf(comp='T')
import numpy as np
import matplotlib.pyplot as plt
m = vmodel.model1d()
m.model_ak135_cps()
m.flat = 1
#
# m.add_perturb_layer(0, 35., 0, 3.2, False)
# m.add_perturb_layer(0, 35., 1, 3.2, False)
# m.add_perturb_layer(0, 35., 2, 5.8, False)
# m.add_perturb_layer(0, 35., 3, 5.8, False)
# m.add_perturb_layer(0, 35., 4, 1.0, False)
# m.add_perturb_layer(0, 35., 5, 2.73, False)
rsolver = ref.ref_solver(m)
rsolver.init_default_2()
rsolver.solve_theo()
rsolver.solve_raysum()
# rsolver.solve_anirec()
rsolver.deconvolve_raysum()
ax = plt.subplot()
ax.plot(rsolver.time, rsolver.rf, 'k-')
tfill = rsolver.time[rsolver.rf > 0]
yfill = rsolver.rf[rsolver.rf > 0]
ax.fill_between(tfill,
0.,
yfill,
color='red',
linestyle='--',
# m.add_perturb_layer(0, 30., 1, 3.2, False)
# m.add_perturb_layer(0, 30., 2, 5.8, False)
# m.add_perturb_layer(0, 30., 3, 5.8, False)
# m.add_perturb_layer(0, 30., 4, 1.0, False)
# m.add_perturb_layer(0, 30., 5, 2.73, False)
# m.add_perturb_layer_love(0, 30., 0, 0.1, True)
# m.add_perturb_layer_love(0, 30., 1, -0.1, True)
m.add_perturb_layer_love(0, 20., 3, -0.05, True)
m.add_perturb_layer_love(0, 20., 4, 0.05, True)
#
m.init_tilt()
m.dipArr[-1] = 90
m.dipArr[-2] = 90 #; m.dipArr[-3] = 90; m.dipArr[-4] = 90
# m.dipArr[-1] = 10; m.dipArr[-2] = 10; m.dipArr[-3] = 10; m.dipArr[-4] = 10
m.strikeArr[-1] = 0
m.strikeArr[-2] = 0
#
#
rsolver1 = ref.ref_solver(m)
rsolver1.init_default_2()
rsolver1.solve_raysum(bazin=[0., 90])
ax = plt.subplot()
for i in xrange(rsolver1.bazArr.size):
ax.plot(rsolver1.time, rsolver1.trROT[0, :, i], '--', lw=3)
plt.xlim([0., 20.])
plt.show()
# m.add_perturb_layer_love(0, 30., 1, -0.1, True) m.add_perturb_layer_love(0, 30., 3, -0.05, True) m.add_perturb_layer_love(0, 30., 4, 0.05, True) # m.init_tilt() m.dipArr[-1] = 90 m.dipArr[-2] = 90 m.dipArr[-3] = 90 m.dipArr[-4] = 90 # m.dipArr[-1] = 10; m.dipArr[-2] = 10; m.dipArr[-3] = 10; m.dipArr[-4] = 10 m.strikeArr[-1] = 0 m.strikeArr[-2] = 0 # # rsolver1 = ref.ref_solver(m) rsolver1.init_default_3() rsolver1.solve_anirec(baz=0.) rsolver1.solve_anirec(baz=90.) rsolver2 = ref.ref_solver(m) rsolver2.init_default_3() rsolver2.solve_raysum(bazin=[0., 90.]) # # m = vmodel.model1d() m.model_ak135_cps() m.flat = 1 #
import numpy as np
import matplotlib.pyplot as plt
m = vmodel.model1d()
m.model_ak135_cps()
m.flat = 1
#
# m.add_perturb_layer(0, 35., 0, 3.2, False)
# m.add_perturb_layer(0, 35., 1, 3.2, False)
# m.add_perturb_layer(0, 35., 2, 5.8, False)
# m.add_perturb_layer(0, 35., 3, 5.8, False)
# m.add_perturb_layer(0, 35., 4, 1.0, False)
# m.add_perturb_layer(0, 35., 5, 2.73, False)
rsolver = ref.ref_solver(m)
rsolver.init_default_2()
rsolver.solve_theo()
rsolver.solve_anirec()
ax = plt.subplot()
ax.plot(rsolver.time, rsolver.rf, 'g-', lw=4, label='theo, T. Shibutani')
time = rsolver.time
# ax.fill_between(time, y2=0., y1=rsolver.rf, where=rsolver.rf>0, color='green', lw=0.01, alpha=0.3, interpolate=True)
ax.plot(rsolver.time,
rsolver.rfr,
'r--',
lw=3,
label='anirec, V.Levin & J. Park')
# ax.fill_between(time, y2=0., y1=rsolver.rfr, where=rsolver.rfr>0, color='red', lw=0.01, alpha=0.3, interpolate=True)