[ncl-talk] error in shapefile_mask_data.ncl

[Rick Brownrigg]

Hi Xiao,

Where/when did you get the file shapefile_mask_data.ncl?   Things may have
been refactored;  the shapefile examples webpage makes reference to
function shapefile_mask_data in a utilities script named
shapefile_utils.ncl.  The version of the function in that script requires
only 3 variables.  It looks as though the svarname and mask_names arguments
in the 5-arg version have been refactored to be parameters on the "opt"
argument (make sense?)   Perhaps you might want to try with
shapefile_util.ncl?

http://ncl.ucar.edu/Applications/shapefiles.shtml

Hope that helps...
Rick

On Wed, Jul 6, 2016 at 1:56 PM, Zhang, Xiao <xiao.zhang at pnnl.gov> wrote:

> Many thanks for pointing this out! I added the delete function and the old
> error message disappeared! But it gave me this
>
>
>
> fatal:syntax error: function shapefile_mask_data expects 5 arguments, got 3
>
> fatal:error at line 623 in file fourier_analysis_HL.ncl
>
>
>
> In the shapefile_mask_data.ncl, it only lists three arguments though, any
> idea where went wrong?
>
>
>
> undef("shapefile_mask_data")
>
> function shapefile_mask_data(data[*][*]:numeric,sfilename[1]:string,\
>
>                              svarname[1]:string,mask_names[*]:string,\
>
>                              keep_opt[1]:logical)
>
>
>
> *From:* Rick Brownrigg [mailto:brownrig at ucar.edu]
> *Sent:* Wednesday, July 06, 2016 12:37 PM
> *To:* Zhang, Xiao
> *Cc:* ncl-talk at ucar.edu
> *Subject:* Re: [ncl-talk] error in shapefile_mask_data.ncl
>
>
>
> Hi Xiao,
>
> 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:
>
>   delete(lat_sub)
>
>   delete(lon_sub)
>
> Hope that helps...
>
> Rick
>
>
>
> On Wed, Jul 6, 2016 at 12:43 PM, Zhang, Xiao <xiao.zhang at pnnl.gov> wrote:
>
>
>
> Hi
>
>
>
> 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!
>
> “fatal:syntax error: line 158 in file ./shapefile_mask_data.ncl before or
> near :
>
>           lat_sub :
>
> ------------------^
>
>
>
> fatal:error in statement
>
> fatal:syntax error: line 159 in file ./shapefile_mask_data.ncl before or
> near :
>
>           lon_sub :
>
> ------------------^
>
>
>
> fatal:error in statement
>
> fatal:Syntax Error in block, block not executed”
>
>
>
> Here are the codes I used.
>
>
>
>
> https://www.ncl.ucar.edu/Support/talk_archives/2014/att-0432/shapefile_mask_data.ncl
>
> part of my code
>
>
>
>                 shpfile   = "USA_adm_shp/USA_adm0.shp"
>
>                 opt             = True
>
>                 opt at return_mask = True
>
>
>
>                  ;load self-define function
>
>                 load "./shapefile_mask_data.ncl"
>
>
>
>                  land_mask_a1 =
> shapefile_mask_data(dif_a1(0,:,:),shpfile,opt)
>
>
>
> Best,
>
> Xiao
>
>
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160706/22b970b0/attachment.html