Exemple #1
0
import viscojapan as vj

f1 = vj.EpochalG('../G_Rake66.h5')
G1 = f1[0]

f2 = vj.EpochalG('../G_Rake90.h5')
G2 = f2[0]

f3 = vj.EpochalG('../G_Rake95.h5')
G3 = f3[0]

fno = 285

bm = vj.MyBasemap(region_code='near')
scale = 2e-2
mplt = vj.MapPlotDisplacement(basemap=bm)
mplt.plot_disp(G1[:, fno], f1.sites, scale=scale)
mplt.plot_disp(G2[:, fno], f2.sites, scale=scale, color='red')
mplt.plot_disp(G3[:, fno], f3.sites, scale=scale, color='blue')

vj.MapPlotFault('../fault_He50km.h5', basemap=bm).plot_fault(fno)

vj.plt.show()
Exemple #2
0
#nth = 34
nth = 15

point_source_file = './outs_He63km_Rake81/day_0000_flt_%04d.out'%nth

tp = np.loadtxt(point_source_file)
disp0 = tp[:,2:5].flatten()

disp1 = G1[:,nth]


#scale = 2e-3
scale = 4e-5

bm = vj.MyBasemap(region_code='I')

mplt = vj.MapPlotFault(fault_file,basemap=bm)
mplt.plot_fault(fno=nth)

mplt = vj.MapPlotDisplacement(basemap=bm)
mplt.plot_disp(disp0, sites, scale=scale,
               X=0.2, Y=0.8, U=1e-6, label='point source')
mplt.plot_disp(disp1, sites,color='r', scale=scale,
               X=0.2, Y=0.9, U=1e-6, label='reg source')



plt.savefig('~compare.pdf')
plt.show()
Exemple #3
0
from pylab import plt
import numpy as np

import viscojapan as vj
##
sites = np.loadtxt('../sites', '4a')
##ep = vj.EpochalDisplacement('../../cumu_post_with_seafloor.h5','../sites')
##slip_co = ep[0]
##
##
##mplt = vj.MapPlotDisplacement()
##
##mplt.plot_disp(slip_co,sites)
##plt.show()

ep = vj.EpochalSlip('../slip0/incr_slip0.h5')
slip_co = ep[0]

mplt = vj.MapPlotFault('../fault_model/fault_bott50km.h5')
mplt.plot_slip(slip_co)

plt.show()
from pylab import plt
import numpy as np

import viscojapan as vj

mplt = vj.MapPlotSlab()
mplt.plot_top()

tp = np.loadtxt(
    '/home/zy/workspace/viscojapan/Pollitz_slip_model/Pollitz-slipmodel.txt')
lons = tp[:, 0].reshape((42, -1))
lats = tp[:, 1].reshape((42, -1))

mplt.basemap.plot(lons, lats, color='green', latlon=True)
mplt.basemap.plot(np.ascontiguousarray(lons.T),
                  np.ascontiguousarray(lats.T),
                  color='green',
                  latlon=True)

mplt = vj.MapPlotFault('fault_bott50km.h5')
mplt.plot_fault(color='red')

plt.show()
Exemple #5
0
 def _plot_fault(self):
     if self.fault_file is not None:
         mplt = vj.MapPlotFault(self.fault_file, basemap=self.basemap)
         mplt.plot_slip(self.fault_slip)