function plotp1c1
t=caltomjd([2000,4,2]):7:caltomjd([2005,6,1]);
s={}; for n=1:31, s={s{:},sprintf('GPS%02d',n)}; end
for n=1:length(s)
    for m=1:length(t)
        bias(m,n)=biasp1c1(s{n},t(m));
    end
end
bias(bias==0)=nan;
tt=(t-t(1))*24;
gut('newfigg','','',[600,400]);
ggt('newaxes','',[0.08,0.08,0.86,0.86],'taxis',t(1),'xlim',[tt(1),tt(end)]);
plot(tt,bias,'-'); plot(tt,bias,'.');
ylabel('P1-C1 Bias (m)'); title('P1-C1 Bias Values')
