R1i{ZdZ ddlmZmZddlZddlZddlZddl Z ddl m Z m Z ddl mZmZddlmZ ddlmZddlZddlmZmZmZmZmZdd lmZdd lm Z m!Z!m"Z"dd l#m$Z$m%Z%dd l&m'Z'm(Z(dd l'm)Z)m*Z*mZm+Z+ddl,m-Z-m.Z.ddl/m0Z0edZ1de1_2e1jfZ3e1jhZ4e1jjZ5e6e1de1jnZ8dZ9dez Z:dZ;dZd8Z?e jdd9d:k(se jdd9d;k(rejne_8d<ZAd=ZBd>ZCd?ZDd@ZEdAZFdBZGdCZHdDZIdEZJdFZKdGZLdHZMdIZNdJZOdKZPdLZQe)jsdMndNZSdOZTdPZUdQZVdRZWdSZXdTZYdUZZdVZ[dWZ\dXZ]dYZ^dZZ_d[Z`d\Zad]Zbd^Zcd_Zdd`ZedaZfdbZgdcZhdxddZideZjdfZkdgZldhZmdiZndjZodkZpdlZqdmZrdnZsdoZtdpZudqZvdrZwdsZxdtZyduZzdvZ{e|dwk(reyy#ddlmZmZY xYw#e$r ddlmZYwxYw)yzTest suite for SGP4.)TestCasemainN) DocTestSuiteELLIPSIS)piisnan)get_data)StringIO)WGS72OLDWGS72WGS84Satrecjday)wgs72)invjdaynewtonnurv2coe) days2mdhms_day_of_year_to_month_day)sgp4sgp4init)api conveniencesioomm) export_omm export_tlesetUpi'assertRaisesRegexgH׊>gf@zE1 00005U 58002B 00179.78495062 .00000023 00000-0 28098-4 0 4753zE2 00005 34.2682 348.7242 1859667 331.7664 19.3264 10.82419157413667zE2 00007 34.2682 348.7242 1859667 331.7664 19.3264 10.82419157413669satnum satnum_str00005classificationU operationmodeiepochyr epochdaysg;Pyf@ jdsatepochg@BA jdsatepochFg`;P?bstarg0yv>ndotgXWKh=nddoteccogǐR?argpog0b)@inclogv“#?mogy?no_kozaig'7u.?nodeognuX@gxrr?r@s rB!test_legacy_built_with_twoline2rvrHQs --ue ,Cc$'rDct}|jtdtdtgt tt |yNr'r )rrr VANGUARD_ATTRSVANGUARD_EPOCH sgp4init_argsr?r@s rB%test_satrec_initialized_with_sgp4initrNUsA (CCLL  x   ~ & crDct}|jtdtdtgt tt |jdy)Nar )rrr rKrLrM assertEqualr&r@s rB3test_satrec_initialized_with_sgp4init_in_afspc_moderR`sI (CCLL  x   ~ & !!3'rDctj}ttdtdt gt t|fzt|dy)Nr'r TrF)model SatelliterrrKrLrMr?r@s rB%test_legacy_initialized_with_sgp4initrVksD // C  sN8,n ~ &# /c$'rDcd}d}tj||}tjddg}tjddg}|j ||\}}}t |ddgk(sJtj |jgdgd gk(sJtj |jgdgd gk(sJy) NzE1 44160U 19006AX 20162.79712247 +.00816806 +19088-3 +34711-2 0 9997zE2 44160 095.2472 272.0808 0216413 032.6694 328.7739 15.58006382062511g@BAgBAgq?r9)FFF)TTT)rr<nparray sgp4_arraylistrtolist)l1l2rAjd0jd1ervs rB9test_whether_array_logic_writes_nan_values_to_correct_rowrevs PB PB   B #C ((Iy) *C ((J + ,CnnS#&GAq! 7q!f   88A;   $9;M#N NN N 88A;   $9;M#N NN NrDc4tdd}t|dy)Ngfffff`@)r! !r/)rrQ)tups rBtest_days2mdhmsrls T9 %C()rDcZtdddddd\}}t|dt|dy) Ni 9gFBAg[?)rrQassertAlmostEqual)jdfrs rB test_jday2rvs/ $Ar2r *FBIb,-rDc Gddtj}tjddddddd|}tj|\}}t |d t |d y) NcHeZdZdZej dZdZdZdZ y)&test_jday_datetime..UTC_plus_4UTC)hoursc|jSr;offsetselfdatetimes rB utcoffsetz0test_jday_datetime..UTC_plus_4.utcoffset ;; rDcy)Nz UTC plus 4rs rBtznamez-test_jday_datetime..UTC_plus_4.tznamesrDc|jSr;r~rs rBdstz*test_jday_datetime..UTC_plus_4.dstrrDN) __name__ __module__ __qualname____doc__dt timedeltarrrrrrDrB UTC_plus_4rys$ A&   rDri r{rr!g@BAgͶJ?)rtzinforr jday_datetimerQrs)rdatetime_localrtrus rBtest_jday_datetimers]RYY[[r1aB:<HN  ' ' 7FBIb,-rDc tjtt}t j |}tj }t|tjddddddd|y)Nr92i1 ) rr<r=r>rsat_epoch_datetimerzrQrr)rArzones rBtest_sat_epoch_datetimersP   E5 )C..s3H   D"++dAr2r2vtLMrDcttfD]f}t|d}tt j ||tt j |dd|t j|hy)ND)r=r>intrQrcompute_checksum fix_checksumverify_checksum)linechecksums rBtest_good_tle_checksumrsYu tBx=B''-x8BOOD"I.5 4 rDctd}t|dtdddz}tttj |tt j |ty)Nr3r7)r=rQ assertRaises ValueErrorrrr)rbads rBtest_bad_tle_checksumrsJRyH# *s CR//5$e,rDcZttd}t|jdj }t gd}t ddg}|j d|D]}|jdst|}|dd }|dd }tj||}tj||t}t|\}} t|\} } |j|vrt||t| ||j|vst| |t| |y) a&Check `export_tle()` round-trip using all the TLEs in the test file. This iterates through the satellites in "SGP4-VER.TLE", generates `Satrec` objects and exports the TLEs. These exported TLEs are then compared to the original TLE, closing the loop (or the round-trip). SGP4-VER.TLEascii)ibeiq5i67rri.1NE)r riterdecode splitlinessetadd startswithnextrr<rrrr rQ) data tle_linesexpected_errs_line1expected_errs_line2line1line2satrec satrec_old actual_line1 actual_line2actual_line1_oldactual_line2_olds rBtest_tle_exportrs Hn -DT[[)4467IABuen-E"$ Ycr cr ""5%0]]5%7 &0%7" l-7 -C** == 3 3  e , (% 0 == 3 3  e , (% 0-rDcdD]d}t}tj}d||<|jtd|dt gt |ttt|fy)N)r2r5r1r3gr'r ) rrKcopyrr rLrMrrr)anglerAwrong_vanguard_attrss rB4test_export_tle_raises_error_for_out_of_range_anglesrsf0h-224&*U# 3,X6 / 0  ZS11rDcdD]A\}}tj||}t|\}}t||t||Cy)N))zE1 44542U 19061A 21180.78220369 -.00000015 00000-0 -66561+1 0 9997zE2 44542 54.7025 244.1098 0007981 318.8601 283.5781 1.86231125 12011)rr<rrQ)rrrAoutline1outline2s rB"test_tle_import_export_round_tripsrsH uu-'_(E8$E8$rDcZttjttt t tjtttttjtttt tjttyr;) assert_wgs72oldrr<r=r>r assert_wgs72r assert_wgs84r rrDrB-test_all_three_gravity_models_with_twoline2rvrs` F%%eUH=>""5%78""5%78""5%01rDcVt}tt}|jtdtdt g|t ||jtdtdt g|t||jtdtdt g|t|yrJ) rrMrKrr rLrr rr r)rAargss rB+test_all_three_gravity_models_with_sgp4initr s (C  (DCLL3x 8.P4PCCLL^H5~MMCLL^H5~MMrDr{rnc|jd\}}}t|ddtt|ddtt|ddty)NBڼZs@rg;TrXgXľ@r6g8o@ sgp4_tsincersGRAVITY_DIGITSrArbrcrds rBrr)sJoo01GAq!ad.?ad-~>ad-~>rDc|jd\}}}t|ddtt|ddtt|ddty)NrrgA>TrXgtXľ@r6g 8o@rrs rBrr/sJoo01GAq!ad/@ad-~>ad-~>rDc|jd\}}}t|ddtt|ddtt|ddty)Nrrg*|TrXg`-Zľ@r6guW:o@rrs rBrr5sJoo01GAq!ad/@ad.?ad-~>rDcd}d}tj||}t|jdt|jdt|j dy)NzE1 4859U 21001A 21007.63955392 .00000000 00000+0 00000+0 0 9990zE2 4859 000.0000 000.0000 0000000 000.0000 000.0000 01.00000000 09ir%21001A)rr<rQr r$intldesg)r^r_rAs rBtest_satnum_leading_spacesr?sM QB PB   B #C D!""C( h'rDc(d}gd}|D]9\}}||}t|j|t|j|;tt}|D]9\}}j t d|tg|t|j|;y)Nc~tjtjd|tjd|S)Nr#)rr<r=replacer>) satnum_strings rBmake_satz-test_satnum_alpha5_encoding..make_satJs2  w !F!&w !FH HrD))iA0000)i DE8493)iJ2931)i"P4018)ihW1928)i0Z9999r')rQr r"rMrKrr rL)rcasesr rrArs rBtest_satnum_alpha5_encodingrIsH  E"' }%CJJ'CNNM2"'  (D!&  UC?$?CJJ'"'rDct}ttd5|jtddt gt tdddy#1swYyxYw)Nzcannot exceed 339999r'i 0)rrrrr rLrMrKr@s rBtest_satnum_that_is_too_larger`sI (C :'= >      > *   ? > >s ,AAcltjtt}t |j dy)N58002B)rr<r=r>rQrr@s rBtest_intldesg_with_6_charactersrks"   E5 )C h'rDc\tjdd}t|jdy)NzE1 39444U 13066AE 20110.89708219 .00000236 00000-0 35029-4 0 9992zE2 39444 97.5597 114.3769 0059573 102.0933 258.6965 14.8209894934469713066AE)rr<rQrr@s rBtest_intldesg_with_7_charactersros,      C  i(rDct}|jtdtdtdz gt tt |jdy)Nr'r gv@c)rrr rKrLrMrQr(r@s rB+test_1990s_satrec_initialized_with_sgp4initrxsL (CCLL  x   ~ &  R rDc8t}d|_|jdk(sJd|_|jdk(sJd|_|jdk(sJd|_|jdk(sJd|_|j dk(sJd|_|j dk(sJy)NSabcdefgr{iabcde)rr$rephtypeelnumrevnumr"r@s rB test_settersr s (CC    $$ $CL <<9 $$ $CK ;;"  CI 99  CJ ::  CN >>W $$ $rDctdd\}}t|ddt|ddtdd\}}t|d dt|d dy) N?g333333@gKy @rh)placesg `8Sh@g?g@gRU~ @g<saA@)rrs)e0ms rBtest_hyperbolic_orbitrsX S# EBb+B7a+B7 S# EBb+B7a*26rDcltjtt}t |j dy)Nr)rr<r=r>rQr(r@s rBtest_correct_epochyrrs$   E5 )C QrDcvtjttt}t |j dy)Nr)rr<r=r>rrQr(r@s rBtest_legacy_epochyrrs$ --ue ,C T"rDctjttt}|j |j k(sJyr;)rr<r=r>rnor4)ss rB!test_support_for_old_no_attributers* eUE*A 441::  rDc:gd}d}t|dD]8\}}td|dzD]!}t|d}t||f||dz }#:d|d<d}t|dD]8\}}td|dzD]!}t|d}t||f||dz }#:y)N) rrrrrrrrrrXFT) enumeraterangerrQ) month_lengths day_of_yearmonthlengthdayrks rBtest_months_and_daysr$sEMK"=!4 vFQJ'C+K?C  c * 1 K(5 M!K"=!4 vFQJ'C+K>C  c * 1 K(5rDc Nd}d}tjddddddd}tj||t}t |j ||jtj }tj||}t tj||y) NzE1 25544U 98067A 19366.82137887 .00016717 00000-0 10270-3 0 9129zE2 25544 51.6392 96.6358 0005156 88.7140 271.4601 15.49497216 6061rgrXr*/i )r) rrrr<rrQepochrrrzrr)rPb correct_epochrAs rBtest_december_32r+s PAOAKKaBB?M --1e $C =)!))1A1A)BM   Aq !C //4mDrDctjdkryttt j d5t jtjddttdddy#1swYyxYw)Nzyour TLE lines are broken because they contain non-ASCII characters: 1 00005U 58002B 00179.78495062 .00000023\xa0 00000-0 28098-4 0 4753 2 00005 34.2682 348.7242 1859667 331.7664 19.3264 10.8241915741366723 u23 ) sys version_inforrreescaperr<r=rr>rrrDrBtest_non_ascii_first_liner4s` $ :ryy2I(J K emmE84eUC  K K K 4A55A>ctjdkryttt j d5t jttjddtdddy#1swYyxYw)Nr-zyour TLE lines are broken because they contain non-ASCII characters: 1 00005U 58002B 00179.78495062 .00000023 00000-0 28098-4 0 4753 2 00005 \xa034.2682\xa0348.7242 1859667 331.7664 19.3264 10.82419157413667 34u 34) r0r1rrr2r3rr<r=r>rrrrDrBtest_non_ascii_second_liner8sZ $ :ryy2(  eU]]5(;UC    r5ctttjd5t j t jddttdddy#1swYyxYw)NaTLE format error The Two-Line Element (TLE) format was designed for punch cards, and so is very strict about the position of every period, space, and digit. Your line does not quite match. Here is the official format for line 1 with an N where each digit should go, followed by the line you provided: 1 NNNNNC NNNNNAAA NNNNN.NNNNNNNN +.NNNNNNNN +NNNNN-N +NNNNN-N N NNNNN 1 00005U 58002B 00179.78495062 .000000234 00000-0 28098-4 0 4753r/234) rrr2r3rr<r=rr>rrrDrBtest_bad_first_liner;sO :ryy2I(J K emmE515%@ K K K 4A!!A*ctttjd5t j t tjddtdddy#1swYyxYw)NaTLE format error The Two-Line Element (TLE) format was designed for punch cards, and so is very strict about the position of every period, space, and digit. Your line does not quite match. Here is the official format for line 2 with an N where each digit should go, followed by the line you provided: 2 NNNNN NNN.NNNN NNN.NNNN NNNNNNN NNN.NNNN NNN.NNNN NN.NNNNNNNNNNNNNN 2 00005 34 .268234 8.7242 1859667 331.7664 19.3264 10.82419157413667r7z34 ) rrr2r3rr<r=r>rrrrDrBtest_bad_second_liner>sO :ryy2I(J K eU]]5%8%@ K K Kr<cd}tttj|5t j t ttdddy#1swYyxYw)Nz,Object numbers in lines 1 and 2 do not match) rrr2r3rr<r=BAD2r)msgs rBtest_mismatched_linesrBs5 8C :ryy~ 6 eT5) 7 6 6s $AAct}|r|j}|d=|d=|d=|d=|jD]\}} tt |||y#t $r#}|j \}d|d|f|_|d}~wwxYw)Nr(r)r*r+zfor attribute z, )rKritemsrQgetattrAssertionErrorr)rArGattrsnamevaluerbmessages rBr?r?s E   )  +  ,  - {{} e  T*E 2% vvHG04g>@AFG sA BA<<Bc N|d|d|d|d|d|d|d|d|d f S) zBGiven a dict of orbital parameters, return them in sgp4init order.r,r-r.r0r1r2r3r4r5r)ds rBrMrM"sD gJ& 1W:qy!G* gJ$:'  <test_satrec_against_tcppver_using_julian_dates..invoke1si &#6x    &   ( *++b"%1aFLL   !QwrDrXrXr9r9r{r.r9run_satellite_against_tcppverrr<rUs rB.test_satrec_against_tcppver_using_julian_datesrZ/s"&"3"3V_MrDcDd}ttj|gdy)NcX|j|\}}}||jk(sJ|||fSr;)rrQ)rrRrbrcrds rBrUz8test_satrec_against_tcppver_using_tsince..invoke=s5$$V,1aFLL   !QwrDrVrWrYs rB(test_satrec_against_tcppver_using_tsincer];s "&"3"3V_MrDc2d}d}gd}t|||y)Nc<tj||t}|Sr;)rr<r)rrrAs rBmake_legacy_satellitez:test_legacy_against_tcppver..make_legacy_satelliteFsmmE5%0 rDcVt||\}}|j|jf||fSr;)rrQ error_message)rrRrcrds rBrun_legacy_sgp4z4test_legacy_against_tcppver..run_legacy_sgp4Js.FF#1 f223Q99rD))rXz;mean eccentricity -0.001329 not within range 0.0 <= e < 1.0)rXz;mean eccentricity -0.001208 not within range 0.0 <= e < 1.0)r9zBmrt 0.996159 is less than 1.0 indicating the satellite has decayed)r9zBmrt 0.996252 is less than 1.0 indicating the satellite has decayed)r{z,semilatus rectum -0.103223 is less than zero)r.zCperturbed eccentricity -122.217193 not within range 0.0 <= e <= 1.0)r9zBmrt 0.830534 is less than 1.0 indicating the satellite has decayed)rX)r`rcerrss rBtest_legacy_against_tcppverreDs": D""7$OrDcttd}|jdd}|jdj d}g}t t |||}t|t|cxk(rdk(sJJd}t||}t|dD]\} \} } | d k(rd |d d z} d | vr| | k(} n| j} | j}| ddDcgc] }t|}}|ddDcgc] }t|}}t|t|k(xr)tdt||Dxr | dd|ddk(} | std| | | fzd | vs| }t||ycc}wcc}w)Nz tcppver.out rDrTirX)start(Use previous data line)z 0.00000000kxxr8c3ZK|]#\}}t ||z cxkxr tknc%ywr;)rQ).0rPrbs rB z0run_satellite_against_tcppver..s/!8AFa!e,u,,!8s)+zc#(Kttd}t|jdj }|D]Q}|j dst |}|||}d|jfzt||||D]}|Syw)zGenerate lines like those in the test file tcppver.out. This iterates through the satellites in "SGP4-VER.TLE", which are each supplemented with a time start/stop/step over which we are supposed to print results. rrrz%ld xx N) r rrrrrrr generate_satellite_output) r<rUrxrrrrrrs rBrprps Hn -DT[[)4467I$ YE5)FMM+++-z3DJ3sBBc#lKtj}||d\}}}t|dr2t|dr$t|dr|j|dyt d||d|ddj D\}} } |} | | krY| |cxk(rdk(r nn| | z } ||| \}}}|dk7r|d k7r|j|yt || |||| | z } | | krY| | z | d z kr;||| \}}}|dk7r|d k7r|j|yt || |||yyw) zAPrint a data line for each time in line2's start/stop/step field.r/rrXr6riNc32K|]}t|ywr;)rt)rnfields rBroz,generate_satellite_output..sH5GE5<5Gsr)rNgư>)rmurappendformat_short_linersformat_long_line) rrUrrxrrbrcrdtstarttendtsteprRs rBrrs\ BVS!GAq! QqT{uQqT{uQqT{!(( CA &&HU23Z5E5E5GHFD% F D. V "s " eOF (1a 6a9n   a vvr1a88% D.}ut|#&1a 6a9n   a vtRA66 $s C*D4-AD4c >d||d|d|d|d|d|dfzS)z?Short line, using the same format string that testcpp.cpp uses.z2 %16.8f %16.8f %16.8f %16.8f %12.9f %12.9f %12.9f rrXr6r)rRrcrds rBrrs? A!adAaD!A$!adD4 44rDc>t|||jd}|j|jz|dz z}t |\}}} } } } t |||\ } }}}}}}}}}}|d|||t z|t z|t z|t z|t z||| | | | f zzS)z>Long line, using the same format string that testcpp.cpp uses. rOzI %14.6f %8.6f %10.5f %10.5f %10.5f %10.5f %10.5f %5i%3i%3i %2i:%2i:%9.6f )rstripr*r+rrrad)rrRrrcrdshortrtyearmonr#hrminutesecprPeccinclnodeargpnurarglattruelonlonpers rBrrs fa + 1 1$ 7E   V// /&6/ AB&-bk#D#sB 1b Q3dD"a&  % 3S$s(DHbf #tS#r63  rDzVANGUARD 1 1 00005U 58002B 20287.20333880 -.00000016 00000-0 -22483-4 0 9998 2 00005 34.2443 225.5254 1845686 162.2516 205.2356 10.84869164218149 zMARIO 1 55123U 98067UQ 23115.44827133 .00787702 29408-3 15680-2 0 9999 2 55123 51.6242 216.2930 0014649 331.8976 28.1241 15.99081912 18396 a
MARIO1998-067UQEARTHTEMEUTCSGP42023-04-25T10:45:30.64291215.99081912.001464951.6242216.2930331.897628.12410U551239991839.1568E-2.787702E-2.29408E-3 avOBJECT_NAME,OBJECT_ID,EPOCH,MEAN_MOTION,ECCENTRICITY,INCLINATION,RA_OF_ASC_NODE,ARG_OF_PERICENTER,MEAN_ANOMALY,EPHEMERIS_TYPE,CLASSIFICATION_TYPE,NORAD_CAT_ID,ELEMENT_SET_NO,REV_AT_EPOCH,BSTAR,MEAN_MOTION_DOT,MEAN_MOTION_DDOT MARIO,1998-067UQ,2023-04-25T10:45:30.642912,15.99081912,.0014649,51.6242,216.2930,331.8976,28.1241,0,U,55123,999,1839,.1568E-2,.787702E-2,.29408E-3 ctj\}}}tj||}t t j tt}t}t j||t||yr;) MARIO_TLErrr<rr parse_xmlr MARIO_XML initializeassert_satellites_matchline0rrsat1fieldssat2s rBtest_omm_xml_matches_old_tler`#..0E5%   UE *D #-- 34 5F 8DNN4 D$'rDctj\}}}tj||}t t j tt}t}t j||t||yr;) rrrr<rr parse_csvr MARIO_CSVrrrs rBtest_omm_csv_matches_old_tlerrrDc t|D]Q}|jdrt||d}|%t|r1t||}t |||d|d|Sy)N_ z != )dirrrEcallablerQ)rrattrvalue1value2s rBrr's_D  ??3  tT* >  F  t$FFT66$JKrDc tj\}}}tj||}t |d}t |idddddddd d d d d ddddddddddddddddddd d!d"dd#d$d%d&d'y)(Nz VANGUARD 1ARG_OF_PERICENTERguq Hd@BSTARg+; CENTER_NAMEEARTHCLASSIFICATION_TYPEr% ECCENTRICITYg1+o?ELEMENT_SET_NOiEPHEMERIS_TYPErEPOCHz2020-10-13T04:52:48.472320 INCLINATIONg48EA@ MEAN_ANOMALYg i@MEAN_ELEMENT_THEORYSGP4 MEAN_MOTIONg`쵇%@MEAN_MOTION_DDOTr/MEAN_MOTION_DOTg:0y NORAD_CAT_IDr! OBJECT_IDz 1958-002B OBJECT_NAMEg0*0l@TEMEi6Urz)RA_OF_ASC_NODE REF_FRAME REV_AT_EPOCH TIME_SYSTEM) VANGUARD_TLErrr<rrQ)rrrrArs rBtest_omm_exportr6s &113E5%   E5 )C \ *FX w s    #  ! - w  v { C 8  [!" |#$#+rDc ddlm}|||ttjdk\rfd}d}t t ||}|jtdi||jtdi||jtd i||S) z>Run our main documentation as a test, plus all test functions.r)add_test_functionsr7ctj|_tjtjj t tj|_ dt_y)NT) osgetcwdolddirchdirpathdirname__file__r acceleratedoldaccelsuites rBrzload_tests..setUp_s:99;EL HHRWW__X. / __EN"COrDcltj|j|jt_yr;)rrrrrrrs rBtearDownzload_tests..tearDownds HHU\\ "#nnCOrD) optionflagsrr)r)zsgp4.conveniences)zsgp4.functions) sgp4.wulfgarrrr0r1dictraddTestsr)loadertestsignorerrroptionss rB load_testsrUs~0vuh/ 6! #  -858L |6g67 |C7CD |@@A LrD__main__)F)}r unittest2rrunittestrrr2rr0doctestrrmathrrpkgutilr rr ImportErrornumpyrYsgp4.apir r r rrsgp4.earth_gravityrsgp4.extrrrsgp4.functionsrrsgp4.propagationrrrrr sgp4.exporterrr sgp4.modelrT _testcasemaxDiffrQrsrrEassertRaisesRegexprrQrr=r>r@rKrLr1rCrHrNrRrVrerlrvrrrrrrrrrrrrrrrrrrrrr rrrrr$r+r4r8r;r>rBr?rMrZr]rerXrprrrrrrrrrrrrrrrDrBrsX(( *"99$..@+++0 W   ## //%% I':%88:  bjOOO a'c  S q): Z  S I  !" #$ %&#'( ), BQ6!S%5%5bq%9V%C!)!4 zFU(''"!!"sH H H H*)H*