if undefined(dltau) then begin ; restore, 'ttau.idlsave' ; restore, 'd3t36g10mm00n02.30-87.avg_z_tau.idlsave' ltau1 = -7.0 & ltau2 = 7.0 & nltau = 201 dltau = (ltau2-ltau1)/float(nltau-1) altau = ltau1 + findgen(nltau)*dltau t1 = 2000.0 & t2 = 18000.0 & nt = 401 dt = (t2-t1)/float(nt-1) at = t1 + findgen(nt)*dt d=hist_2d(ltau,t,BIN1=dltau, BIN2=dt, MAX1=ltau2, MAX2=t2, MIN1=ltau1, MIN2=t1) d = float(smooth(d,4,/edge_truncate)) ; d = float(d) ; for i=0,nltau-1 do begin ; d(i,*) = total(d(i,*),/cumulative)/total(d(i,*)) ; endfor endif if keyword_set(psplot) then begin a4ps, /half, /encapsulated, file='ttaupdf.eps' redsave = red bluesave = blue xmargin = !x.margin ymargin = !y.margin red=1 & blue=5 !x.margin=[5.7,0.5] !y.margin=[3.3,0.5] endif ;levs=[0.01, 0.1, 0.3, 0.5, 0.7, 0.9, 0.99]; findgen(11)/10.0 levs = [31.6, 100.0, 316.0, 1000.0, 3162.0] anno = ['', '100', '', '1000', ''] contour, d, altau, at/1e3, levels=levs, $ xrange = [-5, 5], /xstyle, $ xtitle='log !9t!3!Dross!N', $ yrange = [2,12.0], /ystyle, ytitle='Temperature [10!E3!N K]', $ c_thick=2, c_charthick=2, c_annotation=anno, c_labels=[0,1,0,1,0] oplot, avg_z_tau.tau.tau, aver(avg_z_tau.tau.tavg1,2)*1e-3, thick=4 if keyword_set(psplot) then begin a4ps, /close red = redsave blue = bluesave !x.margin = xmargin !y.margin = ymargin endif end