Software initiation

README

@@ -0,0 +1,55 @@
+Infomation fo relocating seismic scatterers Code
+http://seis.ustc.edu.cn/PKPprecursor
+
+Copyright (c) 2000-2015, L. Wen
+###########################
+The package is free and distributed in the hope that it will be useful. 
+You can use it and/or modify it under the terms of the GNU General Public License 
+for NON-COMMERCIAL purposes.
+
+It does not come with any warranties, nor is it guaranteed to work on your computer.
+The user assumes full responsibility for the use of this system.
+We are NOT responsible for any damage that may follow from correct or incorrect
+use of this package.
+
+Permission is granted to make and distribute verbatim copies of this license provided
+that the copyright notice and these paragraphs are preserved on all copies.
+
+
+Description
+===========
+
+This package is a code developed by Lianxing Wen
+for determining the location of seismic scaterers
+using PKP precursor onset times.
+
+Reference
+=========
+
+Wen, L., 2000. Intense seismic scattering near the Earth's core-mantle boundary beneath
+the Comoros hotspot. Geophys. Res. Lett. 27,3627-3630
+
+Requirement
+===========
+
+#. Seismic Analysis Code (SAC)
+
+   You need SAC installed on your system, and
+   $(SACHOME) set to the SAC installation directory.
+
+#. Generic Mapping Tools
+
+   Some figures are plotted using GMT, it is not necessary.
+
+Installation
+============
+
+Go to src directory, and run make.
+
+Examples
+========
+
+Just change to the example/Run diretory and run:
+
+    perl LocateScatterer.pl
+

example/Run/.gmtcommands4

@@ -0,0 +1,20 @@
+# GMT common arguments shelf
+-B20f10/a20f10wSnE
+-JQ195/30/6.5i
+-JX5/4d
+-Jx0.09d
+-R230/290/20/70
+-X1.5
+-x-1.5
+-Y5
+-y-5
+-jx0.09d
+EOF
+
+
+
+
+OF
+ursors/Test1/patchs2891
+-jx0.12d
+EOF

example/Run/.gmtdefaults4

