function plotsdisp(td,t,rcv,sdisp,range,oload)
if nargin<4, sdisp=[1,1,1,0]; end
if nargin<5, range=[-0.1,0.1]; end
if nargin<6, oload='K:\gps\pos\oload_gsi.blq'; end
posr=readpos(td,0,rcv,'','approx')';
[p,E]=geodtoecef(eceftogeod(posr));
[odisp,ophas]=readoload(rcv,oload);
for n=1:length(t)
    tu=td+(t(n)-13)/86400;
    erp=readerp(tu,'K:\gps\erp\iers');
    [U,a,b,gmst]=ecsftoecef(tu,erp);
    [rsun,rmoon]=sunmoonpos(tu);
    d(n,:)=(E*sitedisp(tu,posr,U*rsun,U*rmoon,odisp,ophas,gmst,erp,sdisp))';
end
gut('newfigg','','',[600,400]);
margin=[0.08,0.03,0.04,0.012];
yl={'East (m)','North (m)','Up (m)'};
for n=1:3
    if n<3, topts='nolabel'; else topts=''; end
    ggt('subplotv','',n,3,'taxis',td,'xlim',[t(1),t(end)]/3600,'ylim',range,...
        'topts',topts,'margin',margin);
    plot(t/3600,d(:,n)); ylabel(yl{n});
    if n==1, title(['Site Displacements : ',rcv]), end
end
