X-Sympa-To: ids.analysis.forum Received: from hector.cls.fr (osis.cls.fr [62.161.32.21]) by hermes.cls.fr (8.9.3/8.9.3) with ESMTP id WAA26826 for ; Fri, 24 Mar 2006 22:11:08 +0100 Received: by hector.cls.fr (Postfix) id 1E82686ABB; Fri, 24 Mar 2006 20:13:03 +0000 (GMT) Delivered-To: ids.analysis.forum@cls.fr Received: from cluster-d.mailcontrol.com (orange-103.cls.fr [62.161.32.14]) by hector.cls.fr (Postfix) with ESMTP id 2D4AD867DC for ; Fri, 24 Mar 2006 20:13:01 +0000 (GMT) Received: from bowie.gsfc.nasa.gov (bowie.gsfc.nasa.gov [128.183.106.180]) by rly38d.srv.mailcontrol.com (MailControl) with ESMTP id k2OKCv4O018044 for ; Fri, 24 Mar 2006 20:12:59 GMT Received: by bowie.gsfc.nasa.gov (Postfix, from userid 247) id 398B82855; Fri, 24 Mar 2006 15:12:56 -0500 (EST) X-UIDL: 5?G"!dT,!!^1$#!_Ej!! X-Original-To: flemoine@bowie.gsfc.nasa.gov Received: from gsfcrelay3.gsfc.nasa.gov (gsfcrelay3.gsfc.nasa.gov [128.183.244.243]) by bowie.gsfc.nasa.gov (Postfix) with ESMTP id 293AD2869 for ; Thu, 23 Mar 2006 18:00:15 -0500 (EST) Received: from gsfcrelay3.gsfc.nasa.gov (localhost [127.0.0.1]) by gsfcrelay3.gsfc.nasa.gov (8.12.11.20060308/8.12.8) with ESMTP id k2NN0Dd5014863 for ; Thu, 23 Mar 2006 18:00:13 -0500 Received: (from root@localhost) by gsfcrelay3.gsfc.nasa.gov (8.12.11.20060308/8.12.9/Submit) with LDAP id k2NN0ChA014850 for flemoine@bowie.gsfc.nasa.gov; Thu, 23 Mar 2006 18:00:12 -0500 Received: from sol.csr.utexas.edu (sol.csr.utexas.edu [129.116.190.1]) by gsfcrelay3.gsfc.nasa.gov (8.12.11.20060308/8.12.8) with ESMTP id k2NN0BaW014845 for ; Thu, 23 Mar 2006 18:00:12 -0500 Received: from localhost (localhost.localdomain [127.0.0.1]) by sol.csr.utexas.edu (8.12.11.20060308/8.12.11) with ESMTP id k2NMxsYK014372; Thu, 23 Mar 2006 16:59:54 -0600 Received: from sol.csr.utexas.edu ([127.0.0.1]) by localhost (sol.csr.utexas.edu [127.0.0.1]) (amavisd-new, port 10024) with LMTP id 14064-10; Thu, 23 Mar 2006 16:59:52 -0600 (CST) Received: from [129.116.190.133] (ries.csr.utexas.edu [129.116.190.133]) by sol.csr.utexas.edu (8.12.11.20060308/8.12.11) with ESMTP id k2NMxl3F014351; Thu, 23 Mar 2006 16:59:48 -0600 Mime-Version: 1.0 X-Sender: ries@sol.csr.utexas.edu Message-Id: Date: Thu, 23 Mar 2006 16:59:42 -0600 To: altamimi@ensg.ign.fr, Bruce.J.Haines@jpl.nasa.gov, eelco@deos.tudelft.nl, flechtne@gfz-potsdam.de, Frank.Lemoine@gsfc.nasa.gov, Jean-Michel.Lemoine@cnes.fr, jean-paul.berthias@cnes.fr, luca.cerri@cnes.fr, nzelensk@bowie.gsfc.nasa.gov, pascal.r.willis@jpl.nasa.gov, pschaeffer@cls.fr, Richard.Biancale@cnes.fr, sabine.houry@cnes.fr, Scott.B.Luthcke@nasa.gov, wib@cobra.jpl.nasa.gov From: "John C. Ries" Subject: SRP test Content-Type: multipart/Related; boundary="============_-1068967708==_mr============" ; type="text/html" X-Virus-Scanned: amavisd-new at csr.utexas.edu Sender: flemoine@bowie.gsfc.nasa.gov X-Modified-HTML: 1 X-Mailcontrol-Inbound: Bl3AkpjjvVpE8ExYxQDHVBwHPXXURTA5THfcd2uWoYQ= X-Spam-Score: 2.4 X-Scanned-By: MailControl A-06-00-05 (www.mailcontrol.com) on 10.68.0.148 X-Validation-by: flemoine@bowie.gsfc.nasa.gov --============_-1068967708==_mr============ Content-Type: multipart/alternative; boundary="============_-1068967708==_ma============" --============_-1068967708==_ma============ Content-Type: text/plain; charset="us-ascii" ; format="flowed" Hi Everyone, Bruce suggested a test where I look at the impact of the SRP model on the SLR/DORIS dynamic orbits. So I re-ran the first 90 cycles with a 3% change to the full SRP (box-wing) effect. It does indeed lead to a 120-day (117 days to be exact) variation in the centering, mainly in Z. There was an effect on X and Y, but it was only 0.5 mm. The effect on Z was about 2 mm (see plot, where I multiplied the z-shift by -30 and compared it to the beta-prime angle...you can see the correlation). The effect on the inertial centering, as measured by the crossover means, was negligible. So I presume a larger scaling of the SRP could lead to a larger Z variation. I've tested the effect of drifts in the Z component of the reference frame and found that the orbits tended to follow the Z-drift, at least about 90% of it. But it looks like there is some wiggle room in that remaining 10%, and one can get some modulation on top of that due to SRP. JR --============_-1068967708==_ma============ Content-Type: text/html; charset="us-ascii" SRP test
Hi Everyone,

