2CPMIO MYC/DSKDRV MYC 0/DTBF MYC<F80 CYM?JFCHAIN MYCFORLIB RYLFORTRAN HYPINIT MYCiQIOINIT MYCL80 CYM7Y%LIB80 CYMLPTDRV MYC2LUNTB MYCM80 CYMj̭PL1 FOR%jsPL2 FOR PLANET FOR5TEN MYC %TTYDRV MYCjvCPMIO.MAC t:W7_v;blnYU߭nBQKܾS'J=y:tROOmdj,' 32HlOg..t';vN\4|x77#fđ*mOQ+@埒r*B8\D:aLSvEP_Qm}x,QA[ڕqM9 jы+A7k-T=7xkj@T&Na:rPU|?tw4c5 1j`D[b\*\EOg~z',2bm98r{$g]j)=f,ĬnPlfgCO!dqqN]֯z4$`ī/`d$)"r{FiN~|䖚ITS+.(\ŘzCx* O;'ͷTNb̕`m]˦N]Qzx64v9Tt r2r7 )#DB6W "Wd#aհU սYGl&_o{uRnO32/ȹ܏v;|n/%_FlhV8G[ p8N#N"\mf֏!p_{LH"ThFie?w>ixP,~}Ÿ>VY.^hpźB%ܗ+sp=9Hӏ?22xWiD\;Aw>Z&b ]|:LMwsrҲJ9rG0ډu~G2=tˎQ"ѦȤcDTA`!Ϟ ==X'W̿\o8^,>î'Kf6xE;B WjtnYcmC=w5Tj2p5U[4x{d7_Sӫᅁ% ,/@SA͏S,Ɔ=!F#OԞ"$R>2xN'^ӟ ?qx{1´dOTy 7R Ůp_Rwuzf`k'pPƠ`SocHMZ A՚0Ncnһ Tae2fͽ{lVBۄNWj^lyI;^Cc*I(-ÿݮ7NnR 77\_{4 ^ei3RF|{lL.7݋nyZ>1*\s++Ǒk"IK/&dFo #YRq]~M2N<2H,?iL$ppǰr:v~ͥyX{9ףk2 rMnU3E<ʴCaqy8aH~珟|Z]llϋ68Y0ZT%N@{޾yeV+_ϱrq+2AG{8o53ό4X[uT$Pݫ8Η%@Ok[͒>e~4|wVwrx\"ۼ7]2elf3ki$騠ި?!O9wޕxw=!ub>;V̝6wӟ^5iO_Q:@,di5a{{B X`*5I^A>P[b4P'A3I;ȮE9 мT,>h:~{=◚kW0|O=nnA~Y |(qR/F|Z[{QU亟7WclJ \W(Wh x\⽸@.xnV8fSrj\G;>}zITr|ke9T _&0&1X*ݱ*%fK_dT ģ'7U<|CɯRB2:A~x",~Pԝ{boM9ćcC69 Fx3Fސ3;78n# %]S{)<>ݽ_9Jsn"+: zv~{!{}ΒQ:Mi D7(XѩgQݬ"?*%jB-pdjAWVkJUDF~ "ZI;u B(+ЬZ zfklGGxmnnۥ\;Klqq~Oڑu1ې9Alh/s33xb/-KS36!*IԖA-~B# zwP@\0a3Y -T'Yt3.&mXGy~3ߗ>_sU"m/zB/R`Tӫ#w+瓗 }^ѷܦ9`U={7 ۚTR7cA$КVw͡zb-H. FKWʖq%5 ~Yڍ^S+_e.g5Lݑ?/BO_OoЪeΓx=r,yiSjL ښ1E;̍ї Y?p'MO%^)3&Ljp&IJ}̄?vY-v;InV=Aŭr ý|Sπu\ش+fw ĬGWnѻ9`ZGl757wЙm~~Q(6dk|ٻ o0VY'.ۻS <.C?Fm1z ^Ee❞Be]e=ws[W[Zn輶m[iK5U.>{-~k(~.kPNcMЧrkj>p>C, %~~J^9zQ!6xE:jyޚo~WC/\8;$rUcĕ1HdTK xI9>yI+'d^s]Ju bvw+g|3[5n, ?ךK^\l+pGas?b矁 YS]^V~IOr$KOsbd ﴊt3ɚ`$V\sffngwBa?wExo:k Y/Hg]<٩f^+d+ѩuh9BiFsZ^}dܭ$9Hlc$^b{ځ~^>ܸ|3^EJͧ)0}^X~/'Z+a x/u9~#!&]RrI#yJd7T{~ri.:U;~W` s/vDTBF.MAC t::3??U.~̾S[Y~{٪MOˑU ûY{FCqځ7s,LƬVng2XN"N e7񷮗Yˁwpo lufM8O:P<=,5dFaa2 jXoQ mp5ίSX,+OS , 'פww_l-w9$ .B- Ti¼Is%cDyI&Y/11W DvF80.COM M(b@0/_ *@2[[7cgǰRYR}lgWZu%QAϗf1̵ًoir'%23RY/wm>%n!}x[)Rb9&P%In'u Kw"eEpJdPDjgXI~ߏ*^y/_] ʳͳY,) &t<s_D"D = ;ܰv} =~\ 4˼iN>ivKWxkzQbuoc49:anbw53jzc=Y܈VBrl\ҢF/jY1(9y?_&2HP=n;Nbx4,#?s~n~̏?2(-[{i;ikr*a"eF QjLavc#W={C߶m%cՍպY/8 G?Vރp_/ݵ}/.>A%0Êx6DUk moI9sBmGFg&3}BE9^T}(qAV)Ql]{,Hy`2v8֘b,#vX1͗jC䵘ۚ~] p_X%TZ! ; 6L@aY^H%/b+npɥ(qO_u*p'7Bņ "&,U< ƺN~4A :McI ZaԽW\])5l k_裩7`Cp,ygCӒ?׊p; _X|qof4>씆g9t-O>hC,l(.7&n7#sOjo9yWU΋]]1/w=>nϣNP~_,++ܦv\ h<__ 2m9:o_ @~/.ɼvrʁ/Vs!.s/CJS8 GwP5%̫chN].HWڦŵ04CFM/TwΙYGM /u2 !|`}ͯQI`bqkޓ<ѥ*z*uHE~ωz[@OΗ;Ɂؽ5ul_ 5k|9 g7> 9Mx>&D.yWڐS:xra1ǠhWA֔dG{:CBT*1}79:,V&f2>ι{YeE]'GIC4r\j%޴ z]]Ue$bw03g*Z7,%L`.¹3"ODa_9rW B5CaR[e?%-&tS}$/QϘ_HQ9ΟFvf9@V|{ʣŕ3ԒAя̰{׏Y h8-B,|o̒] 5Zm|™d{X 4F[T=}F5"Pu eq3XCJ|/H*6:k6!`6 {GaMݾ6i{/UOQD|,4=qDEQo/#''ޗZ &9+LXY?1VHخ_:?}=rՎ$J>߷~=j?.jB/FVpd? (+jweAi.>L_6񻺑"/s֨=yJ̳~ccr-J^|5%!fQ߰sr[1 3<$Unɍ]o>/:amFc?O3L4o KLl )LtwZ$;(#K=Zc$C*\ 0băzVu6uRgmQ.xl1a㖥w4Lت"3y LEU=,_?]ʦME[p#CvHSaBg?n'6EL'X-UÖ֛%#s|oOř)sE6|]z%k;V?k&Sce$OX)y)Uo;ǥX,L:xf㊜jFrDӆ yNFpsY)ē0q©^ڡԲ|Ɋd7gKe`s Z-01c@_ITdH$p#_åaS?4n>O NsJ`A!4w-dOfgFqဎߕEa_ޔ>Fshom ^6)'/fFC}ZCJWe>br,  8|l$߿'"7s̈ q @w[2؄biS?2'rX$PֻT_vMUy)͝hkBR WpOG!҆>y&ۢ*=lK/u>5?NAn9[uT;Z@mCd!WVgVصw8cx&Xp{¨>$ٺc sJeFZ yĿri U+3XljuOsf|g'm+"*!l;K("/|Dn2_1 $7.w^ n-یs*ߍ'(ez"E*]珼wS_^=yxszfb[3?5PAGsm vpϽ*D|0ͺӽ|8'XeٿFft~&8{_sݾ=;@}QR a#0ղnu(kzPὀ!՟f;b#:›S*Z"bW45ƧmoA3[?E5D2e socېq9n #DU9G#a"[ T 8.=Ct]hL|,E.HܸdYT~@\+dQXߕn/e²F(BbYxk3_{; t~0/cp%.XU) n&(:"I$` ɫEoHLF͑3z19<ު"#*\3 0k~_K9EygQU?OJ+' \bR_uS=*&`2,kRlrTw::^nt߄1ӈOij.^]̞ca.o3]}3r˨{-_3?ոO)3v!56ALvmpuYѾQ7ſ;t-wa[|lv_K!δgOߎ~y҈p($Yǻ){8Nd׻}t0e"lv')vOgk[$Wݟ% `D{F\Ŷ;qb3OŲcYI:-tOTdr}_q{G4Iu9xbQYJR6UO3r'ڥ-#p{/+~BKۍnlXqlVaHu &~hqrg_/jM[HV[/ .N3{dm8I Sm'dO!?,dٚxq)z6Xeq&ퟦT>Lvy?XȔbUSłF%r֭pW˪4;gurbwJ`L2w4cCͭ 1.wJ_׶C_G1tw`}?Xkg)EZ/ҬUf"{] qNOSnd{Cm$^->)#u=0~}6 ?*o=@Gg %:{4EUca6AhHڞSu> vru8bBaXہy["ͭ<_ن^k\$iHZ_؊AmOCg|hOxN{Br2>\43_ ̘GaM_o$)l P&ϱ\l30I"(SE` u1b"5|"9`[*?5g sxv{ML`_1 T=e7ڀ}XLUSnF2+8i浐}P-mTCoCZ,)Sj0]5tVlPD=R8d k$tD> g `I*A/ku$[ʹ,^S2! jt hGY;Rm~[VACfqB H~F]Vnt> ynYP$PHkN*NcZcX2 `k֛Y v@C4^)Pdb\9<)?GA`0Wk;#hdMt5-; 5CHfpw$cE?G }nB^$Z#I,-oM$shL- x*;j}t+{OWEO4>i_9o(AoMG jmI9N.~ " Ȼl$i9=aOƮ"̍_N..⍾O) &>W!Za]b[׫ʬIoyc ɤןFsW@]2>{_M \Ju!)q[)eDɯ<![Xsd~f 's4`E.?$ \/8RV:.h X˦] ioJ U(Uܝ(uQaq^R @ҿ~bmhBpTέ fR+'~7=3jQ_|| )A/31VOY,p/qjpgҺ]܎ 0@Fbʰ#Yo;~+ł*~}3^^&vY|ӕh}mxRUTxO6qd+Bv Y2X6l85ؓϦ$L`a!K-T6Ev芷adN@RBWYC1^X= Ex}^g$6lq=,+G{d;)>!p9?q!po}M n\s9옫Ê7}0o(_vEv m` 3qVʆPs9X,fezAR Ut*zvD^˸-SrQP{FdlaC}b~#) Idաۡ'YTouɄ Kk͌5xJ~KM Zi06Vm%^Cl ;P7@%EwLuCH֞hiI.Q\a`I4}Bz.{QgIaVnjPnAźh#S*cA nk/Z'8gRY;Er+ #dtiB"vbZB&}`"B? fVi7:%``<\?>xfMnmwL_}rBp(.'f֓&X}Q:.N&q|5=4F#mO $1RW&2BI EqQro##U)6(yڀ ܉_LT*=/p3,QN_cK\!ѯ ?#OS@J[VFD+x R}e=r MpFE ߘ>vdۺbV߻"świ?m| NPAcAUv8.쨂;FESJN^:} h ,lv]7d"+gjj@yBI/PK>-<2sݨx0tA7յM{>~Ș|.}]}{iѬ)|nC{p}rNyN(٣yWڿפ'O-7~,[dc^?7>9zf s5]g٪z37<._ |]snr|KæFyuز8 ϑ&G >t6vps3ek]ϗH*{޻o6̔Gs4R|L_OI>. 1_{߅Re X- b Q!at@ٝ`%.<4r 6>X=w ђ*)fB,y7u%6ԥJjPvvtط2OSVYwɌ|s]$\$謠V6!&.grQT [eKAa_xW;F#v({Hؽper$.9łLV DWh:j"OhؽWp&ܣ\[#j&bW)ø?y';SWc*-yF=%6FL. %#ŰWKy6$THg~օ@:6Sksճd=NIBB.6\$)3I)gԙcʜ{ZpnU~y%}HS33t.&E6ͺʅ06`lޭ*{B'E$`T;5Hr TU{t:Ar^n)[[(, #Zٸ@fRS`~7aKK,Y`2U=[coMbVHH1.zi2EdR̖uk]ZlMr jOP=y>L ^Ym䐚_FttYx!<`l/}RuHϾF'xqXs9I˩Z:F̫:~3+#_w@D>;h1'gpId-:>V20>sꋬzuQɛcgouq17ՕR>WJl2%''%XϹyJבO(ie{F%d/2)Z\a~__a?"6-#+ s;W^GN*lxzjڤ16͹$˾d* yo!lB;yۢ_AU˿b_}YY؇W}K1vN~M(s)<@2c6ב2?I DiwӆU5{xQY]_[ EğW'V~ͽ췓GS0%PIypL x#S `]>TdQ`_;K.` .vjߺfqueP=؜M3Czr,YJ3CVvj?pCmC g^ /_rFug(1XZS`.a NvI3ɉBeKZ7&~ -| ~{gcGwy,X(.*eD 9/TPp_BTzs5@sMH[X\g\٣rmR0c *g Oz\'c!WVu=\~r N+9o}E`&]@zKY"2iooV|r~~xW}1lͷ(j2L|F֓nO]޺~Vy˅Q I>7^X!mףxYr esg.fC3Iiyb#=G38SQ7"t1du򵤴An4?SsԀ`~:[/5VO~qO |u{3K_Ow5#^|PeKޢu=˝  Bpv:-jt꒒xZj_84L O[Wv=|/c.*ڼZYH'+ 9#c `Fz]=cΑғ`=o}m yfZ5H I:tǽ*7n2@O;Ĝ!QlKTu9=VGSLWUDvyybka= F-䪗%'"`YT&]^gJ $J5IV›W@./ ܁z*y[ָ47z@fmMzs4y|XwE4DŽ~BY@oUj"yhZ=F; y1"Yo=(9mDɸ'[$ΊXçRYcť#BeM믅aBگe"cI=47ci9y6B lHfk~Yiz Ku_H_UYHYŔtӕW/*C v4^#H.ߴr l I{ Za?qZhueIQC+Vvv,1n-0,UK / ೻Bx G=*,G+O܀gDKB30%,FK"@f% k/S`vP$-Hg}23SRFF9\;۪Q@ijH]xׄj(]c,A@G'};a]a$sK] q+m$&T)=|c} a&{bx"&#SjF ͒o??ݰ#|PaĉE,x}2|B{ک룷ʡ}AaYl1^ ݹ|,]lߏ%V9-iTCs,7+Uf}۶cz .TzG%ĒٖS Hp3]UЌoO6IO[9y ĆV>Xe-ړb|kyl jW" 0aRK۽gVAAr!r 뻟KzDrIfطV}GC3!7%FS0"L!0B=%uC@YyO1v(c ]9`^3TK#h|yfέޠo&E\tn AZrRVUgkܜ؁<z{XA1G{ rν+z`1ѷn~9ыo<[\ Ӻ|H񓷰}4%w{D}r+qoukHuuV-QΟy,.8DU%ȸE'fRQWcGDs}Ii~ByAfxJ 8B#4J ~K^/̍t~GMBxP*xrXhcScɣ!K-[eY~Rءr ST}~;N%v_RZDi)(f#ZmsE8aԕփ-j*C޽ _R d&Ec^)2ψ튗7 95+|ߋ ]})?[ͫOl8u\tz<}* .;ί% :Vb;&o1vf"Ӝ f̾6m"m69ax1Wrcu=Ym 9^}WTò`?\!M T^|;Cw=D{Y=4!/:ZPN3egkDI헮jd1Ak |/U8[j"ct#B \Rumc{FJRG>qL)wkߕS"5tM+7ԏt@zWK#5A~+~k2O=ES 5@MȎ0Em$/zrC#jЏ̤66`j~o\&D]sb"cB9C(%nut~wt7x ^{{(?3U8߼8ne2lySx7TH(k &|eTF ak^ƀ0-Oϴ.]š'#[X26h?/b zѕH.ݒ8U\=U,3|3F%=x@"Cg!67t/ AV?gVlZ#|Vxkxg|vok}; u_޲Qq=Oy7ֽOw=ߞJ `6~ ^̲zIQR/ioXډEe[֙ I I|g~%Qx L3h_ml fU tQ25E1JO a4+'Ȧ ]U!3dQdݐ=uB,\{ P]ڜvɢg;B]l>b%-1#Z($3)hiyw 9T6r2Xw"0%X*Jt_J BVzz1 }$kAstQY35e@JP7_| ڻ>ߖcyَ|6<ՊkJymHѯ6dRխMyu~1!J6g7_938q rZIfK4 K9seI#ΊC7)1j~ԋŻPջ'f|1u~#Ӳ$Pq_Ov+tnv~V<ق'x[ZYN\!gJ9D|R*"17oޟjղ72PzU}|7.j*(}[AI>gH` 2^6G+g1#ek6]~cgb΋,11X>JTFc9#â 0S 5bU0EDvxO-NWW4n;;; @j@Mz-W=aJr?YcG_ah aiy~&oeM}3p^`25Qhޤ)7,qR, ezMx9zgSK? I wD[o7p3(;ōSli(5i+ EW6c$ ٝLM_Ouhrmӟ\ם$DD6XiӏU@,^C^;^&1ގ{fʹڥ4>O/U7Rs9TΣ:yяx ??#@LÇ!a ys]1".̀ޥ8JOt oS]dž|ͶnR))i*+{s_C͇ITTLnx2Th,?ÂHQiI)lXԱ ysFJٖ|Ů3gc6e4݅az=p4K-7f>?R̯r19瘝:hjO/^@7%gE1xJy[&e˶g>r.ĺGqNܱ7_IoH~7\Nja&9ξZa ؊i6|h:Oڟbt7;\,0A 8/~?Y&08P \ ^Ò)i;6#{Ed?Nby]/_9r8q?cu6pFCr'ؔ49vOd`x8W ߴ;܃|87`^BDŘxMU]a1e}a78p'Rr9js ~[#vIތ!:3XW-*y圣a{SK8k+u/01nlx(qMNpn2Q/H)ACwW11z '}PoUْhht&RZa& #UwQ$,U~0y.\{|6 >_eǩqoT+W?I~}g"/ԝDJSYm5#߳[N؏5=k nPQPZ7AGEVBHm[pi0 y.-ƜŮr74|oieMG34x,\נ˨9&w#/]T5Of6;#J҇p| |lekt1&«O`1~K'^Fx"LDϪJ{#NWԭBhr^Aжvz~kʒ4߰o߾0$#98gtqyOvk:Z>p ${ӕ}F]\ԏJ*ȽN |_R4ϙ縭}5XZugM+h٠ov @v׳:v^}p<Ѥ# e̐v!Ȏ8q)ut#G-1Eo.o4(/ҿۏ~ҥ)Wmj/Xx5)I;S8%M ډhkRU/9իy5fn7_+D?I\d:a)x3 دXgba$M7n_V`*BrvFCHAIN.MAC {[=wuvm>;.uO,Rޗ.sEyC?[q\8W}iY4Υ HUqс{Z!A<-'n%y3wCڟxgƅ"v[ANCPCP$p dó{l(㈻$-o~uOh4̄ķr6*ťVq241ͅuڝgm˸5ţzU3M(a э5Zo_g޳1O.8E͌;R~a8ۢ/ J [Ӑ5awm36#Z5gy`u^RO}/$3"Y;&tq.71B)\ƪ Hg.rQ&/ٔ^1/Jh誼x& S\x%M{SU; ${x#Z1L6?ngyBc+f+_sKT]]Sw V/[HU:-Ƈ/, .dͻ˻W@[YI=~ dxpR2"snh5T6+ " nսJ|.r.pGVk,c=ݏ1Q)v<Fd%eHצ;Tɩ.r5>}g"A JLe|s][Vr BD^Qf@nV@f}$+-x;?^l?$r8%L ]!5jHm~t5fdj%\\ 86"!\<S8cۈ΢o˂Cs%,o9sGNs˶a 7%yKj4onOd] iR ttck#Dkϣ١s1w߰uc7d,,mJ+Cv TVʙ%ߕW,PNd ;:-k hE -E-ݛ:V41 G^Dg~` U/ߥʹrd] cnQU^"D&{ c vDuN؅I`_w]ϦD N`zw u̷~}Qb#RtX)=Re)$픪Rq '&ʷF[rxߖяt[w̫\^!|+M+5k*&.))6LDBYUR]ʼn&!+1/rexG64amM ONUy)E$*::ءƪb:~{N>-dY]=5O!ea P(Fջ_N>G"Kg0^f?GC̘YdJx_30+΍#k|%7-UB帉ʬU&IH)[;' Zl/ݯ3_ H;G4{(Rb~,DqH9naGPuԐ"=n+TEJƮQ-k3lil?4> Jsi:5>*'"Z#)p$6uGdv[ mle4Ƹucr=[U*2ad=:Nۙ{Kڪc@8ݍ/'u-UV2tsɀ#S%_\ y?6 ݻZwqY\ $mZ"60:%wvx(.=$9u5%`^řĻoC3nPW\o=ZXOu5Mҭwl?#36vzѠ o6IBxF/؉<)ϡ ҥý17`Uxdo&li7f=Kc&$tbnCjY;9E-/V8Z,SHKaotc N%t[5UCWo $Vy p$ /_~] 1iZ 3f.AhiL~XݚfĔȍOæْ㹢ډq s<.Qß6 夝&lkT׏*3G چyo }{ԯq2Klx(b\.> ,7 ŭ,3P۲_EF~n71..%] m0(<:MN< ,K[ vnsu6 j\LYq4Cz`oݠsGR8c8em|vͲ lM_oraI3(1 l)ްX6OG_L*ïk /M/3D&#L*=O@>xS!;#OЖAjdzm'+Jhx;m54|f ҍ2M,*ܻJ[:;gϤo %X$DZO?r!T0ࣄ1,;[,4Cp:?LRHӍ؛vؓļ'_- JK'4:#o$ Lz~bⱖr1^e:0koBmJ]MM,R7x=>1VYku*`Ά4_4ɯ MvzSgo8Lc}m x7^|?0g߹HNby 昭CY8` dC8$8ǵ? FLn_N񸤜D;EC:أ8_+YZf?ϲNc?EzAͻйYUAzRa: $][IKayCL}Idfm/ƥhd _/ggGi`ll7΄J#mm~v~Rx>r=ԑETt)TN:}Ϝ,}CeW8rv2xnכ#{olΊ冃h|lO/HaH6hY'DRaNoHx.c;oCQd74˕KV-YWR69@+8}VJ{@Y=^|Z;C{} Fpi2Đ/ڂOK ;AtJ1ɉ<]}w?'M5#$BaAOtdVy^O9nEBv6c%Y2` aDh晳~ \ a[@_sPCZ\95IN}rZrٲ+oб٘:=ֈ"y[7噻N,+pv +2~{*P :,تO9xɨY{oELDŽ}"{}6H?“o)9FA:*M{F>fe)/I8:\o{$Gj6:53. EY9-z$wr<ŮrCmW֦? so/HniCriop)h|S|,e\>4z#WDZqحeŬfo >wM;,ϨoP@xUWqyu?n^6 ma2~n%,Ӑ1R2tZ2!$0@M{"٫uX[Yu5:yX&4m\b;6*G72SnœXBygehwjȂ,cCޘ-=.n<1հsI%[; [G ~ox r |}>&lCzUkcjzaN8Ĉ/:J[r6O$_*P"d%^ϳ_<`Qe~,ˢ†~'+<* ߁ GG;ZF FK!PR/H""T.םxo)=kƙvvn-i7 - ߅}ð~LE~u >mI'ɸp {2GFO-`ޖgÕ+ rp[9">^g˜;O O~8T Vd1x䃿c,yf| eoW|diBcGi&cm"ݫZJ-~tVKѝ֚?0GslCjku}iCt^p#_*/>:ϊr#R#w҈.RK8WZb#s[~7d/\m69U$T-h ta莔Va+G}EO9:Mg$(Sjz%<3oԠJ/K;^%t(k x6zřrԇ \[g՗Yo^9hz_I=/#e@9)skmfB!V&qHq];sff/3(rtĿџ?螞pENkjsoOD2>H8K]# Oo޷W#7Y@ Bj|T ֝gW"ʇL?>{ d>`6[|ѭOpa(ɿ|o+iOדЦOJ<=F[1Z{-TaÒx:qKa nqrtRF[&`yK\Bݍh9_N_\Q:sIz=)Oh~R>B˕iʼ[JC h%hkujAp1\ׯM}ae.Y:GS90"}M+Sr;IרIP@QۜZ% 6 0F+ed%O;|%=&ٽvuG.}TwWVw5kO͗.pV.N`tzɳ^B*(O&lڽG &~2;rz0ϺZ[ٙy͍?4HDhzIǰ9 }M*q)ab2΄ǼDKG,<*Ds߳i/v_3is@.3vKywYJAC&>0Q0V ؆Fܴ!aE(? F3JVҘ+?5f|p^(v-+eF#7QMуn+ x냌@5㱀w/Iy"hVKsGl 2 J3<=JMn;kO+] 8$NTؾlDcD1Mޭp=@lG*u":a0+W*oBqY`g[K@xT]2@뭾zW^k{y^M]umqɧ,ki fY n״W-`b|ג!c(>B|̈́ru|Ȓf{3`2)F7){Ot_DJF! <d § T3R Tr_{XC(twJ&˄'eq].0S[ͣdCGؗc2_l8 #C& ³7*u#P{^vү b)M/}_/ k@'BX^'6ME:nqGY[aO,錔J*iU%Ȧ+~jҼTXo*Uk Q Tq›dmpawB[pLk'X 7 )GSΈӥ>[Дwg='SƢݸLſM@&S|b)T!WNc{̐ZUxB d.Qw,dJ}O[|7mr*MhDL8)5K!)M)j{aH.6ZxeM,E׶xv1Z*ku~V gwنo e+eT :m|%&|F@~2hm,3 @un`AiAQ.n&O2kZl݆ypfW~6T%Ŕ_:8Dod.&'wjna"!:YK ssNGB2;hnE(LL9ɱCD5mgfs;`%?ZĦ.?_ ~D6m[6I#lK CZAއoIJЌ#s|-W]o1ZIM t֢UD~4d[ Lh+zVDJ3y ~s 6V0? Ԟ^6_ xu61JԊbsW:JMphե{Kt Yh1\3q՚#iyB,82o/(_d9ܪ4O v,f- OoV-4 q' _[?K(tf#fGNYI,Ťɭ@x@q;Jʜ\(FfF4ߡ8F⿧.@!^ՍѢ&]_n*RjpQ?64IQy&@"=ر46l`KF ʰi'֮Ue4sV؍VlS2QOlj`RH6l"I_S $0HG9D{Ҩgi},۳qDyzQGޡQtTWqV/Xh#ߪb1(ĵӀѸwWn̈́fC6ߞaᵈR凘V-˜hOwC . [?t]M%++û~SPETߨO8>{=? }]Ū)ꮫ9+Fܖy^ە5rBY.tGi^HF"ɡ3DhD҅@(&"5eivu&z$x[K5au5?~jŻڕAeR*n2pfC-/hyDmu3XtJ$Lr\85NJ%;뜼XrL+eH3p|4hf;}ڒFJ_VPL)7/)z12P,tZڈEXd-iZ_Y&ÕKOs8YԍGLڶ6U>!\l̢O,K#SMHY-NF;C*95d|M/> w;|v|DO6Ѓ@N|5)>tFbZu|7Z>9|vq_8aJ^+ j)$rYΠ䟉rtKhA@L`o '=0@lVR"k |r߽vP}u:d^kE8; %9Gc|#mV,Ù$`tݼjxyȏxϒYgD =BTW9.f3^HH!{F{F _uc|%ADaO6 6z9' ɩFP̸z|0J[1;jt0V|* oo~мfE]1އՏ1^=_oq1r u{&Ӟ$3<+ sɩΉm.p'/Gٞqf&BOXX=}z?UNop2,z; HULCc=kg3$ooIi*~ b$oL/_:y,2PŻM?vKGXx1JѸf%uW%>CAgЕ> C]^M\嚹:=o\TWf#~qÁW\R4{WH%yYԈ1^#]߀/t>}Jx;70S>̌cQĈ3u= 7 Sl i}U#Uo{}ѮMAgvLm⚇#y38㯲mi>R߃*B.g;ML~'g~,}$ksvzToS٩b8r :NJmh%i@ aڼ(xVVygJ8ЅR,R2sQd!0[IϹ.3y]?#i6;S:ܘ9+ s ?@`}kL7|7 "& @H$6jR \ސP21 enfHsebdP$\3S>@4%YM{/ޣ!El:>;|FQYŢ~l:8|-WMktaSQ?yW+{=SX+6R 517blة>W6xv*`)rK}kd<5OoFeds=AT!h(Qi{Uv\UviUӟy ΆxIlIapw\3s3K19ZNK[iɳa [m5񺖮bjRH\wP?p{M| qj&= u)װv+|9?]17)z ox0N 5:@B{ou&[rNб]i[)k#BPU/%Z\+2/k'զA0/w#"La9iL riDޅ. Uv#nOɷ]?J~SA0'1? )2v;* 2f 菑7/.} Я\օl][H ʖf@h_ra7+U86\ݡ YuT)D,1M_+n/Exvoo[5i8C7PP=8 ̦c[IR<88N]kЮ_SD6l Cօ\VɷSE ZYM :s&&jEX)%ׂF*QS6|['4k: bz33r -N؅hݴmΊր;*: U*UpeP^-[An}6TTv?4\{A c"RwI0z0~-P`s[X,rM]2U!7saVQeDR~OU)~e[$[rT'vn!&aFWL\~U{@ *4Ax0ړwK«|uh?c4ޞXoB< s ӂ7WPGeW2ÀI V<$օ SZ2pCH,>g|%}.ƪX\*[ކ)Q`s7T{3ӛAI@;M3_(a[M'mFFי?oy0ݙ=z$˦$,(Ϝ% 7 r0]p>8#+"P+[hfKXyA"ܼ;OosSUKﮥecy{^Zx F롆<q^d8^%Ih,. {Ÿy߰:[ϒF}Jߡ=x u?6g Z_Ec,3\ ΍EKS}CG$+zt7>-!7hL_IB(*/ÃD8wQ?gl03:۞Zi7NSerT=@-.quS:s٥ozq jy7_#lG{m fE* x3t(&gL՗8l;:-~\kr~Md/Ѿ36hr}xj^y4:"˜8fC-K7Fd9?^_jPOnӗh(R8"wcM)fe K'2h p ;mHx\4g`5&:>5E%54Qyfw^]z]׳U,U?-uQMg>e6v?J IU\ a=& p4 o }FeDz+x/JzowkVayX|ɑ`!Ov~'@^ZHF"-ǜƨݸ0yQ~_a|&^ҫOW&{"D'I5:^piF\.ݛߎ;nM=(w HJR|=+D GwZ|Tq+T[6Tk![@ ³1C0ؗ0(U6E7n:rboegw5~k8 t l7$ujG7&p |?_Ԯٳ|CŅٓKWcdĚ#urqm3_ 6aŚ-'O>O)IDፗ>!jO9'^/^*&j~;wqve-ΎMtЕm<YvUQTQp3Kb,rsZKoyͧkhBCR="bzG~b}w_+N r>CwիҴjf=xزL.@9x{OHuSwrUst0pXG~h,;Ҍ`xi}ɶGUm6(/Tpc9"6 l/g)@pyWM8ꗢ|Q8>S_uRUP$Vy{u@.nyClYQ:Ǽ kgd i xuim_tV~!?wG ç*3=Z O%g,Mts#;U k# vw DǓa*ujHopc[1W237~`/G-ZfXDTsײ&N@3'Ձ۬) p?l@}-I'l4Vl 5AlE|66%Q.Eƪ#U’ashY P1?]\vbL+Iy{GgRYW}_.X)*SETPOb[~0bKa5q~*^3Z3Ign3l?oƲb}y” ZwboSl{Z|Q)|{FmZ, }牄R(,0CtCbws&D6=dNѫ݁oKo.@. WM>ɡa ΥTbkhw$(9b .%cM/s{xzq>O: :@@xSM CZW.#ۮTmkըjAV=3IѶn0dlmv ^k}.pݔdmq{lj-ICAc/ޘٕ%^&䈀dA'+X$o}оymѫ?/w%SG`7,ucVs)[tE}gs4E{Wط~@|?Vơ֡U=v6߬t~~`rGLF0Q׶8Dޫ MV7;~\AUku\N0 갚ec85 5|ti>2#Ni=i+#2 Ȳs\[_9 KGP\~dw)rt8$eSW1~2cȇ.t5!1VNbT#{9|[2|L"6X0S9 .(eonjَN;sMRۻ'nM<6rWwL~!aM4k3t?"oPgA% &]d4;?eïNhg[u!#$cxMsIa+YHBw+ڷL7#JV-&7F:9qmoK׎ӊHV, ԕ06R_ܜoј-ݼulhWw~xjmU(b4Lqv`cxy-ͷ=R2|eq V0N}ٙ&}fZuNJ;Sw27!G2E&‰:1~ҟ^ et86pJ7I^iM+~ެ7K~2}Ԛ<7 Q8 2%]hu^$ e-})*i$ <%܄/lh7Qw;9aE6}JݳYc,of?#9zZ= Ġkr;ܠQ(қGG IHbWg?.ەB7V跅'*x_gt-XHG{wn6;x4>M ྘E8~`[Xh48|\R"U'#gޢef}q,x\Z΋ˏ"uoʦdqO~мb@7PlF>4׹r(UwN^6Kfjy2Iaj7FO#̻+o=)W'uT17ڿ.~mSrO*K6 U>f[/b*C'N_7NO64.-5Kt0^p~$e=g'ZIs02Sc jqUs2x޹!AOk6π*!fiΧ)=fWރ@-{O˂xfrUcWNP. # /K&c/#$goYyw9kQQ,gV3XdAM>OGyX\د/$iϭ_p4'q>fl47>):N}//0#FO 8^.pk)6ޙ |?jMސ*tb6GT@GFUɮŎy#}.eЉYkHk8M?k!#1F˙]XJi ]bi v@.<@Fz ',#x8&5Œvkk+;/<@ t z챦BKS(n~COo$߱ʦgD@uiY,X1/D9]Udxg_wÆ˩wOVr?V.ȚY' ЃRt9?8%O\EzR{ LZ^: Jodɾ1s=T~Ӏm܁:t?y > euP \ v\Rá?p任Uú_vd,/<8SX{XES\QB9' B1@a] (活^2o$“iш>VE 'οP%_xgIslETo[\ "sfKFP&քsۨPk]Ǿ{ ` -Z_C6{go_f?jTZ]?/߄hs~Z0{fc)l/F| ?y^$z \!$2l*o 9WoYEQW 2tl|Yi/LNfX>ߎew>]:)U`^~ BR*EIlQuά.螘@JjV>˰B\|:QqG`( nM2[n5ց9á ʮ7^N&*WZǝ%?{߼ԅo8zyl RW_.J}:lM$OUFԒWTQuv8ȳcy,ep//Nv!`6zUS&,%4圿Jψ3 $/:fڗٶ G 012"-gj[8Ð7.0ᣕ?is8΄Vk|wR[jMID[˺&FGu] C(ѨWv+w >I9Z5] cM:= ]4dҖlT{Է̫ɇ"!^*h2| 1~>̨qcρC>ctPf9 (;d2a4c$/$bCKv9g3^-Nφf-Ȇ" 4-3fpIdS ~-rY̯B3#x]hq`U,tOnl4V8O;!XA4 1^sK؇IKw6Z*d 1h'5#ϡd^4dicuY۝>φDƂF?qyAz&Cߥ{ֳd аO ϑ;۹bl+/Nz. C|SֽG$h$Ib0+xWdK?%}\|n׫'~;b| 쐱BwiӺ2gaGʄ*BCE<~ߏآM^_m(*iI־0r|WEd|4*#wdC*rO3#4YR/HO#OŠ%MCˑ ȡrp"51몬Ȭ 4aFYNc˲UH[kp@ eUx9ã+{iwvrJ%|PvFORTRAN.HLP ??_U;\߅Wsjņec?6}^ߗw u#J;>B4NܾnM(µvH׵&[߫YL)L㢽:؆o7I[שQO1-Oط01LpoțW>mϨDuEU"14sz@ EP p&(>FE1{~m;EXg2k2 rehØO"\BeZ%@\d>bFϹsfTl"kǙ#.>EqGaQ*:HJH= JvH?e]TosujZM]/|G^>?~t׃gƐ @iլ&LMqKK͸ Oek{U7֟HsZx?yՑ *GIb{$&khCF-#U߶_J9p{KX4N<|!+t{]*2(P~wCq~j9%kPY3C3M/C{llBfme{|1Kud{x6S7axM͓f2+#sr~DQiJǣb}`ZGoҼ|$|X\&s6~g s Yp)շ~ ~CNk[n}O<"Jna=l_|Vl=lvEl_־|ڷƛV<x'#jLq;sV#~KJx(^~%/.J\D'0(KtpSa|y]'l0L|(~/J-dV)GS !~)TB?[^&NIKpY %*PREb+7fb,&_r+w)z@^``sǑvT Ea+PUv߇j6f͞ oe~1l u s"(O둣FIX#k€?dȑ.lo668ۇUü;bS?hk{q @9A\53/{;6ʬ뫮g! x;'~4)gvLic{#Hhrxgq _ 1覆@9PSdSLW{p>|uظoȲٵvMI{䌙ר:NCS⮒us]?bX%C<І֋$ و Xxn9*ף  87^9mس{~p~{8a(xzzc.Z"w Zl+mN50uEg7NR1W#WĬ-EٔMss6+;9u& CDAd2"؝Nn`W&k㍬3Cv#vdW5Ҥ-beNGFj26p1ٮC_*_6 )5q:#6~X[ &fLJv>9մQ|K1=]ျ0Uܵc CUܶ1dl^CFQzkjdNePI1Rqn]Lkv8vvIOINIT.MAC x}n~&fbV{ܻm˼nSfף6o3UME&-,b?}  9 wiVړڅ\؎4[o}YZ*Ѓe?kaflFn"n<c2Ih𡱧T%I gXd"V)QaSD) x$r~&Uu=y5&p_<{D]߮0Pf> :k;3h_,ZI\}ַP(|.w2+pJڦhvƛh~_v)'fgzk@~ J?hnjwv滺jW\-Pl%o@ĭZfEغiNJЇG@xl/xr"NaVQ(*d9nL$A܇(\f7|<#6^ΤCHYghl,+ +6KVAV([!OkhQ`mƏ1OJ`.pQE4UORfgʥ')YF߳~vqM/:.UՇ'*PQy"z"Ɓ?A !W^<^،YT]~&Ff%S0J< 䮲6Ι^{0ΰ^ $s7=uZ a.HSЮY@,pƵvt_tnrP Զ0` HdH ;0I )vr3ڌlA,O]Ax^bXD5O}|G/oAȚJrZ"ݗ?4hg(b^F0;*wf]t1$KʻW԰Oavtd9]?Cq`Ƅ>#hPTY(߃Aj63nكʐ9eΠuJfΰ[#1 .rٸܑ=ޱNh?+ ?Mw$ ORJ~Q|bet̷" ZI@H^9rj<"@'J_xDݖ 2}c E&9r섺g8]!H ?6t^y^_8<ﻄM׌fl0|!a% |g[KCȽ5(?IOBexbS`W*QB"Ĥb=N>co9O]itk a:ޙsL,(枷rB.(ǯD,|d2q%2Z|6$xFiie˭B m5{ 8bCxZ*R RdL꜄ sg2m婉k;Dx 'FkU$LlT4Il}vxmleD.qշ0n tdn]Ӝ\F(z(*[F0 -THx9*Sg$@EȻux,XIh-[2` d6Jl&wnlW՘2ǏJF o>GVrl5S ^m{0iWm0Tck/ fL`Qhbo ٜq?^n,uH!ظ6g럨w3wf&r$7SᣯѡAݮ62yH$ؾn [w8lsa@2^VQPBϳcx6 x<)|$?1 įynzeg']$xG̶/Ś;4Z /Fs :1Vv{ {o'4&jC9w AW!>Zٔ3QsLΊKQ4< d\:VH7/SҕiEEW5ʝ xkd]B pBo^`{J9}/gȤ64xps8?Mjμ>o|6b1Ÿ慃꿄W=yTQ& 褙x z,+95;@.$W羸~Pwg*I>Cxj?vs$coTN]UbxMGQ%vO _+ŻD>FހIsKʓBL9ڱQVqwնem 9F꿒qN|؃[l۴c7z»Ǧ.> 6V|ޓ;#xy?;bpy.aNNt5[܃XO WeurR%!7І]I`Yɾ8^PDQvDTs'E$JD m[n"wkHqF-Ο?oPun?Ï`Rn`Ck5#׆5~"][',-5@5<᜷a/J9^(7!@$zrvOx9pabV0/R&1B}HcX`|ܯ9U<'=B\1'fS"z8lWjjL^qDL ڗ `hk)A_4$q|j^j`LApcOF}n-'")OhF XZ 0V 1-N".v ]fX<֩]̼<|}rF "eU>aJdWwtn].xg%+L7|Ϟ'ׂ*}p7W&;;eN/ Sr R6Ne HnLzwWOOօusAÏ]|}ԒNr~=8+օ=Ln՘GK W9-yj{E}czRmV[JRmO6~Pzlm^{.,qKzLxnZ0$z5] }I *ȱ8Yqe,^+aCQ]R|\s2~ j2wQv=Oڋ žeXb2uc\'p Sv~d*'遼-amhgH?~ǼKΪ9MLNk0t}'qHn?UC؏=UH?KL#E8i\PF/_YA^\GUreO>~YB,&Y(<9wE'  Hg*3$-ZǪ[NRtZe[T: RkKZ~MkD+4X="qLg/]\_ ! %"J(Bz<4[yT^?GBR1)&}&lG ᒊA掑+T%2y*4/(۸?nU|_ؠ6J+u'뉠%[tV0_~ 1p9.ǎL]";Q|ZYmu?$L;SiT03cR1cjnkŞlA7*t 7Cø@sX7|M<{$&kJ^,p *6ZB ձlXS'(`Mex?Ns1[0 MtUp]>]hNH73ng$d́m HfO?T{?q Z]wyYSsVrۥb08uWu…qr~&}cGHC?l<̰0@LP+n8gvz k@]@ R6WARL";WPhs/l z `>cfq;qߕ3o\TUVћy{bEw s<22nCooupی;M 6/P-%Ɂ_>= ѐ6O.z(8Y!!kAR &_dxb wCRn$73oN>T9<{ϯm+|iZh惭w YW\a2^P,gv8+f]xsS3m^s: _*[[ӚUB1fm rl]u=Mwǘa-[ü gR֎1_8W\Z27Rn]vz|DNEV/I cF8 =[uN,wGtd3D{Լc=r/ D6 6Ƃ%*'sKu\|\ r;8^lRގ&[oEj@(>◯Zi1|=q8 UW+([,m` $M ȶ.c;zMI)b?_!EHȸʼJڦBV̙Grگ/UM$#s*SÜC/\ Fw؉a. 9@⿌.tpph! wdV Y/[V !d{;@>eʧD"%{cddpiI#T\8oz#fM.!OEC?+hEb&x/8YAx<GT&Z";g n7u@q׳ּޮoFoMV01{`"'D(TSl'<|yyQl5AM<.|!xx{z=ad9]=k _s1̈)QV]'/gn&SELW\4/{QS*Y%lwSHu1t/'{Yb55TH-c\7ec1z\_6b8*:"HjŒuc;EC8R)) $8[ d x*H2* NI:wЍ #զd(8[iOai&q5Jiڣ#epwT> H;^%Z!z!'nÌӯHS^Ҝ3.N`u/^:Ɍy3M2PbT-*StSσa"kNƔ+\֤^!z.p^%˫UIPqi NGDN^;Cyzu pV92~%y{(6"sp{pz$g9aDw?wHt .䏽WgaY 8w$+a=*iQ}@rIҷJxs63p7{JAf:MNJXd"٪LwM6͕hӎ ;X[J ܿ ˾*=z9>-7L)36!}Zs+"k{N ; whu@u&޿lm 4#"*{=2Aw|ZGl9&oMpoϥ.笞-:o\K"7{~?oyG]OpO̫J+Y[@Jd&ʇ x!4@6xX˃k;J 'g׌k0fM*OL=ABsEḭŋ84Ǐ`}2j zr%}n8m (%30ty`⿒-f`ųrܖ"$W7qHjnfZV&p Z^۩ YQ{ƞ` C ^G!AFS(y8NϮfqzFNtડ SНn;Oa\PoJktgǏ A H$'7ۖ! [CF)CU궯 IZB4wg~_LR9j;v<]t ;21ye|/ UR$U6mvNFn,ڹ `ltsD;j9e/A`y7n}QTJiDSxMU?EV/}&V֊g#|2h!`m9rO\hrSՓ!Tˎ_cnKqN_H`N۫GqCSLܿxd!|MݶmN'jx{`IU'qظ~żk )UsiJ]Re{+(،da<?Y\I*ѤP˾0>̉Rtq@?"P>WP;:q1x!/+.#> lfh7gI?DއZƵ 3*5W}Ϡe=FHκŴQ]%dW޳~Q@ $@^`U{ؑ h-/n6sQ^;0Y!q%Jy*FTk{SWGcƷ>NqdCc݃͟mӫۤv ]Xsr PSW_水; 1e݆x.ǩG|9~8-(4o,-Hj9VycE7ۨ5#YH]b¬g{CS$GT^vuHW؈k$@ oޑ\3 %>!{8qmgY4%7jz㨟 e(Ar kճ\fxv"~0V<7{JIRxtykRI_wcKz)At?I}"RV-3=\P6BfgܧRzC:ȍc`'ܻݮ.vYm֤j)70.\dIQu-d畨a6]9xO/f)x[Ux'CaYq34.XaRᲵ=oc)ovTgQ^1 1x0VK~ȟ;6%YRWAӨc 3fF^3L 6{MW "ۗX1\77&%|E57x9NMuh$=zӼzi4SGYȻvΒ",b!%SM92ƈI R gp"*]OL hEɎ!g$zY,r hyF- L,;l˜vLUNTB.MAC x}n~ؽ̿y:7/[]W+)x sٗܓGD} u1Z76XpH(o8׺K6wh 0"L*AbfK%4YdDTð-P]z E}R Džΰ!T aDQ07j.%\ioD9`onr<;<,i3M i;e&rCbQ~A)ex`#&tuI8\<ʍYYkrWW"I~'XPm6 fL8j(IMM& |ܾ3Gs+~ b(Vp|Dx[? woY7w\D? `ݔ~X;s(d-鏎LgrnU.NFytf cڕ+fYjߔ,6|S7K>Sf0IC>+s.7*t{ >Bj+>[u{h6JPMp8γP6΅pmWB߷U3K6zSuPEȓ-628x̸W)b^JmD5-O DE=ܼ32&$(_ZhJW>߷;W]] .ercj[)z >r4:\M~%*{6[_q_Sq9_f 0 7T?_a4B[TEX['%aU%਄hz@؂!6Ji+RJ uhSn{DQE`h]08(zd4fH]5pqxC (h0U]c VɑQ[\ŎSMG6_ //v@. ^efTɑ9!\N끭u+7YǗ[ӿ"]恟[K<8@amD@J?gSLp˟'KY s_u,?k1jу}CY԰U4xzbQO0:L8rc~::(nE0oÓOXXVxYqx辉#78h:?(U_pʯishu #^u-洈8aZ0qJA_=߂NT~* @ߝ,yV呯Ih1|J 5ŗ#5Gdhy=z99+c骲O&-͒tUYrbTj!u5UziRKcd hui8@"&J,D2EXe~󛙃oH=DH=A._ic ֯0ԚĖ7ٴgӞoV Bp <x[]B[X=1m m!gsҗS<e5.eӮuh`̑a3e@G_ KfUfnЖWq^pWdRLsMQXx;gun`pQ_4aYvR䴃8c!$:k?oP:JhY~M{2n[n4<{~|%U>Pz_d %utT])&ιOUfdZY;šϖF=ThtYbZ.-4 }M*&N4nih!}V1=\͈VWӿfVt!vnI3'͐|=E뤙qEH/Y ˵[ ItQ_B71jC@gh Bsy5@ya8^3ou^|Eo%*g1/n%9@^dߠr5rCcӫx (=ˠ|N0ϐg@) ]= ܊A.tX^{5,OAɗc>Җ`_b$.kFe?UR*n&ffsL"ke_m@S]kF5ghu=kH 0KVEo2=8;G]%h*e ^gb`6lଋUA4@F:YX W`4?$jR>V,:Ɯ 2D1"xKAJhE>Mhlٻucݺ$ 5J5NfVz(}57do.2:孜3аRDWM.%Xǔ vALOGhW;_tK/H1tWNv7Rm<3P~BRꙅv_4>4Z֐s[]'J8m{>qNigNGOy=#Sv,}<5>TVP]TI]g&CvuTNo-\ rf}!fN/Snp|t.%ڪ kiUV< G=x TXv{T1؀kz,6nVFؔ[ $ZZw{s D G]pO=  GP+hEXQX=|%ZۚkT <5ie& mr6Šuu&@OW( N~in 9zx:pa| aIFVz d~uАj``f?~*<:g '_ br3A7vd6cC(lgܞ%{kċtΡ#F/)BֈPvԘJZևAg.)Y&@o`cjE,^XPRfh { Fv[lfu5:fތ (3d W#F([|\* ;"_PC.ed=6fsc g ]+lW&+֓IIyL606vYY_K- %աE?]=bEW 퉥/yj6 n/{rL Jn8HҶzS'&g7M/\4&gMvsNqajGҾ!1@Y߉j 6z!(Y;C;[(X_(.r:}Nz,/B7BAϞl8sSC; (L o|A]SgћGpu}I7.vI5L㍅ j\aOϊv̮ٿgI4Een} hI+xn Poڡ( =@lЂD•|maNV9:[B_3*y2<|~b 9ߓ8Fm.I$zeAnb(#h]LmVqWA&@}zY_J5J/n-[=. ߩV:ڠI{Ƌdv6bjíy#+1*AOIdn = ~S0^xF)5+x`amڷu`P~Zx4W E ݌$h4[#͍Y+X'cLx,A5 gInVt u,!X@.g<4G϶+<`_ƏJ<D[@ n 2dʼn}#e.A1 ֩ /OW3D[v''xW 7*h^#`8m%5;lVhT#"RpΛNLtVޅ*s*fv/h%i- ۔R.G^iy6t,]Ȏ  p p]0~K$Vbލ0j 05ܥ*,PЊ5IQjTP-s84ҮRV<xTpd=(鬄x%dp_qKs#mU6c?}gK~+JpF. e<f UuoE86[V$םJ*jġ8lT5I=s:Xy*B*##nzh71!][HÔOvS+ 8ׇ8l-p Zgh0*EU(T*&_S1t,ĉlσ7#*S tQL:fֽIXO4Vֆ#N{9.urV/r{(|Mƾ1P~AZ\%#5KV2p';U7o8ugIuqt5f7ﲩusbela+HD *grVl4Hv[+nxu0?/#IR.llFҐ0#ԥR&n32Fpa ?³'p^{#OpG9:„/.38rxNç{>Vݤ.r剈O)_. صV{ŸvbϋFgC z'R=)b5L&,CQps;}Rb'j>,&.9idAD6-|(!ټѿNr{?x@~y (nhޛwG}+89O]pQGl>q48l?fEyf7$7\y㵽D* 'NN-SOK.4~ u ::P: hWȮFK"O^ql-UWmp, b}|1@^+!g.>%ח:jQ}4h24嶝%x\( Mb[ ݱ R')`Z& ?5ZE aL籲etWlg#QgtwZs},"̺li+J..>p*\;`WRSyQcVfW ۽WXc+WZǠ1qgdGyv{~{Ś[9ؓZ $)a'GYMW+WZbUȠ ϲ"*OGol_9l,*]wYzF"0.+MqKv<9Q? u"/Z=DkQ. G,FZQHdUw{JQYX fKҕ0 1%̺?wθek~į.eA^opM,BӘq]o&6V %L5˕ÍdGݲ6]\toG#9FUlKFN&nr:ڴ U@[ɦI Xʘ*SzCrT5dup*NOſ[gcfZE_еkO/wE.XsF94rjXl#~+OL'] U8';ϕu+UqVHnr9݉.q?"{?9ȡOa>."5U%^t\ŝ_CLDlD$i:'\S<8) aW0~Ͱ[]bsߖTLEˎ&Lf hx'~R\Mտnc<ݟBAz3;Xc bQ؆t^C}F͢omv޶H 7(a|&aAK VO>-n:{VmBUZKϾJB2Ck: adDe!S׃(e8L)Xʝ0h[WLN%?Q0.%Vn zbNkC< GutC3ˊ|Ǝ8ZM&:>oy` ثմ:~l/j~b@~e8*bW\/ءLmy@k>-]ub7#6Ϣ2{ E$sk܄l6zSc$%)&)~ti}7>HCLKWAka_7e?GKe(?8tV0#0ߐMP!xBů>ZAF˛KbXiF}и=r-|Qx1&ξlUGUiTɳ9ʰ>wcV1k{5a$[US{1ˇ7C (B9[wѓ=LPJo F9~qd 9]>Ėqwj=[fWtRfB|9E,>Wi-.޴Dm*ݺ?s.u<PweMn5nΙ>pX:ݦcw }Vf;"Ė[w n(Y~G.?hMkQS"U1o>ZZGNժ*l/lsIM]#{fkoE9d@5qBVs})KܶY])ט-ʭ#?ӟ]?^tjx#"cxȣaO{1o'D\AEN$Ke&B9c+54 A>WO-\t؊Nޱ_c>隵G[F p }ƸܣSûfnC9kƀ] Dhj&\<#4ysaDsIt]l/J+T< [)Lk &Za.FUFd5&7m٠G^(6p^f#ʘVei7x_'>+~X5Ǽ#U2yV(~`;`g9|fWտ7ae@VTJE)ҸWW&&f6ص$V@H˫>YYE}eF=ɝzdA4TQ!E*o{jVSa5>]]m#à oma7Yji 쥡‹3f>sU"Ht+4[g3p ijU1Z{hߨ 8olmId7yn,k_}x~pH~,AS"oZ*Q{>cVoJ3$U"q7Ԩgf-Y28JMwN>@4#+DT2`am_Mooܓdo{/$< plirv&%벋[f-J6AK]Ki,~L(6)&Zt')ܝᯄHP4E6|mӼwVTY|THXoFϱ%뗦_XDMQFRr?}5n'^ AfrL~ʜKœXH{zݾGkO77$qOtNڕ3A~[?bIaXTJ.]@hOZwħ1ʖWY'籐CQIV^"Ma[(כQuIn {K&L=ENbMhkRea^Bj#KyZ/txnr N.E3Y@}ld;#|nxpBH9' *F Q`a"U{x\sO DU~2gz/\pwx0*w8|O6Z Ҥ<޲&/~VυS{>*f3&ʧģ)6Q99}gd BHWᷫ>+JbyMUF>p( k$_pjپL)=&Ϡ#RW6 9MX̕n|py| Ӌ~>ӿ2zz?ofa?Ă d$1=P+ :pwxCO)@ Pr; (@x|%u]EZƪ?%*|<~O*{CK t)&zӿwQ|o u jp(쿁N(Z013SN_M4{sdyv-A!rV(7c{M4Jһi}|4Ndfd 9Zul%jƾwG1zp(݌njp2)_Ԙi]Sr͍:mI/8Th^b[fz [#/fWz#ܔ?1]HL6X k"t0ilAxZ|w@LQ4vb&4mmC27K,[1iI5~voȦo ZS+?p7I r.p^ @uw5`LYh[q"}PgF0{tJH?qD g,J闾Da|nL&ɖGYK^AvJs8)z( 8M@`vTEN.MAC 8q`KN;m=֪ NX܃,ұEG޺ң/ %L{GZɖ$MӬ<o#:^'U =[QMc=:jfڋ̉!) HPNs 8?d Sz!iJ ;5pOh-àa_WT ψ] oa1?|xzt4$z `l\zl-nt8 h!,/eiU^ߪϵWa4|ix&#& :\,c >d|~6s5k8iZcoޔWA?\BC«å^%T[-Gr͕̋҃5 f @uuBVݥ;1=`HKHxk5{])J.*11Gc=6*ZǮֆ9ΖW6tOB8c) f&&{l8zX4}ԋL%cK+yբrrui[2Y׎^&\h( ̛]Z?K3QK4c#]E32FNy9hzՇ5| < 7<+d:,y OW磿v#h s%]"L9ZY@桡E;U˘CJ}:шa7%spj[>{߯h(i5,w82 Ě;sިmOHj{a(dS0<㆗ PVZdu>s؉HBAa.^Vi)տ6B5{]X}"P2w_76f~5.G k.5[xŨf>_Ekػ=& LB1$ia.U}C\jU9ws;7u.8[qbu6Zu._)SQGqcat&Oilğ,if hǦȗŇ4 ADžSu3UH׶Cו! % VظҬQg"7lپecO*BįĄ¿\>VOc^ a -A&XM2! &Pa=4A ּ %e sͻd77fQC#Ǔ^5/ʻ"ˏm`sh9+*8A@ ڃZr,U>p sGs_Y"\K6w#n//#ɗ7ŞQr69`ޕY.dXpDC>1bO!fS9 >bhq%Γi38Wm>YKG J\3BO{}>,$qpHEdmj@wWi\92!p~RE8mjHeT%Zj =,΢,'O FQbDZ³瑀,g L,d9ɝjhQX Lzmrcq3B᡿ RҺ4v4\)4:/jj3)u32]ᦺ]\ʑJJbT kDh?Q>ʺcD2 #$P? W0){q-xl(~Q 8I7ztkTDUD[UhZ߱R- a25.QRn`E:g:{W`t:D`hӒExu"ʻ9g վ>V"|_DLP]\9Qꥼ\oSbԾզyWlXGcW(e=O1s-} G۟4dOPg`'y*\rLUs 뻤um|H{Yʹ]_?_^NuD?’ZBxőҖ}“*`Ln[P%ڵqn/}A{;!v 18~2Y ."͗/B LUw;QgKu-(~qD.9o:yڥ @1 f\~ ]id)b'Wj c:"E-R7Sk7RmX\5*'/3bMxv Q KwVtY QHJP49UѲꦟΓz9I55MM+*RɊmet)+^]zeb5 Ǣ\5QPMUA{?KtdZ8M1rQجhF~9!P0¹3Y &~|ݠK5Ԧx>6*@ z CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C Program PLANETF C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC PROGRAM XPLANET BYTE OPTION, ZOD1(12),CALEND, 1 RETRO(8),BLAH,GLYPH(21,14) LOGICAL HLP INTEGER ID(21,14),IM(21,14),IS(21,14),IH(4,3) INTEGER*4 JD1 DOUBLE PRECISION T,AN1(9),XJD,CONV,GMT,XL,YL,XP,YP,XN, 1 DAY,PI2,WEEK(7),WEEKR(7),TP, 2 YBES,TT DIMENSION EC0(9),EC1(9),AN0(9),P0(9),P1(9),EC(9), 1 TL(9),TB(9),TH0(9),TH1(9),XI(9),R(9),A(9),G(9),GE(9), 2 X(9),Y(9),Z(9),XX(9),YY(9),ZZ(9),H(4), 3 P(21,14), 8 MONTH(12),MOND(12),XNODE(2,10),YNODE(2,10) CCCCC8CCCC C CCCCCCCCC DATA MOND/31,0,31,30,31,30,31,31,30,31,30,31/ DATA MONTH/'JA','FE','MR','AP','MY','JN','JL','AU', 1 'SE','OC','NO','DE'/ DATA WEEK/'Monday ','Tuesday ','Wednsday', 1 'Thursday','Friday ','Saturday','Sunday '/ DATA WEEKR/'Moon) ','Mars) ','Mercury)','Jupiter)', 1 'Venus) ','Saturn) ','Sun) '/ DATA ZOD1/'A','T','G','C','L','V','=','S','/','K','Q','P'/ CCCCCCCCC C C PLUTO ELEMENTS FROM SHARAF (1964) C CCCCCCCCC DATA A/1.00000023,.387098599,.723331619,1.523688395, 1 5.202802875,9.53884320,19.19097811,30.0706724,39.672599/ DATA AN0/358.475833,102.279381,212.603222,319.529425, 1 225.444651,175.758444,74.313628,41.269550,231.002308/ DATA AN1/35999.049750D0,149472.515289D0,58517.803875D0, 1 19139.858500D0,3034.906654D0,1222.116782D0,428.502578D0, 2 218.466783D0,144.072477D0/ DATA EC0/.01675104,.20561421,.00682069,.09331290, 1 .04825382,.05606075,.04704433,.00853341,.24706226/ DATA EC1/-.00004180,.00002046,-.00004774,.00009206, 1 .0,.0,.0,.0,.0/ DATA P0/101.220833,75.899697,130.163833,334.218203, 1 11.907422,90.110981,169.048778,43.755611,221.592475/ DATA P1/1.719175,1.555489,1.408036,1.840758, 1 .0,.0,.0,.0,1.388888/ DATA TH0/0.,47.145944,75.779647,48.786442, 1 98.932822,112.347606,73.490250,130.678889,108.937165/ DATA TH1/0.,1.185208,.899850,.770992, 1 .0,.0,.509667,1.100972,1.358056/ DATA XI/0.,7.002881,3.393630,1.850333, 1 1.311614,2.494239,.7726658,1.779256,17.109816/ CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C Code C CCCCCCCCC ASIN(DUM)=ATAN(DUM/SQRT(-DUM*DUM+1)) TAN(DUMMY)=SIN(DUMMY)/COS(DUMMY) ARG(P,Q)=P-Q*INT(P/Q)-Q*INT(SIGN(.5,P-Q*INT(P/Q))-.5) HLP=.FALSE. LINEP=0 PI=ATAN(1.0) PI2=6.283185307179586476D0 CONV=PI2/360.D0 CONVS=CONV DO 130 I=1,21 DO 130 J=1,14 130 P(I,J)=0. WRITE(5,10000) 10000 FORMAT('$Sample input (Gregorian date,time', 1 'conversion to GMT,latitude,longitude,OPTION)'// 2 ' 1899/12/31 23:59 +08:00 34N59 122W59) : ') CCCCCCCCC C C HOROSCOPE LOOP STARTS HERE C CCCCCCCCC READ(5,11000)NY,CALEND,NM,ND,NH,NMI,NCH,NCM,LATD,LATS , 1 LATM,LONGD,LONGS,LONGM,OPTION 11000 FORMAT(I4,A1,4(I2,1X),I3,2(1X,I2),A1,I2,1X,I3,A1,I2,1X,A1) WRITE(5,11005) 11005 FORMAT('$Do you care for a hardcopy ? ') READ(5, 11006) BLAH 11006 FORMAT(A1) C IF (BLAH.NE.'Y'.AND.BLAH.NE.'y') GO TO 11035 C CALL OPEN(7,'B:PLANET.DAT') C IOUT=7 C GO TO 11008 11035 IOUT=5 11008 LINEP=1 WRITE(IOUT,11500)NY,CALEND,NM,ND,NH,NMI,NCH,NCM,LATD, 1 LATS,LATM,LONGD,LONGS,LONGM,OPTION 11500 FORMAT('1'//' ---Date--- -Time -Zone- -Lat- -Long- ', 1 'Comments--------'/' ',I4,A1,I2,'/',I2,I3,':',I2, 1 I4,':',I2,I3,A1,I2,I4,A1,I2,1X,A1) 144 GMT=((NH+NCH-12)*60+NMI+SIGN(NCM*1.,NCH*1.))/1440.D0 XJD=(ND-32075+1461.D0*(NY+4800+(NM-14)/12)/4+367.D0*(NM-2- 1 (NM-14)/12*12)/12-3.D0*((NY+4900+(NM-14)/12)/100)/4)+GMT IF(CALEND.NE.'B')GO TO 150 NY=1-NY 150 IF(CALEND.NE.'J'.AND.CALEND.NE.'B')GO TO 160 XJD=(ND-32075+1461.D0*(NY+4800+(NM-14)/12)/4 1 +367.D0*(NM-2-(NM-14)/12*12)/12-38)+GMT 160 T=(XJD-2415020.D0)/36525.D0 C JFINAL=2 CCCCCCCCC C C Planetary Ephemerides comxzted below (Newcomb/Hill,1898) C CCCCCCC 170 DO 270 J=1,JFINAL TT=T+(J-1)/36525.D0 TP=TT+18262.D0/36525.D0 W=TT+0.0 DO 174 I=1,9 EC(I)=EC0(I)+EC1(I)*W N=AN1(I)*TT/360.D0 174 G(I)=(AN1(I)*TT-N*360.D0+DBLE(AN0(I)))*CONV G5=G(5) G6=G(6) G7=(220.169542+428.49311*TP)*CONVS G(5)=G(5)+(.6506*SIN(2*G6-2*G5+336.9*CONVS) 1 +(3.9987-.002213*36525./4332.58*TP) 2 *SIN(5*G6-2*G5+(67.15-8197.0/3600.*TP)*CONVS) 3 +.5380*SIN(5*G6-3*G5+176.5*CONVS) 4 +.4112*SIN(2*G6-G5+1.4*CONVS) 4 +.0399278*36525./4332.58*TP*SIN(-G5+227.46*CONVS) 5 +.2763*SIN(3*G6-2*G5+127.4*CONVS) 6 +.2669*SIN(G6-G5+79.2*CONVS) )*299.12837/3600.*CONVS CCCCCCCCC C CCCCCCCCC G(6)=G(6)+(24.153*SIN(5*G6-2*G5+(247.11-2.277*TP)*CONVS) 1 +5.679*SIN(4*G6-2*G5+277.39*CONVS) 2 +3.505*SIN(2*G6-G5+181.43*CONVS) 3 +.657765*36525./10759.20*TP*SIN(G6+238.0*CONVS) 4 +.278*SIN(3*G6-G5+121.2*CONVS) 5 +.266*SIN(2*G6-2*G5+157.0*CONVS) 6 +.238*SIN(6*G6-2*G5-3*G7+6.9*CONVS) 7 +.223*SIN(10*G6-4*G5+(133.6-14814.5/3600.*W)*CONVS) 8 +.234*SIN(3*G7-G6+321.7*CONVS) )*120.455/3600.*CONVS CCCCCCCCC C C CCCCCCCCC G(7)=G(7)+33.086*W*W/3600.*CONVS G(8)=G(8)-22.401/3600.*W*W*CONVS DO 177 I=1,9 GE(I)=G(I)+EC(I)*SIN(G(I)+EC(I)*SIN(G(I)+EC(I)*SIN( 1 G(I)+EC(I)*SIN(G(I)+EC(I)*SIN(G(I)))))) 177 X(I)=2.*ATAN(SQRT((1.+EC(I))/(1.-EC(I)))*TAN(.5*GE(I))) 1 +(P0(I)+P1(I)*W)*CONVS + (I/7-I/9)*W*5025.3/3600.*CONVS GT=(358.415+35998.928*W)*CONVS GJ=(225.209+3034.462*W+.332*SIN((134.4+38.5*W)*CONVS)) 1 *CONVS X(4)=X(4)+(25.384*COS(GJ-G(4)-48.9*CONVS) 1 +52.490*SIN((47.48+19.771*W)*CONVS)-37.05-13.50*W 2 +21.869*COS(2*GJ-G(4)-188.3*CONVS) 3 +16.035*COS(2*GJ-2*G(4)-191.9*CONVS) 4 +13.966*COS(-GT+2*G(4)-20.5*CONVS) 5 +8.559*COS(-GT+G(4)-35.1*CONVS) )/3600.*CONVS CCCCCCCCC C C CCCCCCCCC GN=(225.417+3034.904*W)*CONVS GG=(175.753+1222.113*W)*CONVS G1=(74.412+428.498*W)*CONVS GP=(74.320+428.498*W)*CONVS GPP=(41.339+218.467*W)*CONVS X(7)=X(7)+(142.938*SIN(GG-2*G1)+19.508*COS(GG-2*G1) 1 +75.70*COS(3*G1-GG)-102.30*SIN(3*G1-GG) 2 -48.623*SIN(GN-G1)-21.320*COS(GN-G1) 3 -27.871*COS(GP-GPP)+19.869*SIN(GP-GPP) 4 +28.793*COS(2*GP-2*GPP)+10.035*SIN(2*GP-2*GPP) 5 +(18.37*SIN(3*GP-3*GPP)+8.91*COS(3*GP-3*GPP))*COS(G(7)) 6 +(8.35*SIN(3*GP-3*GPP)-16.44*COS(3*GP-3*GPP))*SIN(G(7)) 7 -18.585*COS(GG-G1)+12.603*SIN(GG-G1) 8 +4.327*COS(3*GP-3*GPP)+14.280*SIN(3*GP-3*GPP) 9 )/3600.*CONVS CCCCCCCCC C C CCCCCCCCC X(7)=X(7)+((112.317*W-1.551*W*W-.516*W*W*W)*SIN(G(7)) 1 -(68.339*W+9.721*W*W)*COS(G(7))+6.605*W*SIN(2*G(7)) 2 -(29.44-.410*W)*SIN((20.45-22.61*W)*CONVS) )/3600.*CONVS CCCCCCCCC C C JUP/SAT TERM & QUADRATURES OF LONG PERIOD UR`WNEP TERMS C CCCCCCCCC X(8)=X(8)+( 33.972*W*COS(G(8)) 1 +18.553*SIN((180.966+1004.034*W+.1403*W*W)*CONVS) 2 +34.138*SIN((153.267+2816.296*W-.0573*W*W)*CONVS) 3 +(26.50*W+3.92*W*W)*SIN(G(8)) )/3600.*CONVS CCCCCCCCC C C CCCCCCCCC DO 185 I=1,9 TH=(TH0(I)+TH1(I)*W)*CONVS DO 180 K=1,2 CCCCCCCCC C C MOTION OF JUP/SAT NODES FROM LEVERRIER C CCCCCCCCC TH2=TH+(I/5*3636.6-I/6*493.1-I/7*3143.5)/3600.*TP*CONVS R(I)=A(I)*(1.-EC(I)*COS(2.*ATAN(TAN((TH2+(K-1)*PI-(P0(I) 1 +P1(I)*W)*CONVS)*.5)*SQRT((1.-EC(I))/(1.+EC(I)))))) XNODE(K,I)=R(I)*COS(TH2+(K-1)*PI) 180 YNODE(K,I)=R(I)*SIN(TH2+(K-1)*PI) TL(I)=TH+ATAN2(SIN(X(I)-TH)*COS(XI(I)*CONVS), 1 COS(X(I)-TH))+(I/5-I/7)*TP*5026.1/3600.*CONVS TB(I)=ASIN(SIN(XI(I)*CONVS)*SIN(X(I)-TH)) 185 R(I)=A(I)*(1.-EC(I)*COS(GE(I))) TB(5)=TB(5)+4.37431*36525./4332.58*SIN(X(5)+23.62*CONVS) 1 *TP/3600.*CONVS TB(6)=TB(6)+24.266*36525./10759.2*SIN(X(6)-13.05*CONVS) 1 *TP/3600.*CONVS CCCCCCCCC C C TB(5) EFFECTS LESS THAN 35"T, TB(6) LESS THAN 84"T C CCCCCCCCC R(5)=R(5)*10.**(.0002303*COS(2*G6-2*G5+336.9*CONVS) 1 +.0001679*COS(5*G6-3*G5+176.4*CONVS) 2 +.0000125634*36525./4332.58*TP*COS(-G5+227.4*CONVS) ) CCCCCCCCC C C 22",16",10"T,(7",5",3") C CCCCCCCCC R(6)=R(6)*10.**(.0007005*COS(4*G6-2*G5+277.3*CONVS) 1 +.0003783*COS(G6-G5+79.8*CONVS) 2 +.000083491*36525./10759.20*TP*COS(G6+58.0*CONVS) 3 +.0002443*COS(2*G6-G5+176.0*CONVS) ) CCCCCCCCC C C CCCCCCCCC DO 190 I=1,9 X(I)=R(I)*COS(TB(I))*COS(TL(I)) Y(I)=R(I)*COS(TB(I))*SIN(TL(I)) 190 Z(I)=R(I)*SIN(TB(I)) DO 210 I=2,9 DO 200 K=1,2 XNODE(K,I)=XNODE(K,I)-X(1) YNODE(K,I)=YNODE(K,I)-Y(1) 200 P(K+17,I)=ARG(ATAN2(YNODE(K,I),XNODE(K,I))/6.2831853,1.) XX(I)=X(I)-X(1) YY(I)=Y(I)-Y(1) 210 ZZ(I)=Z(I)-Z(1) XX(1)=-X(1) YY(1)=-Y(1) ZZ(1)=-Z(1) XNUT=-17.2327/1296000.*SIN((259.18-1934.142*W)*CONVS) DO 230 I=1,9 P(J+19,I)=ARG(ATAN2(YY(I),XX(I))/6.28318531+XNUT,1.) IF(J-1)230,220,230 220 DIST=SQRT(XX(I)*XX(I)+YY(I)*YY(I)+ZZ(I)*ZZ(I)) P(1,I)=ASIN(ZZ(I)/DIST)/CONVS E=(23.4522944-.0130125*W+.002558*COS((259.18-1934.142 1 *W)*CONVS))*CONVS P(2,I)=ASIN(COS(P(1,I)*CONVS)*SIN(P(20,I)*6.28318531) 1 *SIN(E)+SIN(P(1,I)*CONVS)*COS(E))/CONVS P(6,I)=ARG(TL(I)/6.28318531,1.) P(3,I)=TB(I)/CONVS P(4,I)=ASIN(COS(TB(I))*SIN(TL(I))*SIN(E) 1 +SIN(TB(I))*COS(E))/CONVS 230 CONTINUE CCCCCCCCC C C LUNAR EPHEMERIDES COMPUTED TO WITHIN 1' (BROWN/I.L.E.) C CCCCCCCCC DAY=XJD-2415020.D0+(J-1.D0) XL=.751206D0+DAY*.0366011014634D0 YL=.776935D0+DAY*.0027379092649D0 XP=.928693D0+DAY*.0003094557786D0 YP=.781169D0+DAY*.0000001307457D0 XN=.719954D0-DAY*.0001470942283D0 AL=(XL-IDINT(XL))*PI2 BL=(YL-IDINT(YL))*PI2 AP=(XP-IDINT(XP))*PI2 BP=(YP-IDINT(YP))*PI2 AN=(XN-IDINT(XN))*PI2 U=AL-AP V=BL-BP F=AL-AN D=AL-BL DL=22639*SIN(U)-4586*SIN(U-D-D)+2370*SIN(D+D)+769*SIN(U+U) 1 -668*SIN(V)-412*SIN(F+F)-212*SIN(U+U-D-D)-206*SIN(U+V-D-D) 2 +192*SIN(U+D+D)-165*SIN(V-D-D)+148*SIN(U-V)-125*SIN(D) 3 -110*SIN(U+V)-55*SIN(F+F-D-D)-45*SIN(U+F+F)+40*SIN(U-F-F) 4 -38*SIN(U-4*D)+36*SIN(3*U)-31*SIN(U+U-4*D)+28*SIN(U-V-D-D) 5 -24*SIN(V+D+D)+19*SIN(U-D)+18*SIN(V+D)+15*SIN(U-V+D+D) 6 +14*SIN(4*D)+14*SIN(U+U+D+D)-13*SIN(3*U-D-D) P(J+19,10)=ARG(XL-IDINT(XL)+DL/1296000.+XNUT,1.) P(J+19,11)=ARG(AN/6.28318531,1.) IF(J-1)270,260,270 260 P(1,10)=(18461*SIN(F)+1010*SIN(U+F)-1000*SIN(F-U) 1 -624*SIN(F-2*D)-167*SIN(U+F-D-D)+199*SIN(F-U+D+D) 2 +117*SIN(F+D+D)+62*SIN(U+U+F)-33*SIN(F-U-D-D) 3 -32*SIN(F-U-U)-30*SIN(V+F-D-D)-16*SIN(U+U+F-D-D) 4 +15*SIN(U+F+D+D)+12*SIN(F-V-D-D)+9*SIN(F-U-V+D+D))/3600. P(2,10)=ASIN(COS(P(1,10)*CONVS)*SIN(P(20,10)*6.28318531) 1 *SIN(E)+SIN(P(1,10)*CONVS)*COS(E))/CONVS 270 CONTINUE WRITE(5,271) 271 FORMAT(' ','EPH.READY'/) CCCCCCCCC C CCCCCCCCC E=(23.4522944-.0130125*T+.002558*DCOS((259.18-1934.142 1 *T)*CONV))*CONVS DO 277 I=2,9 RETRO(I-1)='D' IF(ARG(P(21,I)-P(20,I),1.)-.5)277,277,275 275 RETRO(I-1)='R' 277 CONTINUE CCCCCCCCC C CCCCCCCCC XLW=(LONGD+LONGM/60.)/360. XLN=(LATD+LATM/60.)*CONVS IF(LATS.EQ.'N')GO TO 280 XLN=-XLN 280 IF(LONGS.EQ.'W')GO TO 290 XLW=-XLW CCCCCCCCC C C H(1)=GMT H(2)=GST H(3)=LMT H(4)=LST C CCCCCCCCC 290 H(1)=ARG(SNGL(GMT)+.5,1.) W=100.00213590D0*T-IDINT(100.D0*T) H(2)=ARG(H(1)+.27691940+W,1.) H(3)=ARG(H(1)-XLW,1.) H(4)=ARG(H(2)-XLW,1.) DO 300 I=1,4 W=24.*H(I)+.00013888888 IH(I,1)=W+0 W=60.*(W-IH(I,1)) IH(I,2)=W+0 300 IH(I,3)=60.*(W-IH(I,2)) WRITE(IOUT,12000),((IH(I,J),J=1,3),I=1,4) 12000 FORMAT(/' GMT=',I2,2(':',I2),' GST=',I2,2(':',I2), 1 ' LMT=',I2,2(':',I2),' LST=',I2,2(':',I2)) JD1=IDINT(XJD) XJD1=XJD-JD1 I=JD1+INT(XJD1+1.5-XLW)-(JD1-1+INT(XJD1+1.5-XLW))/7*7 YBES=1900.D0+(XJD-2415020.31351528D0)/365.24219878125D0 1 -.00000107523D0*T*T C IYBES=YBES+0 C YBES1=YBES-IYBES WRITE(IOUT,13000)WEEK(I),WEEKR(I),XJD,YBES 13000 FORMAT(' Day of the week (from LMT) IS ',A9, 1 ' (ruled by ',A8/' JD=',D14.7,7X,' Besselian year = ', 2 D14.7) CCCCCCCCC C C HOUSE CUSP COMPUTATIONS (ACCURATE TO LAST DIGIT) C CCCCCCCCC 310 ITER=1+LATD/45 CCCCCCCCC C C "ITER" IS THE NUMBER OF PLACIDEAN ITERATIONS C@@@ C@@@ C CCCCCCCCC RA=6.28318531*H(4) DO 340 N=1,6 AA=6.28318531*(N+2)/12. P(7,N)=ARG(ATAN2(SIN(RA)*COS(AA)+COS(RA)*COS(XLN)* 1 SIN(AA),COS(E)*(COS(RA)*COS(AA)-SIN(RA)*COS(XLN)* 2 SIN(AA))-SIN(E)*SIN(XLN)*SIN(AA))/6.28318531,1.) P(8,N)=ARG(ATAN2(SIN(RA+AA)*COS(E),COS(RA+AA)) 1 /6.28318531,1.) P(9,N)=ARG(ATAN2(SIN(RA+AA),COS(E)*COS(RA+AA)-SIN((3.- 1 ABS(1.-N))/3.*ASIN(TAN(E)*TAN(XLN)))*COS(E))/PI*.5,1.) DO 330 J=1,ITER X(1)=P(9,N)*6.28318531 X(2)=X(1)+CONVS*.5 CCCCCCCCC C C INCREMENT "CONVS*.5" IS ARBITRARILY CHOSEN (]W5 DEG) C CCCCCCCCC DO 320 K=1,2 AD=ASIN(TAN(XLN)*TAN(ASIN(SIN(X(K))*SIN(E)))) DRA=ARG(ATAN2(SIN(X(K))*COS(E),COS(X(K)))-RA-(N-1)*PI/6. 1 ,PI)-PI*.5 320 Y(K)=DRA-(4.-N)*AD/3. 330 P(9,N)=P(9,N)-Y(1)/(Y(2)-Y(1))*.5/360. CCCCCCCCC C C THE ".5" CONSTANT HERE MUST AGREE WITH THE ONE ABOVE C CCCCCCCCC P(10,N)=ARG(ATAN2(SIN(RA+AA),COS(E)*COS(RA+AA)-SIN(E) 1 *TAN(XLN)*SIN(AA))/6.28318531,1.) 340 P(12,N)=ARG(P(7,1)+(N-1)/12.,1.) DO 350 N=1,6 P(11,N)=ARG(P(7,1)+(N-1)/12.+(ARG(P(7,4)-P(7,1),1.) 1 -.25)*(3.-ABS(3.-ABS(N-7.)))/3.,1.) 350 P(13,N)=ARG(P(7,4)+(N-4)/12.,1.) W=ARG(ATAN2(SQRT((COS(XLN)/SIN(E))**2-SIN(P(7,4) 1 *6.28318531)**2),SIN(XLN)*SIN(P(7,4)*6.28318531)),PI)/3. DO 351 N=1,3 351 P(14,N)=ARG(ATAN2(COS(RA+(N-1)*W),-SIN(RA+(N-1)*W) 1 *COS(E)-SIN(E)*TAN(XLN))/6.28318531,1.) W=ARG(ATAN2(SIN(P(7,1)*6.28318531)*COS(E),COS(P(7,1) 1 *6.28318531))-RA,PI)/3. DO 352 N=4,6 352 P(14,N)=ARG(ATAN2(SIN(RA+(N-4)*W),COS(RA+(N-4)*W)*COS(E) 1 )/6.28318531+.5,1.) DO 355 J=7,14 DO 355 N=1,6 355 P(J,N+6)=ARG(P(J,N)+.5,1.) CCCCCCCCC C C PARS FORTUNA & PLANETARY NODES FOLLOW C CCCCCCCCC P(20,12)=ARG(P(20,10)-P(20,1)+P(11,1),1.) DO 360 I=2,9 360 P(17,I)=ARG(P(20,I)-P(20,1)+P(11,1),1.) CCCCCCCCC C C ARGUMENTS TRANSFORMED TO SIGNS/DEGREES/MINUTES C CCCCCCCCC DO 390 I=1,4 DO 390 J=1,12 GLYPH(I,J)='N' IF(P(I,J))370,380,380 370 GLYPH(I,J)='S' 380 DEG=ABS(P(I,J))+.0083333333 ID(I,J)=DEG+0 390 IM(I,J)=60*(DEG-ID(I,J)) IFINAL=JFINAL+19 DO 400 I=6,IFINAL DO 400 J=1,12 W=12.*P(I,J)+.00027777777 IS(I,J)=W+0 ISS=W+1 GLYPH(I,J)=ZOD1(ISS) W=30.*(W-IS(I,J)) ID(I,J)=W+0 400 IM(I,J)=60.*(W-ID(I,J)) CCCCCCCCC C C PROGRESSION & EPHEMERIS OUTPUT ROUTINE FOLLOWS C CCCCCCCCC IF(OPTION.NE.'P')GO TO 420 WRITE(IOUT,14000) 14000 FORMAT(/19X,'Progressed (geocentric) longitudes'// 1 ' YEAR SUN MERC VENUS MARS JUPIT SATUR ', 2 ' URAN NEPT PLUTO MOON NODE') DO 410 I=1,21 II=I 410 WRITE(IOUT,15000)II,(ID(I,J),GLYPH(I,J),IM(I,J),J=1,11) 15000 FORMAT(' ',I6,1X,11(1X,I6,A1,I6)) GO TO 9999 420 IF(OPTION.NE.'E')GO TO 460 WRITE(IOUT,16000)IH(1,1),IH(1,2),NY 16000 FORMAT(//19X,'Daily ephemeris for ',I2,':',I2,' GMT'// 1 ' DATE SUN MERC VENUS MARS JUPIT SATUR URAN', 2 ' NEPT PLUTO MOON NODE'/I6) MY=NY+(NM+7)/10 MOND(2)=1461*(MY+4800)/4-1461*(MY+4799)/4-337 IF(CALEND.EQ.'J')GO TO 430 MOND(2)=MOND(2)-(MY+4900)/100*3/4+(MY+4899)/100*3/4 430 DO 460 I=1,21 WRITE(IOUT,17000)MONTH(NM),ND,(ID(I,J),GLYPH(I,J), 1 IM(I,J),J=1,11) 17000 FORMAT(' ',A2,I2,11(1X,I2,A1,I2)) ND=ND+1 IF(ND-MOND(NM))460,460,440 440 ND=ND-MOND(NM) NM=NM+1 IF(NM-12)460,460,450 450 NM=1 NY=NY+1 WRITE(IOUT,18000)NY 18000 FORMAT(I6) 460 CONTINUE GO TO 9999 CCCCCCCCC 9999 IF(LINEP)9997,9997,9996 9996 WRITE(5,96000) 96000 FORMAT(' $DONE') 9997 END  CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C Program PLANETF C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC PROGRAM XPLANET BYTE COMENT(25),FILNAM(25),OPTION, 1 RETRO(8),BLAH LOGICAL HLP DOUBLE PRECISION T,AN1(9),XJD,CONV,GMT,XL,YL,XP,YP,XN, 1 DAY,PI2,RULER(12),WEEK(7),WEEKR(7), 2 YBES,TT DIMENSION EC0(9),EC1(9),AN0(9),P0(9),P1(9),EC(9), 1 TL(9),TB(9),TH0(9),TH1(9),XI(9),R(9),A(9),G(9),GE(9), 2 X(9),Y(9),Z(9),XX(9),YY(9),ZZ(9),IH(4,3),H(4), 3 ZOD1(12),P(120,14),ID(120,14),IM(120,14),IS(120,14), 4 GLYPH(120,14),IASP(14,14) 8 ,MONTH(12),MOND(12),XNODE(2,10),YNODE(2,10) CCCCC8CCCC C CCCCCCCCC DATA MOND/31,0,31,30,31,30,31,31,30,31,30,31/ DATA MONTH/'JA','FE','MR','AP','MY','JN','JL','AU', 1 'SE','OC','NO','DE'/ DATA WEEK/'Monday ','Tuesday ','Wednsday', 1 'Thursday','Friday ','Saturday','Sunday '/ DATA WEEKR/'Moon) ','Mars) ','Mercury)','Jupiter)', 1 'Venus) ','Saturn) ','Sun) '/ DATA ZOD/'Aries ','Taurus ','Gemini ','Cancer ', 1 'Leo ','Virgo ','Libra ','Scorpio ', 2 'Sagtarus','Caprcorn','Aquarius','Pisces '/ DATA ZOD1/'A','T','G','C','L','V','=','S','/','K','Q','P'/ CCCCCCCCC C C PLUTO ELEMENTS FROM SHARAF (1964) C CCCCCCCCC DATA A/1.00000023,.387098599,.723331619,1.523688395, 1 5.202802875,9.53884320,19.19097811,30.0706724,39.672599/ DATA AN0/358.475833,102.279381,212.603222,319.529425, 1 225.444651,175.758444,74.313628,41.269550,231.002308/ DATA AN1/35999.049750D0,149472.515289D0,58517.803875D0, 1 19139.858500D0,3034.906654D0,1222.116782D0,428.502578D0, 2 218.466783D0,144.072477D0/ DATA EC0/.01675104,.20561421,.00682069,.09331290, 1 .04825382,.05606075,.04704433,.00853341,.24706226/ DATA EC1/-.00004180,.00002046,-.00004774,.00009206, 1 .0,.0,.0,.0,.0/ DATA P0/101.220833,75.899697,130.163833,334.218203, 1 11.907422,90.110981,169.048778,43.755611,221.592475/ DATA P1/1.719175,1.555489,1.408036,1.840758, 1 .0,.0,.0,.0,1.388888/ DATA TH0/0.,47.145944,75.779647,48.786442, 1 98.932822,112.347606,73.490250,130.678889,108.937165/ DATA TH1/0.,1.185208,.899850,.770992, 1 .0,.0,.509667,1.100972,1.358056/ DATA XI/0.,7.002881,3.393630,1.850333, 1 1.311614,2.494239,.7726658,1.779256,17.109816/ CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C Code C CCCCCCCCC TAN(DUMMY)=SIN(DUMMY)/COS(DUMMY) ARG(P,Q)=P-Q*INT(P/Q)-Q*INT(SIGN(.5,P-Q*INT(P/Q))-.5) HLP=.FALSE. LINEP=0 PI=ATAN(1.0) PI2=6.283185307179586476D0 CONV=PI2/360.D0 CONVS=CONV DO 130 I=1,120 DO 130 J=1,14 130 P(I,J)=0. WRITE(5,10000) 10000 FORMAT('$Sample input (Gregorian date,time', 1 'conversion to GMT,latitude,longitude,OPTION)'// 2 ' 1899/12/31 23:59 +08:00 34N59 122W59) : ') CCCCCCCCC C C HOROSCOPE LOOP STARTS HERE C CCCCCCCCC READ(5,11000)NY,CALEND,NM,ND,NH,NMI,NCH,NCM,LATD,LATS , 1 LATM,LONGD,LONGS,LONGM,OPTION 11000 FORMAT(I4,A1,4(I2,1X),I3,2(1X,I2),A1,I2,1X,I3,A1,I2,1X,1A) 6011 CONTINUE WRITE(6,11005) 11005 FORMAT('$Do you care for a hardcopy ? ') READ(5, 11006) BLAH 11006 FORMAT(A1) IF (BLAH.NE.'Y'.AND.BLAH.NE.'y') GO TO 11035 CONTINUE CALL OPEN(6,'PLANET.DAT') IOUT=6 GO TO 11008 11035 IOUT=5 11008 LINEP=1 WRITE(IOUT,11500)NY,CALEND,NM,ND,NH,NMI,NCH,NCM,LATD, 1 LATS,LATM,LONGD,LONGS,LONGM 11500 FORMAT('1'//' ---Date--- -Time -Zone- -Lat- -Long- ', 1 'Comments--------'/' ',I4,A1,I2,'/',I2,I3,':',I2, 1 I4,':',I2,I3,A1,I2,I4,A1,I2) 144 GMT=((NH+NCH-12)*60+NMI+SIGN(NCM*1.,NCH*1.))/1440. XJD=(ND-32075+1461*(NY+4800+(NM-14)/12)/4+367*(NM-2- 1 (NM-14)/12*12)/12-3*((NY+4900+(NM-14)/12)/100)/4)+GMT IF(CALEND.NE.'B')GO TO 150 NY=1-NY 150 IF(CALEND.NE.'J'.AND.CALEND.NE.'B')GO TO 160 XJD=(ND-32075+1461*(NY+4800+(NM-14)/12)/4 1 +367*(NM-2-(NM-14)/12*12)/12-38)+GMT 160 T=(XJD-2415020.D0)/36525.D0 C JFINAL=101 CCCCCCCCC C C Planetary Ephemerides comxzted below (Newcomb/Hill,1898) C CCCCCCC 170 DO 270 J=1,JFINAL TT=T+(J-1)/36525.D0 TP=TT+18262.D0/36525.D0 W=TT DO 174 I=1,9 EC(I)=EC0(I)+EC1(I)*W N=AN1(I)*TT/360.D0 174 G(I)=(AN1(I)*TT-N*360.D0+DBLE(AN0(I)))*CONV G5=G(5) G6=G(6) G7=(220.169542+428.49311*TP)*CONVS G(5)=G(5)+(.6506*SIN(2*G6-2*G5+336.9*CONVS) 1 +(3.9987-.002213*36525./4332.58*TP) 2 *SIN(5*G6-2*G5+(67.15-8197.0/3600.*TP)*CONVS) 3 +.5380*SIN(5*G6-3*G5+176.5*CONVS) 4 +.4112*SIN(2*G6-G5+1.4*CONVS) 4 +.0399278*36525./4332.58*TP*SIN(-G5+227.46*CONVS) 5 +.2763*SIN(3*G6-2*G5+127.4*CONVS) 6 +.2669*SIN(G6-G5+79.2*CONVS) )*299.12837/3600.*CONVS CCCCCCCCC C CCCCCCCCC G(6)=G(6)+(24.153*SIN(5*G6-2*G5+(247.11-2.277*TP)*CONVS) 1 +5.679*SIN(4*G6-2*G5+277.39*CONVS) 2 +3.505*SIN(2*G6-G5+181.43*CONVS) 3 +.657765*36525./10759.20*TP*SIN(G6+238.0*CONVS) 4 +.278*SIN(3*G6-G5+121.2*CONVS) 5 +.266*SIN(2*G6-2*G5+157.0*CONVS) 6 +.238*SIN(6*G6-2*G5-3*G7+6.9*CONVS) 7 +.223*SIN(10*G6-4*G5+(133.6-14814.5/3600.*W)*CONVS) 8 +.234*SIN(3*G7-G6+321.7*CONVS) )*120.455/3600.*CONVS CCCCCCCCC C C CCCCCCCCC G(7)=G(7)+33.086*W*W/3600.*CONVS G(8)=G(8)-22.401/3600.*W*W*CONVS DO 177 I=1,9 GE(I)=G(I)+EC(I)*SIN(G(I)+EC(I)*SIN(G(I)+EC(I)*SIN( 1 G(I)+EC(I)*SIN(G(I)+EC(I)*SIN(G(I)))))) 177 X(I)=2.*ATAN(SQRT((1.+EC(I))/(1.-EC(I)))*TAN(.5*GE(I))) 1 +(P0(I)+P1(I)*W)*CONVS + (I/7-I/9)*W*5025.3/3600.*CONVS GT=(358.415+35998.928*W)*CONVS GJ=(225.209+3034.462*W+.332*SIN((134.4+38.5*W)*CONVS)) 1 *CONVS X(4)=X(4)+(25.384*COS(GJ-G(4)-48.9*CONVS) 1 +52.490*SIN((47.48+19.771*W)*CONVS)-37.05-13.50*W 2 +21.869*COS(2*GJ-G(4)-188.3*CONVS) 3 +16.035*COS(2*GJ-2*G(4)-191.9*CONVS) 4 +13.966*COS(-GT+2*G(4)-20.5*CONVS) 5 +8.559*COS(-GT+G(4)-35.1*CONVS) )/3600.*CONVS CCCCCCCCC C C CCCCCCCCC GN=(225.417+3034.904*W)*CONVS GG=(175.753+1222.113*W)*CONVS G1=(74.412+428.498*W)*CONVS GP=(74.320+428.498*W)*CONVS GPP=(41.339+218.467*W)*CONVS X(7)=X(7)+(142.938*SIN(GG-2*G1)+19.508*COS(GG-2*G1) 1 +75.70*COS(3*G1-GG)-102.30*SIN(3*G1-GG) 2 -48.623*SIN(GN-G1)-21.320*COS(GN-G1) 3 -27.871*COS(GP-GPP)+19.869*SIN(GP-GPP) 4 +28.793*COS(2*GP-2*GPP)+10.035*SIN(2*GP-2*GPP) 5 +(18.37*SIN(3*GP-3*GPP)+8.91*COS(3*GP-3*GPP))*COS(G(7)) 6 +(8.35*SIN(3*GP-3*GPP)-16.44*COS(3*GP-3*GPP))*SIN(G(7)) 7 -18.585*COS(GG-G1)+12.603*SIN(GG-G1) 8 +4.327*COS(3*GP-3*GPP)+14.280*SIN(3*GP-3*GPP) 9 )/3600.*CONVS CCCCCCCCC C C CCCCCCCCC X(7)=X(7)+((112.317*W-1.551*W*W-.516*W*W*W)*SIN(G(7)) 1 -(68.339*W+9.721*W*W)*COS(G(7))+6.605*W*SIN(2*G(7)) 2 -(29.44-.410*W)*SIN((20.45-22.61*W)*CONVS) )/3600.*CONVS CCCCCCCCC C C JUP/SAT TERM & QUADRATURES OF LONG PERIOD UR`WNEP TERMS C CCCCCCCCC X(8)=X(8)+( 33.972*W*COS(G(8)) 1 +18.553*SIN((180.966+1004.034*W+.1403*W*W)*CONVS) 2 +34.138*SIN((153.267+2816.296*W-.0573*W*W)*CONVS) 3 +(26.50*W+3.92*W*W)*SIN(G(8)) )/3600.*CONVS CCCCCCCCC C C CCCCCCCCC DO 185 I=1,9 TH=(TH0(I)+TH1(I)*W)*CONVS DO 180 K=1,2 CCCCCCCCC C C MOTION OF JUP/SAT NODES FROM LEVERRIER C CCCCCCCCC TH2=TH+(I/5*3636.6-I/6*493.1-I/7*3143.5)/3600.*TP*CONVS R(I)=A(I)*(1.-EC(I)*COS(2.*ATAN(TAN((TH2+(K-1)*PI-(P0(I) 1 +P1(I)*W)*CONVS)*.5)*SQRT((1.-EC(I))/(1.+EC(I)))))) XNODE(K,I)=R(I)*COS(TH2+(K-1)*PI) 180 YNODE(K,I)=R(I)*SIN(TH2+(K-1)*PI) TL(I)=TH+ATAN2(SIN(X(I)-TH)*COS(XI(I)*CONVS), 1 COS(X(I)-TH))+(I/5-I/7)*TP*5026.1/3600.*CONVS TB(I)=ASIN(SIN(XI(I)*CONVS)*SIN(X(I)-TH)) 185 R(I)=A(I)*(1.-EC(I)*COS(GE(I))) TB(5)=TB(5)+4.37431*36525./4332.58*SIN(X(5)+23.62*CONVS) 1 *TP/3600.*CONVS TB(6)=TB(6)+24.266*36525./10759.2*SIN(X(6)-13.05*CONVS) 1 *TP/3600.*CONVS CCCCCCCCC C C TB(5) EFFECTS LESS THAN 35"T, TB(6) LESS THAN 84"T C CCCCCCCCC R(5)=R(5)*10.**(.0002303*COS(2*G6-2*G5+336.9*CONVS) 1 +.0001679*COS(5*G6-3*G5+176.4*CONVS) 2 +.0000125634*36525./4332.58*TP*COS(-G5+227.4*CONVS) ) CCCCCCCCC C C 22",16",10"T,(7",5",3") C CCCCCCCCC R(6)=R(6)*10.**(.0007005*COS(4*G6-2*G5+277.3*CONVS) 1 +.0003783*COS(G6-G5+79.8*CONVS) 2 +.000083491*36525./10759.20*TP*COS(G6+58.0*CONVS) 3 +.0002443*COS(2*G6-G5+176.0*CONVS) ) CCCCCCCCC C C CCCCCCCCC DO 190 I=1,9 X(I)=R(I)*COS(TB(I))*COS(TL(I)) Y(I)=R(I)*COS(TB(I))*SIN(TL(I)) 190 Z(I)=R(I)*SIN(TB(I)) DO 210 I=2,9 DO 200 K=1,2 XNODE(K,I)=XNODE(K,I)-X(1) YNODE(K,I)=YNODE(K,I)-Y(1) 200 P(K+17,I)=ARG(ATAN2(YNODE(K,I),XNODE(K,I))/6.2831853,1.) XX(I)=X(I)-X(1) YY(I)=Y(I)-Y(1) 210 ZZ(I)=Z(I)-Z(1) XX(1)=-X(1) YY(1)=-Y(1) ZZ(1)=-Z(1) XNUT=-17.2327/1296000.*SIN((259.18-1934.142*W)*CONVS) DO 230 I=1,9 P(J+19,I)=ARG(ATAN2(YY(I),XX(I))/6.28318531+XNUT,1.) IF(J-1)230,220,230 220 DIST=SQRT(XX(I)*XX(I)+YY(I)*YY(I)+ZZ(I)*ZZ(I)) P(1,I)=ASIN(ZZ(I)/DIST)/CONVS E=(23.4522944-.0130125*W+.002558*COS((259.18-1934.142 1 *W)*CONVS))*CONVS P(2,I)=ASIN(COS(P(1,I)*CONVS)*SIN(P(20,I)*6.28318531) 1 *SIN(E)+SIN(P(1,I)*CONVS)*COS(E))/CONVS P(6,I)=ARG(TL(I)/6.28318531,1.) P(3,I)=TB(I)/CONVS P(4,I)=ASIN(COS(TB(I))*SIN(TL(I))*SIN(E) 1 +SIN(TB(I))*COS(E))/CONVS 230 CONTINUE CCCCCCCCC C C LUNAR EPHEMERIDES COMPUTED TO WITHIN 1' (BROWN/I.L.E.) C CCCCCCCCC DAY=XJD-2415020.D0+(J-1.D0) XL=.751206D0+DAY*.0366011014634D0 YL=.776935D0+DAY*.0027379092649D0 XP=.928693D0+DAY*.0003094557786D0 YP=.781169D0+DAY*.0000001307457D0 XN=.719954D0-DAY*.0001470942283D0 AL=(XL-IDINT(XL))*PI2 BL=(YL-IDINT(YL))*PI2 AP=(XP-IDINT(XP))*PI2 BP=(YP-IDINT(YP))*PI2 AN=(XN-IDINT(XN))*PI2 U=AL-AP V=BL-BP F=AL-AN D=AL-BL DL=22639*SIN(U)-4586*SIN(U-D-D)+2370*SIN(D+D)+769*SIN(U+U) 1 -668*SIN(V)-412*SIN(F+F)-212*SIN(U+U-D-D)-206*SIN(U+V-D-D) 2 +192*SIN(U+D+D)-165*SIN(V-D-D)+148*SIN(U-V)-125*SIN(D) 3 -110*SIN(U+V)-55*SIN(F+F-D-D)-45*SIN(U+F+F)+40*SIN(U-F-F) 4 -38*SIN(U-4*D)+36*SIN(3*U)-31*SIN(U+U-4*D)+28*SIN(U-V-D-D) 5 -24*SIN(V+D+D)+19*Sig(U-D)+18*SIN(V+D)+15*SIN(U-V+D+D) 6 +14*SIN(4*D)+14*SIN(U+U+D+D)-13*SIN(3*U-D-D) P(J+19,10)=ARG(XL-IDINT(XL)+DL/1296000.+XNUT,1.) P(J+19,11)=ARG(AN/6.28318531,1.) IF(J-1)270,260,270 260 P(1,10)=(18461*SIN(F)+1010*SIN(U+F)-1000*SIN(F-U) 1 -624*SIN(F-2*D)-167*SIN(U+F-D-D)+199*SIN(F-U+D+D) 2 +117*SIN(F+D+D)+62*SIN(U+U+F)-33*SIN(F-U-D-D) 3 -32*SIN(F-U-U)-30*SIN(V+F-D-D)-16*SIN(U+U+F-D-D) 4 +15*SIN(U+F+D+D)+12*SIN(F-V-D-D)+9*SIN(F-U-V+D+D))/3600. P(2,10)=ASIN(COS(P(1,10)*CONVS)*SIN(P(20,10)*6.28318531) 1 *SIN(E)+SIN(P(1,10)*CONVS)*COS(E))/CONVS 270 CONTINUE CCCCCCCCC C CCCCCCCCC E=(23.4522944-.0130125*T+.002558*DCOS((259.18-1934.142 1 *T)*CONV))*CONVS DO 277 I=2,9 RETRO(I-1)=' DIR ' IF(ARG(P(21,I)-P(20,I),1.)-.5)277,277,275 275 RETRO(I-1)='RETRO' 277 CONTINUE CCCCCCCCC C CCCCCCCCC XLW=(LONGD+LONGM/60.)/360. XLN=(LATD+LATM/60.)*CONVS IF(LATS.EQ.'N')GO TO 280 XLN=-XLN 280 IF(LONGS.EQ.'W')GO TO 290 XLW=-XLW CCCCCCCCC C C H(1)=GMT H(2)=GST H(3)=LMT H(4)=LST C CCCCCCCCC 290 H(1)=ARG(SNGL(GMT)+.5,1.) W=100.00213590D0*T-IDINT(100.D0*T) H(2)=ARG(H(1)+.27691940+W,1.) H(3)=ARG(H(1)-XLW,1.) H(4)=ARG(H(2)-XLW,1.) DO 300 I=1,4 W=24.*H(I)+.00013888888 IH(I,1)=W W=60.*(W-IH(I,1)) IH(I,2)=W 300 IH(I,3)=60.*(W-IH(I,2)) WRITE(IOUT,12000),((IH(I,J),J=1,3),I=1,4) 12000 FORMAT(/' GMT=',I2,2(':',I2),' GST=',I2,2(':',I2), 1 ' LMT=',I2,2(':',I2),' LST=',I2,2(':',I2)) JD1=XJD XJD1=XJD-JD1 I=JD1+INT(XJD1+1.5-XLW)-(JD1-1+INT(XJD1+1.5-XLW))/7*7 YBES=1900.D0+(XJD-2415020.31351528D0)/365.24219878125D0 1 -.00000107523D0*T*T IYBES=YBES YBES1=YBES-IYBES WRITE(IOUT,13000)WEEK(I),WEEKR(I),JD1,XJD1,IYBES,YBES1 13000 FORMAT(' Day of the week (from LMT) IS ',A9, 1 ' (ruled by ',A8/' JD=',I7,F7.6,7X,' Besselian year = ', 2 I5,F10.9) CCCCCCCCC C C HOUSE CUSP COMPUTATIONS (ACCURATE TO LAST DIGIT) C CCCCCCCCC 310 ITER=1+LATD/45 CCCCCCCCC C C "ITER" IS THE NUMBER OF PLACIDEAN ITERATIONS C@@@ C@@@ C CCCCCCCCC RA=6.28318531*H(4) DO 340 N=1,6 AA=6.28318531*(N+2)/12. P(7,N)=ARG(ATAN2(SIN(RA)*COS(AA)+COS(RA)*COS(XLN)* 1 SIN(AA),COS(E)*(COS(RA)*COS(AA)-SIN(RA)*COS(XLN)* 2 SIN(AA))-SIN(E)*SIN(XLN)*SIN(AA))/6.28318531,1.) P(8,N)=ARG(ATAN2(SIN(RA+AA)*COS(E),COS(RA+AA)) 1 /6.28318531,1.) P(9,N)=ARG(ATAN2(SIN(RA+AA),COS(E)*COS(RA+AA)-SIN((3.- 1 ABS(1.-N))/3.*ASIN(TAN(E)*TAN(XLN)))*COS(E))/PI*.5,1.) DO 330 J=1,ITER X(1)=P(9,N)*6.28318531 X(2)=X(1)+CONVS*.5 CCCCCCCCC C C INCREMENT "CONVS*.5" IS ARBITRARILY CHOSEN (]W5 DEG) C CCCCCCCCC DO 320 K=1,2 AD=ASIN(TAN(XLN)*TAN(ASIN(SIN(X(K))*SIN(E)))) DRA=ARG(ATAN2(SIN(X(K))*COS(E),COS(X(K)))-RA-(N-1)*PI/6. 1 ,PI)-PI*.5 320 Y(K)=DRA-(4.-N)*AD/3. 330 P(9,N)=P(9,N)-Y(1)/(Y(2)-Y(1))*.5/360. CCCCCCCCC C C THE ".5" CONSTANT HERE MUST AGREE WITH THE ONE ABOVE C CCCCCCCCC P(10,N)=ARG(ATAN2(SIN(RA+AA),COS(E)*COS(RA+AA)-SIN(E) 1 *TAN(XLN)*SIN(AA))/6.28318531,1.) 340 P(12,N)=ARG(P(7,1)+(N-1)/12.,1.) DO 350 N=1,6 P(11,N)=ARG(P(7,1)+(N-1)/12.+(ARG(P(7,4)-P(7,1),1.) 1 -.25)*(3.-ABS(3.-ABS(N-7.)))/3.,1.) 350 P(13,N)=ARG(P(7,4)+(N-4)/12.,1.) W=ARG(ATAN2(SQRT((COS(XLN)/SIN(E))**2-SIN(P(7,4) 1 *6.28318531)**2),SIN(XLN)*SIN(P(7,4)*6.28318531)),PI)/3. DO 351 N=1,3 351 P(14,N)=ARG(ATAN2(COS(RA+(N-1)*W),-SIN(RA+(N-1)*W) 1 *COS(E)-SIN(E)*TAN(XLN))/6.28318531,1.) W=ARG(ATAN2(SIN(P(7,1)*6.28318531)*COS(E),COS(P(7,1) 1 *6.28318531))-RA,PI)/3. DO 352 N=4,6 352 P(14,N)=ARG(ATAN2(SIN(RA+(N-4)*W),COS(RA+(N-4)*W)*COS(E) 1 )/6.28318531+.5,1.) DO 355 J=7,14 DO 355 N=1,6 355 P(J,N+6)=ARG(P(J,N)+.5,1.) CCCCCCCCC C C PARS FORTUNA & PLANETARY NODES FOLLOW C CCCCCCCCC P(20,12)=ARG(P(20,10)-P(20,1)+P(11,1),1.) DO 360 I=2,9 360 P(17,I)=ARG(P(20,I)-P(20,1)+P(11,1),1.) CCCCCCCCC C C ARGUMENTS TRANSFORMED TO SIGNS/DEGREES/MINUTES C CCCCCCCCC DO 390 I=1,4 DO 390 J=1,12 GLYPH(I,J)='N' IF(P(I,J))370,380,380 370 GLYPH(I,J)='S' 380 DEG=ABS(P(I,J))+.0083333333 ID(I,J)=DEG 390 IM(I,J)=60.*(DEG-ID(I,J)) IFINAL=JFINAL+19 DO 400 I=6,IFINAL DO 400 J=1,12 W=12.*P(I,J)+.00027777777 IS(I,J)=W ISS=W+1. GLYPH(I,J)=ZOD1(ISS) W=30.*(W-IS(I,J)) ID(I,J)=W 400 IM(I,J)=60.*(W-ID(I,J)) CCCCCCCCC C C PROGRESSION & EPHEMERIS OUTPUT ROUTINE FOLLOWS C CCCCCCCCC IF(OPTION.NE.'P')GO TO 420 WRITE(IOUT,14000) 14000 FORMAT(/19X,'Progressed (geocentric) longitudes'// 1 ' YEAR SUN MERC VENUS MARS JUPIT SATUR ', 2 ' URAN NEPT PLUTO MOON NODE') DO 410 I=20,120 II=I-20 410 WRITE(IOUT,15000)II,(ID(I,J),GLYPH(I,J),IM(I,J),J=1,11) 15000 FORMAT(' ',I3,1X,11(1X,I2,A1,I2)) GO TO 9999 420 IF(OPTION.NE.'E')GO TO 460 WRITE(IOUT,16000)IH(1,1),IH(1,2),NY 16000 FORMAT(//19X,'Daily ephemeris for ',I2,':',I2,' GMT'// 1 ' DATE SUN MERC VENUS MARS JUPIT SATUR URAN', 2 ' NEPT PLUTO MOON NODE'/I6) MY=NY+(NM+7)/10 MOND(2)=1461*(MY+4800)/4-1461*(MY+4799)/4-337 IF(CALEND.EQ.'J')GO TO 430 MOND(2)=MOND(2)-(MY+4900)/100*3/4+(MY+4899)/100*3/4 430 DO 460 I=20,120 WRITE(IOUT,17000)MONTH(NM),ND,(ID(I,J),GLYPH(I,J), 1 IM(I,J),J=1,11) 17000 FORMAT(' ',A2,I2,11(1X,I2,A1,I2)) ND=ND+1 IF(ND-MOND(NM))460,460,440 440 ND=ND-MOND(NM) NM=NM+1 IF(NM-12)460,460,450 450 NM=1 NY=NY+1 WRITE(IOUT,18000)NY 18000 FORMAT(I6) 460 CONTINUE GO TO 9999 CCCCCCCCC 9999 IF(LINEP)9997,9997,9996 9996 WRITE(5,96000) 96000 FORMAT(' $DONE') 9997 END