pro planetary_graph, TvFlag = TvFlag If N_Elements(TvFlag) eq 0 then TvFlag = 1 erase If TvFlag eq 1 then begin window, xsize=1024, ysize=768 EndIf else If TvFlag eq 3 then begin set_plot,'z' device, set_resolution = [1600,1200] endif !p.thick=3 !p.charsize=2 !p.charthick=2 a=findgen(4000)/100. plot,a,a^(3./2), color=fsc_color('black'), background=fsc_color('white'),$ xtitle = 'A', ytitle = 'A!U3/2!N', $ yticklen=1,xticklen=1,ygridstyle=2,xgridstyle=2,$ xticklayout=0, yticklayout=0,$ xticks=6,yticks=6, ytickinterval = 10, xtickinterval=2,$ position = [0.1,0.1,0.95,0.95],/normal If TvFlag eq 3 then begin PngImg = tvrd() tvlct,reds,greens,blues,/get write_png,'~/Desktop/PlanetaryYears0to40.png', $ PngImg,reds,greens,blues EndIf stop erase a=findgen(2000)/1000. plot,a,a^(3./2), color=fsc_color('black'), background=fsc_color('white'),$ xtitle = 'A', ytitle = 'A!U3/2!N', $ yticklen=1,xticklen=1,ygridstyle=2,xgridstyle=2,$ xticklayout=0, yticklayout=0,$ xticks=6,yticks=6, ytickinterval = .2, xtickinterval=.1,$ position = [0.1,0.1,0.95,0.95],/normal If TvFlag eq 3 then begin PngImg = tvrd() tvlct,reds,greens,blues,/get write_png,'~/Desktop/PlanetaryYears0to2.png', $ PngImg,reds,greens,blues EndIf end