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)