load poltidm.prn -ascii
djmd=poltidm(:,3);
dgpol1=poltidm(:,4);
dgpol2=poltidm(:,5);
diff=poltidm(:,6);
plot(djmd,diff,'r')
title('gravity poltide EOP97C04 - EOPC04 at BFO Schiltach 1987-1994')
xlabel('Modified Julian date')
ylabel('nm/s**2')