Esempi in Python per specfem2dASDF, esempi in Python per symdata.specfem2dASDF

Esempio n. 1

0

Mostra file

#!/usr/bin/env python
import symdata

# dbase=symdata.specfem2dASDF('synthetics_50.h5')
# dbase=symdata.specfem2dASDF('/home/lili/code/specfem2d/EXAMPLES/LFMembrane_SH_D/specfem2d_570.h5')
dbase = symdata.specfem2dASDF(
    '/lustre/janus_scratch/life9360/specfem2d_working_dir/LFMembrane_SH_D_990/specfem2d_990.h5'
)
# dbase=symdata.specfem2dASDF('/lustre/janus_scratch/life9360/LFMembrane_SH_0.1_20/specfem2d_mp.h5')
# dbase.readtxt('/home/lili/code/specfem2d/EXAMPLES/LFMembrane_SH_D/OUTPUT_FILES/STATIONS',
#         datadir='/home/lili/code/specfem2d/EXAMPLES/LFMembrane_SH_D/OUTPUT_FILES')

dbase.readtxt(
    '/lustre/janus_scratch/life9360/specfem2d_working_dir/LFMembrane_SH_D_990/DATA/STATIONS',
    datadir=
    '/lustre/janus_scratch/life9360/specfem2d_working_dir/LFMembrane_SH_D_990/OUTPUT_FILES'
)
dbase.AddEvent(x=500., z=500.)
inftan = symdata.InputFtanParam()
inftan.ffact = 5.
inftan.pmf = True
try:
    del dbase.auxiliary_data.DISPbasic1
    del dbase.auxiliary_data.DISPbasic2
    del dbase.auxiliary_data.DISPpmf1
    del dbase.auxiliary_data.DISPpmf2
except:
    pass
#
# dbase.aftanMP(outdir='/lustre/janus_scratch/life9360/specfem2d_working_dir/LFMembrane_SH_D_990/DISP', inftan=inftan, basic2=True,
#             pmf1=True, pmf2=True, tb=-11.9563813705)

Esempio n. 2

0

Mostra file

File: dataprocessing_aftan_travelT.py Progetto: NoisyLeon/specfem2dPy

#!/usr/bin/env python
import symdata



dbase=symdata.specfem2dASDF('/lustre/janus_scratch/life9360/specfem2d_working_dir/kernel_travelT_0.001/seismogram.h5')

dbase.readtxt('/lustre/janus_scratch/life9360/specfem2d_working_dir/kernel_travelT_0.001/DATA/STATIONS',
        datadir='/lustre/janus_scratch/life9360/specfem2d_working_dir/kernel_travelT_0.001/OUTPUT_FILES', factor=10)
dbase.AddEvent(x=200., z=1000.)
inftan=symdata.InputFtanParam()
inftan.ffact=5.
inftan.pmf=True
# inftan.vmin=2.5
# inftan.vmax=4.0
try:
    del dbase.auxiliary_data.DISPbasic1
    del dbase.auxiliary_data.DISPbasic2
    del dbase.auxiliary_data.DISPpmf1
    del dbase.auxiliary_data.DISPpmf2
except:
    pass
# 
# dbase.aftanMP(outdir='/lustre/janus_scratch/life9360/specfem2d_working_dir/LFMembrane_SH_D_780/DATA/DISP', inftan=inftan, basic2=True,
#             pmf1=True, pmf2=True, tb=-11.9563813705)
dbase.aftan(inftan=inftan, basic2=True,
            pmf1=True, pmf2=True, tb=-11.9563813705)
# dbase.aftan(tb=-11.9563813705, outdir='/lustre/janus_scratch/life9360/LFMembrane_SH_0.1_20/DISP', inftan=inftan, basic2=True,
#             pmf1=True, pmf2=True)

try:

Esempio n. 3

0

Mostra file

#!/usr/bin/env python
import symdata
import matplotlib.pyplot as plt

dbase = symdata.specfem2dASDF(
    '/lustre/janus_scratch/life9360/specfem2d_working_dir/LFMembrane_SH_healing_004_adj/seismogram_adj.h5'
)
#
dbase.readtxt(
    '/lustre/janus_scratch/life9360/specfem2d_working_dir/LFMembrane_SH_healing_004/DATA/STATIONS',
    datadir=
    '/lustre/janus_scratch/life9360/specfem2d_working_dir/LFMembrane_SH_healing_004/OUTPUT_FILES'
)
dbase.AddEvent(x=200., z=1000.)

