added weighting option for robustness_function. use robustness_function at every grid point with apply_ufunc and vectorize=True. tried using dask and numba, issue with R being dask array, something wrong
