<div dir="ltr"><div><div>HI,<br><br></div>In what sense are the data not consistent? Do you have a plot that demonstrates the issue?<br><br></div>Rick<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, May 27, 2016 at 4:06 AM, franz <span dir="ltr"><<a href="mailto:f.trotta@sincem.unibo.it" target="_blank">f.trotta@sincem.unibo.it</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear ncl-Team<br>
<br>
I'm using NCL - version 6.3.0<br>
I would like to create a land-sea mask using a shapefile.<br>
I’ve created a shapefile from qgis and I tried to use the function<br>
create_mask_from_shapefile<br>
<br>
mask2dtCoast = new((/ny,nx/),float)<br>
mask2dtCoast(0:ny-1,0:nx-1) =<br>
create_mask_from_shapefile(lon2d_surf(0:ny-1,0:nx-1),pathfile_Coastline)<br>
<br>
but I've got a not expect result:<br>
sea and land grid points are not consistente with the coastline<br>
contained in the shape file in same region of my domain :-(<br>
<br>
I also try to use the function shapefile_mask_data in the following<br>
way:<br>
<br>
opt = True<br>
opt@return_mask = True<br>
mask2dtCoast_surf(0:ny-1,0:nx-1) =<br>
shapefile_mask_data(lon2d_surf,pathfile_Coastline,opt)<br>
<br>
but with the same result.<br>
<br>
where is the error?<br>
<br>
thanks<br>
Francesco T.<br>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
</blockquote></div><br></div>