tr = dbase.get_trace(staid='MEM2D.300S0')
specTR = symdata.specfem2dtrace(tr.data, tr.stats)
specTR.get_adjoint_stf(
    outdir=
    '/lustre/janus_scratch/life9360/specfem2d_working_dir/LFMembrane_SH_healing_004_adj',
    kerneltype='amp')
# specTR.get_adjoint_stf(outdir='/lustre/janus_scratch/life9360/specfem2d_working_dir/multipathing_2000km_4000km_adj_000', kerneltype='t', tmin=360, tmax=460)

Esempio n. 4

0

Mostra file

File: dataprocessing_aftan_990.py Progetto: NoisyLeon/specfem2dPy

#!/usr/bin/env python
import symdata

# dbase=symdata.specfem2dASDF('synthetics_50.h5')
# dbase=symdata.specfem2dASDF('/home/lili/code/specfem2d/EXAMPLES/LFMembrane_SH_D/specfem2d_570.h5')
dbase=symdata.specfem2dASDF('/lustre/janus_scratch/life9360/specfem2d_working_dir/LFMembrane_SH_D_990/specfem2d_990.h5')
# dbase=symdata.specfem2dASDF('/lustre/janus_scratch/life9360/LFMembrane_SH_0.1_20/specfem2d_mp.h5')
# dbase.readtxt('/home/lili/code/specfem2d/EXAMPLES/LFMembrane_SH_D/OUTPUT_FILES/STATIONS',
#         datadir='/home/lili/code/specfem2d/EXAMPLES/LFMembrane_SH_D/OUTPUT_FILES')

dbase.readtxt('/lustre/janus_scratch/life9360/specfem2d_working_dir/LFMembrane_SH_D_990/DATA/STATIONS',
        datadir='/lustre/janus_scratch/life9360/specfem2d_working_dir/LFMembrane_SH_D_990/OUTPUT_FILES')
dbase.AddEvent(x=500., z=500.)
inftan=symdata.InputFtanParam()
inftan.ffact=5.
inftan.pmf=True
try:
    del dbase.auxiliary_data.DISPbasic1
    del dbase.auxiliary_data.DISPbasic2
    del dbase.auxiliary_data.DISPpmf1
    del dbase.auxiliary_data.DISPpmf2
except:
    pass
# 
# dbase.aftanMP(outdir='/lustre/janus_scratch/life9360/specfem2d_working_dir/LFMembrane_SH_D_990/DISP', inftan=inftan, basic2=True,
#             pmf1=True, pmf2=True, tb=-11.9563813705)
dbase.aftan(inftan=inftan, basic2=True,
            pmf1=True, pmf2=True, tb=-11.9563813705)
# dbase.aftan(tb=-11.9563813705, outdir='/lustre/janus_scratch/life9360/LFMembrane_SH_0.1_20/DISP', inftan=inftan, basic2=True,
#             pmf1=True, pmf2=True)

Esempio n. 5

0

Mostra file

File: run_adj_stf.py Progetto: NoisyLeon/specfem2dPy

#!/usr/bin/env python
import symdata
import matplotlib.pyplot as plt

dbase=symdata.specfem2dASDF('/lustre/janus_scratch/life9360/specfem2d_working_dir/LFMembrane_SH_healing_004_adj/seismogram_adj.h5')
# 
dbase.readtxt('/lustre/janus_scratch/life9360/specfem2d_working_dir/LFMembrane_SH_healing_004/DATA/STATIONS',
        datadir='/lustre/janus_scratch/life9360/specfem2d_working_dir/LFMembrane_SH_healing_004/OUTPUT_FILES')
dbase.AddEvent(x=200., z=1000.)


tr=dbase.get_trace(staid='MEM2D.300S0')
specTR=symdata.specfem2dtrace(tr.data, tr.stats)
specTR.get_adjoint_stf(outdir='/lustre/janus_scratch/life9360/specfem2d_working_dir/LFMembrane_SH_healing_004_adj', kerneltype='amp')
# specTR.get_adjoint_stf(outdir='/lustre/janus_scratch/life9360/specfem2d_working_dir/multipathing_2000km_4000km_adj_000', kerneltype='t', tmin=360, tmax=460)

