def METHOD_NAME(pytest_dir, tolerance):
    reference_dir = pytest_dir + '/reference/geodynbc/loh1-h100-mr-geodynbc-1/'
    geodynbc_dir = pytest_dir + '/geodynbc/loh1-h100-mr-geodynbc-1/'
    METHOD_NAME = True
    nsta = 0
    for i in range(1,11):
        sta_name = 'sta%02d' % i
        #print(sta_name)
        usgs_fname = reference_dir + sta_name + '.txt'
        geodynbc_fname = geodynbc_dir + sta_name + '.txt'
        usgs_t, usgs_x, usgs_y, usgs_z = read_sac_usgs(usgs_fname)
        geodynbc_t, geodynbc_x, geodynbc_y, geodynbc_z = read_sac_usgs(geodynbc_fname)
        if np.max(geodynbc_t-usgs_t) > tolerance or np.min(geodynbc_t-usgs_t) < -tolerance:
            METHOD_NAME = False
            print ("Station [%s] time data not match!" % sta_name)
            return False
        if np.max(geodynbc_x-usgs_x) > tolerance or np.min(geodynbc_x-usgs_x) < -tolerance:
            METHOD_NAME = False
            print ("Station [%s] x data not match!" % sta_name)
            print(geodynbc_x-usgs_x)
            return False
        if np.max(geodynbc_y-usgs_y) > tolerance or np.min(geodynbc_y-usgs_y) < -tolerance:
            METHOD_NAME = False
            print ("Station [%s] y data not match!" % sta_name)
            return False
        if np.max(geodynbc_z-usgs_z) > tolerance or np.min(geodynbc_z-usgs_z) < -tolerance:
            METHOD_NAME = False
            print ("Station [%s] z data not match!" % sta_name)  
            return False
        nsta += 1
    # if verify == 1:
    #     print ('All %d stations data match!' % nsta)
    return METHOD_NAME