pro plotwith
;+
;plots withbroe tremps for max and solar min holes
;-
rrr=dindgen(1000)/100d
ww1=withbroe(rrr,1)
ww3=withbroe(rrr,3)
plot,rrr,ww1(1,*),xr=[1.,2.],xs=1,ys=1,xtitle='r/R',ytitle='Temperature ( MK )',yr=[0.,1.5],charsiz=1.3,yticks=3.,thick=2.
oplot,rrr,ww3(1,*),linestyle=3,thick=2.
oplot,[1.65,1.7],[.67,.67],linestyle=2 ,thick=2.
xyouts,1.72,.67,'Ne, Solar Max',charsiz=1.1
oplot,[1.65,1.7],[.6,.6],thick=2.,linestyle=4
xyouts,1.72,.6,'Ne, Solar Min',charsiz=1.1
oplot,[1.65,1.7],[.81,.81],thick=2. 
xyouts,1.72,.74,'Te, Solar Min',charsiz=1.1
oplot,[1.65,1.7],[.74,.74],thick=2.,linestyle=3
xyouts,1.72,.81,'Te, Solar Max',charsiz=1.1
oplot,rrr,(alog10(ww1(2,*))-6.)/2.,thick=2.,linestyle=2
oplot,rrr,(alog10(ww3(2,*))-6.)/2.,linestyle=4,thick=2.
xyouts,2.01,0.,'6.0',charsiz=1.3
xyouts,2.03,1.1,'LOG(Electron Density)',orientation=-90.,charsiz=1.2
xyouts,2.01,1.45,'9.0',charsiz=1.3
end