Esempio n. 6

0

Mostra file

File: dataprocessing_aftan.py Progetto: esiwgnahz/specfem2dPy

#!/usr/bin/env python
import symdata

dbase = symdata.specfem2dASDF(
    '/lustre/janus_scratch/life9360/specfem2d_working_dir/homo_lens_vert_10000km_2000km_R_400km_d_100km/seismogram.h5'
)
#
dbase.readtxt(
    '/lustre/janus_scratch/life9360/specfem2d_working_dir/homo_lens_vert_10000km_2000km_R_400km_d_100km/DATA/STATIONS',
    datadir=
    '/lustre/janus_scratch/life9360/specfem2d_working_dir/homo_lens_vert_10000km_2000km_R_400km_d_100km/OUTPUT_FILES',
    factor=10)
# dbase.AddEvent(x=500., z=2000.)
dbase.AddEvent(x=500., z=1000.)
inftan = symdata.InputFtanParam()
inftan.ffact = 1.
inftan.pmf = True
inftan.vmin = 2.1
# inftan.vmax=3.2
try:
    del dbase.auxiliary_data.DISPbasic1
    del dbase.auxiliary_data.DISPbasic2
    del dbase.auxiliary_data.DISPpmf1
    del dbase.auxiliary_data.DISPpmf2
except:
    pass

#
# dbase.aftanMP(outdir='/lustre/janus_scratch/life9360/specfem2d_working_dir/LFMembrane_SH_D_780/DATA/DISP', inftan=inftan, basic2=True,
#             pmf1=True, pmf2=True, tb=-11.9563813705)
dbase.aftan(inftan=inftan,

Esempio n. 7

0

Mostra file

File: dataprocessing_aftan.py Progetto: NoisyLeon/specfem2dPy

#!/usr/bin/env python
import symdata




dbase=symdata.specfem2dASDF('/lustre/janus_scratch/life9360/specfem2d_working_dir/homo_lens_vert_10000km_2000km_R_400km_d_100km/seismogram.h5')
# 
dbase.readtxt('/lustre/janus_scratch/life9360/specfem2d_working_dir/homo_lens_vert_10000km_2000km_R_400km_d_100km/DATA/STATIONS',
        datadir='/lustre/janus_scratch/life9360/specfem2d_working_dir/homo_lens_vert_10000km_2000km_R_400km_d_100km/OUTPUT_FILES', factor=10)
# dbase.AddEvent(x=500., z=2000.)
dbase.AddEvent(x=500., z=1000.)
inftan=symdata.InputFtanParam()
inftan.ffact=1.
inftan.pmf=True
inftan.vmin=2.1
# inftan.vmax=3.2
try:
    del dbase.auxiliary_data.DISPbasic1
    del dbase.auxiliary_data.DISPbasic2
    del dbase.auxiliary_data.DISPpmf1
    del dbase.auxiliary_data.DISPpmf2
except:
    pass

# 
# dbase.aftanMP(outdir='/lustre/janus_scratch/life9360/specfem2d_working_dir/LFMembrane_SH_D_780/DATA/DISP', inftan=inftan, basic2=True,
#             pmf1=True, pmf2=True, tb=-11.9563813705)
dbase.aftan(inftan=inftan, basic2=True,
            pmf1=True, pmf2=True, tb=-11.9563813705)
# dbase.aftan(tb=-11.9563813705, outdir='/lustre/janus_scratch/life9360/LFMembrane_SH_0.1_20/DISP', inftan=inftan, basic2=True,

Esempio n. 8

0

Mostra file

File: test_specasdf.py Progetto: esiwgnahz/specfem2dPy

import symdata

# dbase=symdata.sw4ASDF('/home/lili/sw4synthetics_ak135_EX_z_1km.h5')
dbase = symdata.specfem2dASDF('../specfem2d.h5')
# dbase.Readsac('station.lst', '../syndata_dir_000/sac_dir/', verbose=True)
dbase.PlotStreamsDistance()