<div dir="ltr"><div><div><div><div><div>Suggestion: <br><br><a href="http://www.ncl.ucar.edu/Applications/mask.shtml">http://www.ncl.ucar.edu/Applications/mask.shtml</a><br><br></div>See examples 9 and 13<br><br></div>In particular, once tou create the appropriate mask, it should be saved to a netCDF file.<br></div>Subsequently read the mask from the file.<br><br></div>Good Luck<br></div><div>D<br></div><div><div><div><div><div><br><br></div></div></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Nov 17, 2014 at 7:16 AM, Ufuk Utku Turuncoglu (BE) <span dir="ltr">&lt;<a href="mailto:u.utku.turuncoglu@be.itu.edu.tr" target="_blank">u.utku.turuncoglu@be.itu.edu.tr</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
<br>
You have to define a polygon around the region that you interested and<br>
you could use gc_inout function to assign different number (oe mask<br>
value) for different regions. You could create a mask for every<br>
sub-regions. So, for example<br>
<br>
       do r = 0, nregs-1<br>
         do i = 0, nlon-1<br>
           do j = 0, nlat-1<br>
             if (gc_inout(lat2d(j,i), lon2d(j,i), ybox(r,:), xbox(r,:)))<br>
then<br>
               dmask(k,r,j,i) = r<br>
             end if<br>
           end do<br>
         end do<br>
       end do<br>
<br>
in this case, nregs is the number of sub-region and xbox, ybox store the<br>
polygon information. You could also add land-sea mask to the if clause<br>
to eliminate the land grid points. Of course, this method works for the<br>
data that has a grid bu it could be extended to other types of data. I<br>
hope it helps.<br>
<br>
Regards,<br>
<br>
--ufuk<br>
_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
</blockquote></div><br></div>