Re: Help: Bugs in "gsn_panel" function with graphic array "plots" in NCAR Command Language Version 5.2.0 ?

Hi,

Just real quick: the problem seems to be with this resource:

            res@gsnContourNegLineDashPattern = 2

If I comment this out, your script works. I had the same problem in V5.1.1, so it's not just a V5.2.0 issue.

I'll look into this further.

--Mary

On May 14, 2010, at 4:19 PM, lxf wrote:

> Hi,
> I am not sure if there still have bug in Version 5.2.0 in gsn_panel if I use a "plots" array with multiple calls.
> My problem is the attached t.ncl using t.nc can draw picture one by one, however, can not panel them up together using gsn_panel?
>
> The script can not be excuted like the following status after ">ncl -x t.ncl ":
>
> Copyright (C) 1995-2010 - All Rights Reserved
> University Corporation for Atmospheric Research
> NCAR Command Language Version 5.2.0
> The use of this software is governed by a License Agreement.
> See http://www.ncl.ucar.edu/ for more details.
> Loading file "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> +
> Loading file "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> +
> Loading file "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> +
> Loading file "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
> +
> + undef("sub_draw_hor")
> + function sub_draw_hor(wks,h,opt,ti)
> + local wks,h,res,res1,plot,plot1,ploth,opt,ti
> + begin
> + ;
> + ; draw Plots
> + ;......
> + print (""+ti)
> + ;wks =gsn_open_wks("eps",outname)
> + ;gsn_define_colormap(wks,"nrl_sirkes") ; choose color map
> + res = True
> + res@gsnDraw = False ; don't draw yet
> + res@gsnFrame = False ; don't advance frame yet
> + res@gsnMaximize = False
> +
> + res@cnFillOn = False ;True ; turn on color for contours
> + res@cnLinesOn = True ;False ; turn off contour lines
> + res@cnLineLabelsOn = False
> +
> + ;res@vpXF = 0.2
> + ;res@vpYF = 0.8
> + res@vpWidthF = 0.14
> + res@vpHeightF = 0.08
> + res@gsnPaperOrientation = "portrait"
> +
> + res@tiMainString = "" ; ti
> + res@tiMainOn = False ; ti
> + res@tiXAxisOn = False
> + res@gsnLeftString = ti
> + res@gsnCenterString = ""
> + res@gsnRightString = ""
> + res@tmXBOn = opt
> + res@tmXTOn = False
> + res@tmYROn = False
> +
> + res@cnLineLabelsOn = True ;False ; turn off line labels
> + res@cnInfoLabelOn = False
> +
> + res@tmYLMode = "Explicit" ; label independently
> + res@tmYLValues = h&time(::8)
> + res@tmYLLabels = (/0,0.5,1,1.5,2/)
> + ;res@tmYLLabelFontHeightF = 0.020
> +
> + res@gsnContourLineThicknessesScale = 1.5
> + res@gsnContourZeroLineThicknessF = 1.5 ; doubles thickness of zero contour
> + res@gsnContourNegLineDashPattern = 2 ; sets negative contours to dash pattern 1
> +
> +
> + plot = gsn_csm_contour(wks, h(time
> +
> +
> + return(plot)
> + end
> +
> + ;====================================================================
> + begin
> +
> +
> + outname = "./t"
> +
> +
> + f = addfile ("t.nc" , "r")
> + a0 = f->a0
> + a1 = f->a1
> + a2 = f->a2
> + a3 = f->a3
> +
> +
> + ;
> + ; draw Plots
> + ;......
> + wks =gsn_open_wks("eps",outname)
> + gsn_define_colormap(wks,"nrl_sirkes") ; choose color map
> + plots = new(4,graphic)
> +
> +
> + opt = False
> + ti = "a)"
> + plots(0)=sub_draw_hor(wks,a0,opt,ti)
> +
> + opt = False
> + ti = "c)"
> + plots(1)=sub_draw_hor(wks,a1,opt,ti)
> +
> + opt = False
> + ti = "b)"
> + plots(2)=sub_draw_hor(wks,a2,opt,ti)
> +
> + opt = False
> + ti = "d)"
> + plots(3)=sub_draw_hor(wks,a3,opt,ti)
> +
> +
> +
> + resP = True ; modify the panel plot
> + resP@gsnMaximize = True
> +
> + resP@gsnFrame = False
> + resP@gsnPaperOrientation = "portrait"
> + resP@gsnPanelBottom = 0.05
> + resP@gsnPanelLabelBar = False
> + resP@gsnPanelScalePlotIndex = 1
> +
> +
> +
> +
> + gsn_panel(wks,plots,(/2,2/),resP)
> + ;draw(plots(0))
> + ;draw(plots(1))
> + ;draw(plots(2))
> + ;draw(plots(3))
> + frame(wks)
> +
> +
> + end
> (0) a)
> (0) c)
> (0) b)
> (0) d)
>
> It just stop here, with cpu 100% using and will not come on!
>
> And the result eps has the last figure "plots(3)" been igoned, which can surely be drawn out if not using gsn_panel and not calling the plots(0:2) array before this. It's just so strange!
>
> Any help or sugestion is appriciated.
>
> Sincerely,
> Li<t.nc><t.ncl><t.eps>_______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri May 14 17:19:09 2010