pylot/inputs/pylot_regional.in

%This is a parameter input file for PyLoT/autoPyLoT.
%All main and special settings regarding data handling
%and picking are to be set here!
%Parameters are optimized for %extent data sets!
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
#main settings#
                                                   #rootpath#       %project path                                                                                        
                                                   #datapath#       %data path                                                                                           
                                                   #database#       %name of data base                                                                                   
                                                   #eventID#        %event ID for single event processing (* for all events found in database)                           
                                                   #invdir#         %full path to  inventory or dataless-seed file                                                       
PILOT                                              #datastructure#  %choose data structure                                                                               
True                                               #apverbose#      %choose 'True' or 'False' for terminal output                                                        
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
#NLLoc settings#
None                                               #nllocbin#       %path to NLLoc executable                                                                            
None                                               #nllocroot#      %root of NLLoc-processing directory                                                                  
None                                               #phasefile#      %name of autoPyLoT-output phase file for NLLoc                                                       
None                                               #ctrfile#        %name of autoPyLoT-output control file for NLLoc                                                     
ttime                                              #ttpatter#       %pattern of NLLoc ttimes from grid                                                                   
AUTOLOC_nlloc                                      #outpatter#      %pattern of NLLoc-output file                                                                        
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
#parameters for seismic moment estimation#
3530.0                                             #vp#             %average P-wave velocity                                                                             
2500.0                                             #rho#            %average rock density [kg/m^3]                                                                       
300.0 0.8                                          #Qp#             %quality factor for P waves (Qp*f^a); list(Qp, a)                                                    
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
#settings local magnitude#
1.11 0.0009 -2.0                                   #WAscaling#      %Scaling relation (log(Ao)+Alog(r)+Br+C) of Wood-Anderson amplitude Ao [nm]                           If zeros are set, original Richter magnitude is calculated!
1.0382 -0.447                                      #magscaling#     %Scaling relation for derived local magnitude [a*Ml+b].                            If zeros are set, no scaling of network magnitude is applied!
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
#filter settings#
1.0 1.0                                            #minfreq#        %Lower filter frequency [P, S]                                                                       
10.0 10.0                                          #maxfreq#        %Upper filter frequency [P, S]                                                                       
2 2                                                #filter_order#   %filter order [P, S]                                                                                 
bandpass bandpass                                  #filter_type#    %filter type (bandpass, bandstop, lowpass, highpass) [P, S]                                          
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
#common settings picker#
local                                              #extent#         %extent of array ("local", "regional" or "global")                                                   
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                                         
True                                               #use_taup#       %use estimated traveltimes from TauPy for calculating windows for CF                                 
iasp91                                             #taup_model#     %define TauPy model for traveltime estimation                                                        
2.0 10.0                                           #bpz1#           %lower/upper corner freq. of first band pass filter Z-comp. [Hz]                                     
2.0 12.0                                           #bpz2#           %lower/upper corner freq. of second band pass filter Z-comp. [Hz]                                    
2.0 8.0                                            #bph1#           %lower/upper corner freq. of first band pass filter H-comp. [Hz]                                     
2.0 10.0                                           #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 0.1 0.5 1.0                                    #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)                             
4.0                                                #pickwinS#       %for initial AIC pick, length of S-pick window [s]                                                   
2.0 0.3 1.5 1.0                                    #tsnrh#          %for ARH/AR3, window lengths for SNR-and slope estimation [tnoise, tsafetey, tsignal, tslope] [s]    
1.0                                                #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.0                                                #minFMSNR#       %miniumum required SNR for first-motion determination                                                
0.2                                                #fmpickwin#      %pick window around P onset for calculating zero crossings                                           
#quality assessment#
0.02 0.04 0.08 0.16                                #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                           
0.8                                                #minAICPslope#   %below this slope [counts/s] the initial P pick is rejected                                          
1.1                                                #minAICPSNR#     %below this SNR the initial P pick is rejected                                                       
1.0                                                #minAICSslope#   %below this slope [counts/s] the initial S pick is rejected                                          
1.5                                                #minAICSSNR#     %below this SNR the initial S pick is rejected                                                       
1.0                                                #minsiglength#   %length of signal part for which amplitudes must exceed noiselevel [s]                               
1.0                                                #noisefactor#    %noiselevel*noisefactor=threshold                                                                    
10.0                                               #minpercent#     %required percentage of amplitudes exceeding threshold                                               
1.5                                                #zfac#           %P-amplitude must exceed at least zfac times RMS-S amplitude                                         
6.0                                                #mdttolerance#   %maximum allowed deviation of P picks from median [s]                                                
1.0                                                #wdttolerance#   %maximum allowed deviation from Wadati-diagram                                                       
5.0                                                #jackfactor#     %pick is removed if the variance of the subgroup with the pick removed is larger than the mean variance of all subgroups times safety factor
[add] latest pylot.in file added to inputs! 2017-07-10 15:49:15 +02:00			`%This is a parameter input file for PyLoT/autoPyLoT.`
[new] pylot.in which will substitute other input files later a more general input file might help prevent ambiguity among use cases 2016-09-08 13:45:16 +02:00			`%All main and special settings regarding data handling`
			`%and picking are to be set here!`
