Re: Overlay is not working

Hi Abhik,
There is still a major draw/frame issue with your coding. Without seeing
your code I am constrained in what I can recommend.

However, I (now) notice that your longitude ranges between your xy plot
and your xy timeseries match. You can try overlaying them as Mary
originally suggested. Modifying the coding I sent earlier:

do j =0,5
.....
res = True
res@gsnDraw = False
res@gsnFrame = False
.......
plotz(j) = gsn_csm_contour(wks,z,res)
xyres =True
xyres@gsnDraw = False
xyres@gsnFrame = False
........
plotx(j) = gsn_csm_xy(wks,lon,x,xyres)
overlay(plotz(j),plotx(j))
delete(ntim) ; do not delete anything else other than ntim, x, and z
delete(x)
delete(z) j = j+1
end do
resP = True
resP@gsnMaximize = True
resP@gsnFrame = False
gsn_panel(wks,plotz,(/6,1/),resP) ; create panel
txres = True
txres@txAngleF = 90. ; Rotate string clockwise
txres@txFontHeightF = 0.01
gsn_text_ndc(wks," Pressure (hPa) ",0.34,0.68,txres)
frame(wks)
end

In the end, overlaying can take a lot of trial and error on your part.
You might have to remove subtitles, change font heights, check that your
draw and frame statements are correct, etc. But the code basis I include
above should get you on the right track, especially in terms of your
draw and frame issues.
Adam

