Sophie,
It sounds like something is wrong with your Cygwin installation. I need more information.
What version of NCL are you running?
ncl -V
Does this cause a seg fault:
ng4ex gsun01n -clean
If this works, it will pop up a window on your screen, which you need to click on with your left mouse button.
Did you download a precompiled NCL binary for Cygwin, or did you build it from source code?
How old is your version of Cygwin?
--Mary
On Jul 8, 2012, at 9:37 PM, Sophie HE wrote:
> Hi there,
>
> Whenever I would love to open an workstation using gsn_open_wks, even with the simplest command, say
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> wks = gsn_open_wks("x11","gsun01n")
> cygwin reports Segmentation fault (core dumped) error.
>
> Could you give me any advice, as I really have no idea of what's going on?
>
> Thanks.
>
> Regards,
> Sophie
>
>
> On Sun, Jul 8, 2012 at 2:00 AM, <ncl-talk-request@ucar.edu> wrote:
> Send ncl-talk mailing list submissions to
> ncl-talk@ucar.edu
>
> To subscribe or unsubscribe via the World Wide Web, visit
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> or, via email, send a message with subject or body 'help' to
> ncl-talk-request@ucar.edu
>
> You can reach the person managing the list at
> ncl-talk-owner@ucar.edu
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of ncl-talk digest..."
>
>
> Today's Topics:
>
> 1. Re: Stripping white space at the edges (Paul Scorer)
> 2. Re: Problem with spacing between plots in a panel plot
> (Adam Phillips)
> 3. eof function (163)
> 4. Re: Stripping white space at the edges (Brandon Capasso)
> 5. Dimension mismatch on dim_sum_n the dimension relating to
> file input (Cannavina, Dominique)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Fri, 06 Jul 2012 23:54:26 +0100
> From: Paul Scorer <paul.scorer01@gmail.com>
> Subject: Re: Stripping white space at the edges
> To: Dennis Shea <shea@ucar.edu>
> Cc: NCL Talk <ncl-talk@ucar.edu>
> Message-ID: <1341615266.4491.7.camel@hp-compaq.home>
> Content-Type: text/plain; charset="UTF-8"
>
> On recent-ish Linux systems, there is available a fork of ImageMagick
> known as GraphicsMagick which is significantly faster, particularly on
> multi-core h/w.
>
> Use thus:
> $ gm convert ...
> etc
>
> Options, etc are (essentially) the same.
>
> HTH
>
> PaulS
>
>
> On Fri, 2012-07-06 at 10:39 -0600, Dennis Shea wrote:
> > the following may be overkill but
> >
> > convert -flatten -trim +repage -background white
> > pr_zonalAvg_clim7100.png pr_zonalAvg_clim7100.trim.png
> >
> >
> > On 07/06/2012 10:34 AM, Noel Aloysius wrote:
> > > Hello NCL,
> > >
> > > I created a panel plot with 14 figures. I want to maximize the plot area
> > > by stripping the white space at the edges. I am not sure which plot
> > > resources do this.
> > > The section of the code and the plot is attached.
> > >
> > > Thank you in advance for your help,
> > > Noel
> > >
> > > ;***********************************************************************************************
> > > wks = gsn_open_wks("png", diro + "pr_zonalAvg")
> > > plot = new(14, graphic)
> > > gsn_define_colormap(wks, "precip3_16lev")
> > >
> > > res1 = True ; plot mods desired
> > > res1@gsnDraw = False
> > > res1@gsnFrame = False
> > > res1@gsnMaximize = True
> > > res1@gsnRightStringFontHeightF = 0.04
> > > res1@cnFillOn = True ; turn on color
> > > res1@gsnSpreadColors = True ; use full color table
> > > res1@cnLinesOn = False ; no contour lines
> > > res1@cnInfoLabelOn = False
> > > res1@pmTickMarkDisplayMode = "Always" ; turn on fancy tickmarks
> > >
> > > res1@tmXBMode = "Explicit" ; label independently
> > > res1@tmXBValues = ispan(0,11,1)
> > > res1@tmXBLabels =
> > > (/"J","F","M","A","M","J","J","A","S","O","N","D"/)
> > >
> > > res1@lbLabelBarOn = False
> > > res1@cnLevelSelectionMode = "ManualLevels"
> > > res1@cnLevelSpacingF = 20.
> > > res1@cnMinLevelValF = 20.
> > > res1@cnMaxLevelValF = 300.
> > >
> > > res1@tmYROn = False
> > > res1@tmXBLabelFontHeightF = 0.025
> > > res1@tmYLLabelFontHeightF = 0.025
> > >
> > > do i=0,11
> > > res1@gsnRightString = gcm_names(i)
> > > plot(i) = gsn_csm_lat_time(wks,
> > > PzoneLat_7100(raw_bc|0,GCM|i,lat|:,month|:), res1)
> > > end do
> > >
> > > res1@gsnRightString = "Obs"
> > > plot(12) = gsn_csm_lat_time(wks, PzoneLat_7100obs(lat|:,month|:), res1)
> > >
> > > res1@gsnRightString = "bias corr. GCM avg."
> > > plot(13) = gsn_csm_lat_time(wks, PzoneLat_7100_BCavg(lat|:,month|:),
> > > res1)
> > >
> > > pres1 = True
> > > pres1@gsnPanelLabelBar = True
> > > pres1@gsnFrame = False
> > > pres1@lbLabelFontHeightF = 0.01
> > > pres1@lbOrientation = "horizontal"
> > > pres1@lbLabelAutoStride = True
> > > pres1@gsnPanelYWhiteSpacePercent = 1.5
> > > pres1@gsnBoxMargin = 0.0
> > > pres1@gsnPanelRowSpec = True ; tell panel what order
> > > to plot
> > > pres1@gsnPanelCenter = True ; centered plot
> > > pres1@gsnMaximize = True
> > >
> > > gsn_panel(wks,plot,(/6,6,2/),pres1)
> > > frame(wks)
> > > ;****************************************************************************************************
> > >
> > >
> > >
> > > _______________________________________________
> > > 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
>
>
>
>
> ------------------------------
>
> Message: 2
> Date: Fri, 6 Jul 2012 18:41:19 -0600
> From: Adam Phillips <asphilli@ucar.edu>
> Subject: Re: Problem with spacing between plots in a panel
> plot
> To: Soumik Basu <sbasu@iarc.uaf.edu>
> Cc: Ncl Talk <ncl-talk@ucar.edu>
> Message-ID: <CCB59B1E-4C9D-40E9-A79C-4B02574960E1@ucar.edu>
> Content-Type: text/plain; charset=us-ascii
>
> Hi Soumik,
> When you are paneling plots, it is recommended that you do not set gsnMaximize = True in your individual plot resource lists. (res and sres in your case) It is of course fine to set gsnMaximize in the panel resource list. Try removing gsnMaximize from your res and sres resource lists and see if that helps... Adam
>
>
>
> On Jul 6, 2012, at 3:58 PM, Soumik Basu <sbasu@iarc.uaf.edu> wrote:
>
> > Hi Adam,
> >
> > Thanks for the reply.I made those changes but I am still getting the similar plots. Can you please suggest what else can cause it. here is the complete plotting resources:
> >
> >
> > ;*******************************************************
> > ;Plot resources for eady growth first set (non shaded)
> > ;******************************************************
> >
> > wks = gsn_open_wks("png","eady_eke_var_diff_spring_winter")
> > gsn_define_colormap(wks,"BlAqGrYeOrReVi200")
> >
> >
> > res = True
> > res@cnLinesOn = True
> >
> > res@tmLabelAutoStride = True
> >
> > res@tmXTOn = False
> > res@tmYROn = False
> >
> > res@cnLevelSelectionMode = "ManualLevels"
> > res@cnMinLevelValF = 0
> > res@cnMaxLevelValF = 2.0
> > res@cnLevelSpacingF = 0.2
> > res@cnLineLabelFontHeightF= 0.015
> > res@cnInfoLabelOn = False ; turn off cn info label
> > res@cnLineLabelsOn = False
> >
> > res@lbLabelBarOn = False
> >
> > res@gsnDraw = False
> > res@gsnFrame = False
> > res@gsnMaximize = True
> >
> > res@pmLabelBarDisplayMode = "Never"
> >
> > ;*************************************************
> > ; Plot resources for eke first set (shaded)
> > ;*************************************************
> > plot = new(4,graphic) ; create a plot array
> >
> > sres = True
> >
> > sres@gsnDraw = False
> > sres@gsnFrame = False
> > sres@gsnMaximize = True
> > sres@gsnSpreadColors = True
> >
> > sres@cnFillOn = True
> > sres@cnLinesOn = False
> > sres@cnLevelSelectionMode = "ManualLevels"
> > sres@cnMinLevelValF = 0
> > sres@cnMaxLevelValF = 250
> > sres@cnLevelSpacingF = 5
> > sres@cnInfoLabelOn = False ; turn off cn info label
> > sres@cnLineLabelsOn = False
> >
> > sres@tmLabelAutoStride = True
> > sres@tmXTOn = False
> > sres@tmYROn = False
> >
> > sres@mpLimitMode = "LatLon"
> > sres@mpMinLatF = 10.
> > sres@mpMaxLatF = 80.
> > sres@mpMinLonF = 130.
> > sres@mpMaxLonF = 330.
> > sres@mpGeophysicalLineThicknessF = 1.75
> > sres@mpCenterLonF = 210
> > sres@mpShapeMode = "FreeAspect"
> > sres@lbOrientation = "vertical" ;vertical label bar
> > sres@lbLabelBarOn = False
> > sres@lbLabelAutoStride = True
> >
> >
> > ;*********************************
> > ; Plot first set (four means)
> > ;*********************************
> >
> > plot_eke1 = gsn_csm_contour_map(wks,eke_var(:,:),sres)
> > plot_eady1 = gsn_csm_contour(wks,eady_var(:,:),res)
> > overlay(plot_eke1,plot_eady1)
> > plot(0) = plot_eke1
> >
> > plot_eke2 = gsn_csm_contour_map(wks,eke_var_ens(:,:),sres)
> > plot_eady2 = gsn_csm_contour(wks,eady_var_ens(:,:),res)
> > overlay(plot_eke2,plot_eady2)
> > plot(1) = plot_eke2
> >
> > res@gsnLeftString = " "
> >
> > plot_eke3 = gsn_csm_contour_map(wks,eke_var1(:,:),sres)
> > plot_eady3 = gsn_csm_contour(wks,eady_var1(:,:),res)
> > overlay(plot_eke3,plot_eady3)
> > plot(2) = plot_eke3
> >
> > plot_eke4 = gsn_csm_contour_map(wks,eke_var_ens1(:,:),sres)
> > plot_eady4 = gsn_csm_contour(wks,eady_var_ens1(:,:),res)
> > overlay(plot_eke4,plot_eady4)
> > plot(3) = plot_eke4
> >
> > ;************************************************
> > ; create first set panel plot(four means)
> > ;************************************************
> > resPanel = True
> >
> > resPanel@gsnFrame = False
> >
> > resPanel@gsnSpreadColors = True
> > resPanel@gsnPanelLabelBar = True ; add color bar
> > resPanel@gsnPanelFigureStrings = (/"a","b","c","d"/)
> > resPanel@gsnPanelFigureStringsFontHeightF = 0.02
> > resPanel@gsnPanelFigureStringsPerimOn = True
> > resPanel@gsnPanelDebug = True
> > resPanel@gsnPanelTop = 0.9
> > resPanel@gsnPanelBottom = 0.3 ;move bottom up from 0 to 0.4
> > resPanel@gsnPanelLeft = 0.0
> > resPanel@gsnPanelRight = 0.9
> > resPanel@gsnPanelDebug = True
> >
> > resPanel@lbOrientation = "Vertical"
> > resPanel@lbLabelAutoStride = False
> > resPanel@lbLabelStride = 5
> >
> > gsn_panel(wks,plot,(/2,2/),resPanel)
> >
> > ;*****************************************************************
> > ; create plot resources eady growth diff next set of plots(non shaded)
> > ;*****************************************************************
> > plot1 = new(2,graphic) ; create a plot array
> >
> > dres = True
> >
> > dres@gsnMaximize = True
> >
> >
> > dres@cnLinesOn = True
> > dres@cnLevelSelectionMode = "ManualLevels"
> > dres@cnMinLevelValF = 0
> > dres@cnMaxLevelValF = 3
> > dres@cnLevelSpacingF = 0.2
> > dres@cnInfoLabelOn = False
> > dres@cnLineLabelsOn = False
> > dres@cnLineLabelDensityF = 1.
> > dres@cnLineLabelFontHeightF = 0.015
> > dres@cnLineLabelInterval = 3
> >
> > dres@tmLabelAutoStride = True
> > dres@tmXTOn = False
> > dres@tmYROn = False
> >
> > dres@lbLabelBarOn = False
> >
> > dres@gsnDraw = False
> > dres@gsnFrame = False
> > dres@gsnMaximize = True
> >
> >
> > ;*****************************************************
> > ; create plot resources eke next set of plots(shaded)
> > ;******************************************************
> >
> > dbres = True
> >
> > dbres@gsnMaximize = True
> > dbres@gsnSpreadColors = True
> > dbres@gsnDraw = False
> > dbres@gsnFrame = False
> >
> > dbres@cnFillOn = True
> > dbres@cnLinesOn = False
> > dbres@cnLevelSelectionMode = "ManualLevels"
> > dbres@cnMinLevelValF = -40
> > dbres@cnMaxLevelValF = 40
> > dbres@cnLevelSpacingF = 0.5
> > dbres@cnInfoLabelOn = False
> > dbres@cnLineLabelsOn = False
> >
> >
> > dbres@tmLabelAutoStride = True
> > dbres@tmXTOn = False
> > dbres@tmYROn = False
> >
> > dbres@lbOrientation = "Vertical"
> > dbres@lbLabelBarOn = False
> > dbres@lbLabelStride = 20.
> > dbres@lbLabelAutoStride = False
> >
> > dbres@mpLimitMode = "LatLon"
> > dbres@mpMinLatF = 10.
> > dbres@mpMaxLatF = 80.
> > dbres@mpMinLonF = 130.
> > dbres@mpMaxLonF = 330.
> > dbres@mpCenterLonF = 210
> > dbres@mpGeophysicalLineThicknessF= 1.75
> >
> > getvalues plot(3)
> > "vpWidthF" : vpw
> > "vpHeightF" : vph
> > end getvalues
> >
> >
> > ; By default, you can't skew the width/height of a map plot.
> > ; You must also set mpShapeMode to "FreeAspect" to indicate
> > ; that you really want to skew the map.
> >
> > dbres@mpShapeMode = "FreeAspect"
> > dbres@vpWidthF = vpw
> > dbres@vpHeightF = vph
> >
> >
> > ; ## SHADED
> >
> > plot1_eke1 = gsn_csm_contour_map(wks,eke_diff,dbres)
> > plot1_eke2 = gsn_csm_contour_map(wks,eke_diff1,dbres)
> >
> > ; ## NON-SHADED
> >
> > plot1_eady1 = gsn_csm_contour(wks,eady_diff,dres)
> > plot1_eady2 = gsn_csm_contour(wks,eady_diff1,dres)
> >
> > overlay(plot1_eke1,plot1_eady1)
> > overlay(plot1_eke2,plot1_eady2)
> >
> > plot1(0) = plot1_eke1
> > plot1(1) = plot1_eke2
> >
> >
> >
> > ;************************************************
> > ; create next set panel plot
> > ;************************************************
> > dresPanel = True
> >
> > dresPanel@gsnFrame = False
> > dresPanel@gsnSpreadColors = True
> > dresPanel@gsnMaximize = True ; use full page
> > dresPanel@gsnPanelLabelBar = True ; add color bar
> > dresPanel@gsnPanelTop = 0.45
> > dresPanel@gsnPanelBottom = 0.1
> > dresPanel@gsnPanelLeft = 0.0
> > dresPanel@gsnPanelRight = 0.9
> > dresPanel@gsnPanelFigureStrings = (/"e","f"/)
> > dresPanel@gsnPanelFigureStringsFontHeightF = 0.02
> > dresPanel@gsnPanelFigureStringsPerimOn = True
> > dresPanel@gsnPanelDebug = True
> >
> > dresPanel@lbOrientation = "Vertical"
> > dresPanel@lbLabelAutoStride = False
> > dresPanel@lbLabelStride = 10
> >
> > gsn_panel (wks,plot1,(/1,2/),dresPanel)
> >
> > frame(wks)
> >
> >
> >
> > Thanks,
> > Soumik
> >
> >
> > On 2012-07-05 19:20, Adam Phillips wrote:
> >> Hi Soumik,
> >> Your top 4 plots are not taking up all of the space you are allotting
> >> to them. They are maxed out on the page horizontally, and they are
> >> only taking up roughly 0.5-0.9 of vertical NDC space. To take up more
> >> space, NCL would need to alter the aspect ratio of the plot(s),
> >> something it won't do unless you tell it to. (Setting res@mpShapeMode=
> >> "FreeAspect" would do that. )
> >>
> >> Try setting gsnPanelTop for your bottom 2 plots to 0.45 (instead of 0.3)..
> >> Hope that helps.. Adam
> >>
> >>
> >> On Jul 5, 2012, at 6:21 PM, Soumik Basu <sbasu@iarc.uaf.edu> wrote:
> >>
> >>> Hi,
> >>>
> >>> I am trying to create a panel plot containing 6 plots. But I am having
> >>> problem with the space between the plots.I am plotting the first four
> >>> plots with a particular colorbar and the last two with another colorbar.
> >>> Can anybody please suggest me what is going wrong.
> >>>
> >>> The panel plot resources (I am sending only the two separate resource
> >>> parts of my script) I am using are as follows:
> >>>
> >>> ;************************************************
> >>> ; create first set panel plot(top four of the panel)
> >>> ;************************************************
> >>> resPanel = True
> >>>
> >>> resPanel@gsnFrame = False
> >>>
> >>> resPanel@gsnSpreadColors = True
> >>> resPanel@gsnPanelLabelBar = True ; add
> >>> color bar
> >>> resPanel@gsnPanelFigureStrings = (/"a","b","c","d"/)
> >>> resPanel@gsnPanelFigureStringsFontHeightF = 0.02
> >>> resPanel@gsnPanelFigureStringsPerimOn = True
> >>> resPanel@gsnPanelDebug = True
> >>> resPanel@gsnPanelTop = 0.9
> >>> resPanel@gsnPanelBottom = 0.3 ;move bottom up
> >>> from 0 to 0.4
> >>> resPanel@gsnPanelLeft = 0.0
> >>> resPanel@gsnPanelRight = 0.9
> >>> resPanel@gsnPanelDebug = True
> >>>
> >>> resPanel@lbOrientation = "Vertical"
> >>> resPanel@lbLabelAutoStride = False
> >>> resPanel@lbLabelStride = 5
> >>>
> >>> gsn_panel(wks,plot,(/2,2/),resPanel)
> >>>
> >>> ;************************************************
> >>> ; create last 2 of panel plot
> >>> ;************************************************
> >>> dresPanel = True
> >>>
> >>> dresPanel@gsnFrame = False
> >>> dresPanel@gsnSpreadColors = True
> >>> dresPanel@gsnMaximize = True ; use full
> >>> page
> >>> dresPanel@gsnPanelLabelBar = True ; add
> >>> color bar
> >>> dresPanel@gsnPanelTop = 0.3
> >>> dresPanel@gsnPanelBottom = 0.1
> >>> dresPanel@gsnPanelLeft = 0.0
> >>> dresPanel@gsnPanelRight = 0.9
> >>> dresPanel@gsnPanelFigureStrings = (/"e","f"/)
> >>> dresPanel@gsnPanelFigureStringsFontHeightF = 0.02
> >>> dresPanel@gsnPanelFigureStringsPerimOn = True
> >>> dresPanel@gsnPanelDebug = True
> >>>
> >>> dresPanel@lbOrientation = "Vertical"
> >>> dresPanel@lbLabelAutoStride = False
> >>> dresPanel@lbLabelStride = 10
> >>>
> >>> gsn_panel (wks,plot1,(/1,2/),dresPanel)
> >>> frame(wks)
> >>>
> >>>
> >>> Thanks,
> >>> Soumik
> >>>
> >>>
> >>> --
> >>> "Numbing the pain for a while will make it worse when you finally feel it.
> >>> ~J.K. Rowling"
> >>> *****************************************
> >>> Soumik Basu
> >>> Graduate Student, Research Assistant
> >>> International Arctic Research Center
> >>> University of Alaska Fairbanks
> >>> Fairbanks,Alaska,USA
> >>> *****************************************
> >>> <eady_eke_var_diff_spring_winter.000001.png>
> >>> _______________________________________________
> >>> ncl-talk mailing list
> >>> List instructions, subscriber options, unsubscribe:
> >>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >
> > --
> > "Numbing the pain for a while will make it worse when you finally feel it.
> > ~J.K. Rowling"
> > *****************************************
> > Soumik Basu
> > Graduate Student, Research Assistant
> > International Arctic Research Center
> > University of Alaska Fairbanks
> > Fairbanks,Alaska,USA
> > *****************************************
> > <eady_eke_var_diff_spring_winter.000001.png>
> > _______________________________________________
> > ncl-talk mailing list
> > List instructions, subscriber options, unsubscribe:
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
> ------------------------------
>
> Message: 3
> Date: Sat, 7 Jul 2012 10:53:13 +0800
> From: 163 <guoyipeng163@163.com>
> Subject: eof function
> To: ncl-talk <ncl-talk@ucar.edu>
> Message-ID: <201207071053134695861@163.com>
> Content-Type: text/plain; charset="gb2312"
>
> Hi all
> When using NCL EOF functions, I find that all functions use the anomalies. It
> means that I have no option to use the original data or the standardized data
> and can only use the anomalies. If there are some options in the function that
> I don't know, would you please tell me?
>
> Best Regards
> Guo Yipeng
>
>
>
>
> 163
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20120707/c130f6d0/attachment.html
>
> ------------------------------
>
> Message: 4
> Date: Sat, 7 Jul 2012 04:12:54 +0000
> From: Brandon Capasso <sblueknigt@gmail.com>
> Subject: Re: Stripping white space at the edges
> To: Noel Aloysius <noel.aloysius@gmail.com>
> Cc: NCL Talk <ncl-talk@ucar.edu>
> Message-ID:
> <CAFrH9e1Fp1k_Cih=mVOhbWwazE58Bwo7EEVYLeCq7f_hrX7BGw@mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> If you're using Linux, you can use imagemagick for this if need be:
>
> > convert infile -trim outfile
>
> you can then use -border to create a border of as much white space as you
> need. hope this helps .
>
> Brandon
> On Jul 6, 2012 12:34 PM, "Noel Aloysius" <noel.aloysius@gmail.com> wrote:
>
> > Hello NCL,
> >
> > I created a panel plot with 14 figures. I want to maximize the plot area
> > by stripping the white space at the edges. I am not sure which plot
> > resources do this.
> > The section of the code and the plot is attached.
> >
> > Thank you in advance for your help,
> > Noel
> >
> >
> > ;***********************************************************************************************
> > wks = gsn_open_wks("png", diro + "pr_zonalAvg")
> > plot = new(14, graphic)
> > gsn_define_colormap(wks, "precip3_16lev")
> >
> > res1 = True ; plot mods desired
> > res1@gsnDraw = False
> > res1@gsnFrame = False
> > res1@gsnMaximize = True
> > res1@gsnRightStringFontHeightF = 0.04
> > res1@cnFillOn = True ; turn on color
> > res1@gsnSpreadColors = True ; use full color table
> > res1@cnLinesOn = False ; no contour lines
> > res1@cnInfoLabelOn = False
> > res1@pmTickMarkDisplayMode = "Always" ; turn on fancy tickmarks
> >
> > res1@tmXBMode = "Explicit" ; label independently
> > res1@tmXBValues = ispan(0,11,1)
> > res1@tmXBLabels =
> > (/"J","F","M","A","M","J","J","A","S","O","N","D"/)
> >
> > res1@lbLabelBarOn = False
> > res1@cnLevelSelectionMode = "ManualLevels"
> > res1@cnLevelSpacingF = 20.
> > res1@cnMinLevelValF = 20.
> > res1@cnMaxLevelValF = 300.
> >
> > res1@tmYROn = False
> > res1@tmXBLabelFontHeightF = 0.025
> > res1@tmYLLabelFontHeightF = 0.025
> >
> > do i=0,11
> > res1@gsnRightString = gcm_names(i)
> > plot(i) = gsn_csm_lat_time(wks,
> > PzoneLat_7100(raw_bc|0,GCM|i,lat|:,month|:), res1)
> > end do
> >
> > res1@gsnRightString = "Obs"
> > plot(12) = gsn_csm_lat_time(wks, PzoneLat_7100obs(lat|:,month|:), res1)
> >
> > res1@gsnRightString = "bias corr. GCM avg."
> > plot(13) = gsn_csm_lat_time(wks, PzoneLat_7100_BCavg(lat|:,month|:),
> > res1)
> >
> > pres1 = True
> > pres1@gsnPanelLabelBar = True
> > pres1@gsnFrame = False
> > pres1@lbLabelFontHeightF = 0.01
> > pres1@lbOrientation = "horizontal"
> > pres1@lbLabelAutoStride = True
> > pres1@gsnPanelYWhiteSpacePercent = 1.5
> > pres1@gsnBoxMargin = 0.0
> > pres1@gsnPanelRowSpec = True ; tell panel what order
> > to plot
> > pres1@gsnPanelCenter = True ; centered plot
> > pres1@gsnMaximize = True
> >
> > gsn_panel(wks,plot,(/6,6,2/),pres1)
> > frame(wks)
> >
> > ;****************************************************************************************************
> >
> >
> > _______________________________________________
> > ncl-talk mailing list
> > 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/20120707/7c89885d/attachment.html
>
> ------------------------------
>
> Message: 5
> Date: Sat, 7 Jul 2012 09:00:19 +0000
> From: "Cannavina, Dominique" <d.cannavina11@imperial.ac.uk>
> Subject: Dimension mismatch on dim_sum_n the dimension
> relating to file input
> To: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu>
> Message-ID:
> <D7809B440AA90D499FD431B2E7289C551F338C7A@icexch-m1.ic.ac.uk>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Hi it seems that I was premature in my celebrations yesterday. The file I sent for runs of trcit below zero degress work okay .........except picking up the files. When I ran the script for various numbers of years (year per file) it was clear that my values were around the same regardless of how many years (or files) that I inputted.
>
>
>
> Then I remembered that I had changed the following;
>
>
>
> spell_tot_uk(:,nr,:,:) = dim_sum_n(cold_uk1(:,{spell(nr)},:,:,:),0)
>
>
>
> spell_tot_uk(0,nr,:,:) = dim_sum_n(cold_uk1(0,{spell(nr)},:,:,:),0) as I was trying to get around a dimension #0 mis-match. This is the variable that relates to file that has been inputted.
>
>
>
> So can I ask another NCL-Talk question on how to get around this please?
>
>
>
> The portion of the script it relates to is below and a copy of the full script and summary of the key variables are attached (join relates to number of files, spell the length of the cold run taking the tcrit from all ensembles at all times/steps)
>
>
>
> n1 = 1
>
> n2 = 60 ; ispan to specify all spells
>
> spell = ispan(n1,n2,1) ; eg n1=1, n2=20
>
> ;;spell = (/1, 1, 10/) ; manually specify
>
> nspell = dimsizes(spell)
>
> spell_tot_uk = new( (/njoin,nspell,nlat,mlon/), "float") ; preallocate
>
> dimspell = dimsizes(spell_tot_uk)
>
> njoin = dimspell(0)
>
> nspell = dimspell(1) ; same as "run_length"
>
> nlat = dimspell(2)
>
> mlon = dimspell(3)
>
> do nr=0,nspell-1
>
> spell_tot_uk(:,nr,:,:) = dim_sum_n(cold_uk1(:,{spell(nr)},:,:,:),0)
>
> end do
>
> spell_tot_uk!1 = "spell"
>
> spell_tot_uk&spell = spell
>
> copy_VarCoords(cold_uk1(0,0,0,:,:), spell_tot_uk(0,0,:,:)) ; trick for coords
>
> printVarSummary(spell_tot_uk)
>
>
>
> Any help greatly appreciated
>
>
>
> kind regards
>
>
>
> Dominique
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20120707/2d461e1f/attachment.html
> -------------- next part --------------
> A non-text attachment was scrubbed...
> Name: demeter_v1_models_check.ncl
> Type: application/octet-stream
> Size: 12992 bytes
> Desc: demeter_v1_models_check.ncl
> Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20120707/2d461e1f/attachment.obj
> -------------- next part --------------
> An embedded and charset-unspecified text was scrubbed...
> Name: Variables.o1827871
> Url: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20120707/2d461e1f/attachment.pl
>
> ------------------------------
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk@ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
> End of ncl-talk Digest, Vol 104, Issue 17
> *****************************************
>
>
> _______________________________________________
> 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 Mon Jul 9 09:52:24 2012
This archive was generated by hypermail 2.1.8 : Mon Jul 09 2012 - 10:45:32 MDT