[a,b]=size(jadson);
[c,d]=size(jad);
rat=d/b;
i=1:d;
j=1:b;
jj=floor(j*rat);
dj=j*rat-jj;
lrgm1(1)=lrgm(1);cogm1(1)=cogm(1);sogm1(1)=sogm(1);imum1(1)=imum(1);imdm1(1)=imdm(1);odec1(1)=odec(1);jad1(1)=jad(1);
for k=2:b-1;
   lrgm1(k)=lrgm(jj(k))*(1-dj(k))+lrgm(jj(k)+1)*dj(k);
   jad1(k)=jad(jj(k))*(1-dj(k))+jad(jj(k)+1)*dj(k);
   sogm1(k)=sogm(jj(k))*(1-dj(k))+sogm(jj(k)+1)*dj(k);
   cogm1(k)=cogm(jj(k))*(1-dj(k))+cogm(jj(k)+1)*dj(k);
	imum1(k)=imum(jj(k))*(1-dj(k))+imum(jj(k)+1)*dj(k);
   imdm1(k)=imdm(jj(k))*(1-dj(k))+imdm(jj(k)+1)*dj(k);
   odec1(k)=odec(jj(k))*(1-dj(k))+odec(jj(k)+1)*dj(k);
end;
lrgm1(b)=lrgm(d);cogm1(b)=cogm(d);sogm1(b)=sogm(d);imum1(b)=imum(d);imdm1(b)=imdm(d);odec1(b)=odec(d);jad1(b)=jad(d);
dc=pi/180;
un1=U.*cos(lrgm1*dc)+sogm1.*cos(cogm1*dc)+V.*sin(lrgm1*dc);
ue1=U.*sin(lrgm1*dc)+sogm1.*sin(cogm1*dc)-V.*cos(lrgm1*dc);
un2=(U+odec1).*cos(lrgm1*dc)+V.*sin(lrgm1*dc);
ue2=(U+odec1).*sin(lrgm1*dc)-V.*cos(lrgm1*dc);
s1=sqrt(un1.*un1+ue1.*ue1);
s2=sqrt(un2.*un2+ue2.*ue2);
dir1=mod(atan2(-ue1,-un1)*180/pi+360,360);
dir2=mod(atan2(-ue2,-un2)*180/pi+360,360);
