[ncl-talk] Plotting overlay plot with significance values: SOLVED
Sri Nandini
snandini at marum.de
Fri Feb 2 04:11:30 MST 2018
Thank you
That works well now.
On Feb 1, 2018 5:46:29 PM, Mary Haley wrote:
> The reason you are getting that one plot on top of your panel plot is because you called gsn_panel first, which draws the paneled plots, and then you called:
>
> draw (plots)
> frame(wks)
> > which drew each plot again on top of each other, each in its original large size. You don't need those two lines, because gsn_panel is already drawing the four plots for you.
>
> > --Mary
>
>
>
> On Thu, Feb 1, 2018 at 2:57 AM, Sri Nandini via ncl-talk > <> ncl-talk at ucar.edu> >> wrote:
> > Hello
> >
> > I tried with and without
> > plots(mm) = gsn_csm_contour_map_polar(wks,> > plots(mm)> > ,res)
> > as well as removing the (mm).
> >
> > Theres no error:
> >
> > However the same plot is being plotted (5 instead of 4).
> >
> >
> >
> >
> >
> >
> > On Jan 31, 2018 10:56:28 PM, Mary Haley wrote:
> > > The error message is telling you what the problem is:
> > >
> > > Argument type mismatch on argument (1) of (gsn_csm_contour_map_polar) can not coerce
> > >
> > > > > > Argument #1 refers to the second argument in the gsn_csm_contour_map_polar call, because argument count starts at 0.
> > >
> > > > > > I looked at your > > > gsn_csm_contour_map_polar calls and found this one:
> > >
> > > > > > plots(mm) = gsn_csm_contour_map_polar(wks,> > > plots(mm)> > > ,res)
> > >
> > > > > > You are trying to pass plots(mm) to > > > gsn_csm_contour_map_polar which doesn't make sense.
> > >
> > > > > > You have already called gsn_csm_contour_map_polar earlier in your script and put the plots in plots(mm) via these calls:
> > >
> > > > > > plots(0) = gsn_csm_contour_map_polar(wks,> > > T1diff,res)
> > > plots(1) = gsn_csm_contour_map_polar(wks,> > > T2diff,res)
> > > plots(2) = gsn_csm_contour_map_polar(wks,> > > T3diff,res)
> > > > > > plots(3) = gsn_csm_contour_map_polar(wks,> > > T4diff,res)
> > >
> > > > > > so I'm not sure why you are trying to call gsn_csm_contour_map_polar again. Try removing or commenting out the offending line and see what happens.
> > >
> > > > > > ; > > > plots(mm) = gsn_csm_contour_map_polar(wks,> > > plots(mm),res)
> > >
> > > > > > --Mary
> > >
> > >
> > >
> > >
> > >
> > >
> > >
> > > On Wed, Jan 31, 2018 at 8:47 AM, Sri Nandini via ncl-talk > > > <> > > ncl-talk at ucar.edu> > > >> > > wrote:
> > > > Hello
> > > >
> > > > I am trying to plot seasonal (DJF, MAM, JJA, SON) temperature values with thier ttest values (all already calculated and saved in netcdf file).
> > > >
> > > > I attach my plotting script below:
> > > > And my error is:
> > > > Argument type mismatch on argument (1) of (gsn_csm_contour_map_polar) can not coerce
> > > >
> > > > Would appreciate any advice on why the plot turns out like this, see attached.
> > > >
> > > > ; Reading all seasonal differences and ttest, then overlaying on each plot with the caspian basin shapefile to show the 95% or 99% significance levels
> > > >
> > > >
> > > > latS = 30. ; N. Hemisphere
> > > > latN = 90.
> > > >
> > > > ADD_SHAPEFILE_OUTLINES = True
> > > > sname = "Caspian_Basin.shp"
> > > >
> > > > ; ==============================================================
> > > > ; Open the files:Read all the seasonal differences and ttest values
> > > > ; ==============================================================
> > > > f = addfile ("T2M_BC_DJF_Diff.nc", "r")
> > > > T1diff = f->T1diff
> > > >
> > > > f2 = addfile("T2M_BC_MAM_Diff.nc","r")
> > > > T2diff = f2->T2diff
> > > >
> > > > f3 = addfile("T2M_BC_JJA_Diff.nc","r")
> > > > T3diff = f3->T3diff
> > > >
> > > > f4 = addfile("T2M_BC_SON_Diff.nc","r")
> > > > T4diff = f4->T4diff
> > > >
> > > > ;==================================================================
> > > > f5 = addfile("T2M_BC_DJF_ttest.nc","r")
> > > > probt1 = f5->probt1
> > > >
> > > > f6 = addfile("T2M_BC_MAM_ttest.nc","r")
> > > > probt2 = f6->probt2
> > > >
> > > > f7 = addfile("T2M_BC_JJA_ttest.nc","r")
> > > > probt3 = f7->probt3
> > > >
> > > > f8 = addfile("T2M_BC_SON_ttest.nc","r")
> > > > probt4 = f8->probt4
> > > >
> > > > ;=====================================================================Start the graphics
> > > > wks = gsn_open_wks("png","T2M_BC_seasons_diff_ttest_")
> > > >
> > > > res = True
> > > > res at gsnDraw = False
> > > > res at gsnFrame = False
> > > > res at cnFillPalette = "BlWhRe" ; select a color map with white in the
> > > > res at gsnPolar = "NH"
> > > > res at mpGridAndLimbOn = "False" ; turn off lat/lon lines
> > > > res at mpFillOn = False ; turn off map fill
> > > >
> > > > res at mpMinLatF = latS ; zoom in on map
> > > > res at mpMaxLatF = latN
> > > >
> > > > res at cnLevelSelectionMode = "ManualLevels" ; manual set levels
> > > > res at cnMinLevelValF =-3
> > > > res at cnMaxLevelValF = 3
> > > > res at cnLevelSpacingF = 0.5 ; contour levels
> > > >
> > > > res at cnFillOn = True
> > > > res at cnLinesOn = False
> > > > res at cnLineLabelsOn = False
> > > > res at cnFillDrawOrder = "Predraw" ; fill and lines before map
> > > > res at mpDataBaseVersion = "MediumRes" ; use GMT coastline
> > > >
> > > > res at mpGridAndLimbOn = False ; set to true if you want grid lines to show
> > > > res at lbLabelBarOn = False ; turn off individual label bars
> > > > ;res at cnRasterModeOn = True ; Raster mode shows grid cells
> > > > res at gsnRightString = "" ; Turn off subtitles
> > > > res at gsnLeftString = ""
> > > >
> > > > ;==================================Create arrays to hold series of plots
> > > > plots = new(4,graphic)
> > > >
> > > > plots(0) = gsn_csm_contour_map_polar(wks,T1diff,res)
> > > >
> > > > plots(1) = gsn_csm_contour_map_polar(wks,T2diff,res)
> > > >
> > > > plots(2) = gsn_csm_contour_map_polar(wks,T3diff,res)
> > > >
> > > > plots(3) = gsn_csm_contour_map_polar(wks,T4diff,res)
> > > > printVarSummary(plots)
> > > >
> > > > ;=====================================================================
> > > > ;Second plot resources and overlay onto first
> > > > ;=====================================================================
> > > > res2 = True ; res2 probability plots
> > > > res2 at gsnDraw = False ; Do not draw plot
> > > > res2 at gsnFrame = False ; Do not advance frome
> > > >
> > > > res2 at cnLevelSelectionMode = "ExplicitLevels" ; set explicit cnlev
> > > > res2 at cnLevels = (/99./) ; only have 1 contour level
> > > > res2 at cnInfoLabelOn = False
> > > > res2 at cnLinesOn = False ; do not draw contour lines
> > > > res2 at cnLineLabelsOn = False ; do not draw contour labels
> > > > res2 at cnFillScaleF = 0.6 ; add extra density
> > > >
> > > > res2 at gsnRightString = "" ; Turn off subtitles
> > > > res2 at gsnLeftString = ""
> > > > res at gsnAddCyclic = True ; add cyclic point
> > > >
> > > > ;========================================================================
> > > > opt = True ; set up parameters for pattern fill
> > > > opt at gsnShadeFillType = "pattern" ; specify pattern fill
> > > > opt at gsnShadeHigh = 17 ; stipple pattern
> > > > opt at gsnShadeDotSizeF = 1 ; make dots larger
> > > >
> > > > ;==================================Create arrays to hold series of ttests or make a loop
> > > > plot2 = new(4,graphic)
> > > >
> > > > plot2(0) = gsn_csm_contour(wks,probt1,res2)
> > > >
> > > > plot2(1) = gsn_csm_contour(wks,probt2,res2)
> > > >
> > > > plot2(2) = gsn_csm_contour(wks,probt3,res2)
> > > >
> > > > plot2(3) = gsn_csm_contour(wks,probt4,res2)
> > > > printVarSummary(plot2)
> > > > ;==========================================================================cycle through seasons
> > > > season = (/"DJF","MAM","JJA","SON"/)
> > > > nof=4
> > > >
> > > > do mm=0,nof-1
> > > > res at gsnCenterString = season(mm)
> > > >
> > > > plots(mm) = gsn_csm_contour_map_polar(wks,plots(mm),res)
> > > > plot2(mm) = gsn_contour_shade(plot2(mm),1,99., opt) ; stipple all areas >= 99% contour
> > > >
> > > > overlay (plots (mm), plot2 (mm))
> > > > end do
> > > >
> > > > printVarSummary(plots)
> > > > printVarSummary(plot2)
> > > >
> > > > ;=====================================
> > > > ;Attach the polylines and shape
> > > > ;=====================================
> > > > pres = True
> > > > pres at gsLineColor = "black"
> > > > pres at gsLineThicknessF = 2.0 ; 2x thickness
> > > > shp1 = gsn_add_shapefile_polylines(wks,plots,sname,pres)
> > > > printVarSummary(shp1)
> > > >
> > > > ;==================================
> > > > ; panel resources for the 4 plots
> > > > ;==================================
> > > > pres1 = True ; modify the panel plot
> > > > pres1 at gsnPanelLabelBar = True ; common label bar
> > > > pres1 at gsnFrame = False ; don't advance frame yet
> > > > pres1 at lbOrientation = "Vertical" ; vertical label bar
> > > > pres1 at gsnPanelMainString = "T2m BigCaspian-NoCaspian Seasonal Difference"
> > > > pres1 at lbBoxEndCapStyle = "TriangleBothEnds" ; Added in NCL V6.4.0
> > > > pres1 at gsnMaximize = True ; large format
> > > >
> > > > gsn_panel(wks,plots,(/2,2/),pres1)
> > > >
> > > > draw (plots)
> > > > frame(wks)
> > > >
> > > >
> > > >
> > > > _______________________________________________
> > > >
ncl-talk mailing list
> > > > ncl-talk at ucar.edu
> > > >
List instructions, subscriber options, unsubscribe:
> > > > 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/20180202/0dc69eae/attachment.html>
More information about the ncl-talk
mailing list