*****README*****
This files contain the source code of Matlab, C-shell for Centroid Moment Tensor inversion process.
Decompress example file Nantou0327.7z first, make sure inversion process can run without problem.
Prepare file for inversion:
1. Data file: [LON LAT dE dN dU sigE sigN sigU X name]
2. Grid file: Grid for inversion
Make Green's function: (Greensfunction.7z)
3. Put Data & Grid file into $PATH/INFO
4. Change path in inv_test.sh file
pathfk: path of fk main directory (Source code from Dr. Lupei Zhu)
gridfile: path of grid file ($PATH/INFO)
datafile: path of data ($PATH/INFO)
-LINE24: velocity model
-LINE33: Units of Green's function
5.Run: ./inv_test.sh
Inversion main program:
* Convert rtz-ENU first to speed the process (program sav_G.m)
6. inver_findmin_res.m: change (1)Data, (2)Grid PATH and (3)n_lon,n_lat,n_z
Results:
mtinv_fom: moment tensor (mxx,mxy,mxz,myy,myz,mzz)
m0: total moment
mw: moment magnitude
dhatCw: model (EEEEE...,NNNNN...,UUUUU...)
NOTE:
misfit: sqrt(sum((dhat-d).^2))
unit of Green function: from fk (default: 3.3e20 dyne-cm)
moment tensor to strike,dip,rake:
-tools/mt2dc ->compile fortran code (csh compile.csh)
-input moment tensor (M1~M6) in mt_inp file
-./get_focalm -> output tmp_inf file