;cjc=gm('yacjc/64/1',sc=24,y='all')
;bwb=gm('yabwb/64/1',sc=24,y=cjc.lbyr)
;c3=gm('ctl3/x',sc=24)

n=n_elements(cjc)

!p.multi=[0,1,3]

cols=[!p.color,3,2,4,1,indgen(20)+5]

plot,pp_date(cjc),pp_area_avg(pp_diff(cjc,c3)),th=3,yr=[-0.5,0.5],/nodata,/xs
polyfill,col=16,(makerange(pp_date(cjc)))([0,1,1,0,0]),([0.1,-0.1])([0,0,1,1,0])
oplot,pp_date(cjc),pp_area_avg(pp_diff(cjc,c3)),th=3
oplot,pp_date(cjc),pp_area_avg(pp_diff(bwb,c3)),th=3,line=2

add_key,['yacjc (new)','yabwb (old)'],[0,2],/line,th=3

for h=0,1 do begin

text=strarr(4)
cols=[!p.color,3,2,6,1,5,4,7,8]
for i=0,3 do begin

  la=(90-i*22.5)*(1-2*h)
  sla=(90-(i+1)*22.5)*(1-2*h)

d=plr(la=la,sla=sla,pp_diff(cjc,c3))
print,yc(d(0))

  if (i eq 0) then begin
    plot,pp_date(cjc),pp_area_avg(plr(la=la,sla=sla,pp_diff(cjc,c3))),th=3,yr=[-1,1],/nodata,/xs
    polyfill,col=16,(makerange(pp_date(cjc)))([0,1,1,0,0]),([0.1,-0.1])([0,0,1,1,0])
  endif
  oplot,pp_date(cjc),pp_area_avg(plr(la=la,sla=sla,pp_diff(cjc,c3))),th=3,col=cols(i)
  oplot,pp_date(cjc),pp_area_avg(plr(la=la,sla=sla,pp_diff(bwb,c3))),th=3,line=2,col=cols(i)
  text(i)='['+shtstr(la)+' - '+shtstr(sla)+']'
endfor

add_key,text,cols,/top,th=3

endfor

gettwogifs,out='mean'

end