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.013 h= 12.1 hd= 13.6 ed= 6.2 lt= 68.994 lo= 23.384 2010 1008 2139 02.8 L KTK1 N k= -0.023 h= 12.1 hd= 13.6 ed= 6.2 lt= 68.994 lo= 23.384 2010 1008 2139 02.8 L KTK1 E k= -0.023 h= 12.1 hd= 13.6 ed= 6.2 lt= 68.994 lo= 23.384 2010 1008 2139 02.8 L KTK1 Z k= -0.037 h= 0.1 hd= 40.1 ed= 40.1 lt= 69.334 lo= 22.787 2010 1129 0526 17.1 L KTK1 N k= -0.039 h= 0.1 hd= 40.1 ed= 40.1 lt= 69.334 lo= 22.787 2010 1129 0526 17.1 L KTK1 E k= -0.036 h= 0.1 hd= 40.1 ed= 40.1 lt= 69.334 lo= 22.787 2010 1129 0526 17.1 L KTK1 Z k= -0.016 h= 8.5 hd= 13.3 ed= 10.2 lt= 69.062 lo= 23.450 2011 0107 2346 04.6 L KTK1 N k= -0.003 h= 8.5 hd= 13.3 ed= 10.2 lt= 69.062 lo= 23.450 2011 0107 2346 04.6 L KTK1 E k= -0.012 h= 8.5 hd= 13.3 ed= 10.2 lt= 69.062 lo= 23.450 2011 0107 2346 04.6 L LOF Z k= -0.043 h= 0.0 hd= 11.9 ed= 11.9 lt= 68.098 lo= 13.811 2010 0227 1546 40.8 L LOF N k= -0.035 h= 0.0 hd= 11.9 ed= 11.9 lt= 68.098 lo= 13.811 2010 0227 1546 40.8 L LOF E k= -0.039 h= 0.0 hd= 11.9 ed= 11.9 lt= 68.098 lo= 13.811 2010 0227 1546 40.8 L LOF Z k= -0.047 h= 15.0 hd= 17.0 ed= 8.1 lt= 68.128 lo= 13.733 2010 0804 1620 04.5 L LOF N k= -0.039 h= 15.0 hd= 17.0 ed= 8.1 lt= 68.128 lo= 13.733 2010 0804 1620 04.5 L
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