Huge speedup using pointers for map_local in subset.f90. subset_to_grid_double used to take 40% of processing time, now drops to insignificant number.