%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 local data sets!
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
#main settings#
/DATA/Insheim                       #rootpath#     %project path
EVENT_DATA/LOCAL		    #datapath#     %data path
2016.08_Insheim                     #database#     %name of data base
e0007.224.16                        #eventID#      %event ID for single event processing
/DATA/Insheim/STAT_INFO             #invdir#       %full path to inventory or dataless-seed file
PILOT				    #datastructure#%choose data structure
0                                   #iplot#        %flag for plotting: 0 none, 1 partly, >1 everything
True                                #apverbose#    %choose 'True' or 'False' for terminal output
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
#NLLoc settings#
/home/ludger/NLLOC                  #nllocbin#     %path to NLLoc executable
/home/ludger/NLLOC/Insheim          #nllocroot#    %root of NLLoc-processing directory
AUTOPHASES.obs                      #phasefile#    %name of autoPyLoT-output phase file for NLLoc
                                                   %(in nllocroot/obs)
Insheim_min1d032016_auto.in         #ctrfile#      %name of autoPyLoT-output control file for NLLoc
                                                   %(in nllocroot/run)
ttime                               #ttpatter#     %pattern of NLLoc ttimes from grid
                                                   %(in nllocroot/times)
AUTOLOC_nlloc                       #outpatter#    %pattern of NLLoc-output file
                                                   %(returns 'eventID_outpatter')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
#parameters for seismic moment estimation#
3530                                #vp#           %average P-wave velocity
2500                                #rho#          %average rock density [kg/m^3]
300 0.8                             #Qp#           %quality factor for P waves ([Qp, ap], Qp*f^a)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
AUTOFOCMEC_AIC_HOS4_ARH.in          #focmecin#     %name of focmec input file containing derived polarities
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
#common settings picker#
15.0                                #pstart#       %start time [s] for calculating CF for P-picking
60.0                                #pstop#        %end time [s] for calculating CF for P-picking
-1.0                                #sstart#       %start time [s] relative to P-onset for calculating CF for S-picking
10.0                                #sstop#        %end time [s] after P-onset for calculating CF for S-picking
2 20                                #bpz1#         %lower/upper corner freq. of first band pass filter Z-comp. [Hz]
2 30                                #bpz2#         %lower/upper corner freq. of second band pass filter Z-comp. [Hz]
2 15                                #bph1#         %lower/upper corner freq. of first band pass filter H-comp. [Hz]
2 20                                #bph2#         %lower/upper corner freq. of second band pass filter z-comp. [Hz]
#special settings for calculating CF#
%!!Edit the following only if you know what you are doing!!%
#Z-component#
HOS                                 #algoP#        %choose algorithm for P-onset determination (HOS, ARZ, or AR3)
7.0                                 #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
3 0.1 0.5 0.5                       #tsnrz#        %for HOS/AR, window lengths for SNR-and slope estimation [tnoise,tsafetey,tsignal,tslope] [s]
3.0                                 #pickwinP#     %for initial AIC pick, length of P-pick window [s]
6.0                                 #Precalcwin#   %for HOS/AR, window length [s] for recalculation of CF (relative to 1st pick)
0.2                                 #aictsmooth#   %for HOS/AR, take average of samples for smoothing of AIC-function [s]
0.1                                 #tsmoothP#     %for HOS/AR, take average of samples for smoothing CF [s]
0.001                               #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
5.0                                 #Srecalcwin#   %for AR-picker, window length [s] for recalculation of CF (2nd pick) (H)
3.0                                 #pickwinS#     %for initial AIC pick, length of S-pick window [s]
2 0.2 1.5 0.5                       #tsnrh#        %for ARH/AR3, window lengths for SNR-and slope estimation [tnoise,tsafetey,tsignal,tslope] [s]
0.5                                 #aictsmoothS#  %for AIC-picker, take average of samples for smoothing of AIC-function [s]
0.7                                 #tsmoothS#     %for AR-picker, take average of samples for smoothing CF [s] (S)
0.9                                 #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
0.2                                 #fmpickwin#    %pick window around P onset for calculating zero crossings
%quality assessment%
#inital AIC onset#
0.05 0.10 0.20 0.40                 #timeerrorsP#   %discrete time errors [s] corresponding to picking weights [0 1 2 3] for P
0.10 0.20 0.40 0.80                 #timeerrorsS#   %discrete time errors [s] corresponding to picking weights [0 1 2 3] for S
4                                   #minAICPslope#  %below this slope [counts/s] the initial P pick is rejected
1.2                                 #minAICPSNR#    %below this SNR the initial P pick is rejected
2                                   #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#
3                                   #minsiglength# %minimum required length of signal [s]
1.0                                 #noisefactor#  %noiselevel*noisefactor=threshold
40                                  #minpercent#   %required percentage of samples higher than threshold
#check for spuriously picked S-onsets#
2.0                                 #zfac#         %P-amplitude must exceed at least zfac times RMS-S amplitude 
#check statistics of P onsets#
2.5                                 #mdttolerance#  %maximum allowed deviation of P picks from median [s]
#wadati check#
1.0                                 #wdttolerance# %maximum allowed deviation from Wadati-diagram

