Received: from cluster-jaune.cls.fr (cl1-jaune.cls.fr [91.207.20.90]) by sympa.cls.fr (8.13.1/8.13.1) with ESMTP id o2402mrp026289 for ; Thu, 4 Mar 2010 00:02:49 GMT Received: by cluster-jaune.cls.fr (Postfix) id 886B09D0DE; Thu, 4 Mar 2010 00:02:48 +0000 (GMT) Delivered-To: ids.analysis.forum@cls.fr Received: from cluster-e.mailcontrol.com (cluster-e.mailcontrol.com [85.115.58.190]) by cluster-jaune.cls.fr (Postfix) with ESMTP id 16B1477B7D for ; Thu, 4 Mar 2010 00:02:47 +0000 (GMT) Received: from sol.csr.utexas.edu (sol.csr.utexas.edu [129.116.190.1]) by rly35e.srv.mailcontrol.com (MailControl) with ESMTP id o2402XOX019566 (version=TLSv1/SSLv3 cipher=DHE-RSA-AES256-SHA bits=256 verify=NO) for ; Thu, 4 Mar 2010 00:02:40 GMT Received: from localhost (sol.csr.utexas.edu [127.0.0.1]) by sol.csr.utexas.edu (8.13.8/8.13.8) with ESMTP id o2402SOP031364; Wed, 3 Mar 2010 18:02:28 -0600 X-Virus-Scanned: amavisd-new at csr.utexas.edu 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 Zc9uf8eW8Uf9; Wed, 3 Mar 2010 18:02:24 -0600 (CST) Received: from [129.116.190.133] (ries.csr.utexas.edu [129.116.190.133]) (authenticated bits=0) by sol.csr.utexas.edu (8.13.8/8.13.8) with ESMTP id o2402DOr031344 (version=TLSv1/SSLv3 cipher=RC4-SHA bits=128 verify=NOT); Wed, 3 Mar 2010 18:02:13 -0600 Mime-Version: 1.0 X-Sender: ries@sol.csr.utexas.edu Message-Id: In-Reply-To: <47C8EBDF-9F35-44F1-AAE5-D14137F4C495@noaa.gov> References: <34F95ED2-4510-4985-A36C-788FFFB363E5@tudelft.nl> <6BCEB074-4D3A-47CE-8B9C-401017E385EE@noaa.gov> <7754E570-027D-4F68-9C9D-319FC68E464E@tudelft.nl> <929FA1B0BAB6F444A31C47865D41F6E80580C1A7@cst-xch-002.cnesnet.ad.cnes.fr> <4B8E64E8.3080107@cnes.fr> <47C8EBDF-9F35-44F1-AAE5-D14137F4C495@noaa.gov> Date: Wed, 3 Mar 2010 18:01:55 -0600 To: Remko Scharroo , Nikita Zelensky From: "John C. Ries" Cc: "jean-paul.berthias@cnes.fr" , Cerri Luca , Eelco Doornbos , Marek Ziebart , Frank Lemoine , Scott Luthcke , Michiel Otten , Eric Leuliette , "Mercier Flavien (contact)" , "ids.analysis.forum@cls.fr" Content-Type: multipart/Related; boundary="============_-944461573==_mr============" ; type="text/html" X-Mailcontrol-Inbound: c!PzNVdyE6pAgYrdZN4HGLIOLgW8bYa616COtuSiUsY= X-Spam-Score: 1.7 X-Scanned-By: MailControl A-09-22-10 (www.mailcontrol.com) on 10.69.0.145 X-Validation-by: ries@csr.utexas.edu Subject: [ids.analysis.forum] Re: TOPEX and Jason surface forces, again --============_-944461573==_mr============ Content-Type: multipart/alternative; boundary="============_-944461573==_ma============" --============_-944461573==_ma============ Content-Type: text/plain; charset="us-ascii" ; format="flowed" Hi Everyone, I doubt if I have been able to understand all the information supplied so far in the extensive email thread, but perhaps I can get some clarification on a few points. 1) The SRP signal is beta-prime dependent, so a full cycle is 120 days, at least in Z. It's almost a perfect scaling of the beta prime angle, including the sign. There shouldn't be a frequency at half that...certainly nothing in my test below. So if 60 days is dominant, perhaps we should be looking something other than SRP. Because of the yaw flip, could something associated with the radiators exhibit a half-beta-prime dependency? 2) There was also the mention of S2, whose alias is 58.6 days. How much did T/P influence the estimation of S2 in the tide models being used? Could S2 have absorbed a T/P perturbation that is now showing up in Jason 1 and 2 (and the same in both)? A tide modeling error will show up 100% in the sea level time series but will be diluted in the orbit differences. For what it's worth, one of the most prominent error aliases in the GRACE estimates for C20 is S2. JR --============_-944461573==_ma============ Content-Type: text/html; charset="us-ascii" Re: TOPEX and Jason surface forces, again
Hi Everyone,

I doubt if I have been able to understand all the information supplied so far in the extensive email thread, but perhaps I can get some clarification on a few points.

1) The SRP signal is beta-prime dependent, so a full cycle is 120 days, at least in Z. It's almost a perfect scaling of the beta prime angle, including the sign. There shouldn't be a frequency at half that...certainly nothing in my test below. So if 60 days is dominant, perhaps we should be looking something other than SRP. Because of the yaw flip, could something associated with the radiators exhibit a half-beta-prime dependency?

