[ncl-talk] My script does not work
Alan Brammer
abrammer at albany.edu
Mon Mar 14 15:25:22 MDT 2016
Errors line by line. Starting at the top.
> warning:wkImageFileName is not a resource in the given object
> warning:wkImageFormat is not a resource in the given object
I don’t know what these are and I don’t get them. They haven’t popped up on ncl talk for many years.
> warning:mpProjection is not a valid resource in test_contour at this time
> warning:mpMinLonF is not a valid resource in test_contour at this time
> warning:mpMaxLonF is not a valid resource in test_contour at this time
> warning:mpMinLatF is not a valid resource in test_contour at this time
> warning:mpMaxLatF is not a valid resource in test_contour at this time
> warning:mpShapeMode is not a valid resource in test_contour at this time
> warning:mpFillOn is not a valid resource in test_contour at this time
> warning:mpOutlineOn is not a valid resource in test_contour at this time
> warning:mpDataBaseVersion is not a valid resource in test_contour at this time
You’re passing mp resources to gsn_csm_contour() . As that is no longer a map function it doesn’t understand mp resources. Delete all the lines where res at mp…. and set the below instead. Then pass the new resources to wrf_map() as below.
mpres = True
mpres at mpFillOn = False
mpres at mpMaxLatF = 38.20 ; specify the plot domain
mpres at mpMinLatF = 37.75 ;
mpres at mpMinLonF = -120.00 ;
mpres at mpMaxLonF = -119.10 ;
mpres at mpOutlineOn = True ; turn the map outline on
mpres at mpProjection = "CylindricalEquidistant" ;The default
mpres at pmTickMarkDisplayMode = "Always"
mpres at mpDataBaseVersion = "MediumRes"
mpres at mpOutlineOn =True
mpres at mpShapeMode = “FreeAspect"
map = wrf_map(wks, a, mpres)
> warning:cnLineOn is not a valid resource in test_contour at this time
Still got a typo here. Not sure how you edit these files, but I’d recommend looking at https://www.ncl.ucar.edu/Applications/editor.shtml the coloring will help pick up simple typos.
> (0) wrf_map_overlay: Warning: This procedure is obsolete. Consider
> (0) using wrf_map_overlays instead.
As it says, use wf_map_overlays() instead. Though this still works for me.
> wrf_map_overlay(wrf,map,(/contour,vector/),True)
You’ve changed wks to wrf here? Not sure why. This causes the main error.
Here’s the edited version, with the newer overlay function. I would go through the edits above and then match up to my version.
> On 14 Mar 2016, at 16:36, Kerwyn Texeira <ktish86 at gmail.com> wrote:
>
> Hi Alan Brammer,
>
> I appreciate your help in this matter!
>
> I still got a few errors:
>
> 1. When i use overlay, I get: warning:wkImageFileName is not a resource in the given object
> warning:wkImageFormat is not a resource in the given object
> warning:PlotManagerSetValues: TickMark annotation cannot be added after NhlCreate
> warning:mpProjection is not a valid resource in test_contour at this time
> warning:mpMinLonF is not a valid resource in test_contour at this time
> warning:mpMaxLonF is not a valid resource in test_contour at this time
> warning:mpMinLatF is not a valid resource in test_contour at this time
> warning:mpMaxLatF is not a valid resource in test_contour at this time
> warning:mpShapeMode is not a valid resource in test_contour at this time
> warning:mpFillOn is not a valid resource in test_contour at this time
> warning:mpOutlineOn is not a valid resource in test_contour at this time
> warning:mpDataBaseVersion is not a valid resource in test_contour at this time
> warning:cnLineOn is not a valid resource in test_contour at this time
> (0) wrf_map_overlay: Warning: This procedure is obsolete. Consider
> (0) using wrf_map_overlays instead.
> warning:ContourPlotDraw: data boundary is out of range
>
> 2. I do not understand part of no.5 "you want to transfer the mp resources to a new variable and pass that to wrf_map() instead of True"
>
> The revised script:
>
> ;*************************************************
> ; overlay_1.ncl
> ;
> ; Concepts illustrated:
> ; - Overlaying line contours on filled contours
> ; - Explicitly setting contour levels
> ; - Selecting a different color map
> ;
> ;*************************************************
> ;
> ; These files are loaded by default in NCL V6.2.0 and newer
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
>
> begin
> a = addfile("./wrfout_d03_2014-01-11_21:00:00.nc <http://00.nc/>","r")
>
> it = 0
> hgt = wrf_user_getvar(a, "HGT", it)
> hgt at lat2d = wrf_user_getvar(a, "XLAT", it)
> hgt at lon2d = wrf_user_getvar(a, "XLONG", it)
> u = wrf_user_getvar(a, "ua", it)
> v = wrf_user_getvar(a, "va", it)
> p = wrf_user_getvar(a, "pressure", it)
>
> u_wind = wrf_user_intrp3d(u, p, "h", 650., 0.0, False)
> v_wind = wrf_user_intrp3d(v, p, "h", 650., 0.0, False)
>
>
> spd = (u_wind*u_wind + v_wind*v_wind)^(0.5) ;m/s
> ; spd = spd*1.94384449
> ; spd at units = "Wind Speed"
> ; spd at units = "m/s"
>
>
>
> wks = gsn_open_wks("png","test") ; send graphics to PNG file
> gsn_define_colormap(wks,"BlAqGrYeOrRe")
>
>
> res = True
> res at mpFillOn = False
> res at mpMaxLatF = 38.20 ; specify the plot domain
> res at mpMinLatF = 37.75 ;
> res at mpMinLonF = -120.00 ;
> res at mpMaxLonF = -119.10 ;
> res at mpOutlineOn = True ; turn the map outline on
> res at gsnDraw = False ; do not draw the plot
> res at gsnFrame = False ; do not advance the frame
> ; res at csnLinesOn = True
> ; res at tiMainString = a
>
> res at mpProjection = "CylindricalEquidistant" ;The default
> res at pmTickMarkDisplayMode = "Always"
> res at mpDataBaseVersion = "MediumRes"
> res at mpOutlineOn =True
> res at lbOrientation = "Vertical"
> res at tiMainOffsetYF = -0.03
> res at mpShapeMode = "FreeAspect"
> res at vpWidthF = 0.9
> res at vpHeightF = 0.9
>
> ; res at cnLevelSelectionMode = "ExplicitLevels" ; use explicit levels
> ; res at cnLevels = ispan(215,265,5) ; set the contour levels
> res at cnLineLabelsOn = False ; do not use line labels
> res at cnFillOn = True ; color fill
> res at cnLineOn = False ; do not draw contour lines
> ; res at cnFillPalette = "BlueDarkRed18"
>
> res at gsnAddCyclic = False
>
> map = wrf_map(wks, a, True)
>
> u_wind = u_wind*1.94384449
> v_wind = v_wind*1.94384449
>
> contour = gsn_csm_contour(wks,hgt,res)
>
> ;------wind vectors
> res2 = True
> res2 at gsnDraw = True
> res2 at gsnFrame = True
> res2 at FieldTitle = "Winds"
> res2 at NumVectors = 47
> res2 at vcWindBarbLineThicknessF= 2.5
>
> vector = wrf_vector(a,wks,u_wind,v_wind,res2)
> wrf_map_overlay(wrf,map,(/contour,vector/),True)
> ;overlay(contour,vector)
>
>
>
> end
>
>
>
>
>
> On Mon, Mar 14, 2016 at 9:04 AM, Alan Brammer <abrammer at albany.edu <mailto:abrammer at albany.edu>> wrote:
>
> Not sure where you’re getting those error messages from but there are many problems in the script. Here are the errors I got on running your code on my own file and the step by step to fix them. Make the below edits. If there are still problems send the whole new script to the list so we can see your edits.
>
>
>> warning:attsetvalues: the value associated with (tiMainString) does not have an HLU representation
>> warning:csnLinesOn is not a valid resource in gsnapp_contour at this time
>> (0) gsn_csm_map_ce: Fatal: The resources mpMinLatF/mpLeftCornerLatF must be less than the resources mpMaxLatF/mpRightCornerF.
>> (0) Execution halted.
>
>
> ;====
> 1.
> so the top error:
>> "warning:attsetvalues: the value associated with (tiMainString) does not have an HLU representation"
>
>
> Some thing is wrong with tiMainString. " res1 at tiMainString = a “ a is a pointer to a file. NCL has no idea how to turn that into a string.
> If you upgrade to 6.3.0 you could do the following.
> res1 at tiMainString = getfilepath(a) ;; needs 6.3.0!!!
> or
> res1 at tiMainString = “Some string value for a title"
>
>
> ;====
> 2.
>> warning:csnLinesOn is not a valid resource in gsnapp_contour at this time
>
> Self explanatory, csnLineOn isn’t a resource. This is a typo. Find that resource and fix it.
>
>
> ;====
> 3.
>> (0) gsn_csm_map_ce: Fatal: The resources mpMinLatF/mpLeftCornerLatF must be less than the resources mpMaxLatF/mpRightCornerF.
>> (0) Execution halted.
>
> Again, does what it says on the tin. minLat must be less than maxLat. Yours are not. Switch your min/max lat resources and you’ll be set.
>
>
> ;====
> 4.
> Once those are fixed, you’ll get the following error.
>> (0) wrf_map_overlay: Warning: This procedure is obsolete. Consider
>> (0) using wrf_map_overlays instead.
>> fatal:NhlAddOverlay: plot class mapPlotClass cannot be overlay plot member
>> warning:NhlRemoveOverlay: plot not found in overlay sequence
>
>
> You can’t overlay 2 maps on top of each other. wrf_map() and gsn_csm_contour_map() so just remove the _map part and make a plain contour plot.
>
>
> ;====
> 5.
> Lastly, you want to transfer the mp resources to a new variable and pass that to wrf_map() instead of True.
> and change the order in the overlay function … map, (/contour, vector/)...
>
>
>
>
>
>> On 13 Mar 2016, at 18:12, Kerwyn Texeira <ktish86 at gmail.com <mailto:ktish86 at gmail.com>> wrote:
>>
>> Hi All,
>>
>> For a month now I have been having trouble overlay windbards on zoomed terrain with specified lat and lon in ncl. I do not want the entire domain plotted with winds. If so, I would have used the multiOverlay.ncl script which works perfectly for many overlay for the entire domain. No images are being plotted with the script below. I'm not getting any fatal errors. I'm using NCL version 6.2.1. All I'm getting is the following; warning:wkImageFileName is not a resource in the given object
>> warning:wkImageFormat is not a resource in the given object.
>>
>> I would like some help with this script please. I looked online for a script to help with this (overlay on zoomed in terrain with specified lat and lon) but I’m having trouble finding it. Any advice will be greatly appreciated. Thanks, Kerwyn
>>
>> Script:
>>
>> ;*************************************************
>> ; overlay_1.ncl
>> ;
>> ; Concepts illustrated:
>> ; - Overlaying line contours on filled contours
>> ; - Explicitly setting contour levels
>> ; - Selecting a different color map
>> ;
>> ;*************************************************
>> ;
>> ; These files are loaded by default in NCL V6.2.0 and newer
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
>>
>> begin
>> a = addfile("./wrfout_d03_2014-01-11_21:00:00.nc <http://00.nc/>","r")
>>
>> it = 0
>> hgt = wrf_user_getvar(a, "HGT", it)
>> hgt at lat2d = wrf_user_getvar(a, "XLAT", it)
>> hgt at lon2d = wrf_user_getvar(a, "XLONG", it)
>> u = wrf_user_getvar(a, "ua", it)
>> v = wrf_user_getvar(a, "va", it)
>> p = wrf_user_getvar(a, "pressure", it)
>>
>> u_wind = wrf_user_intrp3d(u, p, "h", 650., 0.0, False)
>> v_wind = wrf_user_intrp3d(v, p, "h", 650., 0.0, False)
>>
>>
>> ; spd = (u_wind*u_wind + v_wind*v_wind)^(0.5) ;m/s
>> ; spd at units = "Wind Speed"
>> ; spd at units = "m/s"
>>
>>
>> wks = gsn_open_wks("png","test") ; send graphics to PNG file
>> gsn_define_colormap(wks,"BlAqGrYeOrRe")
>>
>>
>> res1 = True
>> res1 at mpFillOn = False
>> res1 at mpMaxLatF = 37.75 ; specify the plot domain
>> res1 at mpMinLatF = 38.20 ;
>> res1 at mpMinLonF = -120.00 ;
>> res1 at mpMaxLonF = -119.10 ;
>> res1 at mpOutlineOn = True ; turn the map outline on
>> res1 at gsnDraw = False ; do not draw the plot
>> res1 at gsnFrame = False ; do not advance the frame
>> res1 at csnLinesOn = True
>> res1 at tiMainString = a
>>
>> res1 at mpProjection = "CylindricalEquidistant" ;The default
>> res1 at pmTickMarkDisplayMode = "Always"
>> res1 at mpDataBaseVersion = "MediumRes"
>> res1 at mpOutlineOn =True
>> res1 at lbOrientation = "Vertical"
>> res1 at tiMainOffsetYF = -0.03
>> res1 at mpShapeMode = "FreeAspect"
>> res1 at vpWidthF = 0.9
>> res1 at vpHeightF = 0.9
>>
>> ; res at cnLevelSelectionMode = "ExplicitLevels" ; use explicit levels
>> ; res at cnLevels = ispan(215,265,5) ; set the contour levels
>> res1 at cnLineLabelsOn = False ; do not use line labels
>> res1 at cnFillOn = True ; color fill
>> res1 at cnLinesOn = False ; do not draw contour lines
>> ; res at cnFillPalette = "BlueDarkRed18"
>>
>> res1 at gsnAddCyclic = False
>>
>> map = wrf_map(wks, a, True)
>>
>> contour = gsn_csm_contour_map(wks,hgt,res1)
>>
>> ;------wind vectors
>> res2 = True
>> res2 at FieldTitle = "Winds"
>> res2 at NumVectors = 47
>> res2 at vcWindBarbLineThicknessF= 2.5
>> vector = wrf_vector(a, wks, u_wind, v_wind, res2)
>>
>> wrf_map_overlay(wks,map,(/vector, contour/),True)
>>
>> end
>> _______________________________________________
>> 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 <http://mailman.ucar.edu/mailman/listinfo/ncl-talk>
>
>
> _______________________________________________
> 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/20160314/993a5f9e/attachment.html
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ncl-talk_wrf_script.ncl
Type: application/octet-stream
Size: 2831 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160314/993a5f9e/attachment.obj
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160314/993a5f9e/attachment-0001.html
More information about the ncl-talk
mailing list