/***************************************************************************

MwP - description

-------------------

begin : NOV 2006

copyright : (C) 2006 by

email : reymond.d@labogeo.pf

***************************************************************************/

This program is distributed in the hope that it will be useful,

but WITHOUT ANY WARRANTY; without even the implied warranty of

MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the

GNU General Public License for more details.

/***************************************************************************

* This program is free software; you can redistribute it and/or modify *

* it under the terms of the GNU General Public License as published by *

* the Free Software Foundation; either version 2 of the License, or *

* (at your option) any later version. *

***************************************************************************/

send feed-back and bugs to

Dominique Reymond

LDG

po box 640

Papeete

98713 French Polynesia

reymond.d@labogeo.pf

///////////////////////////////////////////////////////////////////////////////////////



MwP is a part of STK Seismic ToolKit - Utility Programs

MwP is a tool for computing the moment magnitude on P wave,

using the principle given by Tsuboi 1997.

+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

It is strongly recommended to read the differents publications about Mwp

before operationel uses.

+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

The correction of propagation used is given by E.A. Okal

It includes geometrical spreading and receiver function;

both are function of epicentral distance.

For details, see "Okal, E.A., A student's guide to teleseismic

body-wave amplitudes, Seismol. Res. Letts., 63, 169-180, 1992."

++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++



The window on P wave is selected automatically from the coordinates of epicenter

and origin time given in arguments. The propagation times are calculated with

the IASPEI time tables.



Be very carefull to instrumental response of the sensor used :

sensors must be very broad-band sensors like STS1 and KS54000, with flat response

in velocity until 360s.

The file STATIONS_STS1.txt contains a list of stations with such instruments;

only the stations belonguing to this list will be processed.

Notice thet STS2 sensors have a cutting period of only 120s and cannot be used

//////////////////////////////////////////////////////////////////////////////////////

The file MwP_0.1.ini contains some processing parameters and the path to

differents files and directories described here-after:



FMIN is the cutt_off frequency of the HIGH-PASS filter, 0 for no filter

FMAX is the cutt_off frequency of the LOW-PASS filter, 0 for no filter

FC_INTEGRATOR 2E-3 : cutt_off frequency of the integrator in Hz

Unit_Factor 1.0E9 : conversion factor for having nanometer

LENGTH_before_P 180 : length in second before P wave

LENGTH_before_S 40 : length in second before S wave

Vp_m/s 7900 : velocity of P wave arround the source

Ro_kg/m3 3400 : density of medium

DIST_MIN 25 : minimum distance between source and station in deg

DIST_MAX 95 : maximum distance between source and station in deg

RESULTS_DIR /tmp/MwP : directory for results

LOG_DIR /tmp/MwP : directory for LOG files

FINAL_RESULTS_FILE /tmp/MwP/Final_Results.txt : file name of ALL the results

DEBUG_FILE /tmp/MwP/debug_Mwp.log : file name of debugging

STATION_LOC /home/domi/Work/Gtk+/STK_0.49/Data/geoloc.sta : file of STATION COORD

TABLES_IASPI_P /home/domi/Work/Gtk+/STK_0.49/Data/IASPEI/P.asp : table of IASPEI propagation of P wave

TABLES_IASPI_S /home/domi/Work/Gtk+/STK_0.49/Data/IASPEI/S.asp : table of IASPEI propagation of S wave

GEOM_SPREAD /home/domi/Work/Gtk+/STK_0.49/Data/P_geo_spread.txt : model for geometrical spreading

BAD_STATIONS /home/domi/Work/Gtk+/STK_0.49/Data/bad_stations_IRIS.txt : file of BAD stations (noisy, down ...)

BB_SENSORS /home/domi/Work/Gtk+/STK_0.49/Data/STATIONS_STS1.txt : file of station with STS1 or KS54000 seismometers



///////////////////////////////////////////////////////////////////////////////////

The results files:

RESULTS_DIR will contains indidual results per event, referenced by YEAR_MONTH_DAY_HH_MM_RES_MWP.txt

inside are the results for each station (no limit on the number of stations)



example: file 1991_04_22_21_56_RES_MWP.txt

#sta lat_sta lon_sta start of file julian day AZIM DIST ° Mo N.m MwP

PAS +34.148 -118.171 22/04/1991 21:28 112 118.85 040.32 2.26E+20 07.50

ANMO +34.946 -106.457 22/04/1991 21:31 112 134.08 032.99 4.34E+20 07.69

PFO +33.609 -116.455 22/04/1991 21:31 112 120.08 038.82 2.23E+20 07.50

ISA +35.663 -118.473 22/04/1991 21:40 112 120.06 041.28 2.33E+20 07.51

MBO +14.391 -016.955 22/04/1991 21:41 112 274.02 064.66 2.60E+19 06.88

TAM +22.791 +005.527 22/04/1991 21:42 112 278.36 085.04 2.89E+20 07.57

COL +64.900 -147.793 22/04/1991 21:43 112 109.15 070.69 1.56E+21 08.06

COR +44.586 -123.303 22/04/1991 21:44 112 122.61 049.10 2.05E+20 07.47

GSC +35.303 -116.808 22/04/1991 21:47 112 121.47 039.93 2.59E+20 07.54

SSB +45.279 +004.542 22/04/1991 21:54 112 275.16 081.54 3.73E+20 07.65

SBC +34.442 -119.713 22/04/1991 21:56 112 117.58 041.59 2.09E+20 07.48

KIP +21.423 -158.015 22/04/1991 21:56 112 086.24 072.58 1.47E+20 07.38

CCM +38.056 -091.245 22/04/1991 21:57 112 163.32 029.19 9.24E+20 07.91



/////////////////////////////////////////////////////////////////////////////

FINAL_RESULTS_FILE will contain the average results for all the station per event

example for 2 events:

dd/mo/yyyy hh:mm jj Mo N.m MwP Nsta

22/04/1991 21:57 112 2.66E+20 07.55 13

13/01/2007 04:34 013 2.24E+21 08.17 37

////////////////////////////////////////////////////////////////////////////////

the file debug_MwP.log contains header informations on each channels for debugging

example:

dt:1.000000 sens:1.998850E+09

npt : 10202

sta : UNM

location :

rate:LH ax : Z

NET : G

date: 22/4/1991 21/25/29

dt:1.000000 sens:4.212000E+09

npt : 9061

sta : PAS

location :

rate:LH ax : Z

NET : TS

date: 22/4/1991 21/28/59

... etc

//////////////////////////////////////////////////////////////////////////////////////

prog. arguments:

./MwP PATH_to_signal_SAC_ASCII DD/MM/YYYY HH:MM:SS LAT LON DEPTH *CHANNEL_SELECTION* DD/MM/YYYY HH:MM:SS LAT LON DEPTH are origin time and coord of earthquake

example:

./MwP /home/Work/Signal/Sac_ascii/1991_04_22_costarica/dotdot 22/04/1991 21:56:52 9.68 -83.08 15 *.LHZ.*

/////////////////////////////////////////////////////////////////////////////////////////