<div dir="ltr"><div dir="ltr">Dear NCL users </div><div>I use the following simple code to create map (C) in the attached image. As you can see, I am trying to change the land use index inside the mask from 17 to 11, which I succeeded in achieving. However, due to masking I lose the rest of the data.  </div><div dir="ltr"><div><br></div><div><div>    a = addfile("<a href="http://geo_em.nc">geo_em.nc</a>","w")<br></div><div>    lu=wrf_user_getvar(a,"LU_INDEX",-1)<br></div><div>    f   = addfile("<a href="http://wrf_mask.nc">wrf_mask.nc</a>","r")<br></div><div>     marray   = f->WRF_mask</div><div>     lu_mask = mask(lu,marray,1)</div><div>      a->LU_INDEX=where(lu_mask .eq. 17,11,a->LU_INDEX)</div></div><div><br></div><div>To correct this problem I do step D and E in the attached image but I still am unable to achieve the required result. Can someone tell me how to correct this issue ?</div><div><br></div><div><br></div></div></div>