[new] add extent option "active" to pick high frequency active seismic data

makePyLoT.py

@@ -166,7 +166,7 @@ def installPyLoT(verbosity=None):
     if verbosity > 1:
         print('copying input files into destination folder ...')
     ans = input('please specify scope of interest '
-                '([0]=local, 1=regional, 2=global) :') or 0
+                '([0]=local, 1=regional, 2=global, 3=active) :') or 0
     if not isinstance(ans, int):
         ans = int(ans)
     if ans == 0:
@@ -175,6 +175,8 @@ def installPyLoT(verbosity=None):
         ans = 'regional'
     elif ans == 2:
         ans = 'global'
+    elif ans == 3:
+        ans = 'active'
     link_dest = []
     for file, destination in files_to_copy.items():
         link_file = ans in file

pylot/core/io/default_parameters.py

@@ -27,7 +27,7 @@ defaults = {'rootpath': {'type': str,
                         'namestring': 'Event ID'},
 
             'extent': {'type': str,
-                       'tooltip': 'extent of array ("local", "regional" or "global")',
+                       'tooltip': 'extent of array ("active", "local", "regional" or "global")',
                        'value': 'local',
                        'namestring': 'Array extent'},
 

pylot/core/pick/utils.py

@@ -548,13 +548,14 @@ def getResolutionWindow(snr, extent):
     see also Diehl et al. 2009
     :param snr: Signal to noise ration which decides the witdth of the resolution window
     :type snr: float
-    :param extent: can be 'local', 'regional', 'global'
+    :param extent: can be 'active', 'local', 'regional', 'global'
     :type extent: str
     :return: half width of the resolution window
     :rtype: float
     """
 
     res_wins = {
+        'active': {'HRW': .02, 'MRW': .05, 'LRW': .1, 'VLRW': .15},
         'regional': {'HRW': 2., 'MRW': 5., 'LRW': 10., 'VLRW': 15.},
         'local': {'HRW': 2., 'MRW': 5., 'LRW': 10., 'VLRW': 15.},
         'global': {'HRW': 40., 'MRW': 100., 'LRW': 200., 'VLRW': 300.}