function [index] = searchElementRegion(lat, lon, lat_limit, lon_limit)
lat_ind = find(lat >= lat_limit(1) & lat <= lat_limit(2));
lon_ind = find(lon >= lon_limit(1) & lon <= lon_limit(2));
index = intersect(lat_ind, lon_ind);
