Re: Image resize problem

Sahidul,

I forgot about this, but the wrf_map_overlays function removes the contour plot from the map plot before it exits the function.

I'm not sure why it was coded this way, but I think you can work around this by setting:

    pltres@PanelPlot = True
    pltres@gsnDraw = True
    pltres@gsnFrame = True

--Mary

On Jun 20, 2010, at 10:01 AM, Sahidul wrote:

> Dear Mary,
>
> According to advise, I tried to add text function and it is working fine. BUT It is not plotting the image; only "white" image is creating (contour plot of rainfall is not showing).
> So my purpose is not solved. I have shown my script below. Kindly check it.
>
> =================Start scripting======================
> 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"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>
> begin
>
> datadir = "/casvol5/basantas/cas/WRF3-1/output-wrf/2010061600/"
> date1="2010-06-16_00"
> date2="2010-06-16_23"
> date3="2010-06-17_23"
> date4="2010-06-18_23"
>
> FILES1 = systemfunc ("ls -1 " + datadir + "wrfout_d02_"+date1)
> FILES2 = systemfunc ("ls -1 " + datadir + "wrfout_d02_"+date2)
> FILES3 = systemfunc ("ls -1 " + datadir + "wrfout_d02_"+date3)
> FILES4 = systemfunc ("ls -1 " + datadir + "wrfout_d02_"+date4)
>
> wrf_out = addfile(FILES1+".nc","r")
>
> res = True
> opts = res
> opts@InitTime = False
> opts@Footer = False
>
> ; Total Rainfall calculation
> ;***********************************************************************
> rainc1 = wrf_user_getvar(wrf_out,"RAINC",0)
> rainnc1 = wrf_user_getvar(wrf_out,"RAINNC",0)
>
> wrf_out2 = addfile(FILES2+".nc","r")
>
> ; Total Rainfall calculation
> ;***********************************************************************
> rainc2 = wrf_user_getvar(wrf_out2,"RAINC",0)
> rainnc2 = wrf_user_getvar(wrf_out2,"RAINNC",0)
>
> rain_tot2=((rainc2+rainnc2)-(rainnc1+rainc1))/10.
> ;***********************************************************************
>
> wrf_out3 = addfile(FILES3+".nc","r")
>
> ; Total Rainfall calculation
> ;***********************************************************************
> rainc3 = wrf_user_getvar(wrf_out3,"RAINC",0)
> rainnc3 = wrf_user_getvar(wrf_out3,"RAINNC",0)
>
> rain_tot3=((rainc3+rainnc3)-(rainnc2+rainc2))/10.
> ;***********************************************************************
>
> wrf_out4 = addfile(FILES4+".nc","r")
>
> ; Total Rainfall calculation
> ;***********************************************************************
> rainc4 = wrf_user_getvar(wrf_out4,"RAINC",0)
> rainnc4 = wrf_user_getvar(wrf_out4,"RAINNC",0)
>
> rain_tot4=((rainc4+rainnc4)-(rainnc3+rainc3))/10.
>
> rain_tot2@description = ""
> rain_tot3@description = ""
> rain_tot4@description = ""
> ;***********************************************************************
> ;Open Workstation ps,eps,X11,pdf,ncgm
> ;***********************************************************************
> wks2 = gsn_open_wks("ps","precip_surf_day1")
> wks3 = gsn_open_wks("ps","precip_surf_day2")
> wks4 = gsn_open_wks("ps","precip_surf_day3")
> ;gsn_define_colormap(wks,"gui_default")
> res@InitTime = False
> res@Footer = False
> res@gsnMaximize = True
> res@MainTitle = "Total Precipitation"
>
> opts@cnLinesOn = False
> opts@cnFillOn = True
> opts@cnLevelSelectionMode = "ExplicitLevels"
> ;opts@UnitLabel = "cm"
> ; opts@description = "Precipitation Rate"
>
> opts@cnLevels = (/ 1., 2., 5., 10., 15., 20., 25., 30., 35., 40., 50./)
> opts@cnFillColors = (/"white","seagreen1","green4","greenyellow","yellow","goldenrod1","goldenrod2","tomato","orangered",\
> "red","red2","red4"/)
> opts@cnInfoLabelOn = False
> opts@cnConstFLabelOn = False
> res@gsnMaximize = True ; must include w/ Paper Orientation
> res@gsnPanelBottom = 0.2 ; save space at bottom
>
> txres = True ; text mods desired
>
> pltres = True
> pltres@gsnDraw = False ; don't draw plot
> pltres@gsnFrame = False ; don't advance frame
> pltres@gsnMaximize = True ; must include w/ Paper Orientation
>
> mpres = True
> ; mpres@mpOutlineBoundarySets = "National"
> mpres@mpNationalLineThicknessF = 1.5
> mpres@mpNationalLineColor = "Black"
> mpres@mpGeophysicalLineColor = "Black"
> mpres@mpGeophysicalLineThicknessF = 1.5
> ;***********************************************************************
> contour_tot2 = wrf_contour(wrf_out,wks2,rain_tot2,opts)
> plot = wrf_map_overlays(wrf_out,wks2,contour_tot2,pltres,mpres)
>
> ;*********************plot1**************************************************
> txres@txFontHeightF = 0.013 ; font smaller. default big
> txres@txFont = "text-symbols"
> xlon1 = 53
> xlat = 35
> dumid1 = gsn_add_text(wks2,plot,"S",xlon1,xlat,txres) ; need to replace xlon1 and xlat with appropriate lon/lat values
> ;; gsn_text_ndc(wks2,"S",0.12,.19,txres)
> ; page coordinates are normalized. go from 0->1 not 0->8" etc.
> txres@txFontHeightF = 0.010 ; font smaller. default big
> txres@txFont = "courier"
> ;;txres@txFont = "courier-bold"
> xlon2 = 56
> xlat = 35
> dumid2 = gsn_add_text(wks2,plot,"CDAC",xlon2,xlat,txres) ; need to replace xlon2 and xlat with appropriate lon/lat values
> ;; gsn_text_ndc(wks2,"CDAC",0.17,.19,txres)
> ;;txres@txFontHeightF = 0.01 ; font smaller. default big
> ;;txres@txFont = "courier"
> ;;xlon3 = 57
> ;;xlat1 = 35.5
> ;;dumid3 = gsn_add_text(wks2,plot,"(Continent Boundary may differ from Political boundary) ",xlon3,xlat1,txres) ; need to replace xlon3 and xlat with appropriate lon/lat values
> ;; gsn_text_ndc(wks2,"(Continent Boundary may differ from Political boundary) ",0.6,.19,txres)
>
> draw(plot)
> frame(wks2) ; now advance frame
> end
> =================end script======================
>
> Please help me.
>
> Thanking you,
>
>
> On Thu, Jun 17, 2010 at 8:56 PM, Mary Haley <haley@ucar.edu> wrote:
>
> Sahidul,
>
> Are the text strings showing up on the very last page of your PS file by any chance?
>
> I think what may be happening is that when you call wrf_map_overlays, it is calling a frame advance for you.
> This means that anything you draw after this call will be drawn on a new page.
>
> If this is what's happening, then before you call this function, try setting:
>
> pltres@gsnFrame = False
>
> You also mentioned a problem with resizing the image. I'm not sure what you mean.
>
> Note that if you want the text to appear on the map, you'd be better off using "gsn_add_text",
> because this allows you to attach the text to the plot. This means if you resize the plot, the text
> will be resized automatically.
>
> In order to do this, however, you need to use data coordinates of your plot, rather than NDC
> coordinates. This would be lat/lon coordinates, since "plot" is a map.
>
> Also, the text won't get drawn until the plot gets drawn.
>
> Here's what your code changes might look like:
>
> . . .
> pltres = True
> . . .
> pltres@gsnDraw = False ; don't draw plot
> pltres@gsnFrame = False ; don't advance frame
>
> . . .
> contour_tot2 = wrf_contour(wrf_out,wks2,rain_tot2,opts)
> plot = wrf_map_overlays(wrf_out,wks2,(/contour_tot2/),pltres,mpres) ; plot will not get drawn at this point
>
> ;*********************plot1**************************************************
> txres@txFontHeightF = 0.013 ; font smaller. default big
> txres@txFont = "text-symbols"
> dumid1 = gsn_add_text(wks2,plot,"S",xlon1,xlat,txres) ; need to replace xlon1 and xlat with appropriate lon/lat values
> ; page coordinates are normalized. go from 0->1 not 0->8" etc.
> txres@txFontHeightF = 0.015 ; font smaller. default big
> txres@txFont = "courier-bold"
> dumid2 = gsn_add_text(wks2,plot,"CDAC",xlon2,xlat,txres) ; need to replace xlon2 and xlat with appropriate lon/lat values
> txres@txFontHeightF = 0.01 ; font smaller. default big
> txres@txFont = "courier"
> dumid3 = gsn_add_text(wks2,plot,"(Continent Boundary may differ from Political boundary) ",xlon3,xlat1,txres) ; need to replace xlon3 and xlat with appropriate lon/lat values
>
> draw(plot) ; This will draw plot and the three text strings.
> frame(wks2) ; now advance frame
>
> --Mary
>
>
> On Jun 17, 2010, at 6:29 AM, Sahidul wrote:
>
> > Dear NCL users forum,
> >
> > I am facing problem to resize the image created by following script:
> >
> > We are not able to print the "text" using this following script i.e., required text is not writing in image file. Kindly help me to sort out this problem
> >
> > =================================
> > 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"
> > load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> >
> > begin
> >
> > datadir = "/WRF/output/2010061600/"
> > date1="2010-06-16_00"
> > date2="2010-06-16_23"
> > date3="2010-06-17_23"
> > date4="2010-06-18_23"
> >
> > FILES1 = systemfunc ("ls -1 " + datadir + "wrfout_d02_"+date1)
> > FILES2 = systemfunc ("ls -1 " + datadir + "wrfout_d02_"+date2)
> > FILES3 = systemfunc ("ls -1 " + datadir + "wrfout_d02_"+date3)
> > FILES4 = systemfunc ("ls -1 " + datadir + "wrfout_d02_"+date4)
> >
> > wrf_out = addfile(FILES1+".nc","r")
> >
> > res = True
> > opts = res
> > opts@InitTime = False
> > opts@Footer = False
> >
> > ; Total Rainfall calculation
> > ;***********************************************************************
> > rainc1 = wrf_user_getvar(wrf_out,"RAINC",0)
> > rainnc1 = wrf_user_getvar(wrf_out,"RAINNC",0)
> >
> > wrf_out2 = addfile(FILES2+".nc","r")
> >
> > ; Total Rainfall calculation
> > ;***********************************************************************
> > rainc2 = wrf_user_getvar(wrf_out2,"RAINC",0)
> > rainnc2 = wrf_user_getvar(wrf_out2,"RAINNC",0)
> > rain_tot2=((rainc2+rainnc2)-(rainnc1+rainc1))/10.
> > ;***********************************************************************
> > wrf_out3 = addfile(FILES3+".nc","r")
> >
> > ; Total Rainfall calculation
> > ;***********************************************************************
> > rainc3 = wrf_user_getvar(wrf_out3,"RAINC",0)
> > rainnc3 = wrf_user_getvar(wrf_out3,"RAINNC",0)
> > rain_tot3=((rainc3+rainnc3)-(rainnc2+rainc2))/10.
> > ;***********************************************************************
> > wrf_out4 = addfile(FILES4+".nc","r")
> >
> > ; Total Rainfall calculation
> > ;***********************************************************************
> > rainc4 = wrf_user_getvar(wrf_out4,"RAINC",0)
> > rainnc4 = wrf_user_getvar(wrf_out4,"RAINNC",0)
> > rain_tot4=((rainc4+rainnc4)-(rainnc3+rainc3))/10.
> >
> > rain_tot2@description = ""
> > rain_tot3@description = ""
> > rain_tot4@description = ""
> > ;***********************************************************************
> > ;Open Workstation ps,eps,X11,pdf,ncgm
> > ;***********************************************************************
> > wks2 = gsn_open_wks("ps","precip_surf_day1")
> >
> > res@InitTime = False
> > res@Footer = False
> > res@gsnMaximize = True
> > res@MainTitle = "Total Precipitation"
> >
> > opts@cnLinesOn = False
> > opts@cnFillOn = True
> > opts@cnLevelSelectionMode = "ExplicitLevels"
> > ;opts@UnitLabel = "cm"
> > ; opts@description = "Precipitation Rate"
> >
> > opts@cnLevels = (/ 1., 2., 5., 10., 15., 20., 25., 30., 35., 40., 50./)
> > opts@cnFillColors = (/"white","seagreen1","green4","greenyellow","yellow","goldenrod1","goldenrod2","tomato","orangered",\
> > "red","red2","red4"/)
> > opts@cnInfoLabelOn = False
> > opts@cnConstFLabelOn = False
> > res@gsnMaximize = True ; must include w/ Paper Orientation
> > res@gsnPanelBottom = 0.2 ; save space at bottom
> >
> > txres = True ; text mods desired
> >
> > pltres = True
> > mpres = True
> > mpres@mpOutlineBoundarySets = "National"
> > mpres@mpNationalLineThicknessF = 1.5
> > mpres@mpNationalLineColor = "Black"
> > mpres@mpGeophysicalLineColor = "Black"
> > mpres@mpGeophysicalLineThicknessF = 1.5
> > ;***********************************************************************
> > contour_tot2 = wrf_contour(wrf_out,wks2,rain_tot2,opts)
> > plot = wrf_map_overlays(wrf_out,wks2,(/contour_tot2/),pltres,mpres)
> >
> > ;*********************plot1**************************************************
> > txres@txFontHeightF = 0.013 ; font smaller. default big
> > txres@txFont = "text-symbols"
> > gsn_text_ndc(wks2,"S",0.12,.19,txres)
> > ; page coordinates are normalized. go from 0->1 not 0->8" etc.
> > txres@txFontHeightF = 0.015 ; font smaller. default big
> > txres@txFont = "courier-bold"
> > gsn_text_ndc(wks2,"CDAC",0.17,.19,txres)
> > txres@txFontHeightF = 0.01 ; font smaller. default big
> > txres@txFont = "courier"
> > gsn_text_ndc(wks2,"(Continent Boundary may differ from Political boundary) ",0.6,.19,txres)
> > frame(wks2) ; now advance frame
> > end
> > =================================
> >
> > Thanking you in advance
> >
> > Sahidul
> > _______________________________________________
> > 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

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Jun 22 07:50:46 2010