[add] latest pylot.in file added to inputs! 2017-07-10 15:49:15 +02:00			`%Parameters are optimized for %extent data sets!`
[new] pylot.in which will substitute other input files later a more general input file might help prevent ambiguity among use cases 2016-09-08 13:45:16 +02:00			`%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%`
			`#main settings#`
[update] default pylot.in file updates 2017-09-21 15:19:15 +02:00			`#rootpath# %project path`
			`#datapath# %data path`
			`#database# %name of data base`
			`#eventID# %event ID for single event processing (* for all events found in database)`
			`#invdir# %full path to inventory or dataless-seed file`
[add] latest pylot.in file added to inputs! 2017-07-10 15:49:15 +02:00			`PILOT #datastructure# %choose data structure`
			`True #apverbose# %choose 'True' or 'False' for terminal output`
			`%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%`
[new] pylot.in which will substitute other input files later a more general input file might help prevent ambiguity among use cases 2016-09-08 13:45:16 +02:00			`#NLLoc settings#`
[update] default pylot.in file updates 2017-09-21 15:19:15 +02:00			`None #nllocbin# %path to NLLoc executable`
			`None #nllocroot# %root of NLLoc-processing directory`
			`None #phasefile# %name of autoPyLoT-output phase file for NLLoc`
			`None #ctrfile# %name of autoPyLoT-output control file for NLLoc`
[add] latest pylot.in file added to inputs! 2017-07-10 15:49:15 +02:00			`ttime #ttpatter# %pattern of NLLoc ttimes from grid`
			`AUTOLOC_nlloc #outpatter# %pattern of NLLoc-output file`
			`%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%`
[new] pylot.in which will substitute other input files later a more general input file might help prevent ambiguity among use cases 2016-09-08 13:45:16 +02:00			`#parameters for seismic moment estimation#`
[add] latest pylot.in file added to inputs! 2017-07-10 15:49:15 +02:00			`3530.0 #vp# %average P-wave velocity`
			`2500.0 #rho# %average rock density [kg/m^3]`
			`300.0 0.8 #Qp# %quality factor for P waves (Qp*f^a); list(Qp, a)`
			`%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%`
			`#settings local magnitude#`
[update] default pylot.in file updates 2017-09-21 15:19:15 +02:00			`1.11 0.0009 -2.0 #WAscaling# %Scaling relation (log(Ao)+Alog(r)+Br+C) of Wood-Anderson amplitude Ao [nm] If zeros are set, original Richter magnitude is calculated!`
			`1.0382 -0.447 #magscaling# %Scaling relation for derived local magnitude [a*Ml+b]. If zeros are set, no scaling of network magnitude is applied!`
[new] pylot.in which will substitute other input files later a more general input file might help prevent ambiguity among use cases 2016-09-08 13:45:16 +02:00			`%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%`
[add] latest pylot.in file added to inputs! 2017-07-10 15:49:15 +02:00			`#filter settings#`
[update] default pylot.in file updates 2017-09-21 15:19:15 +02:00			`1.0 1.0 #minfreq# %Lower filter frequency [P, S]`
			`10.0 10.0 #maxfreq# %Upper filter frequency [P, S]`
			`2 2 #filter_order# %filter order [P, S]`
[add] latest pylot.in file added to inputs! 2017-07-10 15:49:15 +02:00			`bandpass bandpass #filter_type# %filter type (bandpass, bandstop, lowpass, highpass) [P, S]`
[new] pylot.in which will substitute other input files later a more general input file might help prevent ambiguity among use cases 2016-09-08 13:45:16 +02:00			`%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%`
			`#common settings picker#`
[update] default pylot.in file updates 2017-09-21 15:19:15 +02:00			`local #extent# %extent of array ("local", "regional" or "global")`
			`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`
			`True #use_taup# %use estimated traveltimes from TauPy for calculating windows for CF`
			`iasp91 #taup_model# %define TauPy model for traveltime estimation`
			`2.0 10.0 #bpz1# %lower/upper corner freq. of first band pass filter Z-comp. [Hz]`
			`2.0 12.0 #bpz2# %lower/upper corner freq. of second band pass filter Z-comp. [Hz]`
			`2.0 8.0 #bph1# %lower/upper corner freq. of first band pass filter H-comp. [Hz]`
			`2.0 10.0 #bph2# %lower/upper corner freq. of second band pass filter z-comp. [Hz]`
[new] pylot.in which will substitute other input files later a more general input file might help prevent ambiguity among use cases 2016-09-08 13:45:16 +02:00			`#special settings for calculating CF#`
			`%!!Edit the following only if you know what you are doing!!%`
			`#Z-component#`
[add] latest pylot.in file added to inputs! 2017-07-10 15:49:15 +02:00			`HOS #algoP# %choose algorithm for P-onset determination (HOS, ARZ, or AR3)`
[update] default pylot.in file updates 2017-09-21 15:19:15 +02:00			`7.0 #tlta# %for HOS-/AR-AIC-picker, length of LTA window [s]`
[add] latest pylot.in file added to inputs! 2017-07-10 15:49:15 +02:00			`4 #hosorder# %for HOS-picker, order of Higher Order Statistics`
			`2 #Parorder# %for AR-picker, order of AR process of Z-component`
[update] default pylot.in file updates 2017-09-21 15:19:15 +02:00			`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`
[add] latest pylot.in file added to inputs! 2017-07-10 15:49:15 +02:00			`0.001 #addnoise# %add noise to seismogram for stable AR prediction`
[update] default pylot.in file updates 2017-09-21 15:19:15 +02:00			`3.0 0.1 0.5 1.0 #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]`
[add] latest pylot.in file added to inputs! 2017-07-10 15:49:15 +02:00			`0.001 #ausP# %for HOS/AR, artificial uplift of samples (aus) of CF (P)`
[update] default pylot.in file updates 2017-09-21 15:19:15 +02:00			`1.3 #nfacP# %for HOS/AR, noise factor for noise level determination (P)`
[new] pylot.in which will substitute other input files later a more general input file might help prevent ambiguity among use cases 2016-09-08 13:45:16 +02:00			`#H-components#`
[add] latest pylot.in file added to inputs! 2017-07-10 15:49:15 +02:00			`ARH #algoS# %choose algorithm for S-onset determination (ARH or AR3)`
[update] default pylot.in file updates 2017-09-21 15:19:15 +02:00			`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`
[add] latest pylot.in file added to inputs! 2017-07-10 15:49:15 +02:00			`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)`
[update] default pylot.in file updates 2017-09-21 15:19:15 +02:00			`4.0 #pickwinS# %for initial AIC pick, length of S-pick window [s]`
			`2.0 0.3 1.5 1.0 #tsnrh# %for ARH/AR3, window lengths for SNR-and slope estimation [tnoise, tsafetey, tsignal, tslope] [s]`
			`1.0 #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)`
[add] latest pylot.in file added to inputs! 2017-07-10 15:49:15 +02:00			`0.9 #ausS# %for HOS/AR, artificial uplift of samples (aus) of CF (S)`
[update] default pylot.in file updates 2017-09-21 15:19:15 +02:00			`1.5 #nfacS# %for AR-picker, noise factor for noise level determination (S)`
[add] latest pylot.in file added to inputs! 2017-07-10 15:49:15 +02:00			`#first-motion picker#`
			`1 #minfmweight# %minimum required P weight for first-motion determination`
			`2.0 #minFMSNR# %miniumum required SNR for first-motion determination`
			`0.2 #fmpickwin# %pick window around P onset for calculating zero crossings`
			`#quality assessment#`
[update] default pylot.in file updates 2017-09-21 15:19:15 +02:00			`0.02 0.04 0.08 0.16 #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`
			`0.8 #minAICPslope# %below this slope [counts/s] the initial P pick is rejected`
[add] latest pylot.in file added to inputs! 2017-07-10 15:49:15 +02:00			`1.1 #minAICPSNR# %below this SNR the initial P pick is rejected`
			`1.0 #minAICSslope# %below this slope [counts/s] the initial S pick is rejected`
			`1.5 #minAICSSNR# %below this SNR the initial S pick is rejected`
[update] default pylot.in file updates 2017-09-21 15:19:15 +02:00			`1.0 #minsiglength# %length of signal part for which amplitudes must exceed noiselevel [s]`
[add] latest pylot.in file added to inputs! 2017-07-10 15:49:15 +02:00			`1.0 #noisefactor# %noiselevel*noisefactor=threshold`
			`10.0 #minpercent# %required percentage of amplitudes exceeding threshold`
			`1.5 #zfac# %P-amplitude must exceed at least zfac times RMS-S amplitude`
			`6.0 #mdttolerance# %maximum allowed deviation of P picks from median [s]`
			`1.0 #wdttolerance# %maximum allowed deviation from Wadati-diagram`
[update] default pylot.in file updates 2017-09-21 15:19:15 +02:00			`5.0 #jackfactor# %pick is removed if the variance of the subgroup with the pick removed is larger than the mean variance of all subgroups times safety factor`