;plots and analysis for synthesis paper and presentations

PRO syn_compare
;compare monthly and annual flux from the major methods

  dta = read_ascii(mydocs()+'synthesis/monthlyflux.txt')
  dta = dta.(0)
  bd = where(dta EQ -999,nbd)
  IF nbd GT 0 THEN dta[bd] = nan()
  header = ['Time','Year','Month','LEF','Msite','ED','Helliker','Hell04','Bakwin']

  ann = average_cols(dta,12)
  ann[3:*,*]*=12.

  annsub = ann[[3,5,6],0:4]
  annsubmn = total(annsub,2)/5.
  annsubmn = reform(replicate_arr(annsubmn,5),3,5)
  anndev = annsub-annsubmn

  
  annsum = dta[[3,5,6],5:53:12]+dta[[3,5,6],6:54:12]+dta[[3,5,6],7:55:12]
  annsummn = total(annsum,2)/5.
  annsummn = reform(replicate_arr(annsummn,5),3,5)
  annsumdev = annsum-annsummn

  stop


  


;LEF cumulative per year
 
  ;97 98 99 00 01 xx 03 04

  black = 0
  red = fsc_color('red',200)
  blue = fsc_color('blue',201)
  orange = fsc_Color('brown',202)
  purple = fsc_color('purple',203)
  cyan = fsc_color('cyan',204)
  yellow = fsc_color('orange',206)
  gyellow = fsc_color('green yellow',207)

  wcolors = [black,red,blue,orange,purple,black,cyan,yellow]
  annlef = average_arr(dta[3,*],12)*12.

  plot,[0,0],[0,0],/nodata,xrange=[0,12],yrange=[-100,200],xticks=11,xtickv=numgen(0,11,1),xtickn=['J','F','M','A','M','J','J','A','S','O','N','D'],position=[0.1,0.1,0.8,0.9],title='LEF Cumulative NEE 1997-2004',xtitle='Month',ytitle='gC m-2 yr-1',xticklen=0.5,yticklen=0.5
  FOR i = 0,7 DO oplot,numgen(0.,12.,0.1),spline(findgen(13),total([0,reform(dta[3,(i*12):((i*12)+11)])],/cum),numgen(0.,12.,0.1)),color=wcolors[i],thick=3
; FOR i = 0,7 DO oplot,total([0,reform(dta[3,(i*12):((i*12)+11)])],/cum),color=wcolors[i],thick=3
  FOR i = 0,7 DO xyouts,12.5+(i/4.),annlef[i],string(1997+i,format='(i0)'),color=wcolors[i],charsize=1.0

;  plot,[0,0],[0,0],/nodata,xrange=[0,12],yrange=[-50,50],xticks=11,xtickv=numgen(0,11,1),xtickn=['J','F','M','A','M','J','J','A','S','O','N','D']
;  FOR i = 0,7 DO
;  oplot,reform(dta[3,(i*12):((i*12)+11)]),color=wcolors[i],thick=3


;monthly flux comparison

;LEF = red, Msite = green-yellow, ED = Blue, Helliker = orange, bakwin
;= black

  stop

 black = 0
  red = fsc_color('red',200)
  blue = fsc_color('blue',201)
  orange = fsc_Color('brown',202)
  purple = fsc_color('purple',203)
  cyan = fsc_color('cyan',204)
  yellow = fsc_color('orange',206)
  gyellow = fsc_color('green yellow',207)

  mcolors = [red,gyellow,blue,yellow,black]
  mnames = ['LEF','Multisite','ED Model','ABL-Helliker','ABL-Bakwin']
  dtaplot = dta[[3,4,5,6,8],*]

  plot,[0,0],[0,0],/nodata,xrange=[1997,2005],yrange=[-150,50],title='Monthly NEE 1997-2004',xticks=8,xtickv=numgen(1997,2005,1),xticklen=0.5
  FOR i = 0,4 DO oplot,dta[0,*],dtaplot[i,*],color=mcolors[i],thick=3
  FOR i = 0,4 DO plots,[1997.5,1998],[(-80)-(14*i),(-80)-(14*i)],color=mcolors[i],thick=3
  FOR i = 0,4 DO xyouts,1998.1,(-81)-(14*i),mnames[i]

END