def gaia_dr2_conesearch_query(ra=165.86, dec=34.829694, radius=3.0, max=100000
    ):
    """
    Generate a query string for the TAP servers
    TBD
    
    Parameters
    ----------
    ra, dec : float
        RA, Dec in decimal degrees

    radius : float
        Search radius, in arc-minutes.
    
    Returns
    -------
    query : str
        Query string
        
    """
    query = (
        "SELECT TOP {3} * FROM gaiadr2.gaia_source  WHERE CONTAINS(POINT('ICRS',gaiadr2.gaia_source.ra,gaiadr2.gaia_source.dec),CIRCLE('ICRS',{0},{1},{2:.2f}))=1"
        .format(ra, dec, radius / 60.0, max))
    return query
