%ajitter2

ok_dday = data.dday(okrefli);
ok_refu = table1([ref_dday ref_u],  ok_dday);
ok_refv = table1([ref_dday ref_v],  ok_dday);
ok_smu = table1([sm_dday sm_u],  ok_dday);
ok_smv = table1([sm_dday sm_v],  ok_dday);
ok_lon = data.lon(okrefli);
ok_lat = data.lat(okrefli);


h1= subplot(511);
h=plot(ref_dday, ref_u,'g');lw(h,.1)
hold on
h=plot(sm_dday, sm_u,'k.');ms(h,1)
yl(-100,100)


h2 = subplot(512);
h=plot(ok_dday, ok_refu - ok_smu, 'go'); ms(h,3);hold on
h=plot(ok_dday, ok_refu - ok_smu, 'k.'); ms(h,4)
h=plot(ok_dday(okosi), ok_refu(okosi) - ok_smu(okosi), 'ro'); ms(h,3);
yl(-20,20)

h3 = subplot(513);
h=plot(ref_dday, ref_v,'m');lw(h,.1)
hold on
h=plot(sm_dday, sm_v,'k.');ms(h,1)
yl(-100,100)


h4 = subplot(514);
h=plot(ok_dday, ok_refv - ok_smv, 'mo'); ms(h,3);hold on
h=plot(ok_dday, ok_refv - ok_smv, 'k.'); ms(h,4)
h=plot(ok_dday(okosi), ok_refv(okosi) - ok_smv(okosi), 'bo'); ms(h,3);
yl(-20,20)


h5 = subplot(515)
plot(data.dday, data.clmps, 'c')
hold on 
h=plot(ok_dday, data.clmps(okrefli), 'k.');ms(h,4)
h=plot(ok_dday(okosi), data.clmps(okrefli(okosi)), 'ro');ms(h,3)
yl(0,8)

for ii=1:5
   eval(['hh = h' num2str(ii) ';']);
   pos = get(hh, 'pos');
   set(hh, 'pos', pos + [0 0 0 .03]);
   if (ii < 5)
      set(hh, 'xtick', [])
   end
end


