%This is a parameter input file for autoPyLoT.
%All main and special settings regarding data handling
%and picking are to be set here!
%Parameters are optimized for regional data sets!
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

#main settings#
/DATA/Egelados                      #rootpath#     %project path
EVENT_DATA/LOCAL		    #datapath#     %data path
2006.01_Nisyros                     #database#     %name of data base
e1412.008.06                        #eventID#      %event ID for single event processing
PILOT				    #datastructure#  %choose data structure
0                                   #iplot#        %flag for plotting: 0 none, 1, partly, >1 everything
AUTOPHASES_AIC_HOS4_ARH             #phasefile#    %name of autoPILOT output phase file 
AUTOLOC_AIC_HOS4_ARH                #locfile#      %name of autoPILOT output location file 
AUTOFOCMEC_AIC_HOS4_ARH.in          #focmecin#     %name of focmec input file containing polarities
HYPOSAT                             #locrt#        %location routine used ("HYPOINVERSE" or "HYPOSAT")
6                                   #pmin#         %minimum required P picks for location
4                                   #p0min#        %minimum required P picks for location if at least
                                                   %3 excellent P picks are found
2                                   #smin#         %minimum required S picks for location
/home/ludger/bin/run_HYPOSAT4autoPILOT.csh #cshellp# %path and name of c-shell script to run location routine
7.6 8.5                             #blon#         %longitude bounding for location map
49 49.4                             #blat#         %lattitude bounding for location map
#parameters for moment magnitude estimation#
5000                                #vp#           %average P-wave velocity
2800                                #vs#           %average S-wave velocity
2200                                #rho#          %rock density [kg/m^3]
300                                 #Qp#           %quality factor for P waves
100                                 #Qs#           %quality factor for S waves
#common settings picker#
20                                  #pstart#       %start time [s] for calculating CF for P-picking
160                                 #pstop#        %end time [s] for calculating CF for P-picking
3.0                                 #sstart#       %start time [s] after or before(-) P-onset for calculating CF for S-picking
100                                 #sstop#        %end time [s] after P-onset for calculating CF for S-picking
3 10                                #bpz1#         %lower/upper corner freq. of first band pass filter Z-comp. [Hz]
3 12                                #bpz2#         %lower/upper corner freq. of second band pass filter Z-comp. [Hz]
3 8                                 #bph1#         %lower/upper corner freq. of first band pass filter H-comp. [Hz]
3 6                                 #bph2#         %lower/upper corner freq. of second band pass filter H-comp. [Hz]
#special settings for calculating CF#
%!!Be careful when editing the following!!
#Z-component#
HOS                                 #algoP#        %choose algorithm for P-onset determination (HOS, ARZ, or AR3)
7                                   #tlta#         %for HOS-/AR-AIC-picker, length of LTA window [s] 
4                                   #hosorder#     %for HOS-picker, order of Higher Order Statistics
2                                   #Parorder#     %for AR-picker, order of AR process of Z-component
1.2                                 #tdet1z#       %for AR-picker, length of AR determination window [s] for Z-component, 1st pick
0.4                                 #tpred1z#      %for AR-picker, length of AR prediction window [s] for Z-component, 1st pick
0.6                                 #tdet2z#       %for AR-picker, length of AR determination window [s] for Z-component, 2nd pick
0.2                                 #tpred2z#      %for AR-picker, length of AR prediction window [s] for Z-component, 2nd pick
0.001                               #addnoise#     %add noise to seismogram for stable AR prediction
5 0.2 3.0 1.5                       #tsnrz#        %for HOS/AR, window lengths for SNR-and slope estimation [tnoise,tsafetey,tsignal,tslope] [s]
3                                   #pickwinP#     %for initial AIC and refined pick, length of P-pick window [s]
8                                   #Precalcwin#   %for HOS/AR, window length [s] for recalculation of CF (relative to 1st pick)
1.0                                 #aictsmooth#   %for HOS/AR, take average of samples for smoothing of AIC-function [s]
0.3                                 #tsmoothP#     %for HOS/AR, take average of samples for smoothing CF [s]
0.3                                 #ausP#         %for HOS/AR, artificial uplift of samples (aus) of CF (P)
1.3                                 #nfacP#        %for HOS/AR, noise factor for noise level determination (P)
#H-components#
ARH                                 #algoS#        %choose algorithm for S-onset determination (ARH or AR3)
0.8                                 #tdet1h#       %for HOS/AR, length of AR-determination window [s], H-components, 1st pick
0.4                                 #tpred1h#      %for HOS/AR, length of AR-prediction window [s], H-components, 1st pick
0.6                                 #tdet2h#       %for HOS/AR, length of AR-determinaton window [s], H-components, 2nd pick
0.3                                 #tpred2h#      %for HOS/AR, length of AR-prediction window [s], H-components, 2nd pick
4                                   #Sarorder#     %for AR-picker, order of AR process of H-components
10                                  #Srecalcwin#   %for AR-picker, window length [s] for recalculation of CF (2nd pick) (H)
6                                   #pickwinS#     %for initial AIC and refined pick, length of S-pick window [s]
5 0.2 3.0 3.0                       #tsnrh#        %for ARH/AR3, window lengths for SNR-and slope estimation [tnoise,tsafetey,tsignal,tslope] [s]
3.0                                 #aictsmoothS#  %for AIC-picker, take average of samples for smoothing of AIC-function [s]
1.0                                 #tsmoothS#     %for AR-picker, take average of samples for smoothing CF [s] (S)
0.4                                 #ausS#         %for HOS/AR, artificial uplift of samples (aus) of CF (S)
1.5                                 #nfacS#        %for AR-picker, noise factor for noise level determination (S)
%first-motion picker%
1                                   #minfmweight#  %minimum required p weight for first-motion determination
2                                   #minFMSNR#     %miniumum required SNR for first-motion determination
6.0                                 #fmpickwin#    %pick window around P onset for calculating zero crossings
%quality assessment%
#inital AIC onset#
0.04 0.08 0.16 0.32                 #timeerrorsP#   %discrete time errors [s] corresponding to picking weights [0 1 2 3] for P
0.04 0.08 0.16 0.32                 #timeerrorsS#   %discrete time errors [s] corresponding to picking weights [0 1 2 3] for S
5                                   #minAICPslope#  %below this slope [counts/s] the initial P pick is rejected
1.2                                 #minAICPSNR#    %below this SNR the initial P pick is rejected
8                                   #minAICSslope#  %below this slope [counts/s] the initial S pick is rejected
1.5                                 #minAICSSNR#    %below this SNR the initial S pick is rejected
#check duration of signal using envelope function#
6                                   #minsiglength# %minimum required length of signal [s]
1.5                                 #noisefactor#  %noiselevel*noisefactor=threshold
60                                  #minpercent#   %required percentage of samples higher than threshold
#check for spuriously picked S-onsets#
1.5                                 #zfac#         %P-amplitude must exceed at least zfac times RMS-S amplitude 
#check statistics of P onsets#
35                                  #mdttolerance#  %maximum allowed deviation of P picks from median [s]
#wadati check#
2.0                                 #wdttolerance# %maximum allowed deviation from Wadati-diagram