On 01/07/2013 09:51 PM, Abhik Santra wrote:
> Hi Adam,
>
> It didn't work properly.
> Please check the attached plot.
>
> Regards,
> Abhik
>
>
> ----- Original Message -----
> From: "Adam Phillips" <asphilli@ucar.edu>
> To: "Abhik Santra" <abhiksantra@tropmet.res.in>
> Cc: "ncl-talk" <ncl-talk@ucar.edu>
> Sent: Saturday, January 5, 2013 11:16:36 AM
> Subject: Re: Overlay is not working
>
>
> Hi Abhik,
> I think it might be best to simplify your script a bit. But to answer your question about why you had to use a ps file instead of a eps file: eps files cannot have more that one frame in them.
> The maximize_output documentation states that the plots will be drawn and the frame will be advanced, which is why you are getting 6 frames. ( http://www.ncl.ucar.edu/Document/Graphics/Interfaces/maximize_output.shtml )
>
>
> If you wish to use maximize_output I would only call it once after the two calls to gsn_panel.
>
>
> Try this, setting vpWidthF and vpHeightF =0.6 in both the res and xyres resource lists:
> (... = leave the coding the same as you have it)
> do j =0,5
> .....
> res = True
> .......
> plotz(j) = gsn_csm_contour(wks,z,res)
> xyres =True
> ........
> plotx(j) = gsn_csm_xy(wks,lon,x,xyres)
> delete(ntim)
> delete(x)
> delete(z) j = j+1
> end do
> resP = True
> .... ; do not set gsnMaximize = True in the panel resource list
> gsn_panel(wks,plotz,(/6,1/),resP) ; create panel
> gsn_panel(wks,plotx,(/6,1/),resP) ; create panel
> txres = True
> txres@txAngleF = 90. ; Rotate string clockwise
> txres@txFontHeightF = 0.01
> gsn_text_ndc(wks," Pressure (hPa) ",0.34,0.68,txres)
> maximize_output(wks,False)
>
> end
>
>
> Let ncl-talk know if that doesn't fix things...
> Adam
>
>
> On Jan 4, 2013, at 4:46, Abhik Santra < abhiksantra@tropmet.res.in > wrote:
>
>
>
>
>
> Hi Adam,
>
> My plot doesn't involve any map. It contains the latitude averaged vertical distribution over the longitude belt.
> So the mpFreeShapeMode = "FreeAspect" is not a valid resource for gsn_csm_contour.
> I have made some progress by introducing "maximize_output(wks,res)" inside the do loop.
> It is plotting the output from gsn_csm_xy over the output from gsn_csm_contour.
> But the panels are now not working and the output file type has been changed from eps to ps.
> Presently, 6 plots appear on 6 different pages without having the common labelbar.
> But I want to make a panel plot for all of them.
> How can I solve this issue?
>
> I am attaching the output ps file and my present scrips is as follows:
>
>
>
>
>
> begin
>
>
>
>
>
>
>
>
>
> plotx = new(6,graphic)
>
>
>
>
> plotz = new(6,graphic)
>
>
>
>
> wks = gsn_open_wks("ps","test")
>
>
>
>
> gsn_define_colormap(wks,"BlWhRe")
>
>
>
>
>
>
>
>
>
>
>
>
>
>
> do j=0, 5
>
>
>
>
>
>
>
>
>
> f1 = addfile ("XZ"+j+".nc" , "r") : res: 1.25Dx1.25D
>
>
>
>
> f2 = addfile ("X"+j+".nc" , "r") ; res: 1Dx1D
>
>
>
>
>
>
>
>
>
> z = f1->zz(0,:,:)
>
>
>
>
> x = f2->xx(0,:)
>
>
>
>
>
>
>
>
>
> lon = xx&lon
>
>
>
>
>
>
>
>
>
>
>
>
>
>
> res = True
>
>
>
>
> res@cnFillOn = True ; turn on color fill
>
>
>
>
> res@gsnDraw = False ; do not draw picture
>
>
>
>
> res@gsnFrame = False ; do not advance frame
>
>
>
>
> res@lbLabelBarOn = False ; no individual label bars
>
>
>
>
> res@gsnSpreadColors = True ; use full colormap
>
>
>
>
> res@cnLinesOn = True ; no contour lines
>
>
>
>
> res@cnLineThicknessF = 0.6
>
>
>
>
> res@cnLineColor = "gray"
>
>
>
>
> res@cnInfoLabelOn = False ; turn off cn info label
>
>
>
>
> res@cnLevelSelectionMode = "ManualLevels" ; set manual contour levels
>
>
>
>
> res@cnMinLevelValF = -0.3 ; set min contour level
>
>
>
>
> res@cnMaxLevelValF = 0.3 ; set max contour level
>
>
>
>
> res@cnLevelSpacingF = 0.05 ; set contour spacing
>
>
>
>
> res@cnLineLabelsOn = False
>
>
>
>
> res@tmXBLabelFontHeightF = 0.0
>
>
>
>
> res@tmYLLabelFontHeightF = 0.018
>
>
>
>
> res@tmYROn = False
>
>
>
>
> res@tmYRMinorOn = False
>
>
>
>
> res@trYReverse = True
>
>
>
>
> res@vpHeightF = 0.2 ; default 0.6
>
>
>
>
> res@vpWidthF = 0.8 ; default 0.6
>
>
>
>
>
>
>
>
>
> plotz(j) = gsn_csm_contour(wks,z,res) ; !! This is only coming in the panel !!
>
>
>
>
>
>
>
>
>
> xyres = True
>
>
>
>
> xyres@gsnDraw = False ; don't draw
>
>
>
>
> xyres@gsnFrame = False
>
>
>
>
> xyres@tmXTOn = False
>
>
>
>
> xyres@tmYLOn = False
>
>
>
>
> xyres@tmYROn = True
>
>
>
>
> xyres@tmYRLabelsOn = True
>
>
>
>
> xyres@trYMinF = -5. ; set Y-axis max and mins
>
>
>
>
> xyres@trYMaxF = 5.
>
>
>
>
> xyres@gsnYRefLine = 0. ; Y-value for ref. line
>
>
>
>
> xyres@xyLineThicknesses = (/3./) ; Define line thicknesses
>
>
>
>
> xyres@xyLineColors = (/"black"/) ; Define line color
>
>
>
>
> xyres@trXMaxF = 360.
>
>
>
>
>
>
>
>
>
> plotx(j) = gsn_csm_xy(wks,lon,x,xyres)
>
>
>
>
>
>
>
>
>
> delete(res)
>
>
>
>
> delete(xyres)
>
>
>
>
> delete(ntim)
>
>
>
>
> delete(x)
>
>
>
>
> delete(z)
>
>
>
>
> psres = True
>
>
>
>
> maximize_output(wks,psres)
>
>
>
>
>
>
>
>
>
> j=j+1
>
>
>
>
>
>
>
>
>
> end do
>
>
>
>
> ;*******************************************
>
>
>
>
> ; create panel resources
>
>
>
>
> ;*******************************************
>
>
>
>
> resP = True ; mods desired
>
>
>
>
> resP@gsnFrame = False ; save panel until all ready
>
>
>
>
> ;resP@gsnMaximize = True ; Maximize in frame
>
>
>
>
> resP@gsnPanelBottom = 0.35
>
>
>
>
> resP@gsnPanelLabelBar = True ; commom label bar
>
>
>
>
> resP@lbLabelStride = 2 ; every 2nd label
>
>
>
>
> resP@lbLabelFontHeightF = 0.01
>
>
>
>
> resP@pmLabelBarHeightF = 0.06
>
>
>
>
> resP@gsnPanelFigureStringsFontHeightF = 0.007
>
>
>
>
> resP@pmLabelBarOrthogonalPosF = 0.01
>
>
>
>
> resP@gsnAttachPlotsXAxis = True
>
>
>
>
> resP@gsnScale = True
>
>
>
>
>
>
>
>
>
> gsn_panel(wks,plotz,(/6,1/),resP) ; create panels
>
>
>
>
> txres = True
> txres@txAngleF = 90. ; Rotate string clockwise
> txres@txFontHeightF = 0.01
> gsn_text_ndc(wks," Pressure (hPa) ",0.34,0.68,txres)
>
>
>
>
> frame(wks)
>
>
>
>
> end
>
> Regards,
> Abhik
>
>
>
> ----- Original Message -----
> From: "Adam Phillips" < asphilli@ucar.edu >
> To: ncl-talk@ucar.edu
> Sent: Thursday, January 3, 2013 10:56:45 PM
> Subject: Re: Overlay is not working
>
> Hi Abhik,
> Without having your script and data it is hard to see what could
> possibly be wrong. One guess would be that as your first panel plots
> involve maps, that NCL will not resize the maps so as to not alter the
> aspect ratio of the maps without you explicitly telling it to do so,
> regardless of whether you have vpWidthF and vpHeightF set. See example
> #3 here:
> http://www.ncl.ucar.edu/Applications/resize.shtml#ex3
> (setting mpFreeShapeMode = "FreeAspect" will tell NCL to alter the
> aspect ratio)
>
> Note though, that setting the above can seriously distort your map. If
> you don't mind doing that then setting mpFreeShapeMode is the way to go.
>
> Finally, hopefully you saw Mary Haley's correction to my reply yesterday
> afternoon. You can use overlay to place a XY plot on a map plot,
> provided that the X and Y axis ranges match between the two plots. If
> the ranges do no match though, going with my method #2 is the way to go..
> Adam
>
> On 01/02/2013 11:09 PM, Abhik Santra wrote:
>
>
> Hi Adam,
>
>
>
>
>
> Thank you for your suggestions.
>
>
> I followed your 2nd method and able to overlap the 2 set of plots.
>
>
> But I am unable to control the viewport and as a result the plots are not exactly overlapped.
>
>
> I used vpWidthF and vpHeightF as 0.8 and 0.2 respectively for both the panels.
>
>
> But the plots at 2nd panel are appeared larger than that of the 1st panel.
>
>
> The gsnPanelDebug for 1st panel are as follows:
>
>
> plot #0
>
>
> (0) new x,y = 0.381421,0.993689
>
>
> (0) orig wdt,hgt = 0.8,0.2
>
>
> (0) new wdt,hgt = 0.265712,0.0664279
>
>
> (0) -------Panel viewport values for each plot-------
>
>
> (0) plot #1
>
>
> (0) new x,y = 0.381421,0.914572
>
>
> (0) orig wdt,hgt = 0.8,0.2
>
>
> (0) new wdt,hgt = 0.265712,0.0664279
>
>
> (0) -------Panel viewport values for each plot-------
>
>
> (0) plot #2
>
>
> (0) new x,y = 0.381421,0.835455
>
>
> (0) orig wdt,hgt = 0.8,0.2
>
>
> (0) new wdt,hgt = 0.265712,0.0664279
>
>
> (0) -------Panel viewport values for each plot-------
>
>
> (0) plot #3
>
>
> (0) new x,y = 0.381421,0.756338
>
>
> (0) orig wdt,hgt = 0.8,0.2
>
>
> (0) new wdt,hgt = 0.265712,0.0664279
>
>
> (0) -------Panel viewport values for each plot-------
>
>
> (0) plot #4
>
>
> (0) new x,y = 0.381421,0.677221
>
>
> (0) orig wdt,hgt = 0.8,0.2
>
>
> (0) new wdt,hgt = 0.265712,0.0664279
>
>
> (0) -------Panel viewport values for each plot-------
>
>
> (0) plot #5
>
>
> (0) new x,y = 0.381421,0.598104
>
>
> (0) orig wdt,hgt = 0.8,0.2
>
>
> (0) new wdt,hgt = 0.265712,0.0664279
>
>
>
>
>
> Should I use vpXF and vpYF for the 2nd panel?
>
>
> Or should I use same vpXF and vpYF for both the panels?
>
>
> What will be their values for the plots at the 2nd panels?
>
>
> Please let me know the same.
>
>
>
>
>
> Regards,
>
>
> Abhik
>
>
>
>
>
>
>
>
>
>
>
>
>
>
> ----- Original Message -----
>
>
> From: "Adam Phillips" < asphilli@ucar.edu >
>
>
> To: ncl-talk@ucar.edu
>
>
> Sent: Wednesday, January 2, 2013 11:05:59 PM
>
>
> Subject: Re: Overlay is not working
>
>
>
>
>
> Hi Abhik,
>
>
> It looks like you are trying to overlay a XY plot onto a spatial map
>
>
> plot. I would venture to guess that NCL isn't sure how to overlay the 2,
>
>
> as the plot types are completely different.
>
>
>
>
>
> Someone might have a better alternative, but I see two possibilities:
>
>
> 1) Create all 6 map plots, call gsn_panel (but don't advance the frame
>
>
> by setting gsnFrame = False in the Panel resource list), create all 6 XY
>
>
> plots, call gsn_panel, manually advance the frame. You may have to
>
>
> adjust the vpWidthF/vpHeightF resources for the XY plots, and you will
>
>
> likely have to set gsnPanelXWhiteSpacePercent,
>
>
> gsnPanelYWhiteSpacePercent, gsnPanelBottom,gsnPanelTop,gsnPanelRight and
>
>
> gsnPanelLeft in the Panel resource lists to try to get things lined up
>
>
> correctly.
>
>
>
>
>
> 2) First create all 6 map plots, call gsn_panel (but don't advance the
>
>
> frame by setting gsnFrame = False in the Panel resource list). Draw each
>
>
> XY plot individually, exactly specifying the placement over each of your
>
>
> 6 map plots via the vpXY/vpYF resources. In this scenario you will want
>
>
> to set gsnFrame = False in the xy resource list, and at the end of the
>
>
> script call frame(wks)..
>
>
>
>
>
> I'd recommend trying method 2 first, as even though that method is more
>
>
> tedious than the 1st method, you have a higher chance of lining things
>
>
> up exactly as you want them..
>
>
>
>
>
> If you have further questions about this issue please respond to ncl-talk.
>
>
> Good luck,
>
>
> Adam
>
>
>
>
>
> On 01/02/2013 01:08 AM, Abhik Santra wrote:
>
>
>
>
> Hello ncl-talk,
>
>
>
>
>
>
>
>
>
> I am trying to make 6 overlay plots in a single panel from two different datasets having different resolutions (1.25x1.25 and 1x1 degrees).
>
>
>
>
> But my overlay plots didn't work appropriately as only one set of figures are appeared in the panel.
>
>
>
>
>
>
>
>
>
> I have the script as follow:
>
>
>
>
>
>
>
>
>
>
>
>
>
>
> begin
>
>
>
>
>
>
>
>
>
> plot = new(6,graphic)
>
>
>
>
> wks = gsn_open_wks("eps","test")
>
>
>
>
> gsn_define_colormap(wks,"BlWhRe")
>
>
>
>
>
>
>
>
>
>
>
>
>
>
> do j=0, 5
>
>
>
>
>
>
>
>
>
> f1 = addfile ("XZ"+j+".nc" , "r") : res: 1.25Dx1.25D
>
>
>
>
> f2 = addfile ("X"+j+".nc" , "r") ; res: 1Dx1D
>
>
>
>
>
>
>
>
>
> z = f1->zz(0,:,:)
>
>
>
>
> x = f2->xx(0,:)
>
>
>
>
>
>
>
>
>
> lon = xx&lon
>
>
>
>
>
>
>
>
>
>
>
>
>
>
> res = True
>
>
>
>
> res@cnFillOn = True ; turn on color fill
>
>
>
>
> res@gsnDraw = False ; do not draw picture
>
>
>
>
> res@gsnFrame = False ; do not advance frame
>
>
>
>
> res@lbLabelBarOn = False ; no individual label bars
>
>
>
>
> res@gsnSpreadColors = True ; use full colormap
>
>
>
>
> res@cnLinesOn = True ; no contour lines
>
>
>
>
> res@cnLineThicknessF = 0.6
>
>
>
>
> res@cnLineColor = "gray"
>
>
>
>
> res@cnInfoLabelOn = False ; turn off cn info label
>
>
>
>
> res@cnLevelSelectionMode = "ManualLevels" ; set manual contour levels
>
>
>
>
> res@cnMinLevelValF = -0.3 ; set min contour level
>
>
>
>
> res@cnMaxLevelValF = 0.3 ; set max contour level
>
>
>
>
> res@cnLevelSpacingF = 0.05 ; set contour spacing
>
>
>
>
> res@cnLineLabelsOn = False
>
>
>
>
> res@tmXBLabelFontHeightF = 0.0
>
>
>
>
> res@tmYLLabelFontHeightF = 0.018
>
>
>
>
> res@tmYROn = False
>
>
>
>
> res@tmYRMinorOn = False
>
>
>
>
> res@trYReverse = True
>
>
>
>
> res@vpHeightF = 0.2 ; default 0.6
>
>
>
>
> res@vpWidthF = 0.8 ; default 0.6
>
>
>
>
>
>
>
>
>
> plotz = gsn_csm_contour(wks,z,res) ; !! This is only coming in the panel !!
>
>
>
>
>
>
>
>
>
> xyres = True
>
>
>
>
> xyres@gsnDraw = False ; don't draw
>
>
>
>
> xyres@gsnFrame = False
>
>
>
>
> xyres@tmXTOn = False
>
>
>
>
> xyres@tmYLOn = False
>
>
>
>
> xyres@tmYROn = True
>
>
>
>
> xyres@tmYRLabelsOn = True
>
>
>
>
> xyres@trYMinF = -5. ; set Y-axis max and mins
>
>
>
>
> xyres@trYMaxF = 5.
>
>
>
>
> xyres@gsnYRefLine = 0. ; Y-value for ref. line
>
>
>
>
> xyres@xyLineThicknesses = (/3./) ; Define line thicknesses
>
>
>
>
> xyres@xyLineColors = (/"black"/) ; Define line color
>
>
>
>
> xyres@trXMaxF = 360.
>
>
>
>
>
>
>
>
>
> plotx = gsn_csm_xy(wks,lon,x,xyres)
>
>
>
>
> overlay(plotz,plotx)
>
>
>
>
> plot(j) = plotz
>
>
>
>
>
>
>
>
>
> delete(plotz)
>
>
>
>
> delete(plotx)
>
>
>
>
> delete(res)
>
>
>
>
> delete(xyres)
>
>
>
>
> delete(ntim)
>
>
>
>
> delete(x)
>
>
>
>
> delete(z)
>
>
>
>
>
>
>
>
>
> j=j+1
>
>
>
>
>
>
>
>
>
> end do
>
>
>
>
> ;*******************************************
>
>
>
>
> ; create panel resources
>
>
>
>
> ;*******************************************
>
>
>
>
> resP = True ; mods desired
>
>
>
>
> resP@gsnFrame = False ; save panel until all ready
>
>
>
>
> ;resP@gsnMaximize = True ; Maximize in frame
>
>
>
>
> resP@gsnPanelBottom = 0.35
>
>
>
>
> resP@gsnPanelLabelBar = True ; commom label bar
>
>
>
>
> resP@lbLabelStride = 2 ; every 2nd label
>
>
>
>
> resP@lbLabelFontHeightF = 0.01
>
>
>
>
> resP@pmLabelBarHeightF = 0.06
>
>
>
>
> resP@gsnPanelFigureStringsFontHeightF = 0.007
>
>
>
>
> resP@pmLabelBarOrthogonalPosF = 0.01
>
>
>
>
> resP@gsnAttachPlotsXAxis = True
>
>
>
>
> resP@gsnScale = True
>
>
>
>
>
>
>
>
>
> gsn_panel(wks,plot,(/6,1/),resP) ; create panels
>
>
>
>
>
>
>
>
>
> frame(wks)
>
>
>
>
> end
>
>
>
>
>
>
>
>
>
> Can anyone suggest me the way to solve this problem?
>
>
>
>
> Thanks in advance and wish you a happy new year.
>
>
>
>
>
>
>
>
>
> Regards,
>
>
>
>
> Abhik
>
>
>
>
> ______________________________
>
>
>
>
> Abhik Santra
>
>
>
>
> CSIR Research Fellow,
>
>
>
>
> Indian Institute of Tropical Meteorology,
>
>
>
>
> Pune - 411008.
>
>
>
>
> India
>
>
>
>
> _______________________________________________
>
>
>
>
> ncl-talk mailing list
>
>
>
>
> List instructions, subscriber options, unsubscribe:
>
>
>
>
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>

-- 
______________________________________________________________
Adam Phillips                                asphilli@ucar.edu
NCAR/Climate and Global Dynamics Division       (303) 497-1726
P.O. Box 3000				
Boulder, CO 80307-3000    http://www.cgd.ucar.edu/cas/asphilli
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk