diff --git a/inputs/pylot.in b/inputs/pylot_global.in similarity index 79% rename from inputs/pylot.in rename to inputs/pylot_global.in index e3b811d1..6024b7d8 100644 --- a/inputs/pylot.in +++ b/inputs/pylot_global.in @@ -4,19 +4,19 @@ %Parameters are optimized for %extent data sets! %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% #main settings# -/home/marcel/marcel_scratch #rootpath# %project path -alparray #datapath# %data path -waveforms #database# %name of data base -e0006.036.13 #eventID# %event ID for single event processing (* for all events found in database) -None #invdir# %full path to inventory or dataless-seed file + #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# -/progs/bin #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 -Insheim_min1d032016_auto.in #ctrfile# %name of autoPyLoT-output control file for NLLoc +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 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% @@ -37,62 +37,64 @@ bandpass bandpass #filter_type# %filter type %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% #common settings picker# global #extent# %extent of array ("local", "regional" or "global") -50.0 #pstart# %start time [s] for calculating CF for P-picking -600.0 #pstop# %end time [s] for calculating CF for P-picking +-150.0 #pstart# %start time [s] for calculating CF for P-picking (if TauPy: seconds relative to estimated onset) +600.0 #pstop# %end time [s] for calculating CF for P-picking (if TauPy: seconds relative to estimated onset) 200.0 #sstart# %start time [s] relative to P-onset for calculating CF for S-picking 1150.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 +True #use_taup# %use estimated traveltimes from TauPy for calculating windows for CF +iasp91 #taup_model# %define TauPy model for traveltime estimation. Possible values: 1066a, 1066b, ak135, ak135f, herrin, iasp91, jb, prem, pwdk, sp6 0.05 0.5 #bpz1# %lower/upper corner freq. of first band pass filter Z-comp. [Hz] -0.01 0.5 #bpz2# %lower/upper corner freq. of second band pass filter Z-comp. [Hz] +0.001 0.5 #bpz2# %lower/upper corner freq. of second band pass filter Z-comp. [Hz] 0.05 0.5 #bph1# %lower/upper corner freq. of first band pass filter H-comp. [Hz] -0.01 0.5 #bph2# %lower/upper corner freq. of second band pass filter z-comp. [Hz] +0.001 0.5 #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) -15.0 #tlta# %for HOS-/AR-AIC-picker, length of LTA window [s] +150.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 -6.0 #tdet1z# %for AR-picker, length of AR determination window [s] for Z-component, 1st pick -2.0 #tpred1z# %for AR-picker, length of AR prediction window [s] for Z-component, 1st pick -3.0 #tdet2z# %for AR-picker, length of AR determination window [s] for Z-component, 2nd pick -1.0 #tpred2z# %for AR-picker, length of AR prediction window [s] for Z-component, 2nd pick +16.0 #tdet1z# %for AR-picker, length of AR determination window [s] for Z-component, 1st pick +10.0 #tpred1z# %for AR-picker, length of AR prediction window [s] for Z-component, 1st pick +12.0 #tdet2z# %for AR-picker, length of AR determination window [s] for Z-component, 2nd pick +6.0 #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 -60.0 10.0 150.0 3.0 #tsnrz# %for HOS/AR, window lengths for SNR-and slope estimation [tnoise, tsafetey, tsignal, tslope] [s] -10.0 #pickwinP# %for initial AIC pick, length of P-pick window [s] -20.0 #Precalcwin# %for HOS/AR, window length [s] for recalculation of CF (relative to 1st pick) +60.0 10.0 40.0 10.0 #tsnrz# %for HOS/AR, window lengths for SNR-and slope estimation [tnoise, tsafetey, tsignal, tslope] [s] +150.0 #pickwinP# %for initial AIC pick, length of P-pick window [s] +35.0 #Precalcwin# %for HOS/AR, window length [s] for recalculation of CF (relative to 1st pick) 6.0 #aictsmooth# %for HOS/AR, take average of samples for smoothing of AIC-function [s] 4.0 #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.1 #nfacP# %for HOS/AR, noise factor for noise level determination (P) #H-components# ARH #algoS# %choose algorithm for S-onset determination (ARH or AR3) -6.0 #tdet1h# %for HOS/AR, length of AR-determination window [s], H-components, 1st pick -4.0 #tpred1h# %for HOS/AR, length of AR-prediction window [s], H-components, 1st pick -6.0 #tdet2h# %for HOS/AR, length of AR-determinaton window [s], H-components, 2nd pick -3.0 #tpred2h# %for HOS/AR, length of AR-prediction window [s], H-components, 2nd pick +12.0 #tdet1h# %for HOS/AR, length of AR-determination window [s], H-components, 1st pick +6.0 #tpred1h# %for HOS/AR, length of AR-prediction window [s], H-components, 1st pick +8.0 #tdet2h# %for HOS/AR, length of AR-determinaton window [s], H-components, 2nd pick +4.0 #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) -15.0 #pickwinS# %for initial AIC pick, length of S-pick window [s] -100.0 10.0 40.0 6.0 #tsnrh# %for ARH/AR3, window lengths for SNR-and slope estimation [tnoise, tsafetey, tsignal, tslope] [s] -2.0 #aictsmoothS# %for AIC-picker, take average of samples for smoothing of AIC-function [s] -3.0 #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) +30.0 #Srecalcwin# %for AR-picker, window length [s] for recalculation of CF (2nd pick) (H) +195.0 #pickwinS# %for initial AIC pick, length of S-pick window [s] +100.0 10.0 45.0 10.0 #tsnrh# %for ARH/AR3, window lengths for SNR-and slope estimation [tnoise, tsafetey, tsignal, tslope] [s] +22.0 #aictsmoothS# %for AIC-picker, take average of samples for smoothing of AIC-function [s] +10.0 #tsmoothS# %for AR-picker, take average of samples for smoothing CF [s] (S) +0.001 #ausS# %for HOS/AR, artificial uplift of samples (aus) of CF (S) 1.2 #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 +3.0 #minFMSNR# %miniumum required SNR for first-motion determination +10.0 #fmpickwin# %pick window around P onset for calculating zero crossings #quality assessment# 1.0 2.0 4.0 8.0 #timeerrorsP# %discrete time errors [s] corresponding to picking weights [0 1 2 3] for P 4.0 8.0 16.0 32.0 #timeerrorsS# %discrete time errors [s] corresponding to picking weights [0 1 2 3] for S 0.5 #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.3 #minAICSSNR# %below this SNR the initial S pick is rejected 5.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 +1.2 #zfac# %P-amplitude must exceed at least zfac times RMS-S amplitude +25.0 #mdttolerance# %maximum allowed deviation of P picks from median [s] +50.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 diff --git a/inputs/pylot_local.in b/inputs/pylot_local.in new file mode 100644 index 00000000..fe019440 --- /dev/null +++ b/inputs/pylot_local.in @@ -0,0 +1,100 @@ +%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 diff --git a/inputs/pylot_regional.in b/inputs/pylot_regional.in new file mode 100644 index 00000000..fe019440 --- /dev/null +++ b/inputs/pylot_regional.in @@ -0,0 +1,100 @@ +%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 diff --git a/makePyLoT.py b/makePyLoT.py index cb74a36b..e50d17db 100644 --- a/makePyLoT.py +++ b/makePyLoT.py @@ -158,17 +158,23 @@ def buildPyLoT(verbosity=None): def installPyLoT(verbosity=None): - files_to_copy = {'pylot.in': ['~', '.pylot'], + files_to_copy = {'pylot_local.in': ['~', '.pylot'], + 'pylot_regional.in': ['~', '.pylot'], 'pylot_global.in': ['~', '.pylot']} if verbosity > 0: print('starting installation of PyLoT ...') if verbosity > 1: print('copying input files into destination folder ...') ans = input('please specify scope of interest ' - '([0]=local, 1=regional) :') or 0 + '([0]=local, 1=regional, 2=global) :') or 0 if not isinstance(ans, int): ans = int(ans) - ans = 'local' if ans is 0 else 'regional' + if ans == 0: + ans = 'local' + elif ans == 1: + ans = 'regional' + elif ans == 2: + ans = 'global' link_dest = [] for file, destination in files_to_copy.items(): link_file = ans in file