[ncl-talk] land/sea mask from shapefile
franz
f.trotta at sincem.unibo.it
Fri May 27 04:06:33 MDT 2016
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.
More information about the ncl-talk
mailing list