[ncl-talk] ensamble script

Luis del Castillo luis.delcastillo at cathalac.int
Wed Jan 25 10:07:14 MST 2017


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<mailto: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<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<mailto: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<mailto: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

[Description: http://www.cathalac.org/firma/format_02a.gif]

[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<tel:(507)%20317-3235>
Email: luis.delcastillo at cathalac.int<mailto:luis.delcastillo at cathalac.int>
<http://www.cathalac.org/>[Description: http://www.cathalac.org/firma/facebook.png]<https://www.facebook.com/cathalac> [Description: http://www.cathalac.org/firma/twitter.png] <http://twitter.com/#!/CATHALAC>  [Description: http://www.cathalac.org/firma/linkindin.png] <http://www.linkedin.com/company/2324157>  [Description: http://www.cathalac.org/firma/youtube.png] <http://www.youtube.com/user/CATHALACTV>




_______________________________________________
ncl-talk mailing list
ncl-talk at ucar.edu<mailto: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/89bb0e9d/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image001.png
Type: image/png
Size: 3374 bytes
Desc: image001.png
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170125/89bb0e9d/attachment.png 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image002.png
Type: image/png
Size: 3990 bytes
Desc: image002.png
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170125/89bb0e9d/attachment-0001.png 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image003.png
Type: image/png
Size: 329 bytes
Desc: image003.png
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170125/89bb0e9d/attachment-0002.png 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image004.png
Type: image/png
Size: 471 bytes
Desc: image004.png
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170125/89bb0e9d/attachment-0003.png 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image005.png
Type: image/png
Size: 357 bytes
Desc: image005.png
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170125/89bb0e9d/attachment-0004.png 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image006.png
Type: image/png
Size: 629 bytes
Desc: image006.png
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170125/89bb0e9d/attachment-0005.png 


More information about the ncl-talk mailing list