function [jd1, jd2, inc]=stkdscl(t1,t2)

jd1=floor(t1)-1;
deltat=t2-t1;
inc=1;
if (deltat > 14), inc=7; end;
if (deltat > 60), inc=30; end;
jd2= jd1 + inc + inc*(floor(deltat/inc));

return
