[ncl-talk] Masked data

[ali mughal]

Dear NCL users
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.

    a = addfile("geo_em.nc","w")
    lu=wrf_user_getvar(a,"LU_INDEX",-1)
    f   = addfile("wrf_mask.nc","r")
     marray   = f->WRF_mask
     lu_mask = mask(lu,marray,1)
      a->LU_INDEX=where(lu_mask .eq. 17,11,a->LU_INDEX)

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 ?
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20211129/7295f709/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Capture2.PNG
Type: image/png
Size: 441191 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20211129/7295f709/attachment-0001.png>