Hello all,
I'm working through some issues with saving postscript files using function graphics. When I run the following code, the resulting .ps file only shows the first of the two plots, and if I expand the window it seems its been replicated multiple times in the image.
I've spent the better part of an hour researching this issue with no leads. Is there any way around this besides reverting all of my code to direct graphics? Because that is not ideal at this time.
TRMMlbl=['TRMM Stratiform', 'TRMM convective', 'TRMM shallow non-isolated', 'TRMM shallow isolated', 'TRMM other']
GPMlbl=['GPM Stratiform', 'GPM convective', 'GPM shallow non-isolated', 'GPM shallow isolated', 'GPM other']
p=plot(TRMM_sf_topcount/16, altitude, xtitle='count', ystyle=1, ytickdir=1, ytitle='Storm top height(km)', $
title='Storm top height histogram for TRMM(climo) and GPM(201404-201503)', color='red', layout=[2, 1, 1], $
name=TRMMlbl[0], dimensions=[1400, 700], margin=.08, font_size=13);,/buffer)
p1=plot(TRMM_conv_topcount/16, altitude, linestyle=1, color='red', /overplot, name=TRMMlbl[1])
p2=plot(TRMM_sh_topcount/16, shalaltitude, linestyle=2, color='red', /overplot, name=TRMMlbl[2])
p3=plot(TRMM_si_topcount/16, shalaltitude, linestyle=3, color='red', /overplot, name=TRMMlbl[3])
p4=plot(TRMM_other_topcount/16, altitude, linestyle=4, color='red', /overplot, name=TRMMlbl[4])
p5=plot(GPM_sf_topcount, altitude, color='blue', /overplot, name=GPMlbl[0])
p6=plot(GPM_conv_topcount, altitude, linestyle=1, color='blue', /overplot, name=GPMlbl[1])
p7=plot(GPM_sh_topcount, shalaltitude, linestyle=2, color='blue', /overplot, name=GPMlbl[2])
p8=plot(GPM_si_topcount, shalaltitude, linestyle=3, color='blue', /overplot, name=GPMlbl[3])
p9=plot(GPM_other_topcount, altitude, linestyle=4, color='blue', /overplot, name=GPMlbl[4])
p10=plot(sfrain_TRMM/16, altitude, xtitle='Total rain (mm/yr)', ystyle=1, ytickdir=1, ytitle='Storm top height(km)', $
title='Total rainfall vs storm top heights for TRMM(climo) and GPM(201404-201503)', color='red', layout=[2, 1, 2], /current, $
name=TRMMlbl[0], margin=.08, font_size=13)
p11=plot(convrain_TRMM/16, altitude, linestyle=1, color='red', /overplot, name=TRMMlbl[1])
p12=plot(shrain_TRMM/16, shalaltitude, linestyle=2, color='red', /overplot, name=TRMMlbl[2])
p13=plot(sirain_TRMM/16, shalaltitude, linestyle=3, color='red', /overplot, name=TRMMlbl[3])
p14=plot(otherrain_TRMM/16, altitude, linestyle=4, color='red', /overplot, name=TRMMlbl[4])
p15=plot(sfrain_GPM,altitude, color='blue', /overplot, name=GPMlbl[0])
p16=plot(convrain_GPM,altitude, linestyle=1, color='blue', /overplot, name=GPMlbl[1])
p17=plot(shrain_GPM, shalaltitude, linestyle=2, color='blue', /overplot, name=GPMlbl[2])
p18=plot(sirain_GPM, shalaltitude, linestyle=3, color='blue', /overplot, name=GPMlbl[3])
p19=plot(otherrain_GPM, altitude,linestyle=4, color='blue', /overplot, name=GPMlbl[4])
leg1=legend(target=[p, p1, p2, p3, p4], position=[.25, .87], /auto_text_color, font_size=10)
leg2=legend(target=[p5, p6, p7, p8, p9], position=[.45, .87], /auto_text_color, font_size=10)
p19.save, '/dwn1/kwhite/idl/THESIS/images/TRMM_GPM/TRMM_climo_stormht_ rvsth.ps'
Thanks,
Keith
|