subroutine ephob(dj,xp,vp,x,v) c ******************************************************************** c * calcul d'une ephemeride de phobos dans le repere defini par ref * c * en entree : date julienne dj * c * en sortie : vecteur position x(3) et vecteur vitesse v(3) * c * unite=km et km/jour * c ******************************************************************** implicit double precision (a-h,o-z) save dimension x(3),v(3),xp(3),vp(3),x2000(3),v2000(3), 1 c(3,3),rmat(3,3) double precision a,b data kle/1/ data d2000/2451545.0d0/ data ((rmat(i,j),j=1,3),i=1,3)/ 1 0.9999256791774783d0,-0.0111815116768724d0,-0.0048590038154553d0, 2 0.0111815116959975d0, 0.9999374845751042d0,-0.0000271625775175d0, 3 0.0048590037714450d0,-0.0000271704492210d0, 0.9999881946023742d0/ if(kle.eq.1) then call lec kle=2 endif t=dj-d2000 c..equatorial plane solution call cal(t,xp,vp) c..fk4 solution call refzzz(t,c) do i=1,3 w=0.d0 wp=0.d0 do k=1,3 w=w+c(i,k)*xp(k) wp=wp+c(i,k)*vp(k) enddo x(i)=w v(i)=wp enddo c..fk5 solution do j=1,3 a = 0.0d0 b = 0.0d0 do k=1,3 a = a + rmat(j,k)*x(k) b = b + rmat(j,k)*v(k) enddo x(j) = a v(j) = b enddo return end subroutine edeim(dj,xp,vp,x,v) c ******************************************************************** c * calcul d'une ephemeride de deimos dans le repere defini par ref * c * en entree : date julienne dj * c * en sortie : vecteur position x(3) et vecteur vitesse v(3) * c * unite=km et km/jour * c ******************************************************************** implicit double precision (a-h,o-z) save dimension x(3),v(3),xp(3),vp(3),x2000(3),v2000(3), 1 c(3,3),rmat(3,3) double precision a,b data kle/1/ data d2000/2451545.0d0/ data ((rmat(i,j),j=1,3),i=1,3)/ 1 0.9999256791774783d0,-0.0111815116768724d0,-0.0048590038154553d0, 2 0.0111815116959975d0, 0.9999374845751042d0,-0.0000271625775175d0, 3 0.0048590037714450d0,-0.0000271704492210d0, 0.9999881946023742d0/ if(kle.eq.1) then call lecd kle=2 endif t=dj-d2000 c..equatorial plane solution call cald(t,xp,vp) c..fk4 solution call refzzz(t,c) do i=1,3 w=0.d0 wp=0.d0 do k=1,3 w=w+c(i,k)*xp(k) wp=wp+c(i,k)*vp(k) enddo x(i)=w v(i)=wp enddo c..fk5 solution do j=1,3 a = 0.0d0 b = 0.0d0 do k=1,3 a = a + rmat(j,k)*x(k) b = b + rmat(j,k)*v(k) enddo x(j) = a v(j) = b enddo return end subroutine cal(t,x,v) c *************************************************************** c * phobos * c * calcul des vecteurs position x et vitesse au temps t * c * t est mesure en jours depuis j2000 * c *************************************************************** implicit double precision (a-h,o-z) save dimension x(3),v(3) dimension nb(3),nbv(3) dimension arg(38,3),freq(38,3),cs(38,3),cc(38,3), 1 freq2(38,3),argv(36,3),freqv(36,3),cvs(36,3),cvc(36,3), 2 frev2(36,3) common/ep5/arg,freq,cs,cc,freq2,argv,freqv,cvs,cvc,frev2 common/ep4/nb,nbv common/ep2/accep deg=3600.d0/206264.8062470964d0 t2=t*t*accep*deg do 1 i=1,3 jmax=nb(i) if(i.ne.3) jmax=jmax+4 xw=0 do 3 j=1,jmax a=arg(j,i)+freq(j,i)*t+freq2(j,i)*t2 sa=sin(a) ca=cos(a) xw=xw+cs(j,i)*sa+cc(j,i)*ca 3 continue x(i)=xw 1 continue do 11 i=1,3 jmax=nbv(i) if(i.ne.3) jmax=jmax+4 xw=0 do 13 j=1,jmax a=argv(j,i)+freqv(j,i)*t+frev2(j,i)*t2 sa=sin(a) ca=cos(a) xw=xw+cvs(j,i)*sa+cvc(j,i)*ca 13 continue v(i)=xw 11 continue return end subroutine blocdatas c block data c *************************************************************** c * description des commons * c * ep1 : dnu(deg/jour), dgam = de corrections aux elements * c * moyens metriques de la theorie de phobos * c * lg, h, pi (en deg)= longitude moyenne, longitudes * c * du noeud et du pericentre en j2000 * c * ep2 : terme en t**2 de la longitude moyenne * c * ep3 : angles donnant la position de mars en deg, correction * c * a la precession de mars en deg/jour * c * ep4 : nombre de termes dans les series position et vitesse * c * de phobos (sauf perturbations par deimos et planetes) * c *************************************************************** c *************************************************************** c * description des commons * c * ed1 : dnu(deg/jour), dgam = de corrections aux elements * c * moyens metriques de la theorie de deimos * c * lg, h, pi (en deg)= longitude moyenne, longitudes * c * du noeud et du pericentre en j2000 * c * ed2 : terme en t**2 de la longitude moyenne * c * ed4 : nombre de termes dans les series position et vitesse * c * de deimos * c *************************************************************** implicit double precision(a-h,o-z) save integer nb(3),nbv(3),nbd(3),nbvd(3) dimension pm(3),al(3) dimension dm(3),ald(3) common/ep1/pm,al common/ep2/accep common/ep3/psi0,pipet,alm,dp common/ep4/nb,nbv common/ed1/dm,ald common/ed2/acced common/ed4/nbd,nbvd nb(1)=27 nb(2)=27 nb(3)=10 nbv(1)=32 nbv(2)=32 nbv(3)=12 psi0=208.5619d0 pipet=71.0053d0 alm=19.3730d0 dp=-0.181103d-6 pm(1)=0.4994d-3 pm(2)=-0.341d-3 pm(3)=0.146d-3 al(1)=171.9160d0 al(2)=125.88d0 al(3)=342.91d0 accep=0.9518d-8 nbd(1)=37 nbd(2)=38 nbd(3)=25 nbvd(1)=27 nbvd(2)=27 nbvd(3)=16 dm(1)=-0.333d-4 dm(2)=0.103d-3 dm(3)=-0.221d-3 ald(1)=215.2172d0 ald(2)=11.20d0 ald(3)=224.01d0 acced=-0.377d-9 return end subroutine cald(t,x,v) c************************************************************* c * deimos * c * calcul des vecteurs position x et vitesse au temps t* c * t est mesure en jours depuis j2000 * c ******************************************************* implicit double precision (a-h,o-z) save dimension x(3),v(3) integer nb(3),nbv(3) dimension arg(38,3),freq(38,3),cs(38,3),cc(38,3), 1 freq2(38,3),argv(36,3),freqv(36,3),cvs(36,3),cvc(36,3), 2 frev2(36,3) common/ed5/arg,freq,cs,cc,freq2,argv,freqv,cvs,cvc,frev2 common/ed4/nb,nbv common/ed2/accep deg=3600.d0/206264.8062470964d0 t2=t*t*accep*deg do 1 i=1,3 jmax=nb(i) xw=0.d0 do 3 j=1,jmax a=arg(j,i)+freq(j,i)*t+freq2(j,i)*t2 sa=sin(a) ca=cos(a) xw=xw+cs(j,i)*sa+cc(j,i)*ca 3 continue x(i)=xw 1 continue do 11 i=1,3 jmax=nbv(i) xw=0.d0 do 13 j=1,jmax a=argv(j,i)+freqv(j,i)*t+frev2(j,i)*t2 sa=sin(a) ca=cos(a) xw=xw+cvs(j,i)*sa+cvc(j,i)*ca 13 continue v(i)=xw 11 continue return end subroutine refzzz(t,c) c repere fk4 implicit double precision(a-h,o-z) save dimension ar(3,3),am(3,3),g(3,3),c(3,3),d(3,3) data g/ 0.999925674124d0,-0.011181963465d0,-0.004859004081d0, 2 0.012192051720d0, 0.917413967951d0, 0.397747363640d0, 3 0.000010121726d0,-0.397777041948d0, 0.917482111431d0/ tm=t/365250.d0 rad=1.d0/206264.8062470964d0 deg=3600*rad ah=178409.13618d0*rad+tm*(-0.5149158068948755d-1 1 +tm*(-0.1117775392103901d-2-0.3427183852553758d-4*tm)) gam=0.1614120767052974d-1+tm*(-0.710928404247926d-3 1 +tm*(-0.1968131454406586d-4-0.2505007528551377d-6*tm)) tet=(35.496817571d0+2.507593d-6*t)*deg csh=cos(ah) snh=sin(ah) cst=cos(tet) snt=sin(tet) csth=csh*cst-snh*snt snth=snt*csh+snh*cst gm2=gam*gam rac=sqrt(1-gm2) gm2=2*gm2 am(1,1)=csth+gm2*snh*snt am(1,2)=-snth+gm2*snh*cst am(2,1)=snth-gm2*csh*snt am(2,2)=csth-gm2*csh*cst coef=2*gam*rac am(1,3)=coef*snh am(2,3)=-coef*csh am(3,1)=coef*snt am(3,2)=coef*cst am(3,3)=1-gm2 q=(25.192028020d0+3.269878d-7*t)*deg csq=cos(q) snq=sin(q) ar(1,1)=1 do 1 i=2,3 ar(1,i)=0 1 ar(i,1)=0 ar(2,2)=csq ar(2,3)=-snq ar(3,2)=snq ar(3,3)=csq do 7 i=1,3 do 2 j=1,3 w=0 do 3 k=1,3 3 w=w+am(i,k)*ar(k,j) 2 d(i,j)=w 7 continue do 4 i=1,3 do 5 j=1,3 w=0 do 6 k=1,3 6 w=w+g(i,k)*d(k,j) 5 c(i,j)=w 4 continue return end subroutine lec c *************************************************************** c * lecture des series esapho. introduction des constantes * c * nouvelles dans les moyens mouvements et les coefficients. * c * transformation des arguments en fonctions du temps * c *************************************************************** implicit double precision(a-h,o-z) save dimension f(6),a0(6),ai(6),pf(6),al(3) dimension nb(3),nbv(3),iar(6),ic1(4),ic2(4) dimension arg(38,3),freq(38,3),cs(38,3),cc(38,3), 1 freq2(38,3),argv(36,3),freqv(36,3),cvs(36,3),cvc(36,3), 2 frev2(36,3) dimension asup(4),fsup(4),rx1s(4),rx2s(4) character*57 x1(18),x2(18),x3(10),x1p(9),x2p(9),x(27,3) character*51 y1(18),y2(18),y3(10),y1p(9),y2p(9),y(27,3) character*61 v1(18),v2(18),v3(12),v1p(14),v2p(14),v(32,3) character*56 w1(18),w2(18),w3(12),w1p(14),w2p(14),w(32,3) common/ep5/arg,freq,cs,cc,freq2,argv,freqv,cvs,cvc,frev2 common/ep4/nb,nbv common/ep1/dnu,dgam,de,al common/ep3/psi0,pipet,alm,dp common/ep2/accep equivalence (x(1,1),x1(1)),(y(1,1),y1(1)),(x(19,1),x1p(1)), 1 (y(19,1),y1p(1)),(x(1,2),x2(1)),(y(1,2),y2(1)), 2 (x(19,2),x2p(1)),(y(19,2),y2p(1)), 3 (x(1,3),x3(1)),(y(1,3),y3(1)),(v(1,1),v1(1)), 4 (w(1,1),w1(1)),(v(19,1),v1p(1)),(w(19,1),w1p(1)), 5 (v(1,2),v2(1)),(w(1,2),w2(1)),(v(19,2),v2p(1)), 6 (w(19,2),w2p(1)),(v(1,3),v3(1)),(w(1,3),w3(1)) c **************************************************************** c * coefficients esapho * c **************************************************************** x1(1)=' 0 1 1 0 0 1 9372991.9756-6247057 -181222'// 1 ' -140805' x1(2)=' 0 1 1 0 -1 1 -210991.0229 140430 4098'// 1 '-14066059' x1(3)=' 0 1 1 0 1 1 70361.4467 -46804 -1370'// 1 ' 4689174' x1(4)=' 0 1 -1 0 1 1 -970.2930 -9509 752'// 1 ' -64574' x1(5)=' 0 1 1 -2 0 1 873.8556 -582 180963'// 1 ' -13' x1(6)=' 0 1 1 -1 0 1 0.0000 0 0'// 1 ' 0' x1(7)=' 0 1 1 0 2 1 791.9808 -526 -15'// 1 ' 105574' x1(8)=' 0 1 1 0 0 0 448.5114 -754 -9'// 1 ' -8' x1(9)=' 0 1 1 0 0 2 -448.4011 754 10'// 1 ' 6' x1(10)=' 0 2 2 -1 0 3 -1155.9707 -7995 -32964'// 1 ' -90' x1(11)=' 0 1 3 0 -1 1 324.0074 3139 -250'// 1 ' 21398' x1(12)=' 0 0 0 1 0 -1 1144.9013 7898 31815'// 1 ' 87' x1(13)=' 0 1 -1 0 0 1 289.6238 -451 -1'// 1 ' 319' x1(14)=' 0 3 1 0 0 3 -288.9041 480 30'// 1 ' 4' x1(15)=' 0 1 1 0 -2 1 265.2220 -175 -5'// 1 ' 35366' x1(16)=' 0 1 1 1 0 1 -0.0058 0 0'// 1 ' 0' x1(17)=' 0 2 2 -1 0 2 -144.5293 591 -15329'// 1 ' 10' x1(18)=' 0 0 0 1 0 0 144.4057 -592 15314'// 1 ' -8' c data x1p x1p(1)=' 2 -3 -3 0 0 -3 117.0670 -90 -6'// 1 ' -3' x1p(2)=' 2 -1 -1 0 0 -1 -89.5280 397 26'// 1 ' 1' x1p(3)=' 1 0 0 0 0 0 86.9834 307 -25'// 1 ' 13' x1p(4)=' 0 1 3 0 -2 1 -135.9361 -1940 62'// 1 ' -10809' x1p(5)=' 0 1 -1 0 1 2 -79.1722 -99 16'// 1 ' -5280' x1p(6)=' 0 1 -1 0 2 1 126.1976 1845 -55'// 1 ' 9510' x1p(7)=' 0 3 1 0 0 4 -63.2560 103 0'// 1 ' 0' x1p(8)=' 0 1 -1 0 0 2 62.8875 -101 6'// 1 ' 26' x1p(9)=' 1 -2 -2 0 0 -2 -51.7171 -183 15'// 1 ' 2' c data x2 x2(1)=' 0 1 1 0 0 1 9372992.3416-6247058 -181236'// 1 ' -140791' x2(2)=' 0 1 1 0 -1 1 -211031.5150 140592 4090'// 1 '-14068768' x2(3)=' 0 1 1 0 1 1 70361.4463 -46804 -1370'// 1 ' 4689174' x2(4)=' 0 1 -1 0 1 1 -980.1755 -9494 752'// 1 ' -65240' x2(5)=' 0 1 1 -2 0 1 874.2215 -583 181040'// 1 ' -13' x2(6)=' 0 1 1 -1 0 1 0.0000 0 0'// 1 ' 0' x2(7)=' 0 1 1 0 2 1 791.9808 -526 -15'// 1 ' 105574' x2(8)=' 0 1 1 0 0 0 448.5885 -755 -9'// 1 ' -4' x2(9)=' 0 1 1 0 0 2 -448.1873 755 5'// 1 ' 6' x2(10)=' 0 2 2 -1 0 3 -1155.9709 -7995 -32964'// 1 ' -90' x2(11)=' 0 1 3 0 -1 1 324.0074 3139 -250'// 1 ' 21398' x2(12)=' 0 0 0 1 0 -1 1146.7820 7893 32010'// 1 ' 87' x2(13)=' 0 3 1 0 0 3 -288.8736 480 31'// 1 ' 4' x2(14)=' 0 1 -1 0 0 1 -288.7630 510 59'// 1 ' 328' x2(15)=' 0 1 1 1 0 1 -0.0058 0 0'// 1 ' 0' x2(16)=' 0 1 1 0 -2 1 264.0031 -174 -5'// 1 ' 35204' x2(17)=' 0 0 0 1 0 0 173.0304 -710 18389'// 1 ' -14' x2(18)=' 0 2 2 -1 0 2 -144.5289 591 -15329'// 1 ' 10' c data x2p x2p(1)=' 2 -3 -3 0 0 -3 -117.0670 90 6'// 1 ' 3' x2p(2)=' 2 -1 -1 0 0 -1 -89.5279 397 26'// 1 ' 1' x2p(3)=' 1 0 0 0 0 0 86.9835 307 -25'// 1 ' 13' x2p(4)=' 0 1 3 0 -2 1 -135.9356 -1940 62'// 1 ' -10809' x2p(5)=' 0 1 -1 0 1 2 -81.3238 -96 16'// 1 ' -5426' x2p(6)=' 0 1 -1 0 2 1 126.2383 1845 -55'// 1 ' 9516' x2p(7)=' 0 1 -1 0 0 2 -63.8267 106 8'// 1 ' 26' x2p(8)=' 0 3 1 0 0 4 -63.2563 103 0'// 1 ' 0' x2p(9)=' 1 -2 -2 0 0 -2 51.7171 183 -15'// 1 ' -2' c data x3 x3(1)=' 0 0 0 1 0 0 181017.8207 -120692 18741486'// 1 ' -2723' x3(2)=' 0 0 0 1 -1 0 -4077.2629 2713 -422134'// 1 ' -271818' x3(3)=' 0 1 1 0 0 1 -1481.2817 6135 -651'// 1 ' 114' x3(4)=' 0 0 0 1 1 0 1360.2586 -903 140832'// 1 ' 90654' x3(5)=' 0 1 1 0 0 2 1054.3326 10180 -1477'// 1 ' 293' x3(6)=' 0 1 -1 0 0 1 -457.2263 1074 -94'// 1 ' 14' x3(7)=' 0 0 0 0 1 0 0.0074 0 0'// 1 ' 0' x3(8)=' 0 1 -1 0 0 2 -110.8858 239 -19'// 1 ' 2' x3(9)=' 0 1 1 0 0 0 -97.1704 265 -24'// 1 ' 4' x3(10)=' 0 3 1 0 0 3 55.2317 0 -4'// 1 ' 0' c data v1 v1(1)=' 0 1 1 0 0 1-184667147.923 -61587324 3570444'// 1 ' 2774147' v1(2)=' 0 1 1 0 -1 1 1606.495 2670 -278'// 1 ' 107003' v1(3)=' 0 1 1 0 1 1 -2771993.617 -927369 53890'// 1 '-184737001' v1(4)=' 0 1 1 -2 0 1 17230.055 5772 3568097'// 1 ' -255' v1(5)=' 0 1 1 0 2 1 -46798.873 -15701 884'// 1 ' -6238465' v1(6)=' 0 1 1 0 0 0 -8832.492 6011 177'// 1 ' 157' v1(7)=' 0 1 1 0 0 2 8838.522 -6027 -197'// 1 ' -118' v1(8)=' 0 2 2 -1 0 3 22776.763 163962 649508'// 1 ' 1773' v1(9)=' 0 1 3 0 -1 1 -12763.760 -136428 9848'// 1 ' -842940' v1(10)=' 0 0 0 1 0 -1 -22555.148 -161818 -626771'// 1 ' -1714' v1(11)=' 0 1 -1 0 0 1 5700.884 -3171 -19'// 1 ' 6279' v1(12)=' 0 3 1 0 0 3 5697.288 -3773 -591'// 1 ' -78' v1(13)=' 0 1 1 0 -2 1 5221.378 1770 -97'// 1 ' 696244' v1(14)=' 0 1 1 1 0 1 0.229 0 0'// 1 ' 0' v1(15)=' 0 2 2 -1 0 2 2846.423 -8794 301896'// 1 ' -197' v1(16)=' 0 0 0 1 0 0 -2846.188 8820 -301833'// 1 ' 157' v1(17)=' 2 -3 -3 0 0 -3 5485.491 2702 -281'// 1 ' -140' v1(18)=' 2 -1 -1 0 0 -1 -667.307 1195 193'// 1 ' 7' c data v1p v1p(1)=' 1 0 0 0 0 0 -532.706 -1880 153'// 1 ' -79' v1p(2)=' 0 1 3 0 -2 1 2677.804 39842 -1221'// 1 ' 212926' v1p(3)=' 0 1 -1 0 2 1 -2486.738 -37783 1083'// 1 ' -187395' v1p(4)=' 0 3 1 0 0 4 1248.008 -785 0'// 1 ' 0' v1p(5)=' 0 1 -1 0 0 2 1237.287 -748 118'// 1 ' 511' v1p(6)=' 1 -2 -2 0 0 -2 -1721.139 -8128 499'// 1 ' 66' v1p(7)=' 0 1 1 0 0 -1 -617.237 405 0'// 1 ' 0' v1p(8)=' 0 1 1 0 0 3 616.876 -428 0'// 1 ' 0' v1p(9)=' 0 3 1 0 0 2 -543.548 343 0'// 1 ' 39' v1p(10)=' 0 1 -1 0 0 0 -540.379 345 -19'// 1 ' 177' v1p(11)=' 0 1 3 0 -1 0 -1054.707 -2316 196'// 1 ' -70261' v1p(12)=' 0 2 0 1 0 2 517.078 -410 53482'// 1 ' 0' v1p(13)=' 3 -4 -4 0 0 -4 -845.723 -5333 60'// 1 ' 302' v1p(14)=' 0 1 1 0 3 1 -832.288 -280 0'// 1 ' -166394' c data v2 v2(1)=' 0 1 1 0 0 1 184667155.134 61587311 -3570720'// 1 ' -2773871' v2(2)=' 0 1 1 0 -1 1 -1606.803 -2670 278'// 1 ' -107023' v2(3)=' 0 1 1 0 1 1 2771993.601 927369 -53890'// 1 ' 184737001' v2(4)=' 0 1 1 -2 0 1 -17237.270 -5759 -3569615'// 1 ' 255' v2(5)=' 0 1 1 0 2 1 46798.873 15701 -884'// 1 ' 6238465' v2(6)=' 0 1 1 0 0 0 8834.010 -6030 -177'// 1 ' -78' v2(7)=' 0 1 1 0 0 2 -8834.307 6051 98'// 1 ' 118' v2(8)=' 0 2 2 -1 0 3 -22776.766 -163962 -649508'// 1 ' -1773' v2(9)=' 0 1 3 0 -1 1 12763.760 136428 -9848'// 1 ' 842940' v2(10)=' 0 0 0 1 0 -1 22592.199 161756 630613'// 1 ' 1714' v2(11)=' 0 3 1 0 0 3 -5696.686 3774 611'// 1 ' 78' v2(12)=' 0 1 -1 0 0 1 5683.941 -4349 -1161'// 1 ' -6456' v2(13)=' 0 1 1 1 0 1 -0.229 0 0'// 1 ' 0' v2(14)=' 0 1 1 0 -2 1 -5197.382 -1766 97'// 1 ' -693055' v2(15)=' 0 0 0 1 0 0 3410.371 -10581 362440'// 1 ' -275' v2(16)=' 0 2 2 -1 0 2 -2846.415 8794 -301896'// 1 ' 197' v2(17)=' 2 -3 -3 0 0 -3 5485.491 2702 -281'// 1 ' -140' v2(18)=' 2 -1 -1 0 0 -1 667.306 -1195 -193'// 1 ' -7' c data v2p v2p(1)=' 1 0 0 0 0 0 532.706 1880 -153'// 1 ' 79' v2p(2)=' 0 1 3 0 -2 1 -2677.794 -39842 1221'// 1 ' -212926' v2p(3)=' 0 1 -1 0 2 1 2487.540 37784 -1083'// 1 ' 187513' v2p(4)=' 0 1 -1 0 0 2 1255.765 -827 -157'// 1 ' -511' v2p(5)=' 0 3 1 0 0 4 -1248.014 785 0'// 1 ' 0' v2p(6)=' 1 -2 -2 0 0 -2 -1721.139 -8128 499'// 1 ' 66' v2p(7)=' 0 1 1 0 0 -1 617.168 -405 0'// 1 ' 0' v2p(8)=' 0 1 1 0 0 3 -617.663 408 0'// 1 ' 0' v2p(9)=' 0 1 -1 0 0 0 -547.417 358 39'// 1 ' -118' v2p(10)=' 0 3 1 0 0 2 543.495 -343 0'// 1 ' -39' v2p(11)=' 0 2 0 1 0 2 -538.076 389 -55652'// 1 ' 0' v2p(12)=' 0 1 3 0 -1 0 1054.707 2316 -196'// 1 ' 70261' v2p(13)=' 3 -4 -4 0 0 -4 -845.723 -5333 60'// 1 ' 302' v2p(14)=' 0 1 1 0 3 1 832.288 280 0'// 1 ' 166394' c data v3 v3(1)=' 0 0 0 1 0 0 3567800.464 1190831369388343'// 1 ' -53586' v3(2)=' 0 1 1 0 0 1 -29184.285 91687 -12826'// 1 ' 2246' v3(3)=' 0 0 0 1 1 0 53599.765 18017 5549359'// 1 ' 3572140' v3(4)=' 0 1 1 0 0 2 20782.155 221432 -29113'// 1 ' 5775' v3(5)=' 0 1 -1 0 0 1 8999.931 -12132 1850'// 1 ' -275' v3(6)=' 0 0 0 0 1 0 0.146 0 0'// 1 ' 0' v3(7)=' 0 1 -1 0 0 2 2181.635 -2517 373'// 1 ' -39' v3(8)=' 0 1 1 0 0 0 -1913.567 3304 -472'// 1 ' 78' v3(9)=' 0 3 1 0 0 3 1089.188 1088 -78'// 1 ' 0' v3(10)=' 2 -2 -2 0 0 -2 13.355 19 0'// 1 ' 0' v3(11)=' 0 1 -1 0 0 0 -664.567 1146 -177'// 1 ' 19' v3(12)=' 0 0 0 1 2 0 905.431 314 93730'// 1 ' 120442' c data y1 y1(1)=' -0.0001 0 0 0 0 -56' y1(2)=' -0.0002 0 0 0 0 1' y1(3)=' 0.0000 0 0 0 0 0' y1(4)=' -0.0028 0 0 0 0 -11' y1(5)=' 0.0000 0 0 0 0 0' y1(6)=' -794.4492 3-76760 -34 0 0' y1(7)=' 0.0000 0 0 0 0 0' y1(8)=' 0.0011 0 0 0 0 -24' y1(9)=' 0.0011 0 0 0 0 23' y1(10)=' 0.1201 0 0 0 0 -78' y1(11)=' -0.0009 0 0 0 0 3' y1(12)=' 0.1199 0 0 0 0 77' y1(13)=' 0.0001 0 0 0 0 -15' y1(14)=' -0.0001 0 0 0 0 15' y1(15)=' 0.0000 0 0 0 0 0' y1(16)=' -264.4599 1-25560 16 0 0' y1(17)=' -0.0079 0 0 0 0 -6956' y1(18)=' -0.0104 0 0 0 0 6953' c data y1p y1p(1)=' -0.0475 0 0 0 0 0' y1p(2)=' 6.3708 0 0 0 0 0' y1p(3)=' -8.6457 -22 1 -1 0 0' y1p(4)=' 0.0103 0 0 0 0 -1' y1p(5)=' 0.0000 0 0 0 0 -1' y1p(6)=' 0.0103 0 0 0 0 1' y1p(7)=' 0.0000 0 0 0 0 2' y1p(8)=' 0.0000 0 0 0 0 -2' y1p(9)=' 4.9764 13 -1 0 0 0' c data y2 y2(1)=' 0.0001 0 0 0 0 -73' y2(2)=' -0.0001 0 0 0 0 2' y2(3)=' 0.0000 0 0 0 0 0' y2(4)=' 0.0029 0 0 0 0 -11' y2(5)=' 0.0000 0 0 0 0 0' y2(6)=' 794.4491 -3 76760 34 0 0' y2(7)=' 0.0000 0 0 0 0 0' y2(8)=' -0.0011 0 0 0 0 -25' y2(9)=' -0.0011 0 0 0 0 35' y2(10)=' -0.1201 0 0 0 0 -78' y2(11)=' 0.0009 0 0 0 0 3' y2(12)=' -0.1199 0 0 0 0 77' y2(13)=' 0.0001 0 0 0 0 16' y2(14)=' 0.0000 0 0 0 0 10' y2(15)=' 264.4599 -1 25560 -16 0 0' y2(16)=' 0.0000 0 0 0 0 0' y2(17)=' 0.0079 0 0 0 0 8947' y2(18)=' 0.0079 0 0 0 0 -6956' c data y2p y2p(1)=' -0.0475 0 0 0 0 0' y2p(2)=' -6.3708 0 0 0 0 0' y2p(3)=' 8.6457 22 -1 1 0 0' y2p(4)=' -0.0103 0 0 0 0 -1' y2p(5)=' 0.0000 0 0 0 0 -1' y2p(6)=' -0.0103 0 0 0 0 1' y2p(7)=' 0.0000 0 0 0 0 1' y2p(8)=' 0.0000 0 0 0 0 2' y2p(9)=' 4.9764 13 -1 0 0 0' c data y3 y3(1)=' 0.0000 0 0 0 0 -3' y3(2)=' 0.0000 0 0 0 0 0' y3(3)=' 0.1294 0 0 0 0 -103296' y3(4)=' 0.0000 0 0 0 0 0' y3(5)=' 0.0184 0 0 0 0 106' y3(6)=' 0.0000 0 0 0 0 5' y3(7)=' 411.3091 0 -40 25602 0 0' y3(8)=' 0.0000 0 0 0 0 1' y3(9)=' 0.0000 0 0 0 0 -2' y3(10)=' 0.0000 0 0 0 0 2' c data w1 w1(1)=' -0.002 0 0 0 0 1103' w1(2)=' 0.000 0 0 0 0 -37' w1(3)=' 0.000 0 0 0 0 -12' w1(4)=' 0.000 0 0 0 0 0' w1(5)=' 0.000 0 0 0 0 0' w1(6)=' 0.022 0 0 0 0 472' w1(7)=' 0.022 0 0 0 0 -453' w1(8)=' 2.366 0 0 0 0 1536' w1(9)=' -0.035 0 0 0 0 -118' w1(10)=' 2.362 0 0 0 0 -1516' w1(11)=' -0.002 0 0 0 0 -295' w1(12)=' -0.002 0 0 0 0 -295' w1(13)=' 0.000 0 0 0 0 0' w1(14)=' -10422.818-10386-1007363 630 0 0' w1(15)=' -0.156 0 0 0 0 136994' w1(16)=' -0.205 0 0 0 0 -137041' w1(17)=' 2.226 2 0 0 0 0' w1(18)=' -47.485 -125 0 0 0 0' c data w1p w1p(1)=' -52.948 -134 6 -6 0 0' w1p(2)=' 0.203 0 0 0 0 19' w1p(3)=' 0.203 0 0 0 0 -19' w1p(4)=' 0.000 0 0 0 0 -39' w1p(5)=' 0.000 0 0 0 0 -39' w1p(6)=' -165.614 -628 33 0 0 0' w1p(7)=' 0.000 0 0 0 0 0' w1p(8)=' 0.000 0 0 0 0 0' w1p(9)=' 0.000 0 0 0 0 59' w1p(10)=' 0.000 0 0 0 0 59' w1p(11)=' 0.000 0 0 0 0 0' w1p(12)=' 0.000 0 0 0 0 0' w1p(13)=' 135.835 841 0 0 0 0' w1p(14)=' 0.000 0 0 0 0 0' c data w2 w2(1)=' -0.002 0 0 0 0 -1438' w2(2)=' 0.000 0 0 0 0 37' w2(3)=' 0.000 0 0 0 0 12' w2(4)=' 0.000 0 0 0 0 0' w2(5)=' 0.000 0 0 0 0 0' w2(6)=' 0.022 0 0 0 0 -492' w2(7)=' 0.022 0 0 0 0 689' w2(8)=' 2.366 0 0 0 0 -1536' w2(9)=' -0.035 0 0 0 0 118' w2(10)=' 2.362 0 0 0 0 1516' w2(11)=' -0.002 0 0 0 0 315' w2(12)=' 0.000 0 0 0 0 -196' w2(13)=' -10422.818-10386-1007363 630 0 0' w2(14)=' 0.000 0 0 0 0 0' w2(15)=' -0.156 0 0 0 0 176342' w2(16)=' -0.156 0 0 0 0 -136994' w2(17)=' -2.226 -2 0 0 0 0' w2(18)=' -47.485 -125 0 0 0 0' c data w2p w2p(1)=' -52.948 -134 6 -6 0 0' w2p(2)=' 0.203 0 0 0 0 -19' w2p(3)=' 0.203 0 0 0 0 19' w2p(4)=' 0.000 0 0 0 0 -19' w2p(5)=' 0.000 0 0 0 0 39' w2p(6)=' 165.614 628 -33 0 0 0' w2p(7)=' 0.000 0 0 0 0 0' w2p(8)=' 0.000 0 0 0 0 0' w2p(9)=' 0.000 0 0 0 0 39' w2p(10)=' 0.000 0 0 0 0 -59' w2p(11)=' 0.000 0 0 0 0 0' w2p(12)=' 0.000 0 0 0 0 0' w2p(13)=' -135.835 -841 0 0 0 0' w2p(14)=' 0.000 0 0 0 0 0' c data w3 w3(1)=' 0.000 0 0 0 0 -26' w3(2)=' -2.549 -2 0 0 0 -2035142' w3(3)=' 0.000 0 0 0 0 0' w3(4)=' -0.363 0 0 0 0 2089' w3(5)=' 0.000 0 0 0 0 -98' w3(6)=' -8100.500 -8096 787-504217 0 0' w3(7)=' 0.000 0 0 0 0 -19' w3(8)=' 0.000 0 0 0 0 -39' w3(9)=' 0.000 0 0 0 0 39' w3(10)=' -1306.279 -1135 271 -108 0 0' w3(11)=' 0.000 0 0 0 0 0' w3(12)=' 0.000 0 0 0 0 0' c **************************************************************** c * perturbations par deimos et les planetes * c **************************************************************** data asup/183.1130d0,199.2812d0,255.194d0,271.3608d0/ data fsup/1128.826758d0,-1128.862761d0,1128.849665d0, 1 -1128.839842d0/ data rx1s/51.71d0,-51.71d0,63.305d0,-63.208d0/ data rx2s/51.71d0,51.71d0,63.305d0,63.209d0/ c **************************************************************** c * moyens mouvements esapho * c **************************************************************** data fipsi/350.8919885d0/ data fipi/0.1772311d-4/ data fid/1128.3207210d0/ data fif/1129.280784d0/ data fil/1128.409439d0/ data filp/0.5240207d0/ data fnu/1128.84426d0/ c **************************************************************** c * arguments en j2000 * c **************************************************************** ai(1)=psi0 ai(2)=pipet ai(3)=al(1)-pipet-alm ai(4)=al(1)-al(2) ai(5)=al(1)-al(3) ai(6)=alm c **************************************************************** dpn=dp*1.d5/fnu dnun=dnu/fnu pf(1)=fipsi pf(2)=fipi pf(3)=fid pf(4)=fif pf(5)=fil pf(6)=filp c **************************************************************** c * conversion en radians * c **************************************************************** deg=3600.d0/206264.8062470964d0 do 100 k=1,6 f(k)=pf(k)*deg 100 a0(k)=ai(k)*deg c **************************************************************** c * lecture et transformation des series position esapho. * c **************************************************************** do 1 i=1,3 jmax=nb(i) do 2 j=1,jmax read(x(j,i),1001)(iar(k),k=1,6),r1,(ic1(k),k=1,3) read(y(j,i),1002)r2,ic2,ic1(4) aw=0 fw=0 do 3 k=1,6 aw=aw+iar(k)*a0(k) 3 fw=fw+iar(k)*f(k) arg(j,i)=aw freq(j,i)=fw freq2(j,i)=iar(3)+iar(4)+iar(5) r1=r1+ic1(1)*dnun+ic1(2)*dgam+ic1(3)*de+ic1(4)*dpn r2=r2+ic2(1)*dnun+ic2(2)*dgam+ic2(3)*de+ic2(4)*dpn if(i.eq.1) then cs(j,i)=r2*1.d-3 cc(j,i)=r1*1.d-3 else cs(j,i)=r1*1.d-3 cc(j,i)=r2*1.d-3 endif 2 continue if(i.ne.3) then do 4 j=1,4 jp=jmax+j arg(jp,i)=asup(j)*deg freq(jp,i)=fsup(j)*deg freq2(jp,i)=0 if(i.eq.1) then cs(jp,i)=0 cc(jp,i)=rx1s(j)*1.d-3 else cs(jp,i)=rx2s(j)*1.d-3 cc(jp,i)=0 endif 4 continue endif 1 continue 1001 format(6i3,f13.4,i8,2i9) 1002 format(f14.4,4i6,i13) c **************************************************************** c * lecture et transformation des series vitesse esapho. * c **************************************************************** do 11 i=1,3 jmax=nbv(i) do 12 j=1,jmax read(v(j,i),2001)(iar(k),k=1,6),r1,(ic1(k),k=1,3) read(w(j,i),2002) r2,ic2,ic1(4) aw=0 fw=0 do 13 k=1,6 aw=aw+iar(k)*a0(k) 13 fw=fw+iar(k)*f(k) argv(j,i)=aw freqv(j,i)=fw frev2(j,i)=iar(3)+iar(4)+iar(5) r1=r1+ic1(1)*dnun+ic1(2)*dgam+ic1(3)*de+ic1(4)*dpn r2=r2+ic2(1)*dnun+ic2(2)*dgam+ic2(3)*de+ic2(4)*dpn if(i.ne.1) then cvs(j,i)=r2*1.d-3 cvc(j,i)=r1*1.d-3 else cvs(j,i)=r1*1.d-3 cvc(j,i)=r2*1.d-3 endif 12 continue if(i.ne.3) then do 14 j=1,4 jp=jmax+j argv(jp,i)=asup(j)*deg freqv(jp,i)=fsup(j)*deg frev2(jp,i)=0 if(i.eq.1) then cvs(jp,i)=-rx1s(j)*1.d-3*argv(jp,i) cvc(j,i)=0 else cvs(jp,i)=0 cvc(jp,i)=rx2s(j)*1.d-3*argv(jp,i) endif 14 continue endif 11 continue 2001 format(6i3,f14.3,i10,i9,i10) 2002 format(f15.3,i6,i8,i7,i6,i14) return end subroutine lecd c *************************************************************** c * lecture des series deimos. introduction des constantes * c * nouvelles dans les moyens mouvements et les coefficients. * c * transformation des arguments en fonctions du temps * c *************************************************************** implicit double precision(a-h,o-z) save parameter(id=38,idv=27) dimension f(6),a0(6),ai(6),pf(6),al(3) integer nb(3),nbv(3),iar(6),ic1(3) dimension arg(38,3),freq(38,3),cs(38,3),cc(38,3), 1 freq2(38,3),argv(36,3),freqv(36,3),cvs(36,3),cvc(36,3), 2 frev2(36,3) character x1(18)*59,x2(19)*59,x3(18)*59,y1(18)*10, 1 y2(19)*10,y3(18)*10,v1(18)*63,v2(18)*63,v3(16)*63, 2 w1(18)*11,w2(18)*11,w3(16)*11, 3 x1p(19)*59,x2p(19)*59,x3p(7)*59,y1p(19)*10,y2p(19)*10, 4 y3p(7)*10,v1p(9)*63,v2p(9)*63,w1p(9)*11,w2p(9)*11, 5 x(38,3)*59,v(27,3)*63,y(38,3)*10,w(27,3)*11 equivalence (x(1,1),x1(1)),(y(1,1),y1(1)),(x(19,1),x1p(1)), 1 (y(19,1),y1p(1)),(x(1,2),x2(1)),(y(1,2),y2(1)), 2 (x(20,2),x2p(1)),(y(20,2),y2p(1)),(x(19,3),x3p(1)), 3 (x(1,3),x3(1)),(y(1,3),y3(1)),(y(19,3),y3p(1)),(v(1,1),v1(1)), 4 (w(1,1),w1(1)),(v(19,1),v1p(1)),(w(19,1),w1p(1)), 5 (v(1,2),v2(1)),(w(1,2),w2(1)),(v(19,2),v2p(1)), 6 (w(19,2),w2p(1)),(v(1,3),v3(1)),(w(1,3),w3(1)) common/ed5/arg,freq,cs,cc,freq2,argv,freqv,cvs,cvc,frev2 common/ed4/nb,nbv common/ed1/dnu,dgam,de,al common/ep3/psi0,pipet,alm,dp common/ed2/accep c ***************************************************************** c * coefficients deimos * c **************************************************************** c data x1 x1(1)=' 0 1 1 0 0 1 23451762.1311-15634508 -727543'// 1 ' -9779' x1(2)=' 0 2 2 -1 0 2 -55151.4362 36767 0'// 1 ' 0' x1(3)=' 0 0 0 1 0 0 55133.9035 -36755 0'// 1 ' 0' x1(4)=' 0 1 1 0 -1 1 -14544.5209 9696 455'// 1 '-35177119' x1(5)=' 0 1 1 -2 0 1 5602.8085 -3735 727543'// 1 ' -2' x1(6)=' 0 1 1 0 1 1 4890.3869 -3260 -151'// 1 ' 11725701' x1(7)=' 0 1 1 0 0 0 4551.1609 -3034 0'// 1 ' 0' x1(8)=' 0 1 1 0 0 2 -4545.4296 3030 0'// 1 ' 0' x1(9)=' 2 -3 -3 0 0 -3 -3251.3962 2167 0'// 1 ' 0' x1(10)=' 0 3 1 0 0 3 -2659.0343 1772 0'// 1 ' 0' x1(11)=' 0 1 -1 0 0 1 2592.0639 -1728 0'// 1 ' 0' x1(12)=' 2 -1 -1 0 0 -1 1868.9622 -1245 0'// 1 ' 0' x1(13)=' 0 1 1 -1 0 1 0.0000 0 0'// 1 ' 0' x1(14)=' 0 3 1 0 0 4 -581.0963 387 0'// 1 ' 0' x1(15)=' 0 1 -1 0 0 2 558.8118 -372 0'// 1 ' 0' x1(16)=' 0 2 2 0 0 2 0.0000 0 0'// 1 ' 0' x1(17)=' 0 2 0 1 0 2 -390.0656 260 0'// 1 ' 0' x1(18)=' 0 1 1 1 0 1 0.0000 0 0'// 1 ' 0' c data x1p x1p(1)=' 3 -4 -4 0 0 -4 280.7530 -187 0'// 1 ' 0' x1p(2)=' 0 1 1 0 0 -1 318.5942 -212 0'// 1 ' 0' x1p(3)=' 0 1 1 0 0 3 -316.7712 211 0'// 1 ' 0' x1p(4)=' 0 0 2 -1 0 0 288.9456 -192 0'// 1 ' 0' x1p(5)=' 0 3 1 0 0 2 255.1685 -170 0'// 1 ' 0' x1p(6)=' 0 1 -1 0 0 0 -252.9696 168 0'// 1 ' 0' x1p(7)=' 1 -2 -2 0 0 -2 157.2017 -104 0'// 1 ' 0' x1p(8)=' 0 3 3 -2 0 3 148.1044 -98 0'// 1 ' 0' x1p(9)=' 3 -2 -2 0 0 -2 -126.2662 84 0'// 1 ' 0' x1p(10)=' 1 0 0 0 0 0 -119.1380 79 0'// 1 ' 0' x1p(11)=' 0 0 0 1 0 1 125.3480 -83 0'// 1 ' 0' x1p(12)=' 0 0 0 1 0 -1 -118.4993 78 0'// 1 ' 0' x1p(13)=' 0 3 1 0 0 5 -98.7740 65 0'// 1 ' 0' x1p(14)=' 0 2 2 -1 0 1 -94.0167 62 0'// 1 ' 0' x1p(15)=' 0 1 -1 0 0 3 93.6378 -62 0'// 1 ' 0' x1p(16)=' 0 2 2 -1 0 3 86.4433 -57 0'// 1 ' 0' x1p(17)=' 0 2 0 -1 0 2 85.1187 -56 0'// 1 ' 0' x1p(18)=' 0 2 0 1 0 3 -84.9226 56 0'// 1 ' 0' x1p(19)=' 0 0 2 -1 0 -1 62.9106 -41 0'// 1 ' 0' c data x2 x2(1)=' 0 1 1 0 0 1 23448856.2674-15632570 -727543'// 1 ' -9779' x2(2)=' 0 0 0 1 0 0 66507.3861 -44338 0'// 1 ' 0' x2(3)=' 0 2 2 -1 0 2 -55124.5490 36749 0'// 1 ' 0' x2(4)=' 0 1 1 0 -1 1 -14794.1622 9862 455'// 1 '-35177119' x2(5)=' 0 1 1 -2 0 1 5652.0343 -3768 727543'// 1 ' -2' x2(6)=' 0 1 1 0 1 1 4889.7715 -3259 -151'// 1 ' 11725701' x2(7)=' 0 1 1 0 0 0 4526.5150 -3017 0'// 1 ' 0' x2(8)=' 0 1 1 0 0 2 -4514.4268 3009 0'// 1 ' 0' x2(9)=' 2 -3 -3 0 0 -3 3250.9885 -2167 0'// 1 ' 0' x2(10)=' 0 1 -1 0 0 1 -2843.1994 1895 0'// 1 ' 0' x2(11)=' 0 3 1 0 0 3 -2653.7375 1769 0'// 1 ' 0' x2(12)=' 2 -1 -1 0 0 -1 1868.6828 -1245 0'// 1 ' 0' x2(13)=' 0 0 0 0 0 0 0.0000 0 0'// 1 ' 0' x2(14)=' 0 1 1 -1 0 1 0.0000 0 0'// 1 ' 0' x2(15)=' 0 1 -1 0 0 2 -632.9468 421 0'// 1 ' 0' x2(16)=' 0 3 1 0 0 4 -579.9308 386 0'// 1 ' 0' x2(17)=' 0 2 0 1 0 2 -399.2736 266 0'// 1 ' 0' x2(18)=' 0 2 2 0 0 2 0.0000 0 0'// 1 ' 0' x2(19)=' 0 1 1 1 0 1 0.0000 0 0'// 1 ' 0' c data x2p x2p(1)=' 3 -4 -4 0 0 -4 -280.7182 187 0'// 1 ' 0' x2p(2)=' 0 1 1 0 0 -1 316.9349 -211 0'// 1 ' 0' x2p(3)=' 0 1 1 0 0 3 -314.5569 209 0'// 1 ' 0' x2p(4)=' 0 0 2 -1 0 0 294.7831 -196 0'// 1 ' 0' x2p(5)=' 0 1 -1 0 0 0 269.2981 -179 0'// 1 ' 0' x2p(6)=' 0 3 1 0 0 2 254.6917 -169 0'// 1 ' 0' x2p(7)=' 1 -2 -2 0 0 -2 -157.1732 104 0'// 1 ' 0' x2p(8)=' 0 3 3 -2 0 3 148.0042 -98 0'// 1 ' 0' x2p(9)=' 3 -2 -2 0 0 -2 -126.2468 84 0'// 1 ' 0' x2p(10)=' 1 0 0 0 0 0 -119.1139 79 0'// 1 ' 0' x2p(11)=' 0 1 -1 0 0 3 -109.5100 73 0'// 1 ' 0' x2p(12)=' 0 2 0 -1 0 2 106.0245 -70 0'// 1 ' 0' x2p(13)=' 0 3 1 0 0 5 -98.5724 65 0'// 1 ' 0' x2p(14)=' 0 2 2 -1 0 1 -93.6819 62 0'// 1 ' 0' x2p(15)=' 0 2 0 1 0 3 -86.9127 57 0'// 1 ' 0' x2p(16)=' 0 2 2 -1 0 3 86.0623 -57 0'// 1 ' 0' x2p(17)=' 0 0 0 1 0 1 70.1733 -46 0'// 1 ' 0' x2p(18)=' 0 0 0 1 0 -1 -66.2919 44 0'// 1 ' 0' x2p(19)=' 0 0 2 -1 0 -1 64.1928 -42 0'// 1 ' 0' c data x3 x3(1)=' 0 0 0 1 0 0 726494.5032 -484329 46897178'// 1 ' -303' x3(2)=' 0 1 1 0 0 1 -367660.3813 245106 0'// 1 ' 0' x3(3)=' 0 1 -1 0 0 1 -6081.1262 4054 0'// 1 ' 0' x3(4)=' 0 2 2 -1 0 2 2149.8425 -1433 0'// 1 ' 0' x3(5)=' 0 1 1 0 0 2 1844.8540 -1229 0'// 1 ' 0' x3(6)=' 0 1 1 -2 0 1 -1799.3963 1199 0'// 1 ' 0' x3(7)=' 0 1 1 0 0 0 -1729.1650 1152 0'// 1 ' 0' x3(8)=' 0 1 -1 0 0 2 -1339.1617 892 0'// 1 ' 0' x3(9)=' 0 1 -1 0 0 0 574.0342 -382 0'// 1 ' 0' x3(10)=' 0 0 0 1 -1 0 -454.8726 303 -29334'// 1 ' -1091184' x3(11)=' 0 0 2 -1 0 0 417.5738 -278 0'// 1 ' 0' x3(12)=' 0 3 1 0 0 3 333.6614 -222 0'// 1 ' 0' x3(13)=' 0 1 1 0 -1 1 230.5734 -153 0'// 1 ' 0' x3(14)=' 0 1 -1 0 0 3 -228.7458 152 0'// 1 ' 0' x3(15)=' 0 0 0 1 1 0 151.3386 -100 9778'// 1 ' 363728' x3(16)=' 0 1 1 0 0 3 126.3330 -84 0'// 1 ' 0' x3(17)=' 0 1 1 0 0 -1 -123.0024 82 0'// 1 ' 0' x3(18)=' 2 -2 -2 -1 0 -2 100.4703 -66 0'// 1 ' 0' c data x3p x3p(1)=' 0 0 2 -1 0 -1 92.7574 -61 0'// 1 ' 0' x3p(2)=' 0 1 1 0 1 1 -77.0433 51 0'// 1 ' 0' x3p(3)=' 0 3 1 0 0 4 72.3397 -48 0'// 1 ' 0' x3p(4)=' 2 -2 -2 1 0 -2 63.1447 -42 0'// 1 ' 0' x3p(5)=' 0 0 0 0 0 0 0.0000 0 0'// 1 ' 0' x3p(6)=' 0 1 1 -1 0 1 0.0000 0 0'// 1 ' 0' x3p(7)=' 2 -3 -3 0 0 -3 -51.4568 34 0'// 1 ' 0' c data v1 v1(1)=' 0 1 1 0 0 1 -116719753.045 -38906584 3620991'// 1 ' 48670' v1(2)=' 0 2 2 -1 0 2 274472.052 91511 0'// 1 ' 0' v1(3)=' 0 0 0 1 0 0 -274419.644 -91460 0'// 1 ' 0' v1(4)=' 0 1 1 -2 0 1 27888.799 9293 3621451'// 1 ' -9' v1(5)=' 0 1 1 0 1 1 -48677.506 -16229 1503'// 1 '-116714258' v1(6)=' 0 1 1 0 0 0 -22609.568 -7578 0'// 1 ' 0' v1(7)=' 0 1 1 0 0 2 22664.240 7514 0'// 1 ' 0' v1(8)=' 2 -3 -3 0 0 -3 -8722.200 -42733 0'// 1 ' 0' v1(9)=' 0 3 1 0 0 3 13282.692 4382 0'// 1 ' 0' v1(10)=' 0 1 -1 0 0 1 12853.324 4332 0'// 1 ' 0' v1(11)=' 2 -1 -1 0 0 -1 -13590.022 18354 0'// 1 ' 0' v1(12)=' 0 3 1 0 0 4 2908.069 955 0'// 1 ' 0' v1(13)=' 0 1 -1 0 0 2 2765.881 939 0'// 1 ' 0' v1(14)=' 0 2 2 0 0 2 0.000 0 0'// 1 ' 0' v1(15)=' 0 2 0 1 0 2 1948.621 642 0'// 1 ' 0' v1(16)=' 0 1 1 1 0 1 0.000 0 0'// 1 ' 0' v1(17)=' 0 1 1 0 0 -1 -1579.820 -534 0'// 1 ' 0' v1(18)=' 0 1 1 0 0 3 1582.369 522 0'// 1 ' 0' c data v1p v1p(1)=' 0 0 2 -1 0 0 -1432.710 -486 0'// 1 ' 0' v1p(2)=' 0 3 1 0 0 2 -1272.311 -422 0'// 1 ' 0' v1p(3)=' 0 1 -1 0 0 0 -1256.720 -424 0'// 1 ' 0' v1p(4)=' 1 -2 -2 0 0 -2 602.052 1166 0'// 1 ' 0' v1p(5)=' 0 3 3 -2 0 3 -737.024 -249 0'// 1 ' 0' v1p(6)=' 3 -2 -2 0 0 -2 1062.989 -1964 0'// 1 ' 0' v1p(7)=' 1 0 0 0 0 0 729.627 -483 0'// 1 ' 0' v1p(8)=' 0 0 0 1 0 1 -625.045 -209 0'// 1 ' 0' v1p(9)=' 0 0 0 1 0 -1 588.726 202 0'// 1 ' 0' c data v2 v2(1)=' 0 1 1 0 0 1 116705290.520 38901767 -3620991'// 1 ' -48670' v2(2)=' 0 0 0 1 0 0 331029.222 110323 0'// 1 ' 0' v2(3)=' 0 2 2 -1 0 2 -274338.243 -91466 0'// 1 ' 0' v2(4)=' 0 1 1 -2 0 1 -28133.827 -9374 -3621451'// 1 ' 9' v2(5)=' 0 1 1 0 1 1 48671.381 16233 -1503'// 1 ' 116714258' v2(6)=' 0 1 1 0 0 0 22487.130 7540 0'// 1 ' 0' v2(7)=' 0 1 1 0 0 2 -22509.655 -7465 0'// 1 ' 0' v2(8)=' 2 -3 -3 0 0 -3 -8721.107 -42727 0'// 1 ' 0' v2(9)=' 0 1 -1 0 0 1 14098.635 4753 0'// 1 ' 0' v2(10)=' 0 3 1 0 0 3 -13256.232 -4370 0'// 1 ' 0' v2(11)=' 2 -1 -1 0 0 -1 13587.990 -18353 0'// 1 ' 0' v2(12)=' 0 1 -1 0 0 2 3132.818 1066 0'// 1 ' 0' v2(13)=' 0 3 1 0 0 4 -2902.236 -954 0'// 1 ' 0' v2(14)=' 0 2 0 1 0 2 -1994.620 -658 0'// 1 ' 0' v2(15)=' 0 2 2 0 0 2 0.000 0 0'// 1 ' 0' v2(16)=' 0 1 1 1 0 1 0.000 0 0'// 1 ' 0' v2(17)=' 0 1 1 0 0 -1 1571.592 531 0'// 1 ' 0' v2(18)=' 0 1 1 0 0 3 -1571.308 -521 0'// 1 ' 0' c data v2p v2p(1)=' 0 0 2 -1 0 0 1461.654 495 0'// 1 ' 0' v2p(2)=' 0 1 -1 0 0 0 -1337.837 -451 0'// 1 ' 0' v2p(3)=' 0 3 1 0 0 2 1269.934 424 0'// 1 ' 0' v2p(4)=' 1 -2 -2 0 0 -2 601.943 1166 0'// 1 ' 0' v2p(5)=' 0 3 3 -2 0 3 736.525 248 0'// 1 ' 0' v2p(6)=' 3 -2 -2 0 0 -2 -1062.825 1963 0'// 1 ' 0' v2p(7)=' 1 0 0 0 0 0 -729.480 483 0'// 1 ' 0' v2p(8)=' 0 1 -1 0 0 3 541.026 184 0'// 1 ' 0' v2p(9)=' 0 2 0 -1 0 2 -525.780 -180 0'// 1 ' 0' c data v3 v3(1)=' 0 0 0 1 0 0 3616003.041 1205108 23342274'// 1 '1 -1508' v3(2)=' 0 1 1 0 0 1 -1829850.937 -609954 0'// 1 ' 0' v3(3)=' 0 1 -1 0 0 1 30154.613 10163 0'// 1 ' 0' v3(4)=' 0 2 2 -1 0 2 10699.117 3568 0'// 1 ' 0' v3(5)=' 0 1 1 0 0 2 9198.737 3053 0'// 1 ' 0' v3(6)=' 0 1 1 -2 0 1 8956.758 2987 0'// 1 ' 0' v3(7)=' 0 1 1 0 0 0 -8590.264 -2883 0'// 1 ' 0' v3(8)=' 0 1 -1 0 0 2 6628.283 2250 0'// 1 ' 0' v3(9)=' 0 1 -1 0 0 0 -2851.726 -959 0'// 1 ' 0' v3(10)=' 0 0 2 -1 0 0 2070.500 699 0'// 1 ' 0' v3(11)=' 0 3 1 0 0 3 1666.741 551 0'// 1 ' 0' v3(12)=' 0 1 -1 0 0 3 1130.103 387 0'// 1 ' 0' v3(13)=' 0 0 0 1 1 0 1506.429 511 97330'// 1 ' 3620558' v3(14)=' 0 1 1 0 0 3 631.072 209 0'// 1 ' 0' v3(15)=' 0 1 1 0 0 -1 -609.935 -205 0'// 1 ' 0' v3(16)=' 0 1 1 0 1 1 -766.867 -259 0'// 1 ' 0' c data y1 y1(1)=' 0.0000' y1(2)=' 0.0000' y1(3)=' 0.0000' y1(4)=' 0.0000' y1(5)=' 0.0000' y1(6)=' 0.0000' y1(7)=' 0.0000' y1(8)=' 0.0000' y1(9)=' -0.0020' y1(10)=' 0.0000' y1(11)=' 0.0000' y1(12)=' 0.0007' y1(13)='-1157.0034' y1(14)=' 0.0000' y1(15)=' 0.0000' y1(16)=' 395.4635' y1(17)=' 0.0000' y1(18)=' -383.9143' c data y1p y1p(1)=' 45.6488' y1p(2)=' 0.0000' y1p(3)=' 0.0000' y1p(4)=' 0.0000' y1p(5)=' 0.0000' y1p(6)=' 0.0000' y1p(7)=' -11.6027' y1p(8)=' 0.0000' y1p(9)=' -20.5254' y1p(10)=' 8.8034' y1p(11)=' 0.0000' y1p(12)=' 0.0000' y1p(13)=' 0.0000' y1p(14)=' 0.0000' y1p(15)=' 0.0000' y1p(16)=' 0.0000' y1p(17)=' 0.0000' y1p(18)=' 0.0000' y1p(19)=' 0.0000' c data y2 y2(1)=' 0.0000' y2(2)=' 0.0000' y2(3)=' 0.0000' y2(4)=' 0.0000' y2(5)=' 0.0000' y2(6)=' 0.0000' y2(7)=' 0.0000' y2(8)=' 0.0000' y2(9)=' -0.0020' y2(10)=' 0.0000' y2(11)=' 0.0000' y2(12)=' -0.0007' y2(13)='-1180.2676' y2(14)=' 1156.5679' y2(15)=' 0.0000' y2(16)=' 0.0000' y2(17)=' 0.0000' y2(18)=' -395.4136' y2(19)=' 383.6747' c data y2p y2p(1)=' 45.6431' y2p(2)=' 0.0000' y2p(3)=' 0.0000' y2p(4)=' 0.0000' y2p(5)=' 0.0000' y2p(6)=' 0.0000' y2p(7)=' -11.6006' y2p(8)=' 0.0000' y2p(9)=' 20.5223' y2p(10)=' -8.8016' y2p(11)=' 0.0000' y2p(12)=' 0.0000' y2p(13)=' 0.0000' y2p(14)=' 0.0000' y2p(15)=' 0.0000' y2p(16)=' 0.0000' y2p(17)=' 0.0000' y2p(18)=' 0.0000' y2p(19)=' 0.0000' c data y3 y3(1)=' 0.0000' y3(2)=' 0.0000' y3(3)=' 0.0000' y3(4)=' 0.0000' y3(5)=' 0.0000' y3(6)=' 0.0000' y3(7)=' 0.0000' y3(8)=' 0.0000' y3(9)=' 0.0000' y3(10)=' 0.0000' y3(11)=' 0.0000' y3(12)=' 0.0000' y3(13)=' 0.0000' y3(14)=' 0.0000' y3(15)=' 0.0000' y3(16)=' 0.0000' y3(17)=' 0.0000' y3(18)=' 0.0000' c data y3p y3p(1)=' 0.0000' y3p(2)=' 0.0000' y3p(3)=' 0.0000' y3p(4)=' 0.0000' y3p(5)=' 57.4788' y3p(6)=' -54.7313' y3p(7)=' 0.0000' c data w1 w1(1)=' 0.000' w1(2)=' 0.000' w1(3)=' 0.000' w1(4)=' 0.000' w1(5)=' 0.000' w1(6)=' 0.000' w1(7)=' 0.000' w1(8)=' 0.005' w1(9)=' 0.000' w1(10)=' 0.000' w1(11)=' 0.005' w1(12)=' 0.000' w1(13)=' 0.000' w1(14)=' 3936.455' w1(15)=' 0.000' w1(16)=' -3821.615' w1(17)=' 0.000' w1(18)=' 0.000' c data w1p w1p(1)=' 0.000' w1p(2)=' 0.000' w1p(3)=' 0.000' w1p(4)=' 44.436' w1p(5)=' 0.000' w1p(6)=' -172.796' w1p(7)=' 53.914' w1p(8)=' 0.000' w1p(9)=' 0.000' c data w2 w2(1)=' 0.000' w2(2)=' 0.000' w2(3)=' 0.000' w2(4)=' 0.000' w2(5)=' 0.000' w2(6)=' 0.000' w2(7)=' 0.000' w2(8)=' -0.005' w2(9)=' 0.000' w2(10)=' 0.000' w2(11)=' 0.005' w2(12)=' 0.000' w2(13)=' 0.000' w2(14)=' 0.000' w2(15)=' 3935.958' w2(16)=' -3819.230' w2(17)=' 0.000' w2(18)=' 0.000' c data w2p w2p(1)=' 0.000' w2p(2)=' 0.000' w2p(3)=' 0.000' w2p(4)=' -44.428' w2p(5)=' 0.000' w2p(6)=' -172.770' w2p(7)=' 53.903' w2p(8)=' 0.000' w2p(9)=' 0.000' c data w3 w3(1)=' 0.000' w3(2)=' 0.000' w3(3)=' 0.000' w3(4)=' 0.000' w3(5)=' 0.000' w3(6)=' 0.000' w3(7)=' 0.000' w3(8)=' 0.000' w3(9)=' 0.000' w3(10)=' 0.000' w3(11)=' 0.000' w3(12)=' 0.000' w3(13)=' 0.000' w3(14)=' 0.000' w3(15)=' 0.000' w3(16)=' 0.000' c **************************************************************** c * moyens mouvements deimos * c **************************************************************** data fipsi/350.8919885d0/ data fipi/0.1772311d-4/ data fid/284.6378363d0/ data fif/285.179876d0/ data fil/285.143868d0/ data filp/0.5240207d0/ data fnu/285.161908d0/ c **************************************************************** c * arguments en j2000 * c **************************************************************** ai(1)=psi0 ai(2)=pipet ai(3)=al(1)-pipet-alm ai(4)=al(1)-al(2) ai(5)=al(1)-al(3) ai(6)=alm c **************************************************************** dnun=dnu/fnu pf(1)=fipsi pf(2)=fipi pf(3)=fid pf(4)=fif pf(5)=fil pf(6)=filp c **************************************************************** c * conversion en radians * c **************************************************************** deg=3600.d0/206264.8062470964d0 do 100 k=1,6 f(k)=pf(k)*deg 100 a0(k)=ai(k)*deg c **************************************************************** c * lecture et transformation des series position deimos. * c **************************************************************** do 1 i=1,3 jmax=nb(i) do 2 j=1,jmax read(x(j,i),1001)(iar(k),k=1,6),r1,(ic1(k),k=1,3) read(y(j,i),1002)r2 aw=0.d0 fw=0.d0 do 3 k=1,6 aw=aw+iar(k)*a0(k) 3 fw=fw+iar(k)*f(k) arg(j,i)=aw freq(j,i)=fw freq2(j,i)=iar(3)+iar(4)+iar(5) r1=r1+ic1(1)*dnun+ic1(2)*dgam+ic1(3)*de r2=r2 if(i.eq.1) then cs(j,i)=r2*1.d-3 cc(j,i)=r1*1.d-3 else cs(j,i)=r1*1.d-3 cc(j,i)=r2*1.d-3 endif 2 continue 1 continue 1001 format(6i3,f14.4,3i9) 1002 format(f10.4) c **************************************************************** c * lecture et transformation des series vitesse deimos. * c **************************************************************** do 11 i=1,3 jmax=nbv(i) do 12 j=1,jmax read(v(j,i),2001) (iar(k),k=1,6),r1,(ic1(k),k=1,3) read(w(j,i),2002) r2 aw=0.d0 fw=0.d0 do 13 k=1,6 aw=aw+iar(k)*a0(k) 13 fw=fw+iar(k)*f(k) argv(j,i)=aw freqv(j,i)=fw frev2(j,i)=iar(3)+iar(4)+iar(5) r1=r1+ic1(1)*dnun+ic1(2)*dgam+ic1(3)*de r2=r2 if(i.ne.1) then cvs(j,i)=r2*1.d-3 cvc(j,i)=r1*1.d-3 else cvs(j,i)=r1*1.d-3 cvc(j,i)=r2*1.d-3 endif 12 continue 11 continue 2001 format(6i3,f15.3,3i10) 2002 format(f11.3) return end