Dear Sir/Madam, I like to lay on geographical region on my the image plot of my study like 5S to 25N and 20W to 60E. How do I do? Here is the script I tried but not wrok for me. Thank you for your help. library(ncdf) library(grid) library(maps) library(fields) nc=open.ncdf('wrfout_mon') T2=get.var.ncdf(nc,'T2') lat=get.var.ncdf(nc,'x_2') lon=get.var.ncdf(nc,'y_2') time=get.var.ncdf(nc,'time') # dim(T2)=c(147,121,8,125) # T2.daily=apply(T2, c(1,2,4),mean) T2.clim=apply(T2,c(1,2),mean) image.plot(lat,lon,T2.clim) map(interior=TRUE,add=TRUE) Bedassa -- --------------------------------------------------- *Bedassa Regassa Cheneka* *Adama Science and Technology Univesity* *Natural Resource Management* *P.O.BOX---193* *Assela * ----------------------------------------------------- [[alternative HTML version deleted]]