def test_get_surround_elem_uids():
'''
CubeGridRemap.get_surround_elem_uids(): ne=30
'''
from cube_remap import CubeGridRemap
from pkg.convert_coord.cs_ll import abp2latlon
ne, ngq = 30, 4
rotated = False
cube = CubeGridRemap(ne, ngq, rotated)
td = (np.pi/2)/ne/3/2 # tiny delta
ij = (1,2,1,1,1)
gid = cube.ij2gid[ij]
alpha, beta = cube.alpha_betas[gid]
lat, lon = abp2latlon(alpha+td, beta+td, ij[0])
ret_uids = cube.get_surround_elem_uids(lat, lon)
a_equal(ret_uids, [3,16,17,18,7,19,20,21,11,22,23,24,15,25,26,27])
