<div dir="ltr"><div><div><div><div>Hi Xiao,<br><br>Its complaining about the reassignment operator := I don't recall exactly when it was introduced in NCL, but it was relatively recently, and I strongly suspect it was after ver. 6.1.0. You can work around it by adding these lines at the end of the inner do-loop in that section of code3:<br><br></div> delete(lat_sub)<br></div> delete(lon_sub)<br><br></div>Hope that helps...<br></div>Rick<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jul 6, 2016 at 12:43 PM, Zhang, Xiao <span dir="ltr"><<a href="mailto:xiao.zhang@pnnl.gov" target="_blank">xiao.zhang@pnnl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div link="blue" vlink="purple" lang="EN-US">
<div>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Hi <u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">I use ncl 6.1.0 and would like to use function shapefile_mask_data. Since it has not been built in, I use the code provided from the website but got error message. Could you please help with this? Thank you!
<u></u><u></u></p>
<p class="MsoNormal">“fatal:syntax error: line 158 in file ./shapefile_mask_data.ncl before or near :<u></u><u></u></p>
<p class="MsoNormal"> lat_sub :<u></u><u></u></p>
<p class="MsoNormal">------------------^<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">fatal:error in statement<u></u><u></u></p>
<p class="MsoNormal">fatal:syntax error: line 159 in file ./shapefile_mask_data.ncl before or near :<u></u><u></u></p>
<p class="MsoNormal"> lon_sub :<u></u><u></u></p>
<p class="MsoNormal">------------------^<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">fatal:error in statement<u></u><u></u></p>
<p class="MsoNormal">fatal:Syntax Error in block, block not executed”<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Here are the codes I used.<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><a href="https://www.ncl.ucar.edu/Support/talk_archives/2014/att-0432/shapefile_mask_data.ncl" target="_blank">https://www.ncl.ucar.edu/Support/talk_archives/2014/att-0432/shapefile_mask_data.ncl</a><u></u><u></u></p>
<p class="MsoNormal">part of my code <u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> shpfile = "USA_adm_shp/USA_adm0.shp"<u></u><u></u></p>
<p class="MsoNormal"> opt = True<u></u><u></u></p>
<p class="MsoNormal"> opt@return_mask = True<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal"> ;load self-define function<u></u><u></u></p>
<p class="MsoNormal"> load "./shapefile_mask_data.ncl"<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal"> land_mask_a1 = shapefile_mask_data(dif_a1(0,:,:),shpfile,opt)<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Best,<u></u><u></u></p>
<p class="MsoNormal">Xiao<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
<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>
<br></blockquote></div><br></div>