[ncl-talk] land/sea mask from shapefile

[Rick Brownrigg]

HI,

In what sense are the data not consistent? Do you have a plot that
demonstrates the issue?

Rick

On Fri, May 27, 2016 at 4:06 AM, franz <f.trotta at sincem.unibo.it> wrote:

> Dear ncl-Team
>
>    I'm using NCL - version 6.3.0
>    I would like to create a land-sea mask using a shapefile.
>    I’ve created a shapefile from qgis and I tried to use the function
> create_mask_from_shapefile
>
>    mask2dtCoast = new((/ny,nx/),float)
>    mask2dtCoast(0:ny-1,0:nx-1) =
> create_mask_from_shapefile(lon2d_surf(0:ny-1,0:nx-1),pathfile_Coastline)
>
>    but I've got a not expect  result:
>    sea and land grid points are not consistente with the coastline
> contained in the shape file in same region of my domain :-(
>
>    I also try to use the function shapefile_mask_data in the following
> way:
>
>    opt = True
>    opt at return_mask = True
>    mask2dtCoast_surf(0:ny-1,0:nx-1) =
> shapefile_mask_data(lon2d_surf,pathfile_Coastline,opt)
>
>    but with the same result.
>
>    where is the error?
>
>    thanks
> Francesco T.
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160527/d4739116/attachment.html