Next: 27.3 Program QSTAT Up: 27. Calculation of coda Previous: 27.1 Program CODAQ_AREA   Contents   Index

# 27.2 Program CODAQ_POLYGON

Using EPIMAP, it is possible to select out events in a polygon when plotting the codaQ midpoints (or other epicentres). It is then possible to calculate the average Q in that area using the program CODAQ_POLYGON. The input polygon coordinates are assumed to be present in file epimap.cor, the one created by EPIMAP when selecting in a polygon. The codaQ values in the midpoints are supposed to be present in file codaq.area and no questions are asked. File coda1.out is also used as input and must be present.

The data in the polygon file is written out in a file POLYGON.MAP so it can be plotted again with EPIMAP. The POLYGON.MAP file is appended to so if CODAQ_AREA is used several times with different polygon areas, they will all be stored in POLYGON.MAP read for plotting with EPIMAP. It is up to the user to ensure the POLYGON.MAP is empty if a new series of polygons are stored.

The section of the codaq.area corresponding to the polygon area is output in codaq_polygon.area and it can be used with QSTAT for statistics of the data in the polygon area.

Example run

codaq_polygon

c:\jens\Q-port\ALL>codaq_polygon
Frequencies    1.00    2.00    4.00    8.00   16.00
Number of points defining polygon:     7
f=   1.00000000

f=   2.00000000
216.   269.   308.
f=   4.00000000
696.   450.   602.   776.   677.   494.   614.   489.   707.   643.
764.   687.   649.   573.   816.   841.   834.   518.   536.   524.
786.
f=   8.00000000
1366.   878.   968.  1295.   803.   999.   979.   967.  1155.  1080.
940.  1110.  1042.   994.   992.   999.  1023.  1020.  1009.   853.
1177.  1007.  1316.  1186.   777.   898.   959.  1052.  1097.   945.
1030.  1206.  1123.  1008.  1140.  1062.  1455.  1029.  1125.   809.
1298.   871.  1386.  1231.  1175.  1100.   771.  1171.   891.   897.
740.  1286.  1352.  1108.  1037.  1240.
f=   16.0000000
1667.  1514.  1648.  1638.  1898.  1777.  1381.  1667.  1620.  1810.
1220.  1251.  1543.  1344.  1436.  1520.  1772.  1575.  1586.  1455.
1666.  1534.  1491.  1696.  1514.  1575.  1690.  1597.  1352.  1450.
1358.  1672.  1674.  1514.  1633.  1543.  1392.  1552.  1619.  1474.
1647.  1425.  1510.  1456.  1374.  1391.  1389.  1366.  1657.  1449.
1342.  1689.  1238.  1667.  2297.  1573.  1972.  1463.  1405.  1449.
Number of Q-midpoints within polygon:      140
f= 2.0  av= 259.0  sd=  46.9  n=    3
f= 4.0  av= 628.9  sd= 119.7  n=   21
f= 8.0  av=1035.7  sd= 163.2  n=   56
f=16.0  av=1531.8  sd= 167.4  n=   60

q0=  172.7  sd 38.8    qalpha= 0.81  sd 0.10   q10=1107.8  sd=182.1   corr= 0.98  rms= 0.00  nf= 4

Use all data, no weight
q0=  166.0  sd 36.8    qalpha= 0.84  sd 0.12   q10=1151.4  sd=188.4  corr= 0.98  rms= 0.00  nf= 4

Append to POLYGON.MAP (n/y=default)

Accumulated polygon coordinates in POLYGON.MAP
Area file selected is codaq_polygon.area

For each frequency, all the Q-values are printed so it is possible to get an idea of the variation of the values. The averages are calculated for each frequency if 2 or more values are available. The Q-realtion is then

calculated if averages are available for at least 3 frequecies. Finally, a Q-relation is calculated if also one value is available, but then there is no weighting with the standard deviations of the averages.

How to calculate coda Q in polygon areas

Calculate coda Q in large area of interest.

indexCoda Q in a polygon Plot all the midpoints between station end events with EPIMAP using command map codaq.area.nor . The input epicentre file is codaq.area.nor which is output from CODAQ.

Do the following for each polygon area:

With EPIMAP, use polygon function to select an area. Terminate EPIMAP.

Run codaq_polygon to calculate coda Q for that area, the polygon outline is stored in file

POLYGON.MAP for later plotting.

After e.g. 5 runs, the 5 areas are stored in POLYGON.MAP and the Q-relations have been printed out.

Next: 27.3 Program QSTAT Up: 27. Calculation of coda Previous: 27.1 Program CODAQ_AREA   Contents   Index
Peter Voss : Fri Jun 7 07:35:13 UTC 2019