80 likes | 91 Views
This progress report discusses the ongoing work on 3-D inversions and multiple ridge fitting for seismic analysis. The code is running under IDL and has the ability to use different size tiles. Kernels are derived from a model with a polytropic atmosphere and are scalable. Future plans include taking errors into account, determining kernels from a better model, and porting the code to Fortran.
E N D
3-D inversions and multiple-ridge fitting - PROGRESS REPORT Deborah Haber, Brad Hindman, Mike Thompson
3-D inversions • Code is running under IDL and has ability to use different size tiles • Kernels are from model with polytropic atmosphere, and are scalable
early 3D inversion 30 From 16o tiles n= 1, 2 and 2o tiles n= 1 -30 -30 30
Future Work – Short Term • Add ability to take errors into account • Determine kernels from better model • Port code to Fortran
Multiple Ridge Fits- Motivation • Tiles of smaller sizes or at lower frequencies generally have a large amount of overlap in the peaks leading to poor determinations of the amplitude, widths, and frequencies of the waves. • By fitting all the n-value peaks for a given k we are able to get more consistent results that allow us to fit for more n values than was previously possible.
Future Work • Finish widget program for getting better starting parameters from azimuthally averaged power spectra. • Update Fortran program to carry out multiple ridge fits on both azimuthally averaged data and for full 3d fits. • Work on better characterizing background and what tolerances to allow for deciding on fit.