Re: WRF_lc_5.ncl contributed script gives incorrect masking result

Hi Kamal,

I believe the issue is that the data for this particular example had missing values in it. Your data does not, so that's why there's no masking of land.

You can read the LANDMASK variable and your file and apply it to sst before plotting:

  nt = 0
  sst = f->SST(nt,:,:) ; (Time, south_north, west_east )
  landmask = f->LANDMASK(nt,:,:)
  sst = mask(sst,landmask,1)

No need to set sst@_FillValue. It will be set by the "mask" function.

--Mary

On Apr 2, 2013, at 7:38 PM, mmkamal@uwaterloo.ca wrote:

> Hi,
>
> I applied the following contributed script (http://www.ncl.ucar.edu/Applications/Scripts/WRF_lc_5.ncl) to my WRF output. But It provides result with no masking
>
> I have changed nothing except the data path. I am attaching the data and plot for your convenience.
>
> https://www.dropbox.com/s/2t7uyyk3jm0bhb9/new_wrfout_d01_2002-07-01.rar
>
>
> Thanks
> Kamal
>
>
>
>
>
> <WRF_lc.png>_______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Apr 3 09:16:05 2013