@@ -0,0 +1,107 @@
+#
+#	GMT-SYSTEM 4.5.5 [64-bit] Defaults file
+#
+#-------- Plot Media Parameters -------------
+PAGE_COLOR		= white
+PAGE_ORIENTATION	= landscape
+PAPER_MEDIA		= a4
+#-------- Basemap Annotation Parameters ------
+ANNOT_MIN_ANGLE		= 20
+ANNOT_MIN_SPACING	= 0
+ANNOT_FONT_PRIMARY	= Helvetica
+ANNOT_FONT_SIZE_PRIMARY	= 14p
+ANNOT_OFFSET_PRIMARY	= 0.2c
+ANNOT_FONT_SECONDARY	= Helvetica
+ANNOT_FONT_SIZE_SECONDARY	= 16p
+ANNOT_OFFSET_SECONDARY	= 0.2c
+DEGREE_SYMBOL		= ring
+HEADER_FONT		= Helvetica
+HEADER_FONT_SIZE	= 36p
+HEADER_OFFSET		= 0.5c
+LABEL_FONT		= Helvetica
+LABEL_FONT_SIZE		= 24p
+LABEL_OFFSET		= 0.3c
+OBLIQUE_ANNOTATION	= 1
+PLOT_CLOCK_FORMAT	= hh:mm:ss
+PLOT_DATE_FORMAT	= yyyy-mm-dd
+PLOT_DEGREE_FORMAT	= ddd:mm:ss
+Y_AXIS_TYPE		= hor_text
+#-------- Basemap Layout Parameters ---------
+BASEMAP_AXES		= WESN
+BASEMAP_FRAME_RGB	= black
+BASEMAP_TYPE		= fancy
+FRAME_PEN		= 1.25p
+FRAME_WIDTH		= 0.2c
+GRID_CROSS_SIZE_PRIMARY	= 0c
+GRID_PEN_PRIMARY	= 0.25p
+GRID_CROSS_SIZE_SECONDARY	= 0c
+GRID_PEN_SECONDARY	= 0.5p
+MAP_SCALE_HEIGHT	= 0.2c
+POLAR_CAP		= 85/90
+TICK_LENGTH		= 0.2c
+TICK_PEN		= 0.5p
+X_AXIS_LENGTH		= 25c
+Y_AXIS_LENGTH		= 15c
+X_ORIGIN		= 2.5c
+Y_ORIGIN		= 2.5c
+UNIX_TIME		= FALSE
+UNIX_TIME_POS		= BL/-2c/-2c
+UNIX_TIME_FORMAT	= %Y %b %d %H:%M:%S
+#-------- Color System Parameters -----------
+COLOR_BACKGROUND	= black
+COLOR_FOREGROUND	= white
+COLOR_NAN		= 128
+COLOR_IMAGE		= adobe
+COLOR_MODEL		= rgb
+HSV_MIN_SATURATION	= 1
+HSV_MAX_SATURATION	= 0.1
+HSV_MIN_VALUE		= 0.3
+HSV_MAX_VALUE		= 1
+#-------- PostScript Parameters -------------
+CHAR_ENCODING		= ISOLatin1+
+DOTS_PR_INCH		= 300
+GLOBAL_X_SCALE		= 1
+GLOBAL_Y_SCALE		= 1
+N_COPIES		= 1
+PS_COLOR		= rgb
+PS_IMAGE_COMPRESS	= lzw
+PS_IMAGE_FORMAT		= ascii
+PS_LINE_CAP		= butt
+PS_LINE_JOIN		= miter
+PS_MITER_LIMIT		= 35
+PS_VERBOSE		= FALSE
+TRANSPARENCY		= 0
+#-------- I/O Format Parameters -------------
+D_FORMAT		= %.12lg
+FIELD_DELIMITER		= tab
+GRIDFILE_FORMAT		= nf
+GRIDFILE_SHORTHAND	= FALSE
+INPUT_CLOCK_FORMAT	= hh:mm:ss
+INPUT_DATE_FORMAT	= yyyy-mm-dd
+IO_HEADER		= FALSE
+N_HEADER_RECS		= 1
+NAN_RECORDS		= pass
+OUTPUT_CLOCK_FORMAT	= hh:mm:ss
+OUTPUT_DATE_FORMAT	= yyyy-mm-dd
+OUTPUT_DEGREE_FORMAT	= D
+XY_TOGGLE		= FALSE
+#-------- Projection Parameters -------------
+ELLIPSOID		= WGS-84
+MAP_SCALE_FACTOR	= default
+MEASURE_UNIT		= cm
+#-------- Calendar/Time Parameters ----------
+TIME_FORMAT_PRIMARY	= full
+TIME_FORMAT_SECONDARY	= full
+TIME_EPOCH		= 2000-01-01T12:00:00
+TIME_IS_INTERVAL	= OFF
+TIME_INTERVAL_FRACTION	= 0.5
+TIME_LANGUAGE		= us
+TIME_UNIT		= d
+TIME_WEEK_START		= Sunday
+Y2K_OFFSET_YEAR		= 1950
+#-------- Miscellaneous Parameters ----------
+HISTORY			= TRUE
+INTERPOLANT		= akima
+LINE_STEP		= 0.025c
+VECTOR_SHAPE		= 0
+VERBOSE			= FALSE

example/Run/CMBout2091

