Add function and wrapper script to subset global data into regional sets
bugfix