2) There was also the mention of S2, whose alias is 58.6 days. How much did T/P influence the estimation of S2 in the tide models being used? Could S2 have absorbed a T/P perturbation that is now showing up in Jason 1 and 2 (and the same in both)? A tide modeling error will show up 100% in the sea level time series but will be diluted in the orbit differences. For what it's worth, one of the most prominent error aliases in the GRACE estimates for C20 is S2.

JR






 
--============_-944461573==_ma============-- --============_-944461573==_mr============ Content-Id: 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 iVBORw0KGgoAAAANSUhEUgAAAkgAAAGNCAMAAADdDX4HAAADAFBMVEUAAAAZ AAAAABoiAAA1AAEmHQAwIw0AACYAADc3ADcAKDItLTZNAQBoAwFbMQl2JgBM MTE5QjJESTtvTi0AAEkAAFULL1U2O1UAAGYAG2YAAHMAOXg6U3ZcW1tsbFV4 Vnxram2JEgCIJACQKQCLNQCTMwDdCAbdCBjhCwbfLxHpLQbdDifhCzjoNDqO SQCWTACfVgCpWgCvZAKRXiaNajm5ezboTQvtbRjiRivucyvdCVXgCEjdK1bn OErdCG/gGnfdOXDhJmmNdluye0eKZ2avfGrhXlvCe0TpdFPeQ2zlXXbkb23H gybLizjSkTz0ijS6jFGXjXa5kXrOkUPyllHao1jJmGTap2DJpHjWqnXrtnjq wX4AFIcAMI1AN48AUpomWZwvY5MAUqgAZ6wAarIAdb4Sc7svd7NTcYVodZgB ANUAANQ0ANYALdozMNUACeEyDeEAJeM+LuBLBdZoAdRNJthnMNhJNuIaQN44 R9kdfMAreMMCTOYoVugMa+sydulORtZrTddFdtFrbNxLaefdD4XdLZDhPo6B f4HdToffYpbgbJrdUKPecLHhZ6WHCtSPLNSVS9SkUdSSc9qwbNbTfs5/hIZx k7Vqrb4Sgccnhcgyicg8kssNgPAuh/BOk8tpnctYotV5rcpiqtt2qNdqsd9Q ju5olOxZpPRwtPJ7wNh9w/KMjJC/m4e5poaMm6mhnqWNqLyrq6rplZbasYrG pZHMuJbTvJDuqZTnmbHOuKTvt7fsxob5zIn00Y/nx5r6zZL205bXy7bpyqP4 2KnkybHh2rz03LP75rizhtuRtc+ptMmUleSrleKUsPK0s+vLjNfinc/YqNzl rtXRseTjuuGPwd2rzteJyfmDyOqSxeyc0O+FyfiWzveU0PGp0PehyOWn2e24 1umk2POz2/O25fzMx8jYzcnIzNfM0trY0tDt1tHW6Nzt4Mn96sP/88jo4tb3 5tj++NfP0+/v0OfH6PnW6PfI8f7b+v7p6en66Ov+++nq6vn16vXr+v7+/v5N fMOcAAAgAElEQVR4nO2dbbQdVXnHjyY1oIlQDUJTXyINYxZQMFSppIVlre1S MYgNqYilqZqCzopYmmgrUqhFsa5WCIKlK7iIQhVoUlsJXWCBpby0XhElJljv t/spgqbrXrm5ISEnubvzcmbOvOyX53n23jOzz9l/WCfnzMzezzOzf/Pfe/bM ObfHvLwMqNd2Al6jIQ+SlxF5kLyMyIPUunpuSLUXjRwrL4ncaAIPUuflRhN4 kDovN5rAg9R5GWiCo5x3huVB6rzATbArjPWN4efb49e5j0cr2NyBZNGea6qF tkyyp+Jtsk/520j9jUk122uRnr7iUz+KoPzMFXexuathWXqQWhe8CY4ciZiZ yj+mILH9bMt29uWUhs17q2W2TPb/It4m+5S/jbUpqWayWmYunNodTrIt1+7f dA/70iQoSw9S68I0wVx4T/Lv98PwLrY7/NvIn+Y+GznVXBjG7b07so9PXhde +3QY8fbJg2z3NyJ0bgqvnvssu/nJjTEbk9HbozeFf72XHd4UXhfGde36u01R ZZv2DrhiKaGxd82z3R9nP7kalKUHqXWlTRAqlG67OTWhufDA0U17I3eaCw9G fdWWe9gtCWARAf1N33hh89XP33JNf+PB+HMExF8eire5Zn53eCjp2rZcG71l m+86cnMK0pVRNZNRHXuuzHOauzk8NHflfNwnzoWFLFV74dWiECDdFLUtu/m6 uw6H1z51JDGOGyfnrow7rHSgE9HAIn6if6Lmz0Dq/9V8ug3bPBmBdOXRTXc9 9VR4NDzEDqdjpKjQl+6OKIpLs7nPXBcNj57+TNS7xSCFEZmFLFV74dWi+E2Q m1BBu6LGZ+yHT/6IHf3hJ8KpQQ+UgpSMZKLB0uEYpEk2d2W/AlLUdW2Zijfv b7rrySd/tC+uqwAS2zi16WD06fCTT049HY3cd10TFUpG6RslWar2wqtBgZvg J2E2LJ77m8hetmdDmQSkpGu75fbUkSYTR9rLbopBYhvn423uZofDg9HmIds0 GQ2Q2MbJqMK40K5oEBQtuiXMr/huid7ddPtcuJd9/5qBbXmQui9wE2xJ+rhk lLT5Y58ID+yK2vvG7REkN25nuxPI9lwdXdLvjT9HS78UXnHF7fH7zVdHn24K PxVem26+O7wiGrPviobqKUjRx2viXiy/epvb+LEwomhXuDE8mE0deJA6L9Tl f/R/ov3755N/j0ZL03XJ56h3GmwQfdx/JN0gWjTPbto7N598il+TNXP5RMCR +N1c0m+mOro/WTcXB/ny3aAsPUity2QTDKaWOLqxNsNU1tMb7+YuPxymxHqQ Oi+jTbBfuGJeUfAp/vKjz6f/epA6LzeawIPUebnRBB6kzsuNJjAM0uzropfH jj9xMnv10tY4gvTom5YyNr3wjvt6UxPJKzkvr1xjCNLs218bgbTuMsZWbU1f NTLzGmgMQWJs5ljG1kT4rL0seb2KnJdXrrEF6azbIoQ+HL+u2UDOyyvX2IIU e9EFW9NXcl5eucYWpHXLDvVX3Jm+FtdlD874/+T/scpnV0DKMy48aFdcj1IM Uv+MRQvWz6evxXX152eMy36IFnbCFZBK0gUp1TPPDl8ldXuBNMYg8eUdiRYC 3ATfDSL9Q/bpQe42D3+OlpRSowXSSArcBBEkR6dX7hl8Op27TX/OREocjRZI 4+1ID18fvZy9hz2+MvKli4LV7PG3BO9MwJl51zuCtx6dPff86x+5of+O81e+ 9cGVEXDTZwfvMsbVaIE0kkqbIFAp6tpOO//8t6xm08Ge/jlfnD13jp36Y3ZR DBebCX7c/+3rZ1a+cy5yrVO/8Nw5q+cuOa9/6hfZJeeZzTKX4yCNsCNBQFr9 4IMPnH3Dd0/b8+BF57Fzj7D+vp1nf6F//vnfnvkdxiZWz5x+JHKt/qlz7D03 sJnTZ1Z++8EH3mAq69ECaSQFBinp2h4577un7dz5wJ7+uax/6u8+eNH1szt3 7pmJjGfivOnV8TgqAWlXBNJ0sHPnzm+bzTKX4yCNsCMBlIB0zj9OnMbYAzew 3zwyExzJuraVc+ySzyeO9HkWg/TjCKTZaOHE75nKerRAGknBQYqNafUcO2fl ytPm+mevjt6c/odBPJyeWXnaW4I90xFhH/hcP+LnPV9mE6ezS95warBHWSst S8dBGm9HYof7Rw/H/+7bl3yK3hxlyYKZz7N4UfrjSIfTjaMPz+0z93NJowXS SMpAE0wbuzgTarRAGnNHEqq/T78OhUYLpJGUv9dWkXckWggPUkW8p1S81HIP JAPPI8nkHYkWou2f9AdKcZzcAsmrK3IcpJF0JBcjuA6SV1fkOEijeTI7GMF1 kLy6IsdBcu1kTp7+sBuCJ+9II6aAC9JIyHGQHDuZ+SA5thPAEG6B5JiC0bUk x0Fy7GTmg+TYTgBDuAWSY/KOZKZu8wqtj1+9IxFDuAWSY+e4W9mi5DhIofWm 8Y5EDOHY80hOneP5F9JGTSPwPJJTjiQAyTsStm4LcuoUj1N1KF2UHAfJfmdh 2JF4IHlHwtZtXm6NOrwjGarbvOw3jXckYggPkj25lS1OHQGJ2j855UjBINtq ut6RsHWLRR3pOHWOi0AaCXmQFDLrSIz3kKR3JGzdQpGvvZw6x0UgjYQ8SAp5 RyKGaBEkAg5JIVfOce9IpuoWigpSYB8k70jEEG2BRMKhAZDEYYmlvCPp1y2U BkisBUcyDJJ3JGzdQnUYJH5UerLekQh1Qx9sI8/U5SDZaxvOHpgGyX1H6sqD bV0GiR/VO1JVXQGJkY5wA21T34mAeI05TLZS1H1H4oVwCaSglZPcOEgjIQ+S QlxH0rHPelHvSNi6RdIAqfRvMzIOkmU1Y4HdAgm5w000jfOORLz5hFUXQAq6 DJIgKqF1WnIk4jUmVh0CCX+E23UkYyA14EhjA1LpX7BadCRSst6RDNUtULSf ycwodn8baRqBI+mBVC5q1ZGIcxVojQpIzZ3kAR2kUgUNiTrrhZXjIKW3eGwe qNpOaINUK2rbkRpBtwMgBVogZSShs6NK98qg4UESdbICq46AxOgghczqgeI5 UuEfuNpxpGC8QGJUkJInGLwjSaOOAEiw55F0QAppJRGSOBIB+8o7QQiDygd0 VknqxvNIRZBQu5tu3rAjBUSQgtYcidkHiXWua8O2TfOOFKRR0TElINl3pHAc QApzkAhtw2gIaqg454Asx3lrXdQRHVrdAQltLGHmSFZ/bLuyE2E+52AOJIuO FIydIxF6qOzoNDpICrO5CjccKRg7RyLgkI1W7ILEcaTcPpHZFt6WStp1pCSC B0mioAVHCgPqpYEEJIvKIo0DSEEJJMT+Jt7AmnWkHCQK9tz3tRBGNU6OpAUS Y807EinmcJDOGhwkBQVHshyyEyAlG6MHsC04UqgLUlpbY44UFBzJBZD+aPny 5VPsseNPnJTXzVXeNMi2CSuO1MRJrgfS4Gxp1pFq7yzJBEgzi2+99z8OTC+8 477elLRubgLUHiqsOJKlA2XQkQZnS7OONIjgBEgTfxq/rruMsVVbpXVzE9AC KSvYzEkeVkGCxyyMBL0jifRQr9f74KE1EURrr5LWzU2AClIwdCSbg6TiThRB wk2nhzlIA+qLJe06UhzREUd6/KPPTx/3lbNuY2zNBmnd3AS0QCKVJGsIEjZm oVyTtwbz2VP7IU1dta3dEDvSBeWuLRWT/jcwljA7yeVbF/7LDhK+JOq/sPQ+ JGZbKmcx2/oxCvOIgbU4oaHnkd69gfVXfGXdskP9FXeWQIIUDjQcKYvTnCNR J6/Cgn02OO1VmHNwwZEeX/zShb92oH/GogXr56V181Q6wuimCdOS+OME/iOp 5TESFftQhr3VMVJmFy6AxNi+n8evzzyrqJun0hFG9OTSpgEUp7SeJZCsKenP wmYua1uf2dYDiepI8D/bXHOk7C0qZtCKI+UgOeNIwLp521BPVWnTqMNS/gB4 aAkka0p7tmauFB0HScORgCAVNisO4vAgFUJXWtWmI4X25/5TuQxSXgf+OBEd iQpS6TEO7KPBpOFcoiJItl2wbZDInj+YKU4LYp/a5k58CDYtlqqBBIxZBYlh HMkESPa7006BhLj+lzeNsnRIaZ7CRRsSXq1sKeaZKgOpiXFZB0Aabo9om7Dg SESQYLFKjkT1T01HIoIUeJAAdWs7UhBgx59hMahBkFRRcTO1QyUgZXWMGUjg 3dVypPx6GFCm6kiUZMvnB9KRBnf4aCCFOUi2H5J0HCR8wcL2WEuyBZIqagoB noIhSGzkQSLzoOlIcJBMOVItOjcEP1ltkOjdI1iOg1QqiDjJB5vrOxKsfAUk 3N3pENoLc6IOHcn2aNsqSOrrDaOOBJ4eLIKkKqRwJAJIrAqS0pHwvXAWddiz WQWp/d9HqoAEPlXLPT7yOA22xk9jBhWQMNmW41McSQOkBi7bugQS4vo/fwY6 L4gDKRjYPsaRQrMgFc4DQLYUDAqOhISXIIdBKn7UAAlhDdV+2BBIkGwpN10L jjTyIAVVkIC7m3+bA1uQDW0eBlLFkUq1QGMGcuylxym7cNcCafQdSRckdMHC xth74mHNPg2BpMo2sU88B5Vh5MiDVCzQkCNhQBI5EhIkXgL1EPxsKRcGrHw/ fNRBIhsLvyBw5iAHCTvFzAMJUlwBkipbOkjFs4DSOSLkOEhlR8KBBAunciQS SCHGkdILA7yjFK/+PUgCcUYd8LnM7JjqOxI8W04GiKgBBQRqR0yT4yAVHIkM kryU1JG0QBpOScmzzXLVACnUGNEB1TGQgKdq9afsMMOdYRD0nRUiSLUf3sNQ MQQJ28KVbMEu6CJI1cOJBml4CwABEiuApCxlwJF4IJX9QpruACTsIKnieeDy xKFUF0FS74eqaaRlCyDhJsQDSyDJsx3+6i+uievDyNEGqVwCCxLNkfIRFqBU wZGqIEFj6jhSsR9GgzSsBVGeOCRvFSTyqarpSBiQCuXq9tmEI5UHdLBcYxHN frxAyg8J2ZGKpaTFCrc9qSDVQtAcCWlJdUcCZwuOUZBVkFQPtumDhC3I0CAV yhkECf7gVd6h6oCEzBYcY1h5vaVdcaRyCCpI6mJKR0JkW8khWyp3JFI/zEbL kZQbaINUKQg5BEU3wDSOPZBgUbGzXk06EpOB1EtFqVRUd20Dg46EaBoUSENH qtYPvf6v7WU1rtKRBu9QfVtpU4Qj4S4NC4lWFwzIyQnSQIkCEuhgKZtGVpYD EqxLtAaSPF2SffKm6GAcGgapV1tiou7aBrwuCgVSYQQDBomhQCo5Uq0m247E +W0vYDOXNgwZuDslciQAqVdfZKLuquo7BuNBz5EKYVH36KggcZoPkW1AA4l7 jgKzBQWoBawuaHCwzfkWMRYkkiMVL5iUxfQdSQmS1JFo1Fe/VVGLKcsWFKAW sLqgYZCqixpyJHQxwYAOYhGcLXCOVClFy7Y9kHral20NgER1JHgxsSNBG5YL UohwJGQ8TrZtOpIBZ2rBkdQHgQsS6PqJCxIgW94W0KntGkhk/xxPkBT7Uhg3 4Aqy8qgDclGcz1TxB3QmQII7ErUjHltHAoNUKAg8yJUfCkLdwieBxNlLoiPp +2crIOmTRAIJsLsh72RGgVT8qChWdKRaVUwPpKqtqrIlg5Q6EgjeERlsg3gQ NQ1odpD6NUeLIIGzhQ+SajcG85iqkiPStWFBQjtS9Wi250j1G4a8bCv9MGna YBABBC/0qrAesbrAJEjy55E4o44mHKkKEvGLK2wIks1sAxz2pYL1mLZAkj2P ZN+RTICk7UhKw5c6EmTMIgApBDtSOVnSpUHmSOrCVEOSjZFoFcrqrqymnqrc tgMORGu/7qrhSLC+jZ8ULGwdJES2gitieUHzIPWsD7ZlIEn3prC+5kjGQZI4 kimQ4I4EHaTXvnmaZduGI+nLIkic2mBXNKaaBhyTu5dUR9LOdrxAUu1NYXXR kWB3V+pNwx2rDddLstUGKSiF4BeuggSbLfOOhAOpUpIAEtUboMkCQJIHrWJP ASnPltYPQ9TihKTwckZ1P57rSDCQqvORTAlSOAxqGKS8pPA4ceilgTRwJDW8 FkBq5vKfsxQFEqog47YD0JEsgiTJ1hBIxbIEs4eoPZC4nQUWJLwjYUEqOBKn NnVMQdPZdaTatEEtpjhbVd2ikNUFzYEkvC5W7A6gacRBRSApx1bEZMXZAojg giSwuPJW3KgtgWR7QtIISEYcSXrZJnMkyGSFaG0hE4wjiYy8slElqTYdyfpg mzx4EKyFjzoovQUVJGGzgxyJUzcYJG5I2vATona7Nt5iFEh4R8KCpHQkxcS4 HKSgEIJfmAsSFnuoI40gSJL9EToSYLBDfk7HJkiYZGkg5SEVyY4VSOVHvcol LYAU5kGNg5SVFB0n0SQo3ZFU9xVtgGTi7r/0eSQhSIrzBtI0kpRojmQVJEyy sMs24g0+MkjS55FaGmwrj7DMkdRNo+FI/PoUMYUeSXUkUGvzHUnp2XYcSV9S kIRJ23Qk7gAD6kiSmOLSEJCk2QKXgqIqQNL6S9x8kBr4EQkzIKEdCQ2S1JGU IIkbvTBgMe9I1ZSGjiTdVfMgNfCzNsKkrTtSfaG6aWQgSdtGvE5thBKQFOnS Dq1Oz9beD22pQBIe/9I6iiPxlwrj5VFFyWqDRHAkxSVFdQOEI0krlsWsLmjq p/9UxgIDCVEw3cRU0xRi2nIkgVMChnTibGVlYbfxhDGrCxobbMt4kO1uaVXZ kZQTSYSmac+RBHxTQMoiKEDSuNXWSZCkBixvGuVJTmsaST+s6C20HEkAktI4 aDf4Rg0kBQ9iRwLeS6I4Erk/VYIUMOFxEgFDAKnsSOIT2EGQJEfDniOROwtl slogybJFLB9Khj1pHUAGQXrs+BMn5XUX1xFB4n2pHVRQEhTmSOJkwdnWStId SdrgtfUFR1I8INUBkKYX3nFfb0pad3EdFSTZ6U8GSdU0ZJCk2RLvfIFAElUJ etKOIuXdfzBY6y5jbNVWad3FdWZAMuJI0qbRcyRZk+ctxz9OwkG1usVrq4eO JCyrOUQyCNKaCKK1V0nrLq6TNJ0KJGGdyidQqE2jA5IkG/kMvugQEUAClAU9 FyiR6BZJD333/6zbIpg2SOsurjMDUsWRiHe3ZU2j7UiSbBSOZA6kPELjIBHG SrEjXVDu2lIx7n9R0oI1YTZm4a4Tl0tAkq8P+GsHh1dYTp6tbK0i2zQqE2cr OD7YbNmwLBMfA2m2sv9C4fNIPQJI65Yd6q+4swSSZGtFFyUyFv4XtgoFlfd7 KY4Eu/hSZ8uLGkgdiVAt79hWHEl8DEw7EgWk/hmLXrR+Xlp3cZ16rMNfpWga 0k1JnScHDYCEDYoHqRhTkK0dkChjJMaeeVZRd2EVFaTKzqIcSQGSoNUU2UpB UozYBiVJjiRp8vraoiPxs9W9aBOBZPvhfwhI9W4X3DTIoMpBumSt7FLRTrYE kApFBbzojrXFIBmQDCS5qUtBEoeggyQ+HUNFthKQFF1QHpV7nCS0qMyjvnIY wYM0WKNoGvI5rhykiyuVDDsUYxn5vXiDIKnLWgKpZ/tbJIpRaHaRWl+jcCT5 49cUkDJHkmSrC5LYkUTl5BWrHale2GVHkqwdgFSrQOsct+NIBkBCJqusmHKp qHvR1lmQ+H2bwpFsgER3JOUTaFnjYR0JD1LJkbjZ2gRJt2sLQuI5jgCpXlB6 xUJzJNV0kHDYoWgaaePZcyTeeu2rfwlI2tbUBEhcR5I1DcWR1EN8QW+hBim/ ScIvLSwoq5gTtehI3F3VHiLZnEeKQCKeUtl9H6IjGQWJKS4xFSBJkpXPOsin HJAg1bLlLHUTJOXOckGqd051R5I91UW5Dgoh2Qp6C1XTFG7bcpMlg1QrUItZ z6S7IDFdkGrX/+DOQrBS0TSybKVBRW0DBAmbrOI4UIaRFkEy8bM2Gk0DBKnm SLL7+CSQQmW24rbRdSRxOSRIrTlSj3TTtlY3FaQMiGrfBu4sBOtkIMku9yyD hE1WsbZrjqQv8Q9tGQOp7lmSpiGBRHUk9TcK6I6EBqnkSJyymL/fzc+n3tKF 7/4nzkStm3XNkRRzTDQny7KtMwOZ4ROzrRipazgSt19MDrTNwTblAbdC3eK2 UZ/jOo4kaBvpdZDYkSD4ikCSFst2iPeUA+CLAaKVckcSgqTEXi4pSD3WKZDA F9QmQVJeYgpKI0BC5iMrF4vQEbsKkrppGB+kWjmOIwmtRX4dJHEkCkiwew78 i1NpPoVy/HW8NWVHEoKkyFYqx0Gq1SwhQvFUkbBupQ3yrBLWNNSvHGFBUlVt GSTOz0ni6qY3DROCVKmpVjUNJEmTQrKlg8REjiTlQdbJ81aUHYlzGC2ClF6x 6V616YFUIUk9RNIBSXy5RwAJ/IdRZU+gSEojQVKU1b36T+qoLjA5j2QfJLgj KS+oReWIIDEYSCHXkZQDLA1H4t0f0L/67zBI9b5Ny5GUow6akwm2gTaNoHLT INXL1rLVvGiTg6QJlWGQuJ17rWrR9T8VJPhf2isXgjVNO45U3WbsQALsrHGQ QI5U+7WdgSF11JEEIClyVdRaXWAOJMkzPsZA4jmSaPxq25GCLD5i0ME/OdST UHqOxOuIXQQJ1DQCkNRBaU9Qi6fEgWHTjdK7l6E2SOp7+HRH4o/orIKkJ8ls GxSkEknc29acyokgCdYDf8a8AFI4sCUYSLza1Td8dRypVtjE1b/weaTyZ1rV DYDEq5x/KGGPe/NWgIahOUgs/7Eg0E4SQZKMDvCOZGKsLfymrYEHJdNfzOKu gGVt0pEAow6Sk+WbJRtlHZsRkJTxuFtwl6scSRckyA9t6dwiEQADa5oqSNA+ xgZI6lzTfSoMtYFN0wxInJiVa0xdQxJ/9z//M0kadRsGibO3DTkSLFlTIEGe HBCmxd+JysdSbkaGSBKQBgMlSyBBygNA4gVtB6RkZrFo8MBiXJDUpXEgyQrD LzHldVYXFOnRfYxENNsGBykgOBK3euVQRzTAxYFUWIAAqTo5ZhwkyZOmmAGd NKHqgpINaT5qy288MkiIzoLTW1gGqTbtBSzWCEjCwvko2TJIxcs3Qt2aIJXa ho8C9FEedb8oaj3rINGeDjLjSMOrLUsgFT5qgaTXNBWQwJ0FrwsEgKSZbRVf uyAJB5ooRxqSbx0krbolwxVQBQWQRN+L5kUlgSTYBNqhGgQJdhUlyIu/tPZo +2Azx0GCNk0ZJPjwtR2QqqU1QILMHaBAEhQeXmWCd1JcZXVBx0AajLZFf2aE 60j8K2oNkEDJVkpDd5IMEj8xfrmaI2UT8dhsJflUF1gHCdc0wXCOj940wCm+ boGELsdJQl44NNizdRuktG/Lnsvguhs3LG/4qmxY8yDBijUBUv2xrRQkRTGM 7INU72VQILH8uQxoMQ8SQGMIEsuey0A5Uq1tAIMOp0ASXP8DHakwapAUw8gy SNwLKCRIw7tBCEdqGyT46LWaLfQeqrYjBSWOOg4S9wIKCxITg8R3pGoE0Bfx BT9OAz/GpS11QAJ+b4C3kaAg5zAZBUn6+0j64v8iDKJpSjsIL0V8Jtk0SLBS tWxtgMRR1T9hpSTpVBfYBgnTNGqQeI5EB6m+EcY+mVGQ0OVkYfmOpCyGUTsg gato1pFaAan+5UrzIHFE9E9xOtUFXQWJ72NddSTMoKM1RyJlK0ynusA8SNUL KJMg8eMKxq/KYvogEbJtyZGo2AvTqS7oHEhB/o4bgrewElXjN2Y0QIKWqoAE f4Jaz5FoowZJNtUFhkHiTLeZBEkQuNo2oJNOGyRGB2m4NfyuolFHgpYSZ1Nd YPbvtemCxBQgcR1JC6TK4Go0QHLekXjjXqMgCQJzQVKXqnXE4McyKzmiBh3l bBEPXtU3w8R1HiRc2gYdSRmrnlsjIDH7II2kI9FA0jjHGwUpzxZcyDBI4KjE SwNhNtUF3QMpKP5bD8Fbag4kFMDNg1TZTlRQ4Egk/xRlU11gHSRc2gWQEIFb Aonpg4T4/jQCJK7GC6TB7mI6//rwVQ8kdLIUkLLtEY/w1U4Xcbbcw0S6jhFn U13QBEiYStKDjKSP11k0BJIMe1FcGkj1mUz0BTEFe1Ey1QWGn0eqTzJTQeKH EER2C6TQDEjIw6QYfmJk/XmkWkeOHWuz9OBg97UUFX7PoZocCXuGP8VNgET4 cRpzIDHrXRtn3EsAKRDuq8SRjICEzZbSD9ccCVq+5rtI46Zcx4iTqS7oKki4 yGF7IOH9kxV+ItEOSHwNO2JEIWEy1QWdA4nJQBI5EhGkWkeMZSIDCVOGFX5r E9VB1XwX70gOg4TvkRtzJM7zJ+jIFP9kRZCYBkgE7BHhFMlUF3QPJNnxAToS MGjVL9FjbWJHPIxK/yaoFCTBYSKN6ETJVBdYBQl/GUSM3BpIjGBIzABIpF8U dQmkmv+aBQngSLC/LjvYtBWQhlFxPwwKn38SOxKJe1AIt0AShg5KJ3mDIJGU ZUv4SnH+npEcySGQ6r2MuQiiU80MSI11xBk7yO+mwycyRYfJHEdGQPre8uXL b2Xs4pe9+kC9brsgiVSIqgNSQ9lmCab3GeAxdR2payCtXf/Nb/2APXTMTy88 6VCt7pYcSQOkxrHPEkT/5HXFd/GOZGqkzQtBAWkyfl11D5tdMlmruwKSucyl MgRSQ9kmCWZ/nYsEEvhvn1qTAZD6p/R6i55gJ09G77bX6rYLEsyRoEErzDUF UhqV8Nvpw71UDK6EjmROuiDdv3z51B880b946f8t2R6B9PVa3ZUr8YZOGipI YQ0kO/lVFRT+kIkOSFaSg+VSXYAE6bEdO+IR9uzivZEjsVMEXVupTYmZciV1 pPwYwy+D6iA1czIHQf4DaggXBIPkgCNFmlmwlz209OCqDWx68VSp7nD4aFuY fu+W5e/t/leOmo5AwOWGn5rJdjDYCdNPASXbcFBSHMHmf/kfzdQDif1+71d7 X40oemlvW2l5OMCpFeYZDGUAAAkOSURBVEfKoiImtlt0JDaMhDhCpZ1035Gi 4fYzz8evP9vPq9sqSEIVjjGibUoJNjeiY8UUaSAhS5qX9ZntdhxJC6RSts05 0uCTDZBccSRp3SWQGpvqMAaSrQQrIoHESiC1Oo3kOkjmHYk55EhgkEbMkZoc dBgDyVaCFQ0p0AHJRmbQVKoLrIJkfGdBjgS3wVYdiQhSfpNOUnDkHKm5s8YI SE2OOvRAwhY0L/sgsS44ErS64tYpSE2dzFlYFLyDvVSCNBKOZBUkceyC66Mu qIct2eQ0kmWQGlBDIOHNARWCt2bo+jSQwuwOiXVVHQlRcuwcyRZIkuAhESTW RrZaIKEfiLOhJkAadBHIXgYTgrfGEEgOOJL6OaYRcaQcpCbPmiFImFGHo46E L2lcjoMEcSQ9kBp0pIAhx9osf/5kUIM0gl1ZBSntvAsgNTgxUwIJUWr4jY5G p5H0QWo022oa9n9oi9kFSeFIARmkMMu2sZM5BwlVdPjIlOzYuu5I2b8BxRx0 g2uBRPoqvZ4MgGQjLXAm1QVugSRzJEYHqXDzqllHwlo2EKRRcSTScEU3eHqM 0SCld2pbmJhJGLIEUgNqDiQb00iyU20AEvqCmpWbpklHovzYW15ibByp6VGH DkiEKR19WQSpATkOksqRAjxIYRmk5k5m70iqunOQGp3qyEDCBS38fn7TEzNB KlSZKvbtyXGQpI4U0kBiJZAadSQNkKRFR8mRrIy1pdGHIKFKsQJIllITiAAS a80/a4lUF1gAiZGGvbgQnFWDY4wPWgKpwZOZwFEJJHUEi2oOJNb8OW4CpAZl D6QG1BhIhPl/TAiuKFN8rF1HQpcdO0eyA5JcdJDayJakMXOksDVHopzlJZDa OJkxZcfCkYZPqbTmSKSLGYccqRv+2czzSLHacSTSBTUrZ9txR4KB5LojDd+2 dNaMgyN1JNsmQbIwZwZxJHytTjkSBKSRciTSBZSuNEBqe6oYrHFzJPytbVwI rvRBcsSR5Ds6Oo5EekaiNRFnxFtSN/zTcZDsnGpFkNxxJHsRIHIcJDty05Ha zaJBkGy4r3ekMXSk1rtxsGxdGtjRODqSubo5IcwpKGDvHQmmpkCizg22I++f aDkOkqVTrcCRE46kOrjekVqSS4bUEf9sFCRzVXNDmJNTjgTA3nVHKj2l0vo5 g1AXTnG42s+2ueeRbMk7EhsHRzJXVdNq/RzHqH1HYs6DZCvEsGUc2AnvSN1V B05xuLwjdTeEU46k7oi9I7Wm9k9xhDpgSK6DZD+ECzuhBMk7khdE3pEcCOHC TnhH8jIi70gOhHBhJ7wjeRlRBwzJdZC8I3UkgusgeXVFjoM0miezgxFcB8mr K2ruwTYrGs2T2bkI7j/Y5tUVOQ6S+ydzIyH8GMnLFTkO0miezA5GcB0kr67I cZBG82R2MILrIHl1RY6DNJons4MRXAfJqyvSBen9/xK9XPyyVx/IXmV1m9do nswORtAE6blvLryNsYeO+emFJxxKX6V1m9dotoGDETRBum/5y7cxtuoeNrtk Mn2V1m1eo9kGDkbQ7tou2Mr6J0f4nPKvyet2ad3kFS3W5HeCWAAM0n3Ll08x tmZr5EPbWf+Uf05ev66VZacORPuh3d4JMEiP7dhxIAGp/xuRF63439iRVsC6 NicORPuh3d4JZNcWgcRWbWATi6fOSl5LdXuNkTRBevdVjE0vfmlv2+DVyysV ZUKy/7P9+auXVyKDM9te4ywPkpcReZC8jMiD5GVEHiQvIzIG0mPHnzip3oqo 7xz/sq9mIWwFer/tCBPHnXjrvNUQ9x9/4hNWd2L2daxc//QrT7wjXWUKpOmF d9zXmzJUWVUTvSfuX/C1NIStQBPx3JjNCLNLtv3P4ntshnik961v9iYtRnj0 TUtZ6SD1T7780cVfS9aZAmndZYyt2mqoslrlH2JszYfSEJYCzS5Zvo1ZjfBw tBP79tsM8fBJ8/3Xb7cXYfbtr13KSgfp8WXz7H1/kqw0BVJ872TtVYYqq+rw fNTQX0lCXGYp0FlbH/k0sxph7ate3vvgIZshZpe8+pUn2Y1wLCsdpIcunWcP mwXprNuiEBsMVVbX/S9aP5+E+LCdQJdcOv/wNmYzAlu7dGr25NtWWQwxsfAj b+tN2tyJmWNZ6SC9LwZpWbLKpCNdYKtr67/xmKkshJVA/de/+BUv6f2ZxQhR +n8ewbTeZoizNsTOajPCzLHlZng4AukRs460btmh/oo7DVVW1SXJM71pCDuB Xtj3/Psvf8FmBPaBkw5GzWwzxJr1yejFYoQYpGL9j/f2Zr5nCqT+GYui3sdQ ZVWt6UW6LQ1hLdAHPs2sRpg9o/crJ+y1GeLxhYtetHTKZoSJY1n5IK3tveSE g8kqcxOSzzxrrKqaXjg8/0IewmYgqxH2PTtvN0T/Fz+3HIENa05f9/1sQKuf 2fYyIg+SlxF5kLyMyIPkZUQeJC8j8iB5GZEHiah/+6cdO3b8/bdsTZ05Jw8S Tf1TeolOmGo7k47Ig0TUmcfsjXC6sLdoqu1MuiEPEk39M9NbA+/r3cnYvvcu X/7EPHs0/lI7++9/P8TuX/6ajx5Q1DBi8iARdebSBKSJhZez7/R6r+j1PsQu ib963F+x9OCa3qJX9NL1YyMPEk39AUgziy9la+JnWt+89ODskkvn2cSCy2eX LJuP6NrWboYNy4NEUwbS7OL1bHo/6//yt6LPb160N+rrJmeXLHriCDt8pO0c G5UHiagBSNMRSM+9bUF8/XaQPRK50Iplh9gfRx9fNWYzAx4kmjJHisZF/ZN7 H/nBLy6MPvePu3RiwbZ47X9FbF3WboYNy4NEU3TVFj+1OXtcb3JmcfyM4Btj sNYuek00Xpo4fjJ+Ev/StnNsVB4kos7s3Xrvve9dHHVmM4tfNfXLi3tLp6Jr uAW99fPRuOnXD7DHFn+47RQblQeJpv6KdGb7q9H7C+MR0nt7L459Kp5WYmvj NfH3FcZIHiSiDr8Q6Rfpldm+e38+ePT0zUsPpAvu/cF4jbU9SEb1vYWXt51C W/IgmVPU3S3a23YSbcmDZFD/uWOq7RRakwfJy4g8SF5G5EHyMiIPkpcR/T/3 QuAe4aBcoQAAAABJRU5ErkJggg== --============_-944461573==_mr============--