Bruce suggested a test where I look at the impact of the SRP model on the SLR/DORIS dynamic orbits. So I re-ran the first 90 cycles with a 3% change to the full SRP (box-wing) effect. It does indeed lead to a 120-day (117 days to be exact) variation in the centering, mainly in Z. There was an effect on X and Y, but it was only 0.5 mm. The effect on Z was about 2 mm (see plot, where I multiplied the z-shift by -30 and compared it to the beta-prime angle...you can see the correlation). The effect on the inertial centering, as measured by the crossover means, was negligible.

So I presume a larger scaling of the SRP could lead to a larger Z variation.

I've tested the effect of drifts in the Z component of the reference frame and found that the orbits tended to follow the Z-drift, at least about 90% of it. But it looks like there is some wiggle room in that remaining 10%, and one can get some modulation on top of that due to SRP.

JR



Click here to report this email as spam.

 --============_-1068967708==_ma============-- --============_-1068967708==_mr============ Content-Type: image/png; name="srp_test_Z.png" ; x-mac-type="504E4766" ; x-mac-creator="514B5054" Content-Disposition: attachment; filename="srp_test_Z.png" Content-Transfer-Encoding: base64 iVBORw0KGgoAAAANSUhEUgAAAokAAAG5CAMAAADVm14jAAADAFBMVEUAAAAd CAAdFgAAABkoAwA3AAAaJgA+KwgAACkAADkhACErADw6ADoNKzUhKyNEAABV AABCAB1oDABHLw14NwBEADxAPjtvQhZURiplXDwAAUkAAFUHJUcAAGMABXoA EncAM3VUAEA3SVI4VGhVUVJ1SkVubVJ4d3GEDgCPKwCQKwCLMwDdCAbeHQbd CBHhCAbeLQbmLgbdCCndCDbgCDqXSACYUwqZVhekXACtZRKocSaudTO4ejDp TwrvbhXiTCXdCETdCFTdCnXdI3eaeVOIcmHhU1bmZlHeSHLkaWLMiybLjDvz izewil63j3fGikXMk0zSm0/KlVTUnVv0kEjWoFj3pF3nh3fXpWDNqXvXsn7g r27zo2jjtHX4uXgADIAAGogAKY4AN4sAK5AARpcAUZcpZpcAUaIIZ68WZqUA YLMSdbs3capneY9Ff7FierAAANUAANQxANUGKtkrKNQABuEAI+NMANR7E9R5 MtQ6TdwASuYwSuAHbutTT9hoXddKY9xqZ9t5W+FVbOVmc+jdC4XdJIvdMJLd Q5TjXZLjdpPjdquIANSOMdSPTNSkU9SWe9q0ctXdfcV1lq4wicoviPBGi8VF k81VkcRFmNVbntN3ktNsocxipdh6s9xRk+9sle9ysvB8xe+Zh4KJjZiUl5in mZG6p4+8o5WPmKmpmqusrq/smpfEp5Lsu4LjibnKvavtqKukx7nqw4P7y4np yZT3y5T10pbZxLb5za/4z63526Xny7b62rX84q365buZgdyygNmOscuwtsqW leaomuSSsO6vr+rWjM/ilsrZqd3pt9bAk+HNt+Xlu+C0yNWLyfqHx+yWyeeS z/SW0fCt0vm70uep2vaw3va04/215frAwMPdzsfb08zNztDZ1tf0zc3o1cL2 083zytTo3tTz0dbe5df+7cL/9Mrs7Nj56dj+9tXM0O/xz+Tnx+z02uXr0vHU 4e/I7f3Y6PjH8f/U9v/p6On46+v9/Ovs6/n16vbq/f7+/v4AAAAAAAAAAAAq ndOhAAAgAElEQVR4nO2dbbRdR3nft2PUSDGtsFPShEsMi7g2wTsFTE1cp6TQ phAiUrBMUicUloNF29Rp1l4KtkiIlXqVUNrUBCemSZvGhNQ4JiR2QLZlfbFW 8IsUvMBaXet+uEsmgRKWZOveD8j3SufqzIfu97Nf5uV55pmZvfec+cu+59yz 98wzs+d3/vOy55wbsaCgMSgaugBBQbkCiUHjUCAxaBwKJE5I0XQFqJz96xdk StNtrECiX5puYwUS/dJ0GyuQ6Jem21iBRL803cYKJPolXGNtHLzzzoN3fmxr 8cJ/Lh4+NmfHk9ULxYtf+aN+wvSEC/nP+vfFL7l+Y2uRXUdPPbWaPRx7Mn04 /llEyQOJExKusc4/+eSTx/Ynm/ULG7flOG18jLGPfGa7OHBh/9lewo3bGEs2 s9Pq3xe/ZJrt3yxP6+rCgV85mHyGsd9O7kw+y87vr4IHEv0SvrG+lpzIH7ef fHKdbdw+P/Z0CuM6W0+OfXV/9pQdv4NdWN8+tsrWn1pPn7Hs/5TY7Gj2+/ZT x+ZZujQNY199OvvJjh3bPpCTePv6U0+nuRWJCj2Twn789vlGivlXki1296fh JQ8kTkjoxtrI/Cl/PPipZG0jSQ7uvyP9ZX538mt3Jr+edrGzAydY9nLy4fTf Vmaa6f8peekJX02fJMkf7P9o9nt25FdTn0uzO5D8RpL7aXr0P6VHf/sOxu7+ aBUw65M/mv2XG+LXbpuDSx5InJCwjXXhwB3Fk8yfnlrdSNZK0FhqarmvsZzE tRSl9LUT2ZEMu/yE7EmWIDlbvPS1FL+NZCt1vDS7gsSz+dE07wMn6pjHk9SG j3+6gnwTXPJA4oRUNlYCUH7iJ1NIMjQ+vbE/+bXVjKMCrNSvzhZjuI0Crxyd uxYkZieUSM5yQNP/jicHDx5M1rIO93xB4u3ZiPFERnAWJw+UddSfSnHN3gEH 1hjbv9YqOaByQVMQksS7k4yD7fX1F9MfdyefydBJ0cv4SjFckJi+XJGYjvSK rnh/eloxxTmwlvnobfPjt6+vrx+b330HK60uy4cdWE0Hg3+ZW++FNNDd2Wx5 /+rx2wtcz1fzoUCiXxI1VuWBbX0lnb1ub29nT38nReWTn84tLyXxtnk6+Z0V dpV1oenLJYkpup/MDTA9+q2sd15N5zxbhSc+k3fGm8dT//vLsnfOjmaemtSd c9Z3H8967FV2PBt0JoFEL4VrrLsKd8xY2Eh+ZX+ymk9IchKzXvMjOU53fTZ3 vqw7/e9/nCcpPPEjyWqSzVg+nE5SCk9M+/oPJ+lc+VeT/b9ezVg+nGQeeNft 9WLjxoE0hz/K7DhJYazWjQKJvgnZWNu58qcX0n41fYHl/xf/5QszaWec/3Kh PLq+zS5UR7czt1yfs0WC/Jds2aZYFZ+ztN9PdddnGkHXi9eKh7uqpe1Aol+y 0VgfOSE+VluaTBc+1Vg876RPqvSBRL9ko7E2npYc698J5Jx0cFV06Kv1kUCi X5puYwUS/dJ0GyuQ6Jem21iBRL803cayROJj30l/zB64b3XxEORCgcS2ju7I Jv7X7boyurd+CHKhQGJL10ZRSuLhlbPsppWt8kGjbEF4BRJbOnfmspTEG/cx duqSzRuKB42yBeEVSOzodSmJe+5h7PTuE+WDRiZBeAUS25r9cEbirSmCl66V D/hMgjQUSGyrIPEWxh7ftVk8dD6UU++RC/9I/3pXcMokNmqVcHex6fbOhy5e ZXtfNc8fbu7cJ+dtlgsyoCmT2JYZEmevzfb6XB9dtHO1flCFsSe32A9btUBi R+fyn8+dfLHxoAgTZEK4xjrzzne+8+3/7UL9+3nBeUd+i1AkoCyRqFLwREvB cI11+vUPPfhnb/nn9a/XCM478hCuWDpaBhKXSUgS35T+ePxNqSl+8YsbjD16 zTPs/JEvfrs4eGbj2w8+w9gL33rhmdlG9pCecuSh9LTiZNNaBhKDJ4p0+h9f YOf/zcfZmTfG74z/ir0j/skz8U++Iz6eH3w4fuu/jj+RPsSfePg3s4e3x2+4 5o1vq042rWUgcZmE7Z3jVL/FPvA2xg6/lZ15EzvyXxn72Y/nBx+9ZpsdedOF R+MN9vDH899+/BPszD/aKE+2XXIfSVxCT4whKj3xyOufeW/89rf/s3j7VMrX kXf8SGqE74rj4yl+Kaobj6bDyPTpw+nDu46z2Rs3ypNNV2UZSFwmIUnMhojs Lcff+y++/Y1vPJP9+oEfeej8Bz6ejgq/cf7wb6ZjxddvZ0BmJKYP70r75JTE 4mRLJa/lI4nBE6WemHa/D771Anv4bezUNRfem+L3lk/keWTd8oNvY00Sc08s TzatZSBxmaQzTnyInXlL/Ib4mRSztz6a/v6OYj6SPn1DeiydrWT/Zw8//juM vX6jPNl2yX0kcQk9EarzqfIn/+9b2c98teYCeyFf6374v8zyF7Nfzuf/zYpn 5cmGtQwkLpMMNlbmgg61DCQGT9SSixsrDS0DicuksAMCqeCJloIFEpEKnmhJ gUSk+HtyLSl44hTUKrm5PdsqBU+0JE9IZH72zkvlidOVqmoekBg0RflI4jJ5 oj/RfCQxaIrykcTlMg5fovlIYtAU5SOJPhhHvp3QVTCRhm61yZPogWIBiR7L RxI9MA4RiR5UDRFs8iR6oHj5TNFHEj0wDhGJHlQNEWzyJHqg4ImZJk9iIpx6 WopmXHHwxEyTJ5E5JdGGhCR6LB9JTJx2bsETrQWbPInTH2YFT8zlwZ5tPzyR UwNfPdHbPdtTtxQhiR7LRxK98ET+pMtXT+QHmzyJ07cUIYkey0cSE6ckBk+0 FmzyJE7fUiZfAQ35SGLitCGDJ1oLNgYSaZ3r1C0lnnoFdOQjiZP3xIrEfg2C J1oIIxNxxjF1S8nLPu3ZP14+kii2FBuy44mMv4sjeKKFMDJNiUQLEpPosTwl kblryOCJ1oINT2JMW5qevKVMvgI6CiRSFTzRWrBRkEghafKWMvkK6CiQSBW/ auQZFz+L4IkWwkg0KRJFRbBCosca5Z5tGomOG5JbLwPjXO4V8NUTR7pnO54U iaIiWCHRY/lJYuvRsnhVs0air57IDzYGEtl0SBQVQR/FEVRgAI2aRL2GGIcn 2iExeKKFMGJNi0RhEYInojReErUbYnhPpA104+CJpQKJVJkjcThTdB58hCTG NBJllmJDfE8kxK/TcWBw5lKkyb+W/CZxIEuhDXQXyQZfEXUZcZwkth6RklmK Dck80TiJrjyRtiCqpXGTqHMpRmApZFPvPXOtQGKm7ApktyEnQmK/avThRe+Z OJgdlQuijqLlGieJxQ3xaZDILYEHnuh6nO0liUn9TCM9VlxPpMQfgSfGgURW XIVij5A+iUn1TCM9XbSBbjwCT6QtiOppnCQyfU9cgDwGTyST2MnAnSeypScx KepPIpFRene6mutIhNQDroj6TCJ4z7YBEhNWe6qDS8n1RMLwoJloaTxxjHu2 2yTiL0WZxh2JvBJYI9GNYtrSvJ58JDEp07tpyF7V4nrONFVPDCTmZy5I1GrJ hicOZCnFQLVRFnxy3nOHIi5D6WlsJCYGSCw80dWUheeJ1kh05on88DY1XhL1 SEqqNMNNnutVqOCJGI2OxLjRvWqSmBDS49WtWvFeoHh68MRSkyYxHtwTKxI1 p1xdEodydefhx0hi8USfxKYn2r+SfE/Ujt9K0svAiSfGgUTWJxF7JeoUQy0o Uge67SSDeOLisgUSGYnE2hOdXMlO1SoSNU09H+cubjgM5YmLwrgIWATrvzQe EjVasuWJQ1iKIRKdbmzraBF0qiQ+d/LkyXXGZg/ctwoIIyiPARKH9MSEOOVq Lqey4Il6ml0dRdG+Obtu15XRveowgvIsSNRpyWRoTzRG4oBbLFue6G4TRP8l XRJPv3wzezi8cpbdtLKlDMMvThslPRIH90RWFUCbRCYg0YUnxi1PnCSJL/3S /c8yduM+xk5dsqkMwy+OCRLL5EMsKCY9EnEFaDkqG8YTm4WeJomHold+f/SL bM89KZS7TyjD8ItDJDGuPdHZ0narai0SNT6e2CCxunPdysCRJ5bvhGSiJL7w BGOP7Fzbc2tK4qVryjD84vRIxF2JRoJBFhQT4lupM84dZEV0QeJUPfH8PEVw x+qeWxh7fNfZTphSTPEvKe6RFM8T1vgN9q95vk566r+SxCp+go3fTJ1frfr3 RZ62/9XrSEVlHERMGutWJkj8qcvnbO/K1qGLV9nem+cdEoF5NMbIWlOOxvmD TJ6TRodG80TKxjKacvidfigtl0ESj74s+q7oy4xdH120s7ugiCGxeqpN4mKc 6OJKdsaJsTkSORVwNU4sSNT+Hg4dmVzZnv3NyRezx+eKB0UYfnE6LYG8EIkR T8T8wY5eUmMkDnaXqCKxDD9FEnFh+Od1WkKLRKon4khsntycsFgh0YEnFssP Sf6fy3eCjyQukmtdSNQfMeomZT0SUQWIh/fEuJixOB+njovEhEhibMITcX9P q+OJlkl044l5aNffXuAjiYvkLkhsJ+2QqFGBZjkG8sRyXdutJ4+fRMyVKM6u PFFzadukJyJbsr35hUOiumCE6VahBonL7IldU9MiUT99kWw4T1SSCCoCTXH1 DVeBxPoX9JVImp6o27vUg3Xo+e20pGUoI55IRDGuv/XP7f3S8ZHY+E2TxEZy XRJjvRZIeKaOSU83dSqJcSAxP80IiSRPrG4XgxO2PdEyifY9sUWiy+7ZRxIb yfXGWSWJOtPu9uQX6yn0ChSjXAI+9YSFaQ5UdeUjiURPLPABI9TzxF5W8OCx ikTldaSMLaoiJAsSHe5QHBWJCZFEZUOCylDmY8YTaRXQuMtkhMRGAQKJTMNT YronLlofGrlRtcQAie3CDOSJdWaBxPJXPImt3HovQYpQxdRozU7DWSARUgTt UW5VhAaJHo4TQVO6zmfJ0C1ZnkvzxJpBHU+0TSLEE5lpEh2gOLbv2e6RiHtP dj4orrUIUabRc4P2Ig6exO4n3fU8kemOcnPF9R9WqntnN6Y4OhJbv+uRSOhc kjaJkLQtT7RCYvPNqcxiQaIeQF0SnXXPPpLYSY5aWk4Wa7v4xuSMc1G5mKuA GRKdDhR9JJEwzGmSCAzd9cTmr8gC8Orf2dWhzKLqT4MngsPwzjLvicil8Sqi 7lokbfbepVZjD0iTRC2CgifmZ9F6J54nYhuStUkEpBV7IpnEPgmq67gY5mqS 2BnoLimJ3YbUJbGZXJtEnYGiARL75UFXQHvqXxaBsqJL0OhJRFxRM55YXHl4 WrOeyC0PNxi/MOipf78Izc55WUnsj5OYDomayZkeia3ko/FE3e65RaLmcFlP IycR1RJVcrInotI2LSS2TSLAE+NqnKtLYjM3hwPFsZHYeUWLxGZq+yQ2UwtI RFSAWx7E2nibRDxCtNULCrQekbj45sI6NXY9r0+iOm3HE9v5QTMpTleTCPNE wlSj48CBxOoVLRK7yREk1uej748wEYnUCuDuF9YkanarwRPzkwx7IrohWyRC 06o8kUZi5wJAPVHXFNtJsJ5ImmZPgkRY9eiW0rzuGgNF8ySi7zKZJRGbjS8k 9hpSl0SqJ+LSSj0RVQAAiYrruFi80eyeOzMuPInIeK1gPY2MRHD97Hgiaubd PR1XAE40tCculrK0TJFcAVy4Vqj+S0Pt2eb89Q/UhaiTG/VEZeKGJ/ZON04i yBObKUkkJs5IHNeebS6JCE8ggsxaloId4+VJ6CRySoSvgEZoURFwnfwkPBFy DrEl6vMontggEZi47Ym9DBHOJK5//TrGE3VMseuJuLdCILEQr2/DtUZC+qov yTjXGInKIjTej3QSWSBx8RrJEykNCSax4cB8EglvpV7vqPTEuEUilo3OtUJ6 In5Y2pQ/JIqSI1ZR7HgipQLIgSKPRAwc/KH2GEiMCulnLwnTP8cUiW1PxJDI OCSqEhvzRAiJ8uvY+dMH+O65c/bCE0GZkECUkRg1RIggCtM7hVdhx57YJhHZ ksN7YsvTdbpn3lB7eBKb/JFZdEqiGU+EJjbmidwGR3oir/yI8QXfExEVAIfi BO+/FDV+dl81GKZ3iqQhIDUUgWybxEVy/ijLKImKEnRIxHbPgtuNA5NoVvok grfBN8bqjeQM3ZD6nsgjEeFL3LOQntgKhZ9z8TwR/lawSKKxYaI2iQhPEIFs m8RFcokpU0mEFUJAIrkGg5NobsLilsRGMGRDMo2GbHkiJ0tGeSt1OwWAJ3aC U0hM2FhIJOSrDtM7xYInEklErovbIBE1UOR9GwqGD+KcMZBYSOiJiIZkXE9R pKqf0Ejk17+7tCTPwvS3ofTiq5JDI3HkI4md5JZJXJztgERcEZAkUhc0bY4T CRkrw3TP4PeD4AuRkD2x963SOE+kkigYByA8MemTiLtX2D11NJ7odO5MJVGc HE5iJxO3nggiUVGC3u06Gom4LDBzI170/kvDzJ3lnZO6ivwJpgMSa08Ujffh FRCUCu6JAhKh7yWJJ+q/laByNU5UfnpASCJwaVucXNtScJNnJyTKS2CaxCIL 6L0FCoiKTw84nbEISAS3ROOcZjD7Ddn0RE6e4O5ZXv+4HYxfEk4FMN1zr6pJ owCg5LA4fPlIIic55BrxSYRPXMUDVUge4vqD73f2vpcHOVCkj9RhcfiSz1go OavCdM+w6okWSZR6Ip3E9gVQeCKZxE6GqCxseuIIZiwmPBF4jcbpiQiYuNvS EH/vjDh5DyQWkngipiHxw6yReWI/tQeeaFBKEgXLieALIUuuTSJuvkMjUTiE gK/ti0mEITJeT3S5sq0iUVXJ5ikUT8QHT4Sp6/jYCvDK0AomOpFEIqf6jfj6 byWgxtI7qzonQEMKk+tbCma+I5p5lgXQrQDGE7kkgrtn4r0FqyRSMlaG6Zwg RgnUEs0z2p5IaUgIx5UnCloC2D0L64/xRF4RKCSWnggeH0GiiMP3XxoXiUBP sGMpmL5dNlA1QKK6DFQS6SN1SBShfCSxFQzRkDwS1cEbnijIlhFJbH5bjywL 0V0ebRKTMofhSXS4sm3bEzUtBbGYZ4tEcBmEcyY4iZT4Nkl0PXcWHKF7Irwh NUi07YmtMuA9kUIixhOhUcQZ9F8a19wZdiHErocjUSu4MHUdn1IBI54IgYTm iXZJJOWsCtM+Thyxt04geKJG8Mo4hCSCCgAkUXIdRQNCAoktT5RnAQ0iyaH/ 0iAkyk0NQqIoZ9skSlIjCiA5AeiJokkyePKs8MRAImg9q3W844mUzg2QWuGJ MBJlfRvYE/nXEDp55t0qrDKgVQAm+YyFlnczjGLPNrQlRKJbinARhDRIY4ZJ VJWAd5Y+iYgCEElU7dl2N2MxSmI7GLghNUl05YnxIpjgNOFIFUQJb+WgVQBZ YvAEXZJF/6XJk8hNrpz78S3FnCfKW9JQBXBFa0k5UtcKDpe0dzYnyyS2k3M9 0RqJUE8kkLgog11P7KWsMggkFsewJPaTQ9bDxLu6KIM04DKOlFTggp6URDUm 8uUHxbvRJoluv8mT2Du1D/c9kUKiahEkqVNLOjdTJGp7ohaJbU+UZIHZGS7M o//SAN9uLH1LoUnsJydYCuwiyyoAKIDcUUCeKLYlICayCqi224JXzyWCzVjI cxYViUSU2ke7nkjq3BCeKMy5V8LeKVASxddRDIM+iVU05RW0TaK7v4KhIFF5 IeiW4oEnikkEDeIII3UTw8SRzFiILdFJ3glmgkRZcFnqKr4jT0QfkhcB6okm holekohNXpwk4E1J4iK5JL6KZkgF5GWgkqgaH02RxNkD960CwrQOO/BEQEMK SFQEr5MLs1bPPaEk2vNEXgGn7onX7boyulcdpnXYgSdqWgrsw3HKyb88D+Ua CeB+IZ1EaQEUA/kxknh45Sy7aWVLGaZ1mIZS5yDHEwkNqbrMSZVan0TFzLMx ZbPqif2EdQZGtvLKZZ7EG/cxduqSTWWY1mESicqGZDZJrFPTSJRnD7jdZtUT pYNl2FBaJfMk7rmHsdO7TyjDNI+qSZT2Dp1jPU8EvGcNeKIkb2aKROF1lN/k IXqi/K2A+F49maQr270nEO25NSXx0jVlmOZRmqkZsBQTniiNb98TFSSqQNEf qZvpnCX3nXXvsOy5hbHHd51VhmkeBY+ThMllwWidm+odPyJPlB3TIRHqie5I ROV46OJVtvfmeSdMKcb9l19FwbH8eJI3lvB45xjjpFdGiIVnyI5V0eQ1KMsv PINaf1UeiTK1Kn/ZFYSUTv6vZENAItMeMF4fXbSzu6Co9kTpcem7rmcHWp4o 69zkhsYrArcA4AoI0jNNT4QMFLnFG4cnEvTcyRcBYZpHaSSq+h7AgqLiHon+ KLMuAOatxEsPeStpFo8xyIrmcL2zw8/2AVtCdCW6hwSeaInERL6EUhXAwFuJ 2fVETrpOAbh5IP8QllDmZyzQMM2jQBJ5Ywl1asiCIs0TVVMCOomwImgdLKS6 POJpG+gWFkDi3nlkJDI4iTxPtEdiQiURwEkNgug6qpbWiZ4o64Htk+iudwaN k4Qkqi8zYIOjYgcDxZCUBcCQKC2BnER5CP2RqqHO2cqMBRamcVDd0gmGxN5Z ABLld/ilzUydcsBIZGpPlCTXIxHjidLcIZKME5193hniaiWJnHwwDSktgb6l QEkULZ0bqgChgGP2xDGSyGAk6niinESpY7okUc8TdUkchycaVZxILUXZd+BI FCaXnCDv3AjznboA4pZUj/ete6J69i4qgKlFnGmRyB0owj1RMZCSRZd6ImRK I/MUBIm2PJFfhKYniuZMoD13EElJNNY7q0mUp5eQCPlcmVUSAZ6p6N3U7aha UFQvratJVJVAQqKJRRz5Ks6YSCzu4vNI7F1GrifKJ8dETwQPL7iHAO1YFUHg ieoFTS0SG9GEFXBDIjn3SlZJhHVuchLpjqNLIuyzwooFRUsktgsgJVGeGqIR kAjrXhEkcjxRfrnU3a+koRNYBSS9G4JEmSdKkyveK2RPlOQNlDMSdQdp5SmJ lERVfCWJAEshVUA6zjJUAWlyLRIhBXBCojkUqSRWnsJZ2uak5g4maSTK9i9A SRS2JIjEaquq4Khy8j5dTzS8F0fyPduYySOIRF5yoqUoxpH6JCaQca6qCMSt OgzkifwYhpYTZXu2DZMoa0hAU1TjrF55uX8Ll5PcHokJaM4hbknYzFPhiXZI bEYTFNTYcqKzHRBkEgULiiaGWUqUVBsk4HMOzcTDe6KMRCOLOKqVbWbgWzzz MAoSlelF+8J4jYD1RMDdPDnH7kgUeaIKBiOeyDvLDYk5g256Z1Am9jyRQiJo RVOUB/grZRRvBiKJ8F4J+LKG5CRWtkgPY4/E/lV07olwErungRuyjML1RMBA VVFGwcG2JwreSk5IZC5IBE1+c0FJ5MUgDbPksMEqwM8D3pCSIoDue6tIBBTA FxJpDZlJSGL3NV4QSQ8EHPALU8MqwB3RYkhkYk+ErMJokNiKNjCJUU0jMYwh Erso4jAgkKgapUGL0HsNRyK+dKBCwj1R/62kzL//UnsvjpnJsz0SeV0eL4j4 goE+wyk+g+aJbGKeqF8BlVT7E5kRS5SSCF8EAJLICyIlUTl3NEMiZw3ehCeC 5t8mPFHwVoK3nzz//kuWVrbtkdjLkRPEJomYOUfnFXhDVl+zJCqAdU8UmrpX JAJzAZLIC6IkUZXcEIlx+xUGb0hhVe2RyIvTe81Q56wk0UTXLL8O5knU8kRp YDmJYFerTqwqgRplCeOA3gp2PVERHKRpkwj+40hGSBQlR7taufEE9FdO1XGA t5kkpQSPtT0hUdKQwGw4JHLXuHhRpF8xRCZRmrp7ZvW9lUVd4KlFd4etkaj0 RIOLOBMksX0LCu5HwltaIBJlq3l4EosmxJLIDwTegSElUa8ABhdxHM1YpF9j jSKRAUgUTDDFlgJZBDFLIlt8oS+4+kISQXnokNi5kHwSDU2dnc2dMxYFR6gk AgswPInlkLjqmsuajJhEZRbGSJTv2e6+QotkyBNjAIloT4QsgkjWlVEk1u+l hSWOlkSeJ/aWoQxZovRzLJ1XSGFcksgPY41EhKs1SWTjJ5GXRdx6wRGJNYtj IbE7eRb4Ed8TSVNPMyQWfy2i3Q8hEvMrC6XBhCe288CtQinL139p8Y3vFYrj JREzyrJFImKglHRNHU1if3oGpYHuiZ22ws64VJn3X6pJrFG0SCKqIYEkojas wBZBZItArkjkLwmCJw3icgoPSDeDomdcqvL1X2rSF5naKSu+QWCeREEJ+Fcc 1rkpSIQVgUwi99amRRIlBWhO/Z2QWH6qih5GSiI4m/aURfgZYn4JaCSKO3dk BXogGiFRM3VVBkESYRaLVRdnJJr6bJ/oOtBIRM08hR/zhE09TZEYy1+QlsEx idws4ioqMrEy8/5L3bHhiEhsdc/Ce8lcTxQOKmE4S4aZVBLBiQf3xLqHT1yS 2HxhrCQiR1l2SMS1hWkS4Ut6BjyxzqMx1nVKookw5kkUt4GGJyoDG/JEEom8 hRjEkp4JT/SXRI2G1NrcJwwF/RKGMXiicxL7WTTvnGMTq/Luv+SWRHxDKjf3 8T2RNswST71xk0cyiZ3TkSRyz4N7Yr2HA5AYK6ckaveNteIKxdoSHQ34LZGI SywgEZiJqKSoMiysoPGCEQ1OIrohi95Btn1A5om83g1IomCYiWsKGon9YBgS BUUVJ+dcyDjurM1PkERxQ+JbEr3LVBhreBIRicUkaqbGl8EHEgU7SZA1iVsr q4L5uLAMXBJB7wXpW0mZeiESiZw/ZW2VRN6FjNt3iZBGIpGPJIoLwbEU6H0a PonYuX83GJ1ERCYmPNESibI922Yl+oJ+HRLj1m/cYIJCUDo3EYn4hpwOiQJP hCXGalgS0ZbSbknkG3IMJLLuO2nMJPLkLYnoqgBIHLMn0klsJQieiA9jiURc Kcgk9ic8WJgAowtJIe7BOqQAAA5ZSURBVHgkgksgKKvL4YVEbknkz1xROQFI FHniWEgkvJW6NbBLosATCaYu0dAkoqtC9cT+OAtMIudM9DiXRmLvvYTajyQo rMvhhUQOSeSu5umQGPefdoKJStHv3JyT2MwE346dcMiPeWJJ5F5I0vBCImck chtNryFVJIpLQejcRkois0kiV16SiF/EYc3Ki9IKPVHUucFI5P/tACKJuLQ2 SJTUX+SJhOGFWAOTqFEVNYmSYhCGWXwS8YWIDVbAMolcxe230nKTGNfPRMFE xfCBxEYaQyQKz+e+GkgsVdde4zKMgUTWfCeh27FPIiIXLInCIminlWh6JDIl iS49UYMmCom9b6ZBk9g7l+iJwNBqOSaxt5pHIxFdjLGRiEzarQH2441IEvki DS/EGphEQgclTir3xFZDmiARlLqhmPJWagXEfmtct7iJaLNedZz7amugPk0S +45kh0R12roAmHtlpknUqkCHRGaVRIGmTiJ3XVmrKjWJkmCiQ1QSOaYOSt0U gcSuJzKcsyIXgQQXkjS8EJet/5KtPdsmSYw1LwNhmGWURPlbSZG4WSQKiZg/ 1sYtgjkSXe7ZNuuJMfrWQHmIQmKvuFp9G4nExCyJTE6i2BO1B7qSsvVfGj2J TEGirByEqSePRMLwQnOQ1UyG9jQ+ibpFMDpMnCyJGm/l/NDESSz2NDXiWyVR dCG9ILEzYLZCoqwcYyEx1u3aGiTiv/rfjCc2SEQnFWuSJDI5iApPbI+zhiJR c3TRLDCepP5qpK4n6k+5EMHskdibumn7u27CVkT4h53L002RSDB1oyS2s0MV QX+gK5JDEjk3TQ1XphlMcChpkYi7a2uaRI2EZj2RKdILL6SXJFoIJisHh0R4 2u7FtzW8UCVdhNcgsVV9zS6WMtAVyUcSJZ7IDJNox9QlZWiUWONPRCFJlHki 4b0EDjZ5EqUFIZHIhiaRNTpF5cJ0X70ugbKqGUiEBRMdo5DYKXAx5ZJFM67G yFbju0xb1QekF1fNPIhGSfzmyWcfW2ds9sB9q9wwE/fEXoEHqUBdZK1hXotE jfTNQoyXxNnVURTtm7Prdl0Z3csLk/DWUDSDyaTyxGZb0El07YllGfT+fiOO REnVzA9LDJJ4+uWb2cPhlbPsppUtTpguiabXRiHqkIi6oOPxxOp7xofzRPMy SOKpl37p/mcZu2Ff+vSSTV6Y7oDZ0nhf6omLIiAXtsfkiQmrNlEN5YnmZZDE h6NXfn/0i2zPPak97j7BC+OIRJlaJGLXQMbjiQmZRMN/nZkuQyTuiaJ7XniC sUd2ru25NSXx0jVeGA6JGrGUAnuiFolx5/dBPLEBIprExhqQh574/F8/953z 8xTBl6zuuYWxx3d1e+ek3rWdsPJf+zeD/5jsaE5iGZ8hS1AOdFt5xdJoNv6V JFbPY0zaPEFde1aYgePyc0qVLEyeSmKmn7p8zvaubB26eJXtvXneIbHkceqe yNqeGA/mieXv2Cu4qH7dOXvmiZmOviz6rmz55seii3Z2FxT7JNpbxJHLHIma nx0gyx2JTmX0HsvXT25nD8+dfJEfJulM3QbwRBqJCeub+hCeWJfBAImKaM7k 8m6fKxLlJWmO2ZGWwCXRuYyQOLqps1sSO8tZtroGpScaJdG9cRBIZAsSAen9 9URHJMpL0nQF7BrIKDzRHYlONTCJFkKpPLHZFiZIHNIT8e/l4Inlw+LSDeUo dBIdvJfkopJYbKDoZDW4fCTRsicuSIxH4InYDMrF+KX3xMYiwjQ9kfVIdK9F WA0Sk5LETk7DazgS7XUNNj2xS+IQxuGMxCXxxIEGKUljpIS1hFahh7rFwhpl MECi2ZIR5COJ8rcyicSkSSIbyBPrMmh0rigSPfZE1iRxoEFK9W7QGajySBxA ZBI1d3xblY8kKjyxcgVDJA7riegMShIZJIMl8kQLkdRFWfRPHngiOi2GRKca gMSia7TYjkBP1Cax/V4axhPj8hGfw7J7YjkyaS7xT9YTq9vWg63GUUmMR+eJ br9nu4jogESrnshapj6MJ5Zl0HonNPeZqzLw0hOrRxckKopigkT817kaFYXE BE6iU/lIon1PbH5N16CeqJFDm0RINDdyTCJrLeZN1RMH3/BcXDvrJDqVjyQC PDHWK0E5OmytCw/kiS5IXAJPHGwnTlWEGP1VJGXaZrkHtBQiic1cxqIhSIwt T1gUb+XCFfRK0Fq7Gc4THZHovSfaJlFVFrMkDqIstN7wJnhi9SSZuid29x8M 5YkxgUTgW8lvT3RAoqosWRm0SSTthDEmbRBZ4xMco1pOHIZEy+vCSk8sp490 EgHRzKoORiAxAZPotSeyqmsc2BNjfUvokDiMDJFotEw0DUMis9o1qN/KJkmc tCcCo7nQICRmj4MOUvQXQdhIPFH/+gVPrJ6Vl9DiZYB7olb2o/BEQhbBE8tn 9klUKy7X4/QTs7FZCljBE2uNwxO1x1ldEoc2DryW2xMbe3JH4ImU9TjS7sBR aISm7n7PdqYReCJh6sm6JA5tHHjF0LfS0FWbPIkAmSHRZIlcapSm7iOJzjwR Gs2czASDtsDQVXNAImENxZAohuCPJw5dkKZ8JBHyVjZH4tDGoaHgiaVi7c11 49AoTJ2i4ImlbJNo+63cJnFo49BQ8MRSlBsco5AHps6CJzL7JFp/K7dIHNo4 NAQdXgxdtcmTaF1Tr0A8RlMfiESrl8GVJzqK1pShYEASh67a5Em0rql74igr MASJtDsc2GAW5IsnOooG03AkWojiSsHULchHEl14Yhw80bB8JNG+giea1yAk Wl5BsP9W9sATIU0wdNVs79lmo1veR2vilji+ge4we7aty5EnOovWkKlgS+yJ FvIcUCNzFLzG5onMTxKDJyoVPNETjc9RkAqe6ERuPNFdtIaCJ1oIM2WNzlCw Gp8lekmig2BeeKKzaCD5SGKQUsETncgt9tOsWvDEoHEoeKILBU9UK3hi0Dg0 Okv0ksTgiVOM5iOJQVOUjyQul3H4Es1HEoOmKB9JXC7j8CXaMHu2g4JaCnu2 pxZtuao2eRKDpigfSVwu4/Almo8kBk1RPpK4XMbhSzQfSQyaonwkcbmMw5do PpIYNEX5SOJyGYcv0XwkMWiK8pHE5TIOX6IZIvGx76Q/Zg/ct7p4UIWxp+Vq Ll+imSHx6O4T6c/rdl0Z3Vs/qMLY03I1ly/RjJB4bRSlJB5eOctuWtkqH5Rh 7Gm5msuXaEZIPP/8ZSmJN+xj7NQlm+WDMgzgkGayKeQ4moKMOUedceIPpyTu uYex07tPlA+UEhCTTSHH0RRkzDmiSNwTRSl5s5zEW1MEL10rHyglICabQo6j KciYc0SR+PxfP/diReItjD2+a7N86IQJClKISGKp16UkHrp4le191Tx/uHmu kUlQUEsaJM5e+8fpzx+LLtq5Wj8EBRGl44nn85/PnXyx8RAURJONu31BQXgF EoPGoUBi0DgUSAwahwKJQeOQeRJ5+8TMa/bA555g0r1phuP96dxZtD//3BPu gs0e+HwjjNVonN2EjXjmSeTtEzOu2dW7fiD60Fy2N82s/ml2I8lNtPdcfNWO W10FO33Zzu+N/oeTaEd39HcTNuIZJ5G7T8y4DqcRDl28JtubZlQPRq/YlO6E M6fHd66xI3UU21U7nNbrxpvnDqLxdhM24xkn8UbePjHjevSXGTsSrd4o2Ztm Uqd2fOnlm8xNtBv2Pf+FZ5mjYOlF/L+z6z8k3eRnSOfOZLsJ29VqxjNOInef mA2duuxVc9neNIOaXf0/z7x0U7oTzpzeE0Wvifa5qhq7Lsq8ykm01/V2Ezbj mSeRt0/MvGY/Hf3bLeneNIN6z+VbR1+65ijae37wLHt055qjqh2Kvvz1n1jZ chGNs5uwGc88ifk+sbOms+1odvXKKuvsTbMWc3Z1lOnEnl9wEe1nfmGe24WT YGzPv0tZ2LHq4kJydhM24xkn0c0+sUM7V9n5ucO9aZknuol2OFpl79+16ahq N/zgJnvEUbTubsJ2PPOrONe72Ce2p3CpWbkpzX7M09+zWYexGy2t09+J/rej YOzoy9Lr+PtOos1e+1nWDdSIZ+Eei5N9YufOZf853Jt2rhHGcrS/bUWxXbWv n/yOo2j5JewEWsQLd/uCxqFAYtA4FEgMGocCiUHjUCAxaBwKJAaNQ4FEizqX /QsfBocpkGhP786X36NXfnnogkxCgUR72hu95oorrvruKPrl4ItqBRLt6cYo 3+/0c1H0e0MXZQIKJNrT3oJE9v5oX/rzsX95xRX/cZXN7r8v23ty5v98fj77 8yu+7zWh6y4VSLSnisTTu1e22Pui6O++JNuU+k+iP0xfe1/072fXRtHfi7IP kQSxQKJNVSSeuewVa7M3R08z9vPRL7FHcod8c3Ti6I5Xb7Gju61/5mciCiTa U0UiuzZ98q9+N31yZMcH52cu27XGTu1+1daj0QdfTF/6fJjO5Aok2lNF4uza V6RDw2/ef9UPRJkfvjvtnt+XdspnLouiH/r8+sCFHI0CifZUkXj0JSubsx/N VhavjPbN2ZHoQ+y6i1cZO/UT6cAxenXonXMFEu2pInFv9MH54ejVz2YfPk9J nF298qXLPjjPd47+7Z98b+Tgc5BTUCDRnsr1xEdekj4ejv5D/sq+/GcU/a/s IdtN//7saVAg0ab2Rj901VVXfV/+OZHD0a4vP3b/jujVZzNnjNJJS/rSyrPn /uIfROHLoXMFEu2p+NhXFP1B+nz25uzZ7/796PK0e/7RdKTIsu8oz/RLYe6c K5BoT+fPpXr+b7aL3x77whNz9sIXsrnyu8ux4Tfv/1yYO1cKJLrX6ZddHubL PQUSXevMP/zu/H5fUFuBRNc6vTsMDXkKJAaNQ4HEoHEokBg0DgUSg8ahQGLQ OPT/AenPnzOu5YyDAAAAAElFTkSuQmCC --============_-1068967708==_mr============--