Python model2d 예제들

예제 #1

파일: plot_model2d_EA_2017.py 프로젝트: iggcas-duyanan/Pygm

# m2d.plot(vmin=2.8, vmax=3.4)
# m2d.save('./10_phv_eigen_EA.npz')

# 
# m2d=model2d.model2d(minlon=minlon, maxlon=maxlon, dlon=0.5, minlat=minlat, maxlat=maxlat, dlat=0.5)
# m2d.read(infname='EA_data/Tph_10.0_lplc.lst')
# # m2d.Zarr=m2d.Zarr+0.025
# # m2d.mask=m2d1.mask
# m2d.np2ma()
# # m2d.plot(vmin=2.8, vmax=3.4)
# m2d.save('./10_T_lplc_EA.npz')
# 
# 
# #
# # 
m2d1=model2d.model2d(infname='./10_phv_eikonal_EA_0.0005.npz')
# m2d1=model2d.model2d(infname='./10_phv_helm_EA.npz')
m2d2=model2d.model2d(infname='./10_phv_eigen_EA.npz')
# m2d2=model2d.model2d(infname='./10_phv_helm_EA.npz')


# # 
# m2d1.np2ma()
# m2d1.minlat=minlat; m2d1.maxlat=maxlat; m2d1.minlon=minlon; m2d1.maxlon=maxlon
# m2d1.plot(label='Helmholtz phase velocity', vmin=2.8, vmax=3.3, geopolygons=basins)

# m2d1.np2ma()
# m2d1.minlat=minlat; m2d1.maxlat=maxlat; m2d1.minlon=minlon; m2d1.maxlon=maxlon
# m2d1.plot(label='Eikonal phase velocity', vmin=2.8, vmax=3.3, geopolygons=basins)
# # #
# m2d2.mask[(m2d2.lonArr>minlon)*(m2d2.lonArr<maxlon)*(m2d2.latArr>minlat)*(m2d2.latArr<maxlat) ]=0

예제 #2

파일: plot_model2d_for_honda.py 프로젝트: iggcas-duyanan/Pygm

import model2d

#
m2d = model2d.model2d(minlon=235, maxlon=295, dlon=0.5, minlat=25, maxlat=50)
# m2d.read(infname='/projects/life9360/US_phase_map/20_ANT.vel.HD_0.05_v1')
m2d.load('US_20sec.npz')
m2d.np2ma()
m2d.plot(vmin=3.2, vmax=3.7)
# m2d.save('US_20sec.npz')

예제 #3

파일: plot_model2d_US_2017.py 프로젝트: iggcas-duyanan/Pygm

minlat = 24. + 3
maxlat = 50. - 3.
minlon = -120.0 + 360. + 3.
maxlon = -80. + 360. - 3.

# m2d=model2d.model2d(infname='/projects/life9360/US_phase_map/10_ANT.vel.HD_0.05_v1')
# m2d.read(infname='/projects/life9360/US_phase_map/10_ANT.vel.HD_0.05_v1')
# m2d.np2ma()
# # m2d.plot(vmin=2.8, vmax=3.4, showfig=True)
# m2d.smooth(10)
# m2d.plot(vmin=2.8, vmax=3.4)
# m2d.save('./10_phv_smooth.npz')
#
m2d = model2d.model2d(minlon=minlon,
                      maxlon=maxlon,
                      dlon=0.5,
                      minlat=minlat,
                      maxlat=maxlat)
