[ncl-talk] ensamble script

Wed Jan 25 22:43:17 MST 2017

It looks like you have a version of NCL that was built without shapefile
support.

Also, you have a rather old version of NCL.  I suggest upgrading to 6.3.0.
If you download one of our precompiled binaries, then shapefile support
will be included.

Please see:

http://www.ncl.ucar.edu/Download/

and follow the instructions for downloading a precompiled NCL binary and
NOT from source code.

If you are familiar with conda, then you can also install NCL using that:

http://www.ncl.ucar.edu/Download/conda.shtml

--Mary

On Wed, Jan 25, 2017 at 10:07 AM, Luis del Castillo <
luis.delcastillo at cathalac.int> wrote:

> Good afternoon Mary Haley I did what you suggest but I am getting this
> error
>
>
>
> [wrf at pantro d2_scripts]$ ncl shape.ncl
>
> Copyright (C) 1995-2013 - All Rights Reserved
>
> University Corporation for Atmospheric Research
>
> NCAR Command Language Version 6.1.2
>
> The use of this software is governed by a License Agreement.
>
> See http://www.ncl.ucar.edu/ for more details.
>
> (0)     Numero de Archivos Dom 2 = 73
>
> (0)     Working on time: 2017-01-25_00:00:00
>
> (0)     Working on time: 2017-01-25_01:00:00
>
> warning:["FileSupport.c":3026]:NCL does not know anything about file
> suffix <shp>.
>
> But NCL will try its best to figure out the file format.
>
>
>
> fatal:An internal error has occurred. The file format requested does not
> appear to be supported, could not open (PAN_adm0.shp)
>
> (0)     Error: gsn_add_shapefile_polylines: Can't open shapefile
> 'PAN_adm0.shp'
>
> (0)            No shapefile information will be added.
>
> (0)     Working on time: 2017-01-25_02:00:00
>
> warning:["FileSupport.c":3026]:NCL does not know anything about file
> suffix <shp>.
>
> But NCL will try its best to figure out the file format.
>
>
>
> fatal:An internal error has occurred. The file format requested does not
> appear to be supported, could not open (PAN_adm0.shp)
>
> (0)     Error: gsn_add_shapefile_polylines: Can't open shapefile
> 'PAN_adm0.shp'
>
> (0)            No shapefile information will be added.
>
>
>
>
>
> *From:* Mary Haley [mailto:haley at ucar.edu]
> *Sent:* martes, 24 de enero de 2017 5:38 p. m.
>
> *To:* Luis del Castillo
> *Cc:* ncl-talk at ucar.edu
> *Subject:* Re: [ncl-talk] ensamble script
>
>
>
> You can only use the High-Res map database for coastal outlines, and not
> for internal outlines.
>
>
>
> If you need more high-res outlines with WRF output files, then I recommend
> using shapefiles.
>
>
>
> See example shapefiles_14.ncl which shows how to add shapefile outlines to
> a plot's been generated with wrf_xxxx routines:
>
>
>
> http://www.ncl.ucar.edu/Applications/shapefiles.shtml#ex14
>
>
>
> The key is that you need to set:
>
>
>
>   pltres = True
>
>   mpres = True
>
>
>
> ; Tells wrf_map_overlays not to remove contours
>
>   pltres at PanelPlot  = True
>
>   plot = wrf_map_overlays(a,wks,contour,pltres,mpres)
>
>
>
> before you call wrf_map_overlays, which keeps the plot intact so you can
> add shapefile outlines with code that looks something like this:
>
>
>
> ;---Attach the shapefile outlines
>
>   lnres                  = True
>
>   lnres at gsLineColor      = "gray25"
>
>   ;lnres at gsLineThicknessF = 0.5
>
>   id = gsn_add_shapefile_polylines(wks,plot,"*xxxx*.shp",lnres)
>
>   draw(plot)
>
>   frame(wks)
>
> --Mary
>
>
>
>
>
>
>
> On Mon, Jan 23, 2017 at 5:43 PM, Luis del Castillo <
> luis.delcastillo at cathalac.int> wrote:
>
> Thanks everything seems to work just fine Mary Haley. Another question how
> to draw in land borders lines when using mpres at mpDataBaseVersion      =
> "HighRes".
>
> I was using
>
>   mpres at mpOutlineBoundarySets  = "National"
>
>   mpres at mpNationalLineColor    = "Black"
>
>   mpres at mpNationalLineThicknessF = 2.0
>
> But It seems to be ignore by the mpDataBaseVersion
>
>
>
> Thanks
>
> *From:* Mary Haley [mailto:haley at ucar.edu]
> *Sent:* lunes, 23 de enero de 2017 1:31 p. m.
> *To:* Luis del Castillo
> *Cc:* ncl-talk at ucar.edu
> *Subject:* Re: [ncl-talk] ensamble script
>
>
>
> Luis,
>
>
>
> What exactly do you mean that each WRF file has "a difference scheme"?
>
>
>
> If all of your WRF plots are defined over the same lat/lon area and the
> same map projection parameters, then you can call wrf_contour with any one
> of your "a" variables:
>
> contour_tend = wrf_contour(a1[0],wks,prep_tot_ens,opts_r) ; total (color)
>
>
> I'm guessing this has to be the case, because you are doing calculations
> on arrays read from different files.  If they weren't on the same lat/lon
> grid, then your calculations would potentially be wrong.
>
> --Mary
>
>
>
>
>
> On Fri, Jan 20, 2017 at 7:09 PM, Luis del Castillo <
> luis.delcastillo at cathalac.int> wrote:
>
> Hi fellows good evening
>
>
>
> I am having some issue trying to create an ensemble script with ncl, the
> idea is to take 3 different wrfout each one has a different scheme and used
> to do a 1hr accum precipitation average. On the attached file you can see I
> got almost everything ready but get stock at the end with the WRF ncl
> function wrf_countour this function needs the  ( file , wks, prep_tot_ens,
> opts_r) but as you see I have 3 a ‘s, a1 is from Schem 1, a2 from Schem 2,
> a3 Schem 3. So I cant not use it with the wrf_countour function so I am a
> little bit lost I am new with NCL and I got this running with GRADS.
>
>
>
> thanks
>
>
>
> [image: Description: http://www.cathalac.org/firma/format_02a.gif]
>
> *[image: Description: http://www.cathalac.org/firma/logocatha.gif]*
> <http://www.cathalac.org/>
> * Luis Alejandro Del Castillo Riley *
> *Gerente de Tecnologia*
> DIVISIÓN DE TECNOLOGÍAS DE LA INFORMACIÓN Y COMUNICACIÓN
> Tel: 507-317-3235 <(507)%20317-3235>
>
> *Email: **luis.delcastillo at cathalac.int <luis.delcastillo at cathalac.int>*
> <http://www.cathalac.org/>[image: Description:
> http://www.cathalac.org/firma/facebook.png]
> <https://www.facebook.com/cathalac> [image: Description:
> http://www.cathalac.org/firma/twitter.png]
> <http://twitter.com/#!/CATHALAC> [image: Description:
> http://www.cathalac.org/firma/linkindin.png]
> <http://www.linkedin.com/company/2324157> [image: Description:
> http://www.cathalac.org/firma/youtube.png]
> <http://www.youtube.com/user/CATHALACTV>
>
>
>
>
>
>
> _______________________________________________
> 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/20170125/d9de31a4/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image006.png
Type: image/png
Size: 629 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170125/d9de31a4/attachment.png 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image003.png
Type: image/png
Size: 329 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170125/d9de31a4/attachment-0001.png 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image001.png
Type: image/png
Size: 3374 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170125/d9de31a4/attachment-0002.png 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image002.png
Type: image/png
Size: 3990 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170125/d9de31a4/attachment-0003.png 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image004.png
Type: image/png
Size: 471 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170125/d9de31a4/attachment-0004.png 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image005.png
Type: image/png
Size: 357 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170125/d9de31a4/attachment-0005.png 