;
; Plots distributions and cutoffs of trends and means of different length
; segments.
;
; Get pre-analysed data
;
restore,filename='trendmean_distrib.idlsave'
;  ntry,tlen,nts,alltit,tren95,mean95,binmean,bintren,freqmean,freqtren,$
;  nbin,binsc
;
; Prepare window/paper for plotting
;
loadct,39
multi_plot,nrow=3,ncol=2
if !d.name eq 'X' then begin
  window, ysize=800
  !p.font=-1
endif else begin
  !p.font=0
  device,/helvetica,/bold,font_size=12
endelse
;
; Select which distributions to show
;
;iuse=[2,6,12]
iuse=[6]
nuse=n_elements(iuse)
;
; Compute distribution totals for scaling by
;
fmtot=total(freqmean,1)/100.
fttot=total(freqtren,1)/100.
;
; Repeat for each time series
;
for its = 0 , nts-1 do begin
  ;
  pause
  plot,binmean*binsc(its),freqmean(*,iuse(0),its,0),/nodata,$
    /xstyle,xtitle='Mean anomaly  (!Uo!NC)',$
    /ystyle,ytitle='Frequency of occurrence!C!C(%)',yrange=[0.,50.]
  oplot,[0.,0.],!y.crange,linestyle=1
  ;
  xyouts,/normal,0.5,0.33+0.33*(!p.multi(0)-1)*0.5,align=0.5,alltit(its)
  ;
  for j = 0 , nuse-1 do begin
    oplot,binmean*binsc(its),freqmean(*,iuse(j),its,0)/fmtot(iuse(j),its,0),$
      thick=4
    oplot,binmean*binsc(its),freqmean(*,iuse(j),its,1)/fmtot(iuse(j),its,1)
  endfor
  ;
;  plot,bintren*binsc(its),freqtren(*,iuse(0),its,0),/nodata,$
;    /xstyle,xtitle='Trend  (!Uo!NC per 30 yr)',$
;    /ystyle,ytitle='Frequency of occurrence!C!C(%)',yrange=[0.,30.]
;  oplot,[0.,0.],!y.crange,linestyle=1
  ;
;  for j = 0 , nuse-1 do begin
;    oplot,bintren*binsc(its),freqtren(*,iuse(j),its,0)/fttot(iuse(j),its,0),$
;      thick=4
;    oplot,bintren*binsc(its),freqtren(*,iuse(j),its,1)/fttot(iuse(j),its,1)
;  endfor
  ;
  if its le 2 then yr=[0.,0.3] else yr=[0.,1.0]
  plot,tlen,mean95(*,its,0),thick=4,psym=-def_sym(10),$
    /xstyle,xtitle='Segment length  (years)',$
    ytitle='95th centile of segment means!C!C(!Uo!NC)',$
    /ystyle,yrange=yr
  oplot,tlen,mean95(*,its,1),psym=-def_sym(15)
  ;
;  if its le 2 then yr=[0.,0.6] else yr=[0.,2.6]
;  plot,tlen,30.*tren95(*,its,0)/tlen,thick=4,psym=-def_sym(10),$
;    /xstyle,xrange=[10,100],xtitle='Segment length  (years)',$
;    ytitle='95th centile of segment trends!C!C(!Uo!NC per 30 yr)',$
;    /ystyle,yrange=yr
;  oplot,tlen,30.*tren95(*,its,1)/tlen,psym=-def_sym(15)
  ;
endfor
;
end