m2d.read(infname='US_data_0.25/Tph_%.1f_appV.lst' % per)
m2d.np2ma()
# m2d.plot(vmin=vmin, vmax=vmax)
m2d.save('US_data_0.25/%d_phv_eikonal_US.npz' % per)
# #
# m2d=model2d.model2d(infname='/projects/life9360/US_phase_map/10_ANT.vel.HD_0.05_v1')
# m2d.read(infname='US_data/helm_%d_phv.lst' %per)
# m2d.np2ma()
# m2d.plot(vmin=vmin, vmax=vmax)
# m2d.save('US_data/%d_phv_helm_US.npz' %per)
# #
m2d = model2d.model2d(minlon=minlon,
                      maxlon=maxlon,

예제 #4

파일: plot_model2d_US_2017_phV_grid_test.py 프로젝트: iggcas-duyanan/Pygm

vmin = 3.7
vmax = 4.1
vArr = [[2.8, 3.4], [3.3, 3.7], [3.5, 4.0], [3.7, 4.1]]

outdir = '/home/life9360/figures_grid_test'
# #
minlat = 24.
maxlat = 50.
minlon = -120.0 + 360.
maxlon = -80. + 360.
for dx in dxArr:
    print 'dx =', dx
    for per, v in zip(perArr, vArr):
        vmin = v[0]
        vmax = v[1]
        m2d = model2d.model2d(infname='US_data_%g/%d_phv_eikonal_US.npz' %
                              (dx, per))
        # # # # #
        plt.figure()
        m2d.np2ma()
        m2d.minlat = minlat
        m2d.maxlat = maxlat
        m2d.minlon = minlon
        m2d.maxlon = maxlon
        m2d.plot(label='Eikonal phase velocity',
                 vmin=vmin,
                 vmax=vmax,
                 dv=0.1,
                 showfig=False)
        phVname = outdir + '/eikonalV_%dsec_%g.pdf' % (per, dx)
        plt.savefig(phVname)

예제 #5

파일: plot_model2d_US_2017_diff_grid_test.py 프로젝트: iggcas-duyanan/Pygm

import model2d
import matplotlib.pyplot as plt
import numpy as np
dxArr = np.array([0.25, 0.75, 1.0, 1.5, 2.0, 2.5, 3.0])
perArr = np.array([10., 20., 30., 40.])

outdir = '/home/life9360/figures_grid_test'
# #
minlat = 24.
maxlat = 50.
minlon = -120.0 + 360.
maxlon = -80. + 360.
for dx in dxArr:
    print 'dx =', dx
    for per in perArr:
        m2d1 = model2d.model2d(infname='US_data_%g/%d_phv_eikonal_US.npz' %
                               (dx, per))
        m2d2 = model2d.model2d(infname='US_data_0.25/%d_phv_eigen_US.npz' %
                               (per))
        # # # # #
        diffm2d = m2d1 - m2d2
        minlat = 24.
        maxlat = 50.
        minlon = -120.0
        maxlon = -80.
        diffm2d.minlat = minlat
        diffm2d.maxlat = maxlat
        diffm2d.minlon = minlon
        diffm2d.maxlon = maxlon
        # # #
        diffm2d.np2ma()
        import numpy as np

예제 #6

파일: run_txt2npy_grid.py 프로젝트: iggcas-duyanan/Pygm

import model2d
import numpy as np
dxArr = np.array([0.5, 0.75, 1.0, 1.5, 2.0, 2.5, 3.0])
perArr = np.array([10., 20., 30., 40.])
for dx in dxArr:
    print 'dx =', dx
    for per in perArr:
        infname = 'US_data_%g/Tph_%.1f_appV.lst' % (dx, per)
        m2d = model2d.model2d(
            infname='/projects/life9360/US_phase_map/10_ANT.vel.HD_0.05_v1')
        m2d.read(infname='US_data_%g/Tph_%.1f_appV.lst' % (dx, per))
        m2d.np2ma()
        m2d.save('US_data_%g/%d_phv_eikonal_US.npz' % (dx, per))

예제 #7

# minlon=86.
# maxlon=132.

minlat = 24.
maxlat = 50.
minlon = -120.0 + 360.
maxlon = -80. + 360.

minlat = 25
maxlat = 50.
minlon = 235
maxlon = 295
# #
m2d = model2d.model2d(minlon=minlon,
                      maxlon=maxlon,
                      dlon=0.05,
                      minlat=minlat,
                      maxlat=maxlat,
                      dlat=0.05)
# m2d.read(infname='/projects/life9360/code/SES3DPy/diffA.lst')
# m2d.read(infname='US_data/phV_100.lst')
#m2d.read(infname='/projects/life9360/code/SES3DPy/MAPS_US/smpkolya_phv_R_40')
# m2d.read(infname='/projects/life9360/US_phase_map/8_ANT.vel.HD_0.05_v1')
m2d.read_numpy_txt(
    infname='/projects/life9360/US_phase_map/8_ANT.vel.HD_0.05_v1')
m2d.np2ma()
m2d.plot(vmin=2.8, vmax=3.4)
# m2d.smooth(10)
#m2d.plot(vmin=3.6, vmax=4.1)
# m2d.plot(vmin=-20, vmax=20, cmap='bwr')
# m2d.save('./10_phv_beam.npz')
# #