@@ -0,0 +1,116 @@
+1.484452e+04 9.595043e+02
+1.485564e+04 9.598629e+02
+1.486676e+04 9.602166e+02
+1.487788e+04 9.605672e+02
+1.488900e+04 9.609133e+02
+1.490012e+04 9.612570e+02
+1.491124e+04 9.615984e+02
+1.492236e+04 9.619393e+02
+1.493348e+04 9.622763e+02
+1.494460e+04 9.626114e+02
+1.495572e+04 9.629446e+02
+1.496684e+04 9.632738e+02
+1.497796e+04 9.636033e+02
+1.498908e+04 9.639311e+02
+1.500020e+04 9.642573e+02
+1.501131e+04 9.645820e+02
+1.502243e+04 9.649052e+02
+1.503355e+04 9.652244e+02
+1.504467e+04 9.655443e+02
+1.505579e+04 9.658629e+02
+1.506691e+04 9.661802e+02
+1.507803e+04 9.664935e+02
+1.508915e+04 9.668077e+02
+1.510027e+04 9.671207e+02
+1.511139e+04 9.674329e+02
+1.512251e+04 9.677403e+02
+1.513363e+04 9.680492e+02
+1.514475e+04 9.683574e+02
+1.515587e+04 9.686608e+02
+1.516699e+04 9.689661e+02
+1.517811e+04 9.692703e+02
+1.518923e+04 9.695699e+02
+1.520035e+04 9.698716e+02
+1.521147e+04 9.701692e+02
+1.522259e+04 9.704682e+02
+1.523370e+04 9.707663e+02
+1.524482e+04 9.710601e+02
+1.525594e+04 9.713558e+02
+1.526706e+04 9.716476e+02
+1.527818e+04 9.719407e+02
+1.528930e+04 9.722332e+02
+1.530042e+04 9.725212e+02
+1.531154e+04 9.728114e+02
+1.532266e+04 9.730974e+02
+1.533378e+04 9.733848e+02
+1.534490e+04 9.736719e+02
+1.535602e+04 9.739541e+02
+1.536714e+04 9.742391e+02
+1.537826e+04 9.745194e+02
+1.538938e+04 9.748017e+02
+1.540050e+04 9.750801e+02
+1.541162e+04 9.753603e+02
+1.542274e+04 9.756402e+02
+1.543386e+04 9.759147e+02
+1.544498e+04 9.761926e+02
+1.545609e+04 9.764653e+02
+1.546721e+04 9.767409e+02
+1.547833e+04 9.770117e+02
+1.548945e+04 9.772852e+02
+1.550057e+04 9.775541e+02
+1.551169e+04 9.778254e+02
+1.552281e+04 9.780924e+02
+1.553393e+04 9.783615e+02
+1.554505e+04 9.786266e+02
+1.555617e+04 9.788938e+02
+1.556729e+04 9.791571e+02
+1.557841e+04 9.794223e+02
+1.558953e+04 9.796836e+02
+1.560065e+04 9.799468e+02
+1.561177e+04 9.802100e+02
+1.562289e+04 9.804676e+02
+1.563401e+04 9.807252e+02
+1.564513e+04 9.809844e+02
+1.565625e+04 9.812436e+02
+1.566736e+04 9.814975e+02
+1.567848e+04 9.817548e+02
+1.568960e+04 9.820068e+02
+1.570072e+04 9.822625e+02
+1.571184e+04 9.825126e+02
+1.572296e+04 9.827662e+02
+1.573408e+04 9.830145e+02
+1.574520e+04 9.832664e+02
+1.575632e+04 9.835128e+02
+1.576744e+04 9.837629e+02
+1.577856e+04 9.840074e+02
+1.578968e+04 9.842557e+02
+1.580080e+04 9.844984e+02
+1.581192e+04 9.847411e+02
+1.582304e+04 9.849860e+02
+1.583416e+04 9.852268e+02
+1.584528e+04 9.854697e+02
+1.585640e+04 9.857087e+02
+1.586752e+04 9.859500e+02
+1.587864e+04 9.861871e+02
+1.588975e+04 9.864268e+02
+1.590087e+04 9.866621e+02
+1.591199e+04 9.869002e+02
+1.592311e+04 9.871337e+02
+1.593423e+04 9.873700e+02
+1.594535e+04 9.876016e+02
+1.595647e+04 9.878364e+02
+1.596759e+04 9.880662e+02
+1.597871e+04 9.882994e+02
+1.598983e+04 9.885274e+02
+1.600095e+04 9.887590e+02
+1.601207e+04 9.889851e+02
+1.602319e+04 9.892151e+02
+1.603431e+04 9.894394e+02
+1.604543e+04 9.896637e+02
+1.605655e+04 9.898902e+02
+1.606767e+04 9.901127e+02
+1.607879e+04 9.903378e+02
+1.608991e+04 9.905584e+02
+1.610103e+04 9.907821e+02
+1.611214e+04 9.910009e+02
+1.612326e+04 9.912230e+02