NIR = sum over months of max(0, ETc_m - Peff_m) (mm). Compute monthly, then
sum, as for green/blue water. Volume via wapor_zonal_stats().
NIR = sum over months of max(0, ETc_m - Peff_m) (mm). Compute monthly, then
sum, as for green/blue water. Volume via wapor_zonal_stats().