<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none"><!-- P { margin-top: 0px; margin-bottom: 0px; }--></style>
</head>
<body dir="ltr" style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">
<div>Dear ncl-Team</div>
<div> I’m using the coastline of NOAA (GSHHS_h_caspSea.shp) (I've upload it in the ftp) and I would like to generate the meshmask from the costline and to mask the bathymetry datasets outside the coastline (landPoint) in the region of Caspian Sea.</div>
<div> I found the ncl Built-in functions '<span style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 16px; background-color: rgb(255, 255, 255);">create_mask_from_shapefile</span>' and I tried to use it in this way:</div>
<div> </div>
<div>//=========================</div>
<div> mask2dtCoast_surf(0:ny(inest)-1,0:nx(inest)-1) = \</div>
<div> create_mask_from_shapefile(bathy(0:ny(inest)-1,0:nx(inest)-1),pathfile_regriddataCoastline)</div>
<div><br>
</div>
<div> ;--- Mashera bathy(dentro o fuori la costline)<br>
</div>
<div> do iy=0,ny(inest)-1</div>
<div> do ix=0,nx(inest)-1</div>
<div> ;if(mask2dtCoast_surf(iy,ix) .eq. 0) then</div>
<div> ; bathy(iy,ix) = 0.</div>
<div> ;end if</div>
<div> if(mask2dtCoast_surf(iy,ix) .eq. 1) then</div>
<div> bathy(iy,ix) = 0.</div>
<div> end if</div>
<div> end do</div>
<div> end do</div>
<div> ;print(bathy)</div>
<div>//============================</div>
<div><br>
</div>
<div>but the bathymetry inside the caspian sea, where there are isle, are not masked (attachment).</div>
<div>Why? Is this function a good choice in order to work with shapfile?</div>
<div><br>
</div>
<div>thanks</div>
<div>Francesco<br>
</div>
<div><br>
</div>
<p><br>
</p>
</body>
</html>