search for: masksim

...10) #Mask: 1=DATA; 10=NODATA vectMask <- c(10,1,10,1,1,10,1,10,1,10,10,1) length(vectSim) length(vectMask) numRow<-3 numCol<-4 numReal<-3 Sim <- array(vectSim, c(numRow,numCol,numReal)) Sim Mask <- array(vectMask, c(numRow,numCol)) Mask SmoothSim <- apply(Sim, c(1,2), ifelse(MaskSim==10,33,99)) SmoothSim ******************** Any help is much appreciated! Thanks, Sander and Alessandro.