[ncl-talk] about shapefiles

[Rick Brownrigg]

Hi,

That certainly looks correct, and I see you are running 6.3.0.  What do you
get if you run:

   grep shapefile_polyline $NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl

or just

   grep shapefile $NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl



Rick



On Sun, Jun 26, 2016 at 7:08 AM, Agud Pique, Enric <eagud.q at ccma.cat> wrote:

> Hello,
>
>
>
> I have a problem with the following code. I can plot the variables Tc and
> the wind, but  I get an error from the shapefile’s command. The error is...
>
>
>
> [Enric at meteo4u em_real]$ ncl plot-shapefiles.ncl
>  Copyright (C) 1995-2015 - All Rights Reserved
>  University Corporation for Atmospheric Research
>  NCAR Command Language Version 6.3.0
>  The use of this software is governed by a License Agreement.
>  See http://www.ncl.ucar.edu/ for more details.
> fatal:Undefined identifier: (gsn_add_shapefile_polylines) is undefined,
> can't continue
> fatal:["Execute.c":8575]:Execute: Error occurred at or near line 49 in
> file plot-shapefiles.ncl
>
>
> Any idea?
>
>
> Best
>
>
>
>
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>
>
>
> begin
>
> a = addfile("./wrfout_d01_2016-06-26_18:00:00","r")
>
> wks = gsn_open_wks("png","shapefiles")
>
>
>
>
>
> T2 = wrf_user_getvar(a,"T2",0)
>
> Tc = T2-273
>
> u10 = wrf_user_getvar(a,"U10",0)    ; u at 10 m, mass point
>
> v10 = wrf_user_getvar(a,"V10",0)    ; v at 10 m, mass point
>
> u10 = u10*1.94386                    ; Turn wind into knots
>
> v10 = v10*1.94386
>
> u10 at units = "kts"
>
> v10 at units = "kts"
>
>
>
> res = True
>
> pltres = True
>
> mpres = True
>
> ;opts = True
>
>
>
> opts = res
>
> opts at cnFillOn = True
>
> opts at cnLevelSelectionMode = "ExplicitLevels"
>
> opts at cnLevels   =(/ -10, -5,0,5,10,15,20,25,30,35,40/)
>
> opts at cnFillColors =(/"Blue","DarkOliveGreen1", \
>
>                       "DarkOliveGreen3","Chartreuse", \
>
>                       "Chartreuse3","Green","ForestGreen", \
>
>                       "Yellow","Orange","Red","Violet"/)
>
> contour_t2 = wrf_contour(a,wks,Tc,opts)
>
> delete(opts)
>
>
>
>
>
> ; Plotting options for Wind Vectors
>
>  opts = res
>
>  opts at NumVectors = 47           ; density of wind barbs
>
> vector = wrf_vector(a,wks,u10,v10,opts)
>
> delete(opts)
>
>
>
> plot = wrf_map_overlays(a,wks,(/contour_t2,vector/),pltres ,mpres)
>
>
>
> shp_name = "/home/Enric/Fronteres/FEsp/ESP_adm2.shp"
>
> lnres = True
>
> lnres at gsLineColor = "gray25"
>
> lnres at gsLineThicknessF = 0.5
>
>
>
> spain_id = gsn_add_shapefile_polylines(wks,plot,shp_name,lnres)
>
> draw(plot)
>
> frame(wks)
> end
>
> _______________________________________________
> 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/20160626/74b185d4/attachment.html