from glob import glob
from shutil import copy

data = glob('data*')
datax = [x.split('_') for x in data]

xlist = []
for i,ilat,ilon in datax:
    lat = float(ilat)
    lon = float(ilon)
    if ((lat >= 46.5) & (lat <= 47.5) & (lon <= -121.75) & (lon >=-123.25)):
        xlist.append(f'{i}_{ilat}_{ilon}')
    #46.5-47.5N, 121.75-123.25W

for x in xlist:
    print(x)
    copy(x, 'thurston')