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