def create_random_src(min_distance=10, nr_sources=100):
    """Create nr_sources random sources in RA and DEC (in degree)
    :param min_distance : create sources with distance > than min_distance
    among each other (in degree)
    :param nr_sources : number of random sources to create
    """
    import astropy.coordinates as coord
    import astropy.units as u
    from astropy.coordinates import SkyCoord
    import numpy as np
    min_distance = min_distance * u.deg
    print('Minimum distance between sources is: ', min_distance,
        '\nNumber of random sources is: ', nr_sources)
    sources_ra = []
    sources_dec = []
    ra_first = np.random.uniform(-180, 180, 1)
    dec_first = np.random.uniform(0, 90, 1)
    sources_ra.append(ra_first)
    sources_dec.append(dec_first)
    while len(sources_dec) < nr_sources:
        ra = np.random.uniform(-180, 180, 1)
        dec = np.random.uniform(0, 90, 1)
        c1 = SkyCoord(ra=ra, dec=dec, frame='icrs', unit='deg')
        c2 = SkyCoord(ra=sources_ra, dec=sources_dec, frame='icrs', unit='deg')
        sep_vect = c1.separation(c2)
        greater_ = [(True) for i in sep_vect if i > min_distance]
        if len(greater_) < len(sources_ra):
            continue
        else:
            sources_ra.append(ra)
            sources_dec.append(dec)
    final_ra = []
    final_dec = []
    for RA, DEC in zip(sources_ra, sources_dec):
        final_ra.append(RA[0])
        final_dec.append(DEC[0])
    return final_ra, final_dec
