Next: 33. Seismic risk related programs Up: 32. Calculating spectra and attenuation Previous: 32.2 Calculating kappa, KAPPA Contents Index
The SPEC and KAPPA programs gives several kappa related outputs. KAPPA_SORT can further sort and analyze the data and prepare outputs that can be used for plotting the results. The main input file to the program is spec_kappa_chan.all (from SPEC) or kappa_chan.all (from KAPPA) that contains all the individual kappa results from each channel. The options are:
1: get all 3 components from same event, but only if all 3 are there the input file is hardwired to be kappa_chan.all or spec_kappa_chan.all. It will first look for kappa_chan.all and if not found look for spec_kappa_chan.all. For some analysis it is important to have only the data if available at all 3 channels. This is done every time the program runs. The output is in kappa_sort.out, example output:
KTK1 Z k= 0.020 h= 9.6 hd= 17.5 ed= 14.6 lt= 68.895 lo= 23.403 baz= 152.7 2013 1026 KTK1 N k= 0.024 h= 9.6 hd= 17.5 ed= 14.6 lt= 68.895 lo= 23.403 baz= 152.7 2013 1026 KTK1 E k= 0.016 h= 9.6 hd= 17.5 ed= 14.6 lt= 68.895 lo= 23.403 baz= 152.7 2013 1026 KTK1 Z k= 0.007 h= 3.0 hd= 8.0 ed= 7.4 lt= 69.076 lo= 23.279 baz= 11.3 2015 0214 KTK1 N k= 0.047 h= 3.0 hd= 8.0 ed= 7.4 lt= 69.076 lo= 23.279 baz= 11.3 2015 0214 KTK1 E k= 0.006 h= 3.0 hd= 8.0 ed= 7.4 lt= 69.076 lo= 23.279 baz= 11.3 2015 0214 KTK1 Z k= 0.044 h= 7.4 hd= 14.3 ed= 12.2 lt= 69.084 lo= 23.467 baz= 48.8 2015 0405 KTK1 N k= 0.048 h= 7.4 hd= 14.3 ed= 12.2 lt= 69.084 lo= 23.467 baz= 48.8 2015 0405 KTK1 E k= 0.024 h= 7.4 hd= 14.3 ed= 12.2 lt= 69.084 lo= 23.467 baz= 48.8 2015 0405 KTK1 Z k= 0.027 h= 6.1 hd= 21.0 ed= 20.1 lt= 68.833 lo= 23.308 baz= 170.2 2015 0501 KTK1 N k= 0.028 h= 6.1 hd= 21.0 ed= 20.1 lt= 68.833 lo= 23.308 baz= 170.2 2015 0501 KTK1 E k= 0.009 h= 6.1 hd= 21.0 ed= 20.1 lt= 68.833 lo= 23.308 baz= 170.2 2015 0501 LOF Z k= 0.016 h= 6.5 hd= 27.5 ed= 26.7 lt= 68.025 lo= 12.967 baz= 243.5 2013 0116 LOF N k= 0.015 h= 6.5 hd= 27.5 ed= 26.7 lt= 68.025 lo= 12.967 baz= 243.5 2013 0116 LOF E k= 0.017 h= 6.5 hd= 27.5 ed= 26.7 lt= 68.025 lo= 12.967 baz= 243.5 2013 0116 LOF Z k= 0.034 h= 3.7 hd= 8.1 ed= 7.2 lt= 68.120 lo= 13.369 baz= 259.1 2013 0124 LOF N k= 0.077 h= 3.7 hd= 8.1 ed= 7.2 lt= 68.120 lo= 13.369 baz= 259.1 2013 0124 LOF E k= 0.063 h= 3.7 hd= 8.1 ed= 7.2 lt= 68.120 lo= 13.369 baz= 259.1 2013 0124
2: optionally make xy files of kappa vs distance and depth for a particular Channel. These xy files can e.g. show if any bias of kappa with distance. Three files are made: kappa vs epi distance, kappa vs hyp distance and kappa vs depth. The files can be used e.g. with program LSQ.
3: optionally make a station file for GMT where kappa has been added to the station code. Input can be from an epimap.are where the spec kappa Nordic files have been use for plotting so kappa is included in the file (spec_kappa.z, spec_kappa.n or spec_kappa.e ) or one of the three kappa Nordic files directly. The average for each componnt as well as the overall average is calculated. This can be to calcualte average kappa in a particular region selecte in EPIMAP using the Kappa Nordic file. Example of input in spec_kappa.z:
2000 0101 0101 0.0 L 69.012 23.237 3 2.2 1 KTK1 Z -0.022 0.013 3 2000 0101 0101 0.0 L 68.132 13.540 2 4.5 1 LOF Z -0.045 0.003 3 2000 0101 0101 0.0 L 67.930 15.242 2 5.1 1 STEI Z -0.051 0.023 3 2000 0101 0101 0.0 L 66.499 13.118 200 5.3 1 KONS Z -0.053 0.012 3 2000 0101 0101 0.0 L 66.333 13.018 47 5.4 1 STOK Z -0.054 0.016 3
Kappa is given on second line and also on header line as a magnitude. Example of output GMT file is kappa_station.gmt:
23.237 69.012 KTK1 22 13.540 68.132 LOF 45 15.242 67.930 STEI 51 13.118 66.499 KONS 53 13.018 66.333 STOK 54