This routine outputs both the industrial areas & the shops within.
@aTarom’s routine searches using only shop nodes & is_in can’t use ways & rels. This is resolved by finding all the nodes with the way/rels shops. using: (._;>;);
There’s no need to map_to_area as landuses are now recognised automatically as areas.
There’s no need to give specific names to all selection sets. ._ is the default name for the previous selection set.
wr[landuse=industrial]({{bbox}});
nwr[shop](area);
out geom qt;
(._;>;);
is_in;
wr._[landuse=industrial];
out geom qt;