import field2d_earth
import GeoPolygon
import raypath
import time
#
# basins=GeoPolygon.GeoPolygonLst()
# basins.ReadGeoPolygonLst('basin1')
#
minlat = 20.
maxlat = 52.
minlon = 80.
maxlon = 134.
#
field = field2d_earth.Field2d(minlon=minlon,
maxlon=maxlon,
dlon=0.1,
minlat=minlat,
maxlat=maxlat,
dlat=0.1,
period=10.)
# field.read_dbase(datadir='./fmst_dbase_0.2')
# dset=raypath.rayASDF('../rays_0.2.h5')
# dset.read_raydat('/projects/life9360/code/Pyfmst/fmm_working_0.2/gmtplot/rays.dat')
# dset.get_pre_amp(field2d=field, outfname='pre_amp_0.2.lst')
field.read_dbase(datadir='./fmst_dbase_0.2')
dset = raypath.rayASDF('../rays_0.2.h5')
dset.read_raydat('/projects/life9360/code/Pyfmst/fmm_working/gmtplot/rays.dat')
# dset.get_pre_amp(field2d=field, outfname='pre_amp_0.2_0.2.lst')
import field2d_earth
minlat = 24.
maxlat = 50.
minlon = -120.0
maxlon = -80.
per = 10.
field = field2d_earth.Field2d(minlon=minlon,
maxlon=maxlon,
dlon=0.5,
minlat=minlat,
maxlat=maxlat,
dlat=0.5,
period=per)
field.read(
fname=
'/lustre/janus_scratch/life9360/ses3d_field_working/stf_100_10sec_US/Tph_%.1f.txt'
% per)
# field.add_noise(sigma=5.)
field.interp_surface(
workingdir=
'/lustre/janus_scratch/life9360/ses3d_field_working/field_working_US',
outfname='Tph_%gsec' % per)
field.check_curvature(
workingdir=
'/lustre/janus_scratch/life9360/ses3d_field_working/field_working_US')
field.gradient_qc(
workingdir=
'/lustre/janus_scratch/life9360/ses3d_field_working/field_working_US',
import field2d_earth
import numpy as np
import matplotlib.pyplot as plt
minlon = 235.
maxlon = 260.
minlat = 31.
maxlat = 50.
field = field2d_earth.Field2d(minlon=minlon,
maxlon=maxlon,
dlon=0.2,
minlat=minlat,
maxlat=maxlat,
dlat=0.2,
period=10.)
field.read(
fname=
'/scratch/summit/life9360/eikonal_working_single/60.0sec/E00979.amp.lst')
field.ZarrIn = field.ZarrIn * 1e9
field.fieldtype = 'amp'
tension = .0
field.interp_surface(workingdir='./eik_working_debug',
outfname='Tph_10sec',
tension=tension)
field.coarse_data()
field.Laplacian('diff2')
field1 = field2d_earth.Field2d(minlon=minlon,
maxlon=maxlon,
dlon=0.1,
basins.ReadGeoPolygonLst('basin1')
minlat = 20.
maxlat = 53.
minlon = 80.
# maxlon=135.
maxlon = 135.
# vmodel=pyfmst.vmodel(minlat=minlat, maxlat=maxlat, minlon=minlon, maxlon=maxlon, dlat=0.5, dlon=0.5, v0=3.0)
# vmodel.read_cv('./China_2015_disp_v1.0/10.phase.map')
# vmodel.smooth(100)
# vmodel.read_cv('./China_2015_disp_v1.0/10.phase.map')
# vmodel.write4field2d('10.phase_extended_map')
#
# minlat=20.2
# maxlat=51.8
# minlon=80.2
# maxlon=133.8
field = field2d_earth.Field2d(minlon=minlon,
maxlon=maxlon,
dlon=0.5,
minlat=minlat,
maxlat=maxlat,
dlat=0.5,
period=10.,
fieldtype='Tph')
field.read(fname='./China_2015_disp_v1.0/10.phase.map')
workingdir = './field_working'
field.interp_surface(workingdir=workingdir, outfname='Tph_10sec')
field.write('10.phase_extended_map', 'txt')
# field.plot_field(contour=False, geopolygons=basins)