d+86FEB09 ACTAN FOR(DEEP FQR@]DRIVER FQRHFMOD2P FOR],SDP4 FQR_%SDP8 FQRzDSGP FQR SGP4 FQR ;SGP8 FQR)SPACETRKDQC'THETAG FORqۊ FUNCTION ACTAN(SINX,COSX) COMMON/C2/DE2RA,PI,PIO2,TWOPI,X3PIO2 ACTAN=0. IF (COSX.EQ.0. ) GO TO 5 IF (COSX.GT.0. ) GO TO 1 ACTAN=PI GO TO 7 1 IF (SINX.EQ.0. ) GO TO 8 IF (SINX.GT.0. ) GO TO 7 ACTAN=TWOPI GO TO 7 5 IF (SINX.EQ.0. ) GO TO 8 IF (SINX.GT.0. ) GO TO 6 ACTAN=X3PIO2 GO TO 8 6 ACTAN=PIO2 GO TO 8 7 TEMP=SINX/COSX ACTAN=ACTAN+ATAN(TEMP) 8 RETURN END vWDEEP.FOR< o !"#$%&'()*+,-./0123456789:;ome@Ń"<@XGu0s ԶySTo[mfnۨw֡9w}W}O߿wU|ҫ^z榪׻ǩ Զ??s|}jC*>3_ws;oַpVjTQ}VnM+۫~nSi*3W >ҨU ė@|c흃+n̻7u|c{hn8a6 |ol13?qzuzCY7_g l?PAU ?8mAӻY&f?gsMf\3g l?J7A 0@|[tu?_o=9N}[X*3AT _:g -:G3AX7 {(=ûì[:?^}|zK',z۷Y\ݯU׻qI۷y$xlx<5S,۷ynx G(T,w_Ywq(4g l>Jq|qw*c*f}|^a@ؾ-|C:gS۷ѕ$7kgkB2ome@u[륍Z8aM_iިޕB=? 8[:M+jSo@|=W񭶹T;@?)LlT:8a a9wHݏVq!T5$:T5g lZ8aM*&}rP P3upV HP)lì#/͖zT_3fK}T_3v:/_2]qQN8|"ޭKP״T u}p.ëNz{h k@إRw:jZPW}pVitwx\w9)_8avK{O~\8Plzy;un3v \yܟ/h k@KUS503vɩb%_@%ݺ%92KpS% 0 usg )L5=˽o:l~_8aMa7۹l 8#^@4 9 |5B3#]/@إR8R8#g :iAf ]*:-B ӂ:h ]*T)8aP3R9aIA@eSirtGg lF@6T gZ*nAfZ*ۖ.LK-™Juxq¶.™-LKPR)|q¶ϴ3- miAfZPWT gZ*J3MZ*3- meJw:'g l;&G g l;P=P{RjV@@J gg lT)TZIg l-T)TڑP)@\zBԦVRPm36im jS];RPm36OuP Wj/BԦ8a=U զ8a3U զ88aⓗ(ʅ3v)2Qz9ʅt%eGyGg]+GOf3vrtҵ" 8p.%Gg[-];R;R @^ZOԙs)Q.3v2^][ 8YRwS~̅t/iګ3@*GGYR~& z q⏏}g kH3֨*̀g ,VWk=Ͱk<ȍc ntt0+ Ɖ3vާ.89 |4d8a⛴+켿h O7O8a"eZ=y'3ҏ*e$3r_dT;I蝗tc$3ݻ1p/12 G)ˍͽ{= a|Əd:<DFU$NE:Ǯ4;acCf+"*@YAT0M8,O]Jduj^d:y].G2 oZ=DZS r"3y9^BFUs2enzN='a 3y'K2 Gߤ7tǽ]ϋ̰Cp5ydT;RLOG/: K̀C($Uq4;O`j8aI6ɏ:^+k$3;+ 3'٤u]W'r$3<\dtQg0M n6tYATp5TJZ᤼߻|˓Gs~'3M(dT;OC).qBfXûяn*&9Qg < Tu]x=''2w2TFF2 D\ S)i-crg/M̰Jxd:<#UqΓ`}7c'{Ld4z:pxFIFUT$$XǦ=7qa7> 8tOW3$*u2z AMrN:RR8ȃ3v'&B ݋ɻ}"AtzJGB=Qpg ^LwD*zXJ @}}*#AKڥl8aq<\O Z%g >I9DDӻU:aU ypA|;WtuJQpO@At!>yItҭJ:8 ė@r%N:R888ai@>TuҭOGB=tAҹz麧J8%J 3vƒMˉP=йJ>U U I{:aJ 36%S0p¶:.S;ߩl7 V%g l4 D<}RrR6DhY63{Hv]JO36!A4S;9an T,3/7 F{@ؾH)!A4ۍ"%EJ@hl7O~ fA(TG@~} ny g l?RBh(i fyHvw@h$fy('OR6 F$%͓l8aI nO<vQ4G84J nOy3R A4ۍi,4Op$[ToTZk58am!A^.իz_} UL?b!:g *x6qF ,GǠ m DtQ^ڨU37?:y@XN3vMk:3@*Թ3v.I5t @ع4M̀rW(w2g .e}ӻRNJ@ؽxܣ'B@+*z8Vw Kiw]P{)Myh/R@>G9 @~vsI:MiY6 Rr.8Mi/4 P;9a{Sԛ^{Y&/ͣyr-1I|?{J<n7֪.՛S[ U{|ӽk|3V1*vu8w/@*{[5 Wu[ŀ̟6gg luu)@ta mLj\pb4[ g ,Fs i g ,Fs@N4o8a?& qF &O*@o VT̾*"g&/Tgw76;Psxǃ36uu!9a6%-nL36L9`kdTh(rfuQS uAD36PA5Pqz n3Vms$P+4gh@kpfp1O 7 683YAT g fл6͸:q}s0+ ]ٸ9\zF7w6ofQ7f}xٸgffQ g lo0ۍ3VpF7׸9@ћ:B8Q|@T?&r 5x@o0[ ;5o8rٷ?kfWGn?fQ/ g l^N菉_ g4GH6p¶z11p ?5 AˡWH{M]{i]G#5n fA5kt⫏!AW-n]4ۍ3Vpj%3k了n]4ۍ3coYY۫|663?5w85ڙ5?fQ$g l7'\uo1ˣM- uS::s׽<g ķM]u8:s׽|g l}!AojO3z[ ė@IG%U O̵^x_7uJYb3vVi O4$PۡSGj g ljӐ MwYm36ZC\.ך36^EHKxB,զ\u ė@RzOÙkR_Hpf.z<8%v}RpYm˵^"@| į]?!Rm˲<g ׮G) u8:/')Bά6Z$@| į]?" u8׬6Z$@| į]?pf.Z4_~K!\tey@>B^SRm˵^g)@| į]/T)UQ5AT?&@YAT৆3Mk r\*RoЫHMr!AˡWzS.$g lfRfԑ3vV菉UU 榪| 98gtq'@1a4;+ 'M!Fsq Qo4o8as)36:g *r5!~P]}-U~f>ZB**@Xo79:To-nTo E4B:Tg lP#$Y ;}Olx:4o@3#9a--n H3Li:h1Hvp¢˞P8 erNòNN{IuXv¹g ,Lh$CSh/1upb kI.Pe 痤]v8 5윁„/IUh/˖Mj˞Sò'e?8awDò쫓g ^&B*p.NHeg 3u.B{ źN.ׅz$a9aq~)/I*.ϯPe WK.*@؛IhCSstjPל:a5g7p_ v8j@+ ;|՜Y dCk6H`f8a~!T;<DR6@'ggo{տ|u돩b gg?JAnkn/@X-nIms[w_9`Aۋ365$uj@YAT KYYn/p3z]!KԿ} l׻ڨU @PgBTo ?SA@ض[B4 ibKF?9aEBgwӡJFBԟOz{n:TT)HB'g l[-ߺw@؟&?M 5B\םRrPA@ؾtZo:T ꠑ36טr]_#g l7*l뽿@\א uw9a6A\zYTo3YT:f8a?& &G;OC).qBfXûяn*&9Qg < Tu]x=''2w2TFF2 D\ S)iv;DRIVER.FOR; o !"#$%&'()*+,-./0123456789:-4ol dRI%[^x,ܨWY`N0r.΢P gutVcV^|HhͭƶT`mE p#72Q* a9w[3E3!uy\f&l6fcz+nHIo>YÂ> aXFodT0r]sGϊLi⺨$XCJBxdHYGl9M0\ל V(:H|stdq$hݖߖ]IBф(:苩`,n!ё}D__R«5}"g![h4iܓ?QT~x ?6RkG?^ <_ǃ|U+IJBx!zxМN_="ʫH.) ۇOqQ:Ƕ㇧q^?\&iW`\p(W0tG ʮꕄ0-4t@ѹ[DoKgo3ͷbͷkWRk+0ް >#]wX($a|VtWjlJIq5vVuQ{}G}q:bӷ$kl;~/x:b+ βa 1U( apS0b*{(G%! {-FW1k`ԍщFn`\&nj-TР6t[}R`pd~{:Y=0s3b¯nK `80 t}zl6b٘ɡVݐ{Cɺ<$W 2qf5ͬ$aa 8;>Y0x#[3E3!uy\f&JR«lLmBڨ$nlb*S$qrAO (z%! nFRJBl6flvLm^Il6An VV5.\㢹5X5ZInkn VV`]׻uy\[U0]|sts4zuF+ a6s+s4ZIFT|و[uT`1m6fV`w[hJ^]Oc56W%[:v]` 0GTRI>q,nQ wv|VmHEg]׾mE֧CIqo0hkv>}\}ۥ,nΧcNSkp}+ a0]`NsW`3xΧ + apS׉onMm=$lу3TV ej[6F%! pyĻ(0i+ apS b_ <$Ma+MQ@¯^=7 F } `RI:g#$JPB$F < qa+JBxZuJ/aa ]婪zP5SuP*^W(عi {%! `J%` M^H.+ aXV* FP8wnZBpy^I RIfBīsӪ#JBJ@u M^H.+ [ n[#6F%! >ͮQIÏ>ͮQI5,ְh[[h~msY&%! Ӱh[[h~msY&%! ? â:bklme0|7N`Y=' kXЧaA? %! pyĻ(\wPN(\0JI |$Y0qSk9Oɗ8K7]b9]-=+ʹfiIӕfBҴ\EJBx-MW:UgYME/K"oie꼥iKNKU-MW{򖔄!9餪+=yK8LZt%{%rZ<髤%}-IIl,͒},%%!zoְ 4> KaA+ pLZttJ7^IB~9Ԋl#6:wc43͑/qԊFX78-&kiȆPA&'H/GVdC-zpO>JO$FjZi%d:~pٻEp"rFOK8J2st\͆#Z paӞ J+ ޫU/d,RU+/n.MFh2jE6_b+Is"ͽ Vdù hZ oxyv8Z !}⚸&o jE68,ҦCez#_v4diө^mVdÀ6=o_b{_Cyv(oU+zTKNjE688|_tpU+l>svqFVB6\qG v @opo?hG З;^v@VdyphWc43kE6$'H R`ِK]~MO=~L+=}jE6WIO8cҊl%.?'OJO$FS+9iE?̯,m浊U`^żN"I)/s/mz_$G7`;WtjE6t4,pՊl`K#kՊlvH(lVd3iE6+"^ܛYtvjE6tkc4KdVdCGdsMȆ>pmR=$ÑdiөMv~`'k:;s_!Z,l{2ܥ,L$ec:;s_U+fZ &daՊlp7r6VdCm/siә/n֊l2OR=$Hp؉Vd7y:Mv澞ey`gyoZ>8IrMRo|M6I؝Vd^8O`Cj|e'v~`'Z"{i ؙJG̓`u]"^/nY+!VnmK8VdCLsaIfIVdCLҦ!}{I%eeR=P@.H! 15yv(ƗN! 1P"t$|v~N`Z"j6Z 0Ԥi1?!Yӹ$|vI+߷?ؽCL7YncsתP/3CLM:Z @.2i,k:?' < CL7YnY+?' U:Z rpzMzӊl㵴!iE6Ckwïm}FYn1IGq}_th2Z LC YL"OV~N`7ۤIzp@#ey4>㱣x|vxt8ґ/n>jE{#db|c8HVZ ٰO1wdnGϝArѤ rsO"jT㸖T@6wЊltZ0 s ;hE6Tc\K8O@6wЊl`s/l lP=(S8Uܵ"brMބM"bfM"b:;^]+R=@#Kztz!kXMVW{#_vCbr"VdgnZ 5XjS80tЊlp`[XGO7qFWȆz(5xtzhE6b;^tU+R7xw"mzh OJ8(.h/ ͂uՊl`}q حZ p( 1I=t^+}CaIGǫVdC=pS="܇/ &ӍzrVB6<|m3^+ȆҾ\#9-q\K{ȆnK;O=od04Z  oYbWnU/&{Hߴ"bzᥣM+!GqKGVdCL?bM+ҷtg^'ng}iә!Ǚ!!Y+!_fH 1+M~H^Z'14HߴJsS]AiwVZ'AG!lj!Z `K+^\?/kE60TtkE6CqVdî8YtVZiN+a7K@}Nj[L 15Ȇu tX!\0Դ"ހݡVnMKNojE6t4,KGUr^dmՊl0 CL,&{Hß]&{nVdþ/&ٯZ &y_IL{hZ ͼ/$^:MZ 16VdCLm"b CLMbJ+!i8wL vGF+ȆVvܦ04;H%4=!&E{ȆVi~CLԟ|ML^+zqMai&VdC;ofOmzCLt4dQ:yٰ[ sS"_f1{Ȇu oewȧ!&Km/t4Ԧ_b[Je8/;Vd~VȆffoVd>odjE6t4y+R۪`R=VȆV:zZ 5r6{ܤɛVdCIG7Ȇ}-5V1SGqfiE6?Jbo'M:ٰULy~ӊl*~;w4717ȆRt4l~iE6|*&IGߴJȆ&sFOZiE6Բd_">nVdïdK+țiwț|UȆ\vyiwț|}hE6e77:y_ZiE6;~h**2:u_:ͅ_":ͅ_" ""!T_:=~_uvpSDP8.FOR8 o !"#$%&'()*+,-./0123456768 kwp~2у \Ì'%( 6$Sg;NhSԻșNJi`]tlSzާ>^'=P:>J 'L홾)}O9"/2ĉS<) S=zhNF+ 6Dz\ äSQl%htg׶ݶJ eӿ:,00A|TI'T^w /&0}0â$ D ^;pC2Gz5\q%A"%XI D 8XO^I%ghAdR/.`5\DQl`>k4_UjRlh6k뮚/%fCUxh6%V[a|MC^4ʋ`jέF}4 <5ܺaG[EI!j6Lឳ!4͆P^ šzέF}Т$;C=֭Q( 6\j"LLJ*`Cq9p!`@PPVlں<( 6J}qXDqchhQlH$ƚHhhQl`m%Q\\Eq3k?,J D\%>k6g&#f%k[|l`y>ں+5$`nzu{F#ޔ3Ena8;``ȹf0}gg%Luy}ǁ|Ɓ|QznL^4B8%~p0S7>+ 6ȭ^3EnH%8"qa}* 'f]C>k6ym?:fs_ȏk?Fn`}~ 3E8][wg&ܗu0S;ηzݡ`2ؿW|l`C(Ji`3wfqwW|l3E|tg&/ &EIE{]Ţ$ؠTM(EQllrl` \ܺܺfEIںU$ؠ5 5%D"5wH:%UQ|Vl9`C̵nl2ηz6l$oU8[sf%x [޺?6l$pi5}gd>SlzM9+ 6\n(+J S /8 w4dQ83۞ij9ɋ`Cju&#fL%ta%+ 6 /7%E$؟0nl2ηn;$~-֡lpwfnCMײ* @l3^ 鎉D״(n\zE~$ؠ#q2yػ]5-J e?X?J ġ³&$+D. g&Dq)OJi`E|#ȏOH|aQD1Dv@8d =ᬧ NfTRInEF%,KZIGu2ynQzo|ki##Vl ?+y=ZɏEqI?#7%QԢwEqI??GnJ 8( 6$s{~"?OIS_\:Ģ%ޭ,KZIA#1(.i%_{o|w[`#7%d0iNG#w8EIS{%DA[_ÅJ k_GQ{[j{EI!JmQl*)/J NyQlx9=&zQlHfu(w * ;m;e%"d}}@woHQ>N&3h%DAibT(0_pBI%׹-L|`Eɾ=j$ #QEq3( 6$$j"1`靯c"ABIڣ(`xlWl[msJ c݆;ؔ[3MITkd( 6 -ю$0n[Ǧ$ ^/ 4Ӧ$ؐ^DB__6%WظYľ( 6U k;H( 6$5 WvőPl.yCG%$ǚ;LbuJ :ثĮ$0OƯڣp`sVH싒``nwhkRl۽^vCC_`Cq~8I=[gpup$˿Su8PW0:nh^+ 6RcwD( 6J۰/.-J*`c5480PAObIBI%J*+qVOsܺY <^o>c,qCQ\$ܭJ`։ľ,nqC>k6M3(KJ~4¯:p6[0(?:f.y:TOa"QgY EqDέ[\) ' N$EIuQ}\GX GL$ eEIUQ\ɘ>^GX GL$ eEIa_$0E;4GEYv9w8f|8=x$qQl"q<'* ccVql7II%lhdž 9҉iJ* 6w?,J I_j"1EIjp~m(acEIYaM0$׆p0$5w86L0(aeQlxuG~X;( 6|!k?I_Xů]_w|ݯ _ۢ$xYAI]_J"0(\ȁqAKI%N^uKI]'|H\H~"1~J H\H:}Ebz) 6#q"1>뭤`ä19%WvSGP.FOR9 o !"#$%&'()*+,-./012345678N/+;RS }4`>]_#zTgSp- Yf}pp-Cw- 9np-/\X׫t7tQ׵d[xzRg33L5L?W[0t ݿ˿2k]2j%|N`!bKt|}k;`lTĢM43_5`%Jimy*AÖ0L0L۹-/%|IXU[8eôJd ݨ'ЍJP&}Ϛ%tvvfb mϡ-a-O%|Tr|,6#7?ln6Q :C`Q : rhTf.!kP -6`ɰ%| ad G[F%|~@  J)?aJ6Taèc 1d G[F%|aӱ-6?1(A3-O~,6`ry.6ϡ-d`T·~,67K{QV$x_VX ѫ zW?pqcspУd8h%|MP:Z!Yt^'[3oQpV%rhG%|aȚR2DpJy¬.0s / yۇe( 9^ Cx pο`#zJ0GpУ%@] D :_ S @Jtޘ-`W{Jҥ yۇnkaF%J,:Oн5NzJ %|>3mJwiTΛGz)A5+A磞-Jq{!4Q :4 mTΛP!9.]\o/% &^##7QϏdz~HC0ul1uJgޘ:4K-L%ڡgS : QQ :o>CҖ)y;Jy -iht%?Fj ̆zr4hR/HZ:4K1u khJyMudXC;ꍩHCSnK0uRoL7N :wڡ-Q :w[;wJ5U;HJ5CLeҖ8\&Z :3T&m{ލzޓ&r5m9tt ȰvSw٨}T>sCЎZK1u6yԧb ԏJXts?"5WR'[_P?Z2+A֧bpop_^rslz%|?%ꏼw%[/ YO0dMKܻM :_ S1Ks6^j4 0u%<D~RiTΓa)0d=J04Q :O͔P iT 34ae. @fRΗd}Ie. %1L6t2, fl'%|Ȱ,'%2@ %Ro6A&KI :JdL6%߁L4X9[wL)QΗN/J[Z'ߔ_il/bZ3D%-CLkht+CLkh/%|\Jei}g}2D%Wwh_ Qt~ձ_Y :9+OCbr"VdgnZ 5XjS80tЊlp`[XGO7qFWȆz(5xtzhE6b;^tUvBSGP4.FOR8 o !"#$%&'()*+,-./01234567.Gz1htA+)f@H9"sht?P>טy~k^cטe4_ݿ< Vd~{`zAv{ e4{ΟV>m^m{}]٫^ELuxWA_׸k|5`40}L+"gAݿ6]'`k#1/ta?"=3E8*v gCZtᴲǘӑ!2ht:|̋1ht\y٩6Vl40K=w..k`х7婖~VWх7媮+o b40ҶI3|̶-,o b4`"sE~b>\^g2,b40kwi$n|fbY@ dyZFCc\>C?sm 1G9m '^1xn64}`U/`T'<B6F @7@-E H17K0 1b]cFu] -<-hRlch)rƘуh3{{tlg$X F 0,Eտ8,z$X=pX9%F QEB<ͳ!\xXX[X*TF po8Xq/=A!tA^٫Ю@`zVe4pegW.$XlyH2~^.7`5]9w_HhtA=k5>?B< W7Wjzy.W">htA=z~'Xw F#b?!!$o b4@*?!]s1]XHE  5] [.$յڷKYu]Fé^٫)ht1Jm ORQ]En4ЦNECd.}/.@V]w2?duY]Fmx/2.@V]eg߁.b4J}+htA1vXF^?mĂhtkc8 (?$c!Q1G.m_d'хS11X,O2Y1ݯy{ Ue U2ui$ht&t͖$IEu]dg"s闩 Ȫ_sK ew$Ȫ_s;UejeYFڽN]é!$s du dȪb4@7 t NF >zҶIR?*c!.4Um.tOwiSR?*)u1FF?01~,x/{=݅I@ht r?مn@v.thta,NZޅn@v2]`߅"cwt&R]dht w}.;Ogr݅n@vڟB7 \FT"Udх1TT߅Igt}BYsWA]de*GS}+4V3{xX^FFT/t[i 'y+X5]3T.o ˓H+F 45>Ҋxh|Lwg&Q1G{*Kȋn|}=69u.}Xͣe9.=Wn40\byHz.B.>! (4.K\:7>aڧ.=j2 eեu]7>>|2<2uiE00ޥuht .ôb4" `h]8WTЎc?-yVkl8aш G?y$bn`4aSe;$bn`40fS3y`R170]`[?p2]}NF T T F101NBj4@_ZR170] []=T F|:Tht!樢Oڱ,`$RF8T'!sхIJc5~8-2b4@SpF Oui'ht=-F#ZdE&hta<ѧ=<>f`3I|:.lUו/ .|Hwt &X#)F+Hc5]F|f+o +o .C%IC1.|>n|2]T>F ;Wvw>m ? z9m4l:֓aJx˘yht6F c9aUFT̐*i1]h Yi1] \E O~gi i4@ӽ1KcO~WE OHֹt |:_|\Nd4O}ҟ:ohrB OtAx%LF0p>f~h3E6<+>w~.*O OO8kc;~ghtZj4 OF k2~ J-.S7u]ib4:n ҜTF|4.t7=,F ձ4XAbǼ,b!.If3?/;_T >хclUWcu\FTӱN]TGiF tS(ROO:.Tj&!htT'tK'^Qx%LFbSF.[]e.R]| d2˲Td,n4VhtOw!i"*h40f`cB$UdQ1OGх11vK O4VEYh4L[]˘"ctD<F40cct{wte1hAxR?ЈWe肊U|хj>n肊9nU<1m1]hVym肊9VUym1]h2olzQ1gޤkcy*nF DZ~q̛4Vht]G-*㶫3oX7хU7Q17х~]̛w4Vbμ.]Uqm Um"}m40uh4]MF uF Tׁi㸯iht{`{_O&?F Lu69i6]Xk6]<_a­a:J2 a»a:J2 :KaaKa͙J\ßa"J4fx:<: <:b: fOvSGP8.FOR: o !"#$%&'()*+,-./0123456789@*e!™E?IEz+Y[q]nMLFM*/[[,䇬bCȏmeeF/?& 28ďw?,0TlV_ZxiG}_[V}FzRl:}lbj~Xa4r4 n,حIEhi˗֯?bѤpq]BI|'zLHBI0h w/܇eJULU 3[ hREqMVhd4x xhR~{xhRr* ܐI(i w ߲:e3X30zSZ_w'IJH+lY :߭F\8&Zx>d4eFzѧ.zO tS_sOZAH=uI=qp t~m7e2T2c.[f|ٲhTnĻ]hpK䪽^?7ԙnW $.RQ|ZdvdwAfQ0됅4dfTfW˻š/iTdnѤpō[f̰<(]1OKI(i}M*Wha<=R> hRߊ@ ziOFSB+]Ѥb)&W&K,G{h4g "sBoÚ@ԅqM= Z?uI=qKn4g M:!/Ѥ"GI+ΰel4hZ_u[o4]  0E %NF(Xb1hRqv+@ avCFn~)>MhR=ۣst{e-gg4c eˌ+`lIE0~@H=u!aÚ@Ѥ"GY0Shԅ1@Ѥ"PZ|'@'IE0GOFNeM=WhѤ"sALj f8;IEk F.$rh&900r8wd4ȞdhROE #pF|f,hRW+p8M*3cx90:d4Vr9+# |f&9|I=';HhR@y 0r`gMH=ۙt~me7o cФ[ /Ѥ▙l&%Ѥs<~J ~z$m}ִv~@ cÚ@FzE.ZO!pN iMH=/F_dH- i"FzRϋ#[a\S?Z?MФK73z^N0y Ni&pC^0`4_l-O~ /d4~-%xIEoeee &1+>:Ft0cDwLJ#J VzvѤSH,# Y4M* dF*`}# ::ft M*r)xIEcv/?ZȌ0cvFOE #d4T#Ν:u0r9|tx8u`n=_&ԟ/V.$fE~2T2eˌ0(3pxFzڅ,qzn^S73p9p&"Fz[-3co[_hR@@f %v'R?=$}OhR@$KHH`4_D'ٯ0Կ~~n4g %ȁ} 80{FBgͧ~&FzR_<}'F|gHH i& nMH ƑWFp 81MFzs 2v^aȌԽO+4MnKFzR[~N $rhhRϿH9t3!e2/Rq M8'hgFpCtl csI="g;C/hRkfƲefF(\ffǖ!$~`˸3M*sQԹcu6T\l[Ɲe238χOd4e.[f|[&KH`2 NF`/7,V;M*_<^ja4M*.Ph"Z4%M* Xx̘a r24ܮv)h _g(&wr^Nn N63YwxFǃjkU~*vD@7&&IEYjU,VvYF@7&&,YeB-K܄Ѥp7ԙN&Z$X&QԹt,V ܄Ѥ":l4es2ҁa2TҁߓѤg([f .- #::#i2T\$1(r ?@ABCDEFGHIJKLMNOPQRQSSTUVWXYZ[\]^_`abcdefghiOKGjklmnopqrstuvwxyz{|}~V_;MN/+1:J>-uF8)LH*lU'&IB,<9]=Y7432i@A_+{+ ;{p~nDNzՃ9X^ڝ\v~$ٻw{yDNzGi9W;wp0G98}I<{;WENDsrJrP OPȉȹDGNC3,@%3l&#Z`gI">CnC'W`jH1| u=ihTЗ09CuE-[`j1t\ihzאBH}¶/# JrbXl >ƶ/# -(4g&cEd6BB cX[*%i` Xlc1t|22P2UƐ+] }Wv&;/ *C%\UWG`_dzL"YCB>rX'@En. 3l=4J򀘂S W%d1 ,`gI*PK#'e) d`ѕȉ';O"B:i2r9TP5(qUC0䉊^.]0ɳ xU W) 3 |D RԸcHzJi`Bpe."| UCrZ;䲷\MO^`xȍ:g=20-ʳP`_d0+gȍ:GK mGΠ-B_YP@ < ĂϠ,Ekȕ'XCG)PSCoҵPw4r"ķ(oʰ+r%|VL- Z[@S*P5N*Ƕpi;OJH/gل jRhjazc(QB_FNFJ|cidtQriXa ؖ6TAMY[J0trCpY(m},r2BU+K#$ 3$&tZx8=C'7T !jK[x2ؖ`NMP)vv=.mbB_ x3\C%c2HTExHT}n}DUDbL-e,FZ50!#\deǶa DBio;OJH/gEQ ۂ:(0t*\#Ԃ`(WfQC3b)j$%O!7T !]74}FN8 ҵPw X<1;OBf$=is2B}Yc4ͻHǡyt-djaːdZ}YC-ArM.C&/YSN2T -r2BU``gIP@ < XȐU"Y2NEdjal >r' A-Q=:rZ$`GKTÕ-ط:d 9Z'l=Bn. U#QzFND +Xv/N.^#9(NcDOd}$ X-BP$Sȕ'XR0t)PSHc!4 `[P`U>nrMKT=ej3HNdWHTE:=xDj- wz_ȉ'BK $`j A֠'T uG-ѣ3/$S&oѣ3/=JJg&cQ1*0;CN*[F ,\ΠET [26\2,E}H瑑ؓK(- ۾p :D. 7 "ɕa]Zb28URU9ZPf*K#'rPL|}$}БyN SУjљKbx)L|љKׂcr94(Dzuۂ'CBZ5 ȉt7Wnv']Gnd~=r"ķ8?}PBz84nh iNEdja[HOK/@-1_`T Uc'i -<-DMr9ۓyti 8Gz>dG6ű B udDP5D0j Ė)na\CQX*#/<a4Tm< חz2 PYc~y-gHsCUI.YLP%\Jg5ط PB_F"&\} O2` A\ΆΞ4o}EQX*#/<УD$ $S-*#/crP`b/ Q!ط$S*Gg^.]T}ܣ:UF^y G-Cg~yx" ɔ}љK#7 4! AqrU-1[WwzdMo(0[P c CGLZ+n}ˠD5e)406)o1`2Ok_>mG:[٬22SЗ!2TmX: bj"g$B[l>-'@ȠI.g{ѣ3/C[ICg'i%l G-Hl^#pը0t+HGKFE,}b3jܨ\nIБk1jABnۂp~$rCpY(9(N&>na^G* 2؂} O&[F ,ɉr C'QdH|lBp[y [x*YX(Y(Adx,l1t*1!nU J x}H|г  XL8Eeecdt#٧%ؗ5طx7TF^I6{jlUI.YLB 1Hde-;OJH/gb,+7 $ }dط_Zk22LFN8]}ܐGz'*2!2LIT`(+e/@|)V!u >d-vv'@E* ࢆGWz[h' 'A%3l0hccK[*%i`B cjK&na׫;=261VK/C}o* Zr%|'1trCpY(9pjzo~eb5889d[ű>%dD\- D9,T|ylA=”ɍ22IQdge\|c[3(mSB `)HGK^-x CgjSB2^ [z`<[`ߢ"ztҵBB uejxԨۛ^&8. U F, V^`<|,zܐGX!Ro^ F,Ŷ/#2Ô+•ApS'4g0ġ[Ŷ`_`Heq |0-ơvYC F,y|2RoR`C3o* 0=_ uRo o1B©Q! m17 0Lط\G`X|ڵ;[:_ez Mlce("$QzR`C3oaĖ\ҞP5!%C{ȉ'BV`܇N YhQ S׫;zu|5e)jo1t*Gg^.] f*Y(0UE-طHA!B4rOQL|DQ[F>ś{Ð2o!t.p~$B[P@ Ҵ 4rceZ⅙+M lyB+nA}#ܨr>',d/'W0=Gf띋7އ10tطҜEgP0RgD)ÔX2 -sr#$[xԨ `[PSBFdȥܲ0ġr -S\:FgPؓ`d%R2y'B Wc2L-LOzͽg:ω@%d2H݂GUZ"7B 9Aq#x•pXj 4KQ} Wi JU3W:S חzZ`ߨp[Tcɔ 띋7v6oRB܇RsУN!jǶ/#HckX2%95{{4r"E -*ԶЗB [cdt[ȍ2)c/C-qL5jC=>%dDDUl{dNA2 2ÔF"+ >C"f)j\F \Mrh,D5طBB RX=r=r"ķ8?} PBz8+kL-BZ%rCpYۂ0lob+N#'rP;ݗq\!,ЗN Yh BAXQ C z0ɳBB ,E  !)f\2RoQ>Cg~埡r2rQ[F&'ޢ XxOUP cH,Yй-{ӎȉA> >X=r=r"ķ8?} PBz8+kL-BZ%rCpYۂ0_:-ȉ'dge\|W%yȴ*@x,eCoSBb^݁ 0tޚW'ς % 3J5d744+rԂKE [ HJ}Fml< חz2Hc>Ua@-"JHdjC4ӎȉA{> Vq`ȉ'`gI$$kP@ < 1=ja Uejl >Ҫ5's"~;S'dge\|W%yȴ*@x,eCoSBb~iu"nVyk~iN켝FND:ɽW0_z9(N&HH֠ x\YcjzX-B|zޚ9WlFN8yWv&;/ *C%\UYc/C} />D |jWvN#'r"'.7{`WW^v>\]Aq2-x ѣ3/E* M׫;zT뒗} /뽽1t jR`ߢ"ztҵph!طH܇2SrPB `L| f*YX(Y(ABı-l`Z-00dj'Cّ"G. Uc[pvIlj\i2ry9 UC}q CgJMT}DQ[F>[`ߢ"ztҵ©Q! M\1_! }CG} z,, 2c)9 ,D{;N8w<nRaFN8caX$[TDμ\^85*d-DyԂPF`b-X@|"Y2,T uz5e^݇- ط dF[G`_d:<؄K*ԜAdu WnHDz5e^݇4Y$$kP@ iȉ'"ug`ߢ"ztҵ©Q! m\74[LiKh-eC -熁djʰP5dAy \Bǡ1_׫X\C4SW `n}4` ! Bʠ2SrPB3 i,ry}љK F,Ŷ/#2o!ܨC"ztҵB-OrPp3طPf*XJHFN8YFZDP5طyt-pjTB[22o1B:^85*d-&&Ϡ,E -J:%C{=-\,0t(W !SOh" `7G?+*Ԝax,3o1t*Gg^.] /XxJ5x,ҜAX=#٧%ؗ5$e \zQUIt9=\=ؙ쾌/$pUe u-pjTB[ 0tZҫ;oo< z,,*RԐ>Bj%S /e3o1tW*'#Y(-edL|,\_- T0VK$(!͒!j Ҹ&N;N#'rOx޹xs}?_\G|48ߛ_Aq2-;oWi6*y-pq]-oQ@. UCbo `[ J{ אeSpia\1}DU>ЙZ u [p e̕ 3W(q9d[w.{dA\p-(4g3Zpu#+!u[:!L|,pUZ>J,TDU D:{ BʠE?+!u d•pXj G HqUZ D W-\i6aeqL4r"ķȍ2)Ҵ t%\ɔ}DOd}$}G-HCciգ`XC*P5dm !3o!k(!:bP!!oCgg" X!SOhׂriD[d[@c )طX\>< 0epc!\N!] [ EE94 Co ` ,\Y2>VBCǕ@Uf2T9(N&zͽ9 9طH$HG*CW%I+Z:n} W%?[p>&|P ΠTBL&-x*]e%Nnfr9d[CP5pjTB[^ !G-oq,2CxB[H݂<-Dþ\@![ad%ΰ-5ط0RoYC/s o1B©Q! mqy"C.=Z5 ,D #+!u[0ġط:Ro:ö`_`Heq |0-ơvY^85*d-ds C%dv-5>\لPþ jo$Q5MH݂![sI{BPdp \G:E #+!4r"Ճ˸Ji WUX [Дas&vLv7qYmbHdsHc!zͽ&vLvӫ;z`ݬҫ;ۛ^J켝FNDΗP:8H.Q-Ȱ%CAXJ>T`"7jȧ =/L(aߢ}  ` Gĥ.{rp>OOb_Aq2-x  Ҵ 2oeSp6o1V& 2X=r"$Eqr}MkDB9Ϡp.A^SBFdȍ2)T CUV_Z "0_;/ڑwWif3(n 2de. Uci`PͬyA\P Uejx"ifgPez D$B9EB5C Ю=-"W`qlд%Y0-G-H´ t.JSep: ikxzJDxP5O\iZ:PBzx peEQ ۂ)\+B՘З1tؗ ogp2$B`UƐ2 ZИ i$r'Z`[0siA>Ҵ@%3 X@Ndq"~;ȉ\^= &{pew79ZИ i<+= aKSZ{u-"Q^0t&y@Q 2RWҴ̕ u [paKS @w ![ ̱ 2A9L[Uȉ'HL-װ%8r"rşMiި3hF(m^YvON;ݏ72-Pl^Z쎞08 nDq3n&:~=7Mō[(n6^?vY/t7i'vĝ2-PlznZv~}ǣă'qgeZ7nM YwҸ;ji7gݦ{e'ŃȉA/~X"x}0Abv8.[wgeiǽA9Y?-;iAbQǟ;+_ڣ[(iӴrEt팚AG7ʴۃvsY?v7#ŋMwFݸFFnӽv{q8 e3ō8mϠ;Ѡ4ȉy'NlA] -˲-;wܙ} z;f;rFy7i{}qYw^? nK3jz9⍿^@>xE`9&Mo%2x4衸;el9(^l^mw: 2 M(iӴGt0(AqLN;.۽~ڍ;Q'E43jQQ3;~5fm"ŋqGetP܈&^3Qܝw&r"'e3viMLNm -A7.P'8ő4];qMLqA?Lq#ФM{ȉ&+xw⃏#'rP-7Pi_zQdдΨQiAfWnō[fn9Jf7@euv܌Y3uiGv9(n2-ăM?{= tq (~n?A:,2mq#'r"["a}$r"ŋ4.hwQ襃nwge{q"8qu(4xQܝiv9(ޘi)G`0BqI{N{P(ΚlQD7fqZƽ^Pܝ5i)^eoڝD7fqZvͨP(Nۃ85M4ȉy'NlA] -˲-;wܙ} z;f;rFy7i{}qAq@g뤃~='xv:N4(nM' q=hG^Q7nu~mw:M7.t(N` zv~\v"~Seqx4t4vN@3hle ڣ^wAߍq.Ӳq eyWv&;/ ДS"JH3j\=6^Y$0EDp Q.yЗD-[ķ/V#-\S8P\ 2-FJ˸Oli1 Vh(:tH ,rCķzL"Ycj[`9 ,DܨcQQ |0tUI.T u- ű'>/'jA}JdH//2r"q=j8 ,peᑫt=9(~*Gqp[{s O,WfCA(9&GVwP܈?,Ϡx5N?Cj"λ(^Ϡx"K(ngPQڤMy^Gqa< ūWӑ9ķ`}0 X[0y ң Ɯ\jl >5(mQS P%\JoPDXSx98Ŷ#(3,HAjcЩK“LA Ml >Bn. UCh&s,4r"'r^௏/;̯GN82d_ * Zs|z)7eP2 BP5 aES}*IH >#BWp2-^/Ci u-ơvYC F,4y>ZxDO%|Q2d/} W%yȴɍ226B xmm ڲ`{N\ ZCܖF\й$Tm uU-BMXi@`}u~ } re8Is}z}gȉ׌ lbnٹ B:+=u]b,\!1LA b2pDq6rM\(_h7CtgNꘒ)-wF4r"XM\I/m@ z_;&-٧Dp>n c\k,LM|gʴ p% U\mYݏ78kڛ8WǔLi륑9boWӝtϾ obo< \2`Ob?xv[V,T uk(bO>\2`zvw=6JDr_h&Io\S2%ΨFN<"5} H#XjߤM<["ɍ:v p>F^'eN1VKo1t(Nv p\2`QӀg7{V2%z:hV-wz4r"CN拏*\sXMHM|~nك[Pn#͒+:O'kҞ5)[ɥdj=JJgbivY(a[Ʌf2r{/JfM\_ sZLtÄ\_ Wcre8W]]dlnRmjnK}ph&>)#(Ú8F=-I齋GXyy&zt-&FarejͯbM.%SZ";:O'dKbMcҞ$RWXTݫmA99 srY(""2طrrX|IS*bmGBB d,Z-  j/2T  T} /'W}4e2 = C0t(^L,r2,g+<3H݂brr SI$30t(^wfmE`E0!͒!u Ji WUQǢD EDz FSx 2B|z)E* M\SC 'ٷ[<.wء\E"O1PSBS SlB_aEN&ZrzL"Ys2¥yxɨi2 r;JVD:yn M{qDjcP0U. UC)c/ =an?A:, FUNCTION THETAG(EP) COMMON /E1/XMO,XNODEO,OMEGAO,EO,XINCL,XNO,XNDT2O,XNDD6O,BSTAR, 1 X,Y,Z,XDOT,YDOT,ZDOT,EPOCH,DS50 DOUBLE PRECISION EPOCH,D,THETA,TWOPI,YR,TEMP,EP,DS50 TWOPI=6.28318530717959D0 YR=(EP+2.D-7)*1.D-3 JY=YR YR=JY D=EP-YR*1.D3 IF(JY.LT.10) JY=JY+80 N=(JY-69)/4 IF(JY.LT.70) N=(JY-72)/4 DS50=7305.D0 + 365.D0*(JY-70) +N + D THETA=1.72944494D0 + 6.3003880987D0*DS50 TEMP=THETA/TWOPI I=TEMP TEMP=I THETAG=THETA-TEMP*TWOPI IF(THETAG.LT.0.D0) THETAG=THETAG+TWOPI RETURN END