# NOT RUN {
BTemp <- raster::raster(ncol=100, nrow=100)
set.seed(2)
raster::values(BTemp) = runif(10000, min=298, max=305)
E <- raster::raster(ncol=100, nrow=100)
set.seed(2)
raster::values(E) = runif(10000, min=0.96, max=0.99)
MWA(BT = BTemp, tau = 0.86, E = E, Ta = 26)
# }
Run the code above in your browser using DataLab