From 513020b6ee0e913eec72c38ba858bc15e970bef4 Mon Sep 17 00:00:00 2001 From: Brian Grefenstette Date: Tue, 25 Sep 2018 09:22:39 -0700 Subject: [PATCH 1/2] Notebook updates --- notebooks/Mosaic Example.ipynb | 250 +------------- notebooks/Planning_20170911.ipynb | 534 ++++++++++++++++++++++++++++++ notebooks/Planning_Example.ipynb | 2 +- 3 files changed, 551 insertions(+), 235 deletions(-) create mode 100644 notebooks/Planning_20170911.ipynb diff --git a/notebooks/Mosaic Example.ipynb b/notebooks/Mosaic Example.ipynb index fdd40d7..2fd6279 100644 --- a/notebooks/Mosaic Example.ipynb +++ b/notebooks/Mosaic Example.ipynb @@ -2,7 +2,7 @@ "cells": [ { "cell_type": "code", - "execution_count": 1, + "execution_count": null, "metadata": { "collapsed": false }, @@ -25,19 +25,11 @@ }, { "cell_type": "code", - "execution_count": 2, + "execution_count": null, "metadata": { "collapsed": false }, - "outputs": [ - { - "name": "stdout", - "output_type": "stream", - "text": [ - "../data/NUSTAR.2018_142.SHADOW_ANALYSIS.txt\n" - ] - } - ], + "outputs": [], "source": [ "fname = io.download_occultation_times(outdir='../data/')\n", "print(fname)" @@ -56,19 +48,11 @@ }, { "cell_type": "code", - "execution_count": 3, + "execution_count": null, "metadata": { "collapsed": false }, - "outputs": [ - { - "name": "stdout", - "output_type": "stream", - "text": [ - "../data/NuSTAR.tle\n" - ] - } - ], + "outputs": [], "source": [ "tlefile = io.download_tle(outdir='../data')\n", "print(tlefile)\n", @@ -86,7 +70,7 @@ }, { "cell_type": "code", - "execution_count": 4, + "execution_count": null, "metadata": { "collapsed": false }, @@ -110,19 +94,11 @@ }, { "cell_type": "code", - "execution_count": 5, + "execution_count": null, "metadata": { "collapsed": false }, - "outputs": [ - { - "name": "stdout", - "output_type": "stream", - "text": [ - "NuSTAR Roll angle for Det0 in NE quadrant: 343.6793280064678 deg\n" - ] - } - ], + "outputs": [], "source": [ "pa = planning.get_nustar_roll(tstart, 0)\n", "print(\"NuSTAR Roll angle for Det0 in NE quadrant: {}\".format(pa))" @@ -141,7 +117,7 @@ }, { "cell_type": "code", - "execution_count": 6, + "execution_count": null, "metadata": { "collapsed": false }, @@ -161,87 +137,12 @@ }, { "cell_type": "code", - "execution_count": 7, + "execution_count": null, "metadata": { "collapsed": false, "scrolled": true }, - "outputs": [ - { - "name": "stderr", - "output_type": "stream", - "text": [ - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ - "Orbit: 0\n", - "Orbit start: 2018-05-29T15:54:40 Orbit end: 2018-05-29T16:56:30\n", - "Aim time: 2018-05-29T16:25:35 RA (deg): 66.237867206401 deg Dec (deg): 21.641611036044274 deg\n", - "\n" - ] - }, - { - "name": "stderr", - "output_type": "stream", - "text": [ - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ - "Orbit: 1\n", - "Orbit start: 2018-05-29T17:31:20 Orbit end: 2018-05-29T18:33:10\n", - "Aim time: 2018-05-29T18:02:15 RA (deg): 66.30622617158947 deg Dec (deg): 21.651913534159284 deg\n", - "\n" - ] - }, - { - "name": "stderr", - "output_type": "stream", - "text": [ - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ - "Orbit: 2\n", - "Orbit start: 2018-05-29T19:08:00 Orbit end: 2018-05-29T20:09:50\n", - "Aim time: 2018-05-29T19:38:55 RA (deg): 66.37459341208549 deg Dec (deg): 21.662187831555375 deg\n", - "\n" - ] - }, - { - "name": "stderr", - "output_type": "stream", - "text": [ - "\r", - "[#################################] 100% deltat.preds\n" - ] - }, - { - "name": "stdout", - "output_type": "stream", - "text": [ - "Orbit: 3\n", - "Orbit start: 2018-05-29T20:44:40 Orbit end: 2018-05-29T21:46:30\n", - "Aim time: 2018-05-29T21:15:35 RA (deg): 66.44296891430218 deg Dec (deg): 21.672433912678994 deg\n", - "\n" - ] - } - ], + "outputs": [], "source": [ "for ind, orbit in enumerate(orbits):\n", " midTime = (0.5*(orbit[1] - orbit[0]) + orbit[0])\n", @@ -264,7 +165,7 @@ }, { "cell_type": "code", - "execution_count": 8, + "execution_count": null, "metadata": { "collapsed": false }, @@ -360,7 +261,7 @@ }, { "cell_type": "code", - "execution_count": 9, + "execution_count": null, "metadata": { "collapsed": true }, @@ -399,131 +300,12 @@ }, { "cell_type": "code", - "execution_count": 11, + "execution_count": null, "metadata": { "collapsed": false, "scrolled": false }, - "outputs": [ - { - "name": "stdout", - "output_type": "stream", - "text": [ - "Step PA 433.6837468416314 deg\n", - "Orbit start: 2018-05-29T15:54:40 Orbit end: 2018-05-29T16:56:30\n", - "Dwell per position: 0:02:28.400000\n", - "\n", - "NuSTAR Roll Angle to get DET0 in top right 343.68 deg\n", - "\n" - ] - }, - { - "name": "stderr", - "output_type": "stream", - "text": [ - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n", - "[#################################] 100% deltat.preds\n" - ] - } - ], + "outputs": [], "source": [ "# Just use the first orbit...or choose one. This may download a ton of deltat.preds, which is a known \n", "# bug to be fixed.\n", @@ -559,7 +341,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.5.2" + "version": "3.5.1" } }, "nbformat": 4, diff --git a/notebooks/Planning_20170911.ipynb b/notebooks/Planning_20170911.ipynb new file mode 100644 index 0000000..2a773f0 --- /dev/null +++ b/notebooks/Planning_20170911.ipynb @@ -0,0 +1,534 @@ +{ + "cells": [ + { + "cell_type": "code", + "execution_count": 1, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + "WARNING: AstropyDeprecationWarning: astropy.utils.compat.odict.OrderedDict is now deprecated - import OrderedDict from the collections module instead [astropy.utils.compat.odict]\n" + ] + } + ], + "source": [ + "from nustar_pysolar import planning, io\n", + "import astropy.units as u" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Download the list of occultation periods from the MOC at Berkeley.\n", + "\n", + "## Note that the occultation periods typically only are stored at Berkeley for the *future* and not for the past. So this is only really useful for observation planning." + ] + }, + { + "cell_type": "code", + "execution_count": 2, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "../data/NUSTAR.2017_250.SHADOW_ANALYSIS.txt\n" + ] + } + ], + "source": [ + "fname = io.download_occultation_times(outdir='../data/')\n", + "print(fname)" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Download the NuSTAR TLE archive.\n", + "\n", + "This contains every two-line element (TLE) that we've received for the whole mission. We'll expand on how to use this later.\n", + "\n", + "The `times`, `line1`, and `line2` elements are now the TLE elements for each epoch." + ] + }, + { + "cell_type": "code", + "execution_count": 3, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "../data/NuSTAR.tle\n" + ] + } + ], + "source": [ + "tlefile = io.download_tle(outdir='../data')\n", + "print(tlefile)\n", + "times, line1, line2 = io.read_tle_file(tlefile)" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Here is where we define the observing window that we want to use.\n", + "\n", + "Note that tstart and tend must be in the future otherwise you won't find any occultation times and sunlight_periods will return an error." + ] + }, + { + "cell_type": "code", + "execution_count": 7, + "metadata": { + "collapsed": false + }, + "outputs": [], + "source": [ + "tstart = '2017-09-11T00:00:00'\n", + "tend = '2017-09-15T00:00:00'\n", + "orbits = planning.sunlight_periods(fname, tstart, tend)" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# We want to know how to orient NuSTAR for the Sun.\n", + "\n", + "We can more or less pick any angle that we want. But this angle has to be specified a little in advance so that the NuSTAR SOC can plan the \"slew in\" maneuvers. Below puts DET0 in the top left corner (north-east with respect to RA/Dec coordinates).\n", + "\n", + "### This is what you tell the SOC you want the \"Sky PA angle\" to be." + ] + }, + { + "cell_type": "code", + "execution_count": 8, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "NuSTAR Roll angle for Det0 in NE quadrant: 23.324561616831126 deg\n" + ] + } + ], + "source": [ + "pa = planning.get_nustar_roll(tstart, 0)\n", + "print(\"NuSTAR Roll angle for Det0 in NE quadrant: {}\".format(pa))" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Set up the offset you want to use here:\n", + "\n", + "The first element is the direction +WEST of the center of the Sun, the second is the offset +NORTH of the center of the Sun.\n", + "\n", + "If you want multiple pointing locations you can either specify an array of offsets or do this \"by hand\" below." + ] + }, + { + "cell_type": "code", + "execution_count": 9, + "metadata": { + "collapsed": false + }, + "outputs": [], + "source": [ + "offset = [1100., -400.]*u.arcsec" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Loop over each orbit and correct the pointing for the same heliocentric pointing position.\n", + "\n", + "Note that you may want to update the pointing for solar rotation. That's up to the user to decide and is not done here." + ] + }, + { + "cell_type": "code", + "execution_count": 11, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Orbit start: 2017-09-10T23:48:20 Orbit end: 2017-09-11T00:49:30\n", + "Aim Time: 2017-09-11T00:18:55\n", + "Aim time: 2017-09-11T00:18:55 RA (deg): 168.8444822500469 deg Dec (deg): 4.675665760939547 deg\n", + "\n", + "Orbit start: 2017-09-11T01:25:00 Orbit end: 2017-09-11T02:26:20\n", + "Aim Time: 2017-09-11T01:55:40\n", + "Aim time: 2017-09-11T01:55:40 RA (deg): 168.90482172609111 deg Dec (deg): 4.650220901438387 deg\n", + "\n", + "Orbit start: 2017-09-11T03:01:40 Orbit end: 2017-09-11T04:03:00\n", + "Aim Time: 2017-09-11T03:32:20\n", + "Aim time: 2017-09-11T03:32:20 RA (deg): 168.96510726129736 deg Dec (deg): 4.624791714579964 deg\n", + "\n", + "Orbit start: 2017-09-11T04:38:20 Orbit end: 2017-09-11T05:39:40\n", + "Aim Time: 2017-09-11T05:09:00\n", + "Aim time: 2017-09-11T05:09:00 RA (deg): 169.02539085542193 deg Dec (deg): 4.599356312960442 deg\n", + "\n", + "Orbit start: 2017-09-11T06:15:00 Orbit end: 2017-09-11T07:16:20\n", + "Aim Time: 2017-09-11T06:45:40\n", + "Aim time: 2017-09-11T06:45:40 RA (deg): 169.0856725353157 deg Dec (deg): 4.573914722834457 deg\n", + "\n", + "Orbit start: 2017-09-11T07:51:40 Orbit end: 2017-09-11T08:53:00\n", + "Aim Time: 2017-09-11T08:22:20\n", + "Aim time: 2017-09-11T08:22:20 RA (deg): 169.14595232776264 deg Dec (deg): 4.548466970481036 deg\n", + "\n", + "Orbit start: 2017-09-11T09:28:30 Orbit end: 2017-09-11T10:29:40\n", + "Aim Time: 2017-09-11T09:59:05\n", + "Aim time: 2017-09-11T09:59:05 RA (deg): 169.2062822224168 deg Dec (deg): 4.522991136560128 deg\n", + "\n", + "Orbit start: 2017-09-11T11:05:10 Orbit end: 2017-09-11T12:06:20\n", + "Aim Time: 2017-09-11T11:35:45\n", + "Aim time: 2017-09-11T11:35:45 RA (deg): 169.26655831847657 deg Dec (deg): 4.49753113344228 deg\n", + "\n", + "Orbit start: 2017-09-11T12:41:50 Orbit end: 2017-09-11T13:43:00\n", + "Aim Time: 2017-09-11T13:12:25\n", + "Aim time: 2017-09-11T13:12:25 RA (deg): 169.32683260705247 deg Dec (deg): 4.472065047118643 deg\n", + "\n", + "Orbit start: 2017-09-11T14:18:30 Orbit end: 2017-09-11T15:19:40\n", + "Aim Time: 2017-09-11T14:49:05\n", + "Aim time: 2017-09-11T14:49:05 RA (deg): 169.38710511464714 deg Dec (deg): 4.446592903978586 deg\n", + "\n", + "Orbit start: 2017-09-11T15:55:10 Orbit end: 2017-09-11T16:56:20\n", + "Aim Time: 2017-09-11T16:25:45\n", + "Aim time: 2017-09-11T16:25:45 RA (deg): 169.44737586769475 deg Dec (deg): 4.421114730438243 deg\n", + "\n", + "Orbit start: 2017-09-11T17:31:50 Orbit end: 2017-09-11T18:33:00\n", + "Aim Time: 2017-09-11T18:02:25\n", + "Aim time: 2017-09-11T18:02:25 RA (deg): 169.50764489255593 deg Dec (deg): 4.395630552942596 deg\n", + "\n", + "Orbit start: 2017-09-11T19:08:30 Orbit end: 2017-09-11T20:09:40\n", + "Aim Time: 2017-09-11T19:39:05\n", + "Aim time: 2017-09-11T19:39:05 RA (deg): 169.5679122155163 deg Dec (deg): 4.370140397965842 deg\n", + "\n", + "Orbit start: 2017-09-11T20:45:10 Orbit end: 2017-09-11T21:46:20\n", + "Aim Time: 2017-09-11T21:15:45\n", + "Aim time: 2017-09-11T21:15:45 RA (deg): 169.62817786278617 deg Dec (deg): 4.344644292013325 deg\n", + "\n", + "Orbit start: 2017-09-11T22:21:50 Orbit end: 2017-09-11T23:23:00\n", + "Aim Time: 2017-09-11T22:52:25\n", + "Aim time: 2017-09-11T22:52:25 RA (deg): 169.68844186050165 deg Dec (deg): 4.319142261618367 deg\n", + "\n", + "Orbit start: 2017-09-11T23:58:30 Orbit end: 2017-09-12T00:59:40\n", + "Aim Time: 2017-09-12T00:29:05\n", + "Aim time: 2017-09-12T00:29:05 RA (deg): 169.74870423471913 deg Dec (deg): 4.293634333346423 deg\n", + "\n", + "Orbit start: 2017-09-12T01:35:10 Orbit end: 2017-09-12T02:36:20\n", + "Aim Time: 2017-09-12T02:05:45\n", + "Aim time: 2017-09-12T02:05:45 RA (deg): 169.8089650114236 deg Dec (deg): 4.26812053379131 deg\n", + "\n", + "Orbit start: 2017-09-12T03:11:50 Orbit end: 2017-09-12T04:13:00\n", + "Aim Time: 2017-09-12T03:42:25\n", + "Aim time: 2017-09-12T03:42:25 RA (deg): 169.8692242165179 deg Dec (deg): 4.24260088957969 deg\n", + "\n", + "Orbit start: 2017-09-12T04:48:30 Orbit end: 2017-09-12T05:49:40\n", + "Aim Time: 2017-09-12T05:19:05\n", + "Aim time: 2017-09-12T05:19:05 RA (deg): 169.92948187583292 deg Dec (deg): 4.217075427366984 deg\n", + "\n", + "Orbit start: 2017-09-12T06:25:10 Orbit end: 2017-09-12T07:26:20\n", + "Aim Time: 2017-09-12T06:55:45\n", + "Aim time: 2017-09-12T06:55:45 RA (deg): 169.9897380151151 deg Dec (deg): 4.19154417384252 deg\n", + "\n", + "Orbit start: 2017-09-12T08:01:50 Orbit end: 2017-09-12T09:03:00\n", + "Aim Time: 2017-09-12T08:32:25\n", + "Aim time: 2017-09-12T08:32:25 RA (deg): 170.04999266003585 deg Dec (deg): 4.166007155725799 deg\n", + "\n", + "Orbit start: 2017-09-12T09:38:30 Orbit end: 2017-09-12T10:39:40\n", + "Aim Time: 2017-09-12T10:09:05\n", + "Aim time: 2017-09-12T10:09:05 RA (deg): 170.11024583618845 deg Dec (deg): 4.140464399767831 deg\n", + "\n", + "Orbit start: 2017-09-12T11:15:10 Orbit end: 2017-09-12T12:16:20\n", + "Aim Time: 2017-09-12T11:45:45\n", + "Aim time: 2017-09-12T11:45:45 RA (deg): 170.17049756908847 deg Dec (deg): 4.114915932751004 deg\n", + "\n", + "Orbit start: 2017-09-12T12:51:50 Orbit end: 2017-09-12T13:53:00\n", + "Aim Time: 2017-09-12T13:22:25\n", + "Aim time: 2017-09-12T13:22:25 RA (deg): 170.23074788416437 deg Dec (deg): 4.089361781493244 deg\n", + "\n", + "Orbit start: 2017-09-12T14:28:30 Orbit end: 2017-09-12T15:29:50\n", + "Aim Time: 2017-09-12T14:59:10\n", + "Aim time: 2017-09-12T14:59:10 RA (deg): 170.29104874490477 deg Dec (deg): 4.063779936090902 deg\n", + "\n", + "Orbit start: 2017-09-12T16:05:10 Orbit end: 2017-09-12T17:06:30\n", + "Aim Time: 2017-09-12T16:35:50\n", + "Aim time: 2017-09-12T16:35:50 RA (deg): 170.35129629914323 deg Dec (deg): 4.038214492087479 deg\n", + "\n", + "Orbit start: 2017-09-12T17:41:50 Orbit end: 2017-09-12T18:43:10\n", + "Aim Time: 2017-09-12T18:12:30\n", + "Aim time: 2017-09-12T18:12:30 RA (deg): 170.4115425113941 deg Dec (deg): 4.0126434445118555 deg\n", + "\n", + "Orbit start: 2017-09-12T19:18:30 Orbit end: 2017-09-12T20:19:50\n", + "Aim Time: 2017-09-12T19:49:10\n", + "Aim time: 2017-09-12T19:49:10 RA (deg): 170.47178740676057 deg Dec (deg): 3.987066820315032 deg\n", + "\n", + "Orbit start: 2017-09-12T20:55:20 Orbit end: 2017-09-12T21:56:30\n", + "Aim Time: 2017-09-12T21:25:55\n", + "Aim time: 2017-09-12T21:25:55 RA (deg): 170.5320829438566 deg Dec (deg): 3.961462590496711 deg\n", + "\n", + "Orbit start: 2017-09-12T22:32:00 Orbit end: 2017-09-12T23:33:10\n", + "Aim Time: 2017-09-12T23:02:35\n", + "Aim time: 2017-09-12T23:02:35 RA (deg): 170.59232527935984 deg Dec (deg): 3.935874889294554 deg\n", + "\n", + "Orbit start: 2017-09-13T00:08:40 Orbit end: 2017-09-13T01:09:50\n", + "Aim Time: 2017-09-13T00:39:15\n", + "Aim time: 2017-09-13T00:39:15 RA (deg): 170.65256637281706 deg Dec (deg): 3.9102816925508765 deg\n", + "\n", + "Orbit start: 2017-09-13T01:45:20 Orbit end: 2017-09-13T02:46:30\n", + "Aim Time: 2017-09-13T02:15:55\n", + "Aim time: 2017-09-13T02:15:55 RA (deg): 170.7128062490058 deg Dec (deg): 3.8846830273527364 deg\n", + "\n", + "Orbit start: 2017-09-13T03:22:00 Orbit end: 2017-09-13T04:23:10\n", + "Aim Time: 2017-09-13T03:52:35\n", + "Aim time: 2017-09-13T03:52:35 RA (deg): 170.7730449326128 deg Dec (deg): 3.859078920823773 deg\n", + "\n", + "Orbit start: 2017-09-13T04:58:40 Orbit end: 2017-09-13T05:59:50\n", + "Aim Time: 2017-09-13T05:29:15\n", + "Aim time: 2017-09-13T05:29:15 RA (deg): 170.83328244824577 deg Dec (deg): 3.8334694001212903 deg\n", + "\n", + "Orbit start: 2017-09-13T06:35:20 Orbit end: 2017-09-13T07:36:30\n", + "Aim Time: 2017-09-13T07:05:55\n", + "Aim time: 2017-09-13T07:05:55 RA (deg): 170.89351882042817 deg Dec (deg): 3.807854492438098 deg\n", + "\n", + "Orbit start: 2017-09-13T08:12:00 Orbit end: 2017-09-13T09:13:10\n", + "Aim Time: 2017-09-13T08:42:35\n", + "Aim time: 2017-09-13T08:42:35 RA (deg): 170.95375407360126 deg Dec (deg): 3.7822342250018606 deg\n", + "\n", + "Orbit start: 2017-09-13T09:48:40 Orbit end: 2017-09-13T10:49:50\n", + "Aim Time: 2017-09-13T10:19:15\n", + "Aim time: 2017-09-13T10:19:15 RA (deg): 171.0139882321224 deg Dec (deg): 3.75660862507583 deg\n", + "\n", + "Orbit start: 2017-09-13T11:25:20 Orbit end: 2017-09-13T12:26:30\n", + "Aim Time: 2017-09-13T11:55:55\n", + "Aim time: 2017-09-13T11:55:55 RA (deg): 171.07422132026522 deg Dec (deg): 3.7309777199590077 deg\n", + "\n", + "Orbit start: 2017-09-13T13:02:00 Orbit end: 2017-09-13T14:03:10\n", + "Aim Time: 2017-09-13T13:32:35\n", + "Aim time: 2017-09-13T13:32:35 RA (deg): 171.13445336222037 deg Dec (deg): 3.7053415369859573 deg\n", + "\n", + "Orbit start: 2017-09-13T14:38:40 Orbit end: 2017-09-13T15:39:50\n", + "Aim Time: 2017-09-13T15:09:15\n", + "Aim time: 2017-09-13T15:09:15 RA (deg): 171.19468438209447 deg Dec (deg): 3.679700103527339 deg\n", + "\n", + "Orbit start: 2017-09-13T16:15:20 Orbit end: 2017-09-13T17:16:30\n", + "Aim Time: 2017-09-13T16:45:55\n", + "Aim time: 2017-09-13T16:45:55 RA (deg): 171.25491440390908 deg Dec (deg): 3.654053446990595 deg\n", + "\n", + "Orbit start: 2017-09-13T17:52:00 Orbit end: 2017-09-13T18:53:10\n", + "Aim Time: 2017-09-13T18:22:35\n", + "Aim time: 2017-09-13T18:22:35 RA (deg): 171.3151434516058 deg Dec (deg): 3.628401594817846 deg\n", + "\n", + "Orbit start: 2017-09-13T19:28:40 Orbit end: 2017-09-13T20:29:50\n", + "Aim Time: 2017-09-13T19:59:15\n", + "Aim time: 2017-09-13T19:59:15 RA (deg): 171.3753715490367 deg Dec (deg): 3.6027445744902167 deg\n", + "\n", + "Orbit start: 2017-09-13T21:05:20 Orbit end: 2017-09-13T22:06:30\n", + "Aim Time: 2017-09-13T21:35:55\n", + "Aim time: 2017-09-13T21:35:55 RA (deg): 171.43559871997257 deg Dec (deg): 3.577082413524341 deg\n", + "\n", + "Orbit start: 2017-09-13T22:42:00 Orbit end: 2017-09-13T23:43:10\n", + "Aim Time: 2017-09-13T23:12:35\n", + "Aim time: 2017-09-13T23:12:35 RA (deg): 171.49582498809883 deg Dec (deg): 3.551415139473964 deg\n", + "\n", + "Orbit start: 2017-09-14T00:18:40 Orbit end: 2017-09-14T01:19:50\n", + "Aim Time: 2017-09-14T00:49:15\n", + "Aim time: 2017-09-14T00:49:15 RA (deg): 171.556050377021 deg Dec (deg): 3.5257427799295926 deg\n", + "\n", + "Orbit start: 2017-09-14T01:55:20 Orbit end: 2017-09-14T02:56:30\n", + "Aim Time: 2017-09-14T02:25:55\n", + "Aim time: 2017-09-14T02:25:55 RA (deg): 171.61627491025928 deg Dec (deg): 3.5000653625181024 deg\n", + "\n", + "Orbit start: 2017-09-14T03:32:00 Orbit end: 2017-09-14T04:33:20\n", + "Aim Time: 2017-09-14T04:02:40\n", + "Aim time: 2017-09-14T04:02:40 RA (deg): 171.6765505278821 deg Dec (deg): 3.474360772700875 deg\n", + "\n", + "Orbit start: 2017-09-14T05:08:40 Orbit end: 2017-09-14T06:10:00\n", + "Aim Time: 2017-09-14T05:39:20\n", + "Aim time: 2017-09-14T05:39:20 RA (deg): 171.7367734192945 deg Dec (deg): 3.4486733182845675 deg\n", + "\n", + "Orbit start: 2017-09-14T06:45:20 Orbit end: 2017-09-14T07:46:40\n", + "Aim Time: 2017-09-14T07:16:00\n", + "Aim time: 2017-09-14T07:16:00 RA (deg): 171.79699552511318 deg Dec (deg): 3.4229808891266287 deg\n", + "\n", + "Orbit start: 2017-09-14T08:22:00 Orbit end: 2017-09-14T09:23:20\n", + "Aim Time: 2017-09-14T08:52:40\n", + "Aim time: 2017-09-14T08:52:40 RA (deg): 171.8572168685288 deg Dec (deg): 3.3972835130024457 deg\n", + "\n", + "Orbit start: 2017-09-14T09:58:40 Orbit end: 2017-09-14T11:00:00\n", + "Aim Time: 2017-09-14T10:29:20\n", + "Aim time: 2017-09-14T10:29:20 RA (deg): 171.91743747265608 deg Dec (deg): 3.3715812177221207 deg\n", + "\n", + "Orbit start: 2017-09-14T11:35:30 Orbit end: 2017-09-14T12:36:40\n", + "Aim Time: 2017-09-14T12:06:05\n", + "Aim time: 2017-09-14T12:06:05 RA (deg): 171.97770927392142 deg Dec (deg): 3.3458518676631885 deg\n", + "\n", + "Orbit start: 2017-09-14T13:12:10 Orbit end: 2017-09-14T14:13:20\n", + "Aim Time: 2017-09-14T13:42:45\n", + "Aim time: 2017-09-14T13:42:45 RA (deg): 172.03792846790213 deg Dec (deg): 3.32013981346965 deg\n", + "\n", + "Orbit start: 2017-09-14T14:48:50 Orbit end: 2017-09-14T15:50:00\n", + "Aim Time: 2017-09-14T15:19:25\n", + "Aim time: 2017-09-14T15:19:25 RA (deg): 172.09814699147296 deg Dec (deg): 3.2944229237971934 deg\n", + "\n", + "Orbit start: 2017-09-14T16:25:30 Orbit end: 2017-09-14T17:26:40\n", + "Aim Time: 2017-09-14T16:56:05\n", + "Aim time: 2017-09-14T16:56:05 RA (deg): 172.15836486742774 deg Dec (deg): 3.268701226603889 deg\n", + "\n", + "Orbit start: 2017-09-14T18:02:10 Orbit end: 2017-09-14T19:03:20\n", + "Aim Time: 2017-09-14T18:32:45\n", + "Aim time: 2017-09-14T18:32:45 RA (deg): 172.21858211848098 deg Dec (deg): 3.2429747498850263 deg\n", + "\n", + "Orbit start: 2017-09-14T19:38:50 Orbit end: 2017-09-14T20:40:00\n", + "Aim Time: 2017-09-14T20:09:25\n", + "Aim time: 2017-09-14T20:09:25 RA (deg): 172.27879876726846 deg Dec (deg): 3.217243521672911 deg\n", + "\n", + "Orbit start: 2017-09-14T21:15:30 Orbit end: 2017-09-14T22:16:40\n", + "Aim Time: 2017-09-14T21:46:05\n", + "Aim time: 2017-09-14T21:46:05 RA (deg): 172.3390148363467 deg Dec (deg): 3.1915075700373032 deg\n", + "\n", + "Orbit start: 2017-09-14T22:52:10 Orbit end: 2017-09-14T23:53:20\n", + "Aim Time: 2017-09-14T23:22:45\n", + "Aim time: 2017-09-14T23:22:45 RA (deg): 172.39923034819873 deg Dec (deg): 3.1657669230832437 deg\n", + "\n" + ] + } + ], + "source": [ + "for ind, orbit in enumerate(orbits):\n", + " midTime = (0.5*(orbit[1] - orbit[0]) + orbit[0])\n", + " sky_pos = planning.get_sky_position(midTime, offset)\n", + "# print(\"Orbit: {}\".format(ind))\n", + " print(\"Orbit start: {} Orbit end: {}\".format(orbit[0].isoformat(), orbit[1].isoformat()))\n", + " print(\"Aim Time: {}\".format(midTime.isoformat()))\n", + " print('Aim time: {} RA (deg): {} Dec (deg): {}'.format(midTime.isoformat(), sky_pos[0], sky_pos[1]))\n", + " print(\"\")\n", + "\n", + " \n" + ] + }, + { + "cell_type": "code", + "execution_count": 12, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "data": { + "text/latex": [ + "$[168.88111,~4.6901735] \\; \\mathrm{{}^{\\circ}}$" + ], + "text/plain": [ + "" + ] + }, + "execution_count": 12, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [] + }, + { + "cell_type": "markdown", + "metadata": { + "collapsed": true + }, + "source": [ + "# Sanity check: Reproduce the 2016-07-26 pointing.\n", + "\n", + "Based on the IDL version of this code, when the aim time was 2016-07-26 19:53:15.00 the NP was 8.86 degrees and the RA and Dec were 126.04630 and 19.33532 deg, respectively for an offset of +1000 arcseconds West and +150 arcseconds North of the center of the Sun." + ] + }, + { + "cell_type": "code", + "execution_count": 8, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "[ 126.04053869 19.33666449] deg\n" + ] + } + ], + "source": [ + "aim_time = '2016-07-26T19:53:15.00'\n", + "offset = [1000, 150]*u.arcsec\n", + "sky_pos = planning.get_sky_position(aim_time, offset)\n", + "print(sky_pos)" + ] + }, + { + "cell_type": "code", + "execution_count": 9, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "8.865884432089102 deg\n" + ] + } + ], + "source": [ + "np = planning.get_nustar_roll(aim_time, 0)\n", + "print(np)" + ] + } + ], + "metadata": { + "anaconda-cloud": {}, + "kernelspec": { + "display_name": "Python [default]", + "language": "python", + "name": "python3" + }, + "language_info": { + "codemirror_mode": { + "name": "ipython", + "version": 3 + }, + "file_extension": ".py", + "mimetype": "text/x-python", + "name": "python", + "nbconvert_exporter": "python", + "pygments_lexer": "ipython3", + "version": "3.5.1" + } + }, + "nbformat": 4, + "nbformat_minor": 2 +} diff --git a/notebooks/Planning_Example.ipynb b/notebooks/Planning_Example.ipynb index c66558b..e5923d7 100644 --- a/notebooks/Planning_Example.ipynb +++ b/notebooks/Planning_Example.ipynb @@ -280,7 +280,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.5.2" + "version": "3.5.1" } }, "nbformat": 4, From 03dc5cd47db78fe7dc54ed3ee73570841dbdc5a0 Mon Sep 17 00:00:00 2001 From: Brian Grefenstette Date: Tue, 25 Sep 2018 09:42:36 -0700 Subject: [PATCH 2/2] Moved the mosaic code into planning.py ...and updated the Mosaic example notebook to use the version of make_mosaic from the planning software. This seems to work (from a ds9 perspective). Note that the doc strings probably aren't correct, and most of my code needs safety checks to make sure the inputs are correctly formatted, but this should be good enough for now. --- notebooks/Mosaic 20180928.ipynb | 257 +++++++++++ notebooks/Mosaic Example.ipynb | 50 ++- notebooks/NuSTAR.tle | 750 ++++++++++++++++++++++++++++++++ nustar_pysolar/planning.py | 132 ++++++ 4 files changed, 1177 insertions(+), 12 deletions(-) create mode 100644 notebooks/Mosaic 20180928.ipynb diff --git a/notebooks/Mosaic 20180928.ipynb b/notebooks/Mosaic 20180928.ipynb new file mode 100644 index 0000000..d338f14 --- /dev/null +++ b/notebooks/Mosaic 20180928.ipynb @@ -0,0 +1,257 @@ +{ + "cells": [ + { + "cell_type": "code", + "execution_count": 1, + "metadata": { + "collapsed": false + }, + "outputs": [], + "source": [ + "from nustar_pysolar import planning, io\n", + "import astropy.units as u\n", + "import warnings\n", + "warnings.filterwarnings('ignore')" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Download the list of occultation periods from the MOC at Berkeley.\n", + "\n", + "## Note that the occultation periods typically only are stored at Berkeley for the *future* and not for the past. So this is only really useful for observation planning." + ] + }, + { + "cell_type": "code", + "execution_count": 2, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "../data/NUSTAR.2018_267.SHADOW_ANALYSIS.txt\n" + ] + } + ], + "source": [ + "fname = io.download_occultation_times(outdir='../data/')\n", + "print(fname)" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Download the NuSTAR TLE archive.\n", + "\n", + "This contains every two-line element (TLE) that we've received for the whole mission. We'll expand on how to use this later.\n", + "\n", + "The `times`, `line1`, and `line2` elements are now the TLE elements for each epoch." + ] + }, + { + "cell_type": "code", + "execution_count": 3, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "../data/NuSTAR.tle\n" + ] + } + ], + "source": [ + "tlefile = io.download_tle(outdir='../data')\n", + "print(tlefile)\n", + "times, line1, line2 = io.read_tle_file(tlefile)" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Here is where we define the observing window that we want to use.\n", + "\n", + "Note that tstart and tend must be in the future otherwise you won't find any occultation times and sunlight_periods will return an error." + ] + }, + { + "cell_type": "code", + "execution_count": 4, + "metadata": { + "collapsed": false + }, + "outputs": [], + "source": [ + "tstart = '2018-09-27T12:00:00'\n", + "tend = '2018-09-29T12:10:00'\n", + "orbits = planning.sunlight_periods(fname, tstart, tend)" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# We want to know how to orient NuSTAR for the Sun.\n", + "\n", + "We can more or less pick any angle that we want. But this angle has to be specified a little in advance so that the NuSTAR SOC can plan the \"slew in\" maneuvers. Below puts DET0 in the top left corner (north-east with respect to RA/Dec coordinates).\n", + "\n", + "### This is what you tell the SOC you want the \"Sky PA angle\" to be." + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": false + }, + "outputs": [], + "source": [ + "pa = planning.get_nustar_roll(tstart, 0)\n", + "print(\"NuSTAR Roll angle for Det0 in NE quadrant: {}\".format(pa))" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Set up the offset you want to use here:\n", + "\n", + "The first element is the direction +WEST of the center of the Sun, the second is the offset +NORTH of the center of the Sun.\n", + "\n", + "If you want multiple pointing locations you can either specify an array of offsets or do this \"by hand\" below." + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": false + }, + "outputs": [], + "source": [ + "offset = [0., 0.]*u.arcsec" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Loop over each orbit and correct the pointing for the same heliocentric pointing position.\n", + "\n", + "Note that you may want to update the pointing for solar rotation. That's up to the user to decide and is not done here." + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": false, + "scrolled": true + }, + "outputs": [], + "source": [ + "for ind, orbit in enumerate(orbits):\n", + " midTime = (0.5*(orbit[1] - orbit[0]) + orbit[0])\n", + " sky_pos = planning.get_skyfield_position(midTime, offset, parallax_correction=True)\n", + " print(\"Orbit: {}\".format(ind))\n", + " print(\"Orbit start: {} Orbit end: {}\".format(orbit[0].isoformat(), orbit[1].isoformat()))\n", + " print('Aim time: {} RA (deg): {} Dec (deg): {}'.format(midTime.isoformat(), sky_pos[0], sky_pos[1]))\n", + " print(\"\")" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# This is where you actually make the Mosaic" + ] + }, + { + "cell_type": "code", + "execution_count": 5, + "metadata": { + "collapsed": false, + "scrolled": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('Step PA', )\n", + "Orbit start: 2018-09-27T11:47:30 Orbit end: 2018-09-27T12:48:50\n", + "('Dwell per position:', datetime.timedelta(0, 147, 200000))\n", + "\n", + "NuSTAR Roll Angle to get DET0 in top right 25.67 deg\n", + "\n" + ] + }, + { + "name": "stderr", + "output_type": "stream", + "text": [ + "[#################################] 100% deltat.data\n", + "[#################################] 100% Leap_Second.dat\n" + ] + } + ], + "source": [ + "# Just use the first orbit...or choose one. This may download a ton of deltat.preds, which is a known \n", + "# bug to be fixed.\n", + "\n", + "orbit = orbits[0]\n", + "planning.make_mosaic(orbits[0], write_output=True, make_regions=True)" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": false + }, + "outputs": [], + "source": [] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], + "source": [] + } + ], + "metadata": { + "anaconda-cloud": {}, + "kernelspec": { + "display_name": "Python [default]", + "language": "python", + "name": "python3" + }, + "language_info": { + "codemirror_mode": { + "name": "ipython", + "version": 3 + }, + "file_extension": ".py", + "mimetype": "text/x-python", + "name": "python", + "nbconvert_exporter": "python", + "pygments_lexer": "ipython3", + "version": "3.5.1" + } + }, + "nbformat": 4, + "nbformat_minor": 2 +} diff --git a/notebooks/Mosaic Example.ipynb b/notebooks/Mosaic Example.ipynb index 2fd6279..604166d 100644 --- a/notebooks/Mosaic Example.ipynb +++ b/notebooks/Mosaic Example.ipynb @@ -2,7 +2,7 @@ "cells": [ { "cell_type": "code", - "execution_count": null, + "execution_count": 1, "metadata": { "collapsed": false }, @@ -25,11 +25,19 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 2, "metadata": { "collapsed": false }, - "outputs": [], + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "../data/NUSTAR.2018_267.SHADOW_ANALYSIS.txt\n" + ] + } + ], "source": [ "fname = io.download_occultation_times(outdir='../data/')\n", "print(fname)" @@ -48,11 +56,19 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 3, "metadata": { "collapsed": false }, - "outputs": [], + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "../data/NuSTAR.tle\n" + ] + } + ], "source": [ "tlefile = io.download_tle(outdir='../data')\n", "print(tlefile)\n", @@ -70,11 +86,20 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 4, "metadata": { "collapsed": false }, - "outputs": [], + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Error in function: sunlight_periods\n", + "No dates found in range. Pick a different occultation file.\n" + ] + } + ], "source": [ "tstart = '2018-05-29T15:37:00'\n", "tend = '2018-05-29T23:10:00'\n", @@ -175,12 +200,14 @@ " ''' \n", " Code to make a mosaic for a 4x4 tiled array on the Sun.\n", " \n", - " Input:\n", + " Input: orbit\n", " \n", - " tstart = '2018-05-28T15:37:00'\n", - " tend = '2018-05-28T23:10:00'\n", + " ...of the type you get from doing:\n", " \n", - " positions = make_mosaic(tstart, tend, write_output=True)\n", + " orbits = planning.sunlight_periods(fname, tstart, tend)\n", + " \n", + " Syntax:\n", + " positions = make_mosaic(orbit, write_output=True)\n", " \n", " Optional flags:\n", " \n", @@ -193,7 +220,6 @@ " make_regions: [False] / True\n", " Make ds9 region files for each tile so that you can see how the FoV moves with each mosaic location.\n", "\n", - "\n", " Output mosaic file has columns of:\n", " \"Arrive By Time\" RA DEC RA_SUN DEC_SUN\n", "\n", diff --git a/notebooks/NuSTAR.tle b/notebooks/NuSTAR.tle index 62e1df2..f327903 100644 --- a/notebooks/NuSTAR.tle +++ b/notebooks/NuSTAR.tle @@ -2364,3 +2364,753 @@ 2 99105 006.0260 141.4355 0012470 267.4386 092.4315 14.87762210275908 1 99105U 12031A 17192.38576808 +.00000920 +00000-0 +22946-4 0 9994 2 99105 006.0261 133.6490 0012454 282.1556 077.7178 14.87762610276062 +1 99105U 12031A 17193.45833028 +.00000913 +00000-0 +22178-4 0 9992 +2 99105 006.0263 125.8635 0012379 296.7657 063.1206 14.87763053276228 +1 99105U 12031A 17194.39682381 +.00000922 +00000-0 +23212-4 0 9990 +2 99105 006.0266 119.0523 0012296 309.9305 049.9744 14.87763478276362 +1 99105U 12031A 17195.06717643 +.00000935 +00000-0 +24537-4 0 9990 +2 99105 006.0266 114.1861 0012222 319.1966 040.7248 14.87763841276461 +1 99105U 12031A 17196.27381280 +.00000962 +00000-0 +27518-4 0 9990 +2 99105 006.0269 105.4288 0012057 336.4496 023.5081 14.87764447276643 +1 99105U 12031A 17197.07823760 +.00000974 +00000-0 +28737-4 0 9990 +2 99105 006.0265 099.5948 0011940 347.7348 012.2490 14.87764869276766 +1 99105U 12031A 17198.01673245 +.00001030 +00000-0 +34740-4 0 9991 +2 99105 006.0268 092.7847 0011809 001.1647 358.8509 14.87765643276906 +1 99105U 12031A 17200.96628195 +.00001022 +00000-0 +33960-4 0 9997 +2 99105 006.0270 071.3796 0011404 044.5053 315.5991 14.87767583277347 +1 99105U 12031A 17202.30698275 +.00000943 +00000-0 +25525-4 0 9992 +2 99105 006.0274 061.6515 0011318 064.5565 295.5734 14.87768069277540 +1 99105U 12031A 17203.11140138 +.00000952 +00000-0 +26433-4 0 9991 +2 99105 006.0274 055.8130 0011263 076.7702 283.3683 14.87768621277660 +1 99105U 12031A 17204.11692308 +.00000944 +00000-0 +25685-4 0 9992 +2 99105 006.0275 048.5166 0011228 092.0292 268.1123 14.87769212277819 +1 99105U 12031A 17205.05540774 +.00000939 +00000-0 +25126-4 0 9999 +2 99105 006.0274 041.7031 0011279 106.2569 253.8801 14.87769710277957 +1 99105U 12031A 17206.46313244 +.00000929 +00000-0 +24056-4 0 9995 +2 99105 006.0273 031.4870 0011392 127.4645 232.6520 14.87770348278165 +1 99105U 12031A 17206.99940770 +.00000922 +00000-0 +23256-4 0 9997 +2 99105 006.0274 027.5958 0011409 135.5193 224.5853 14.87770562278242 +1 99105U 12031A 17207.40161385 +.00000930 +00000-0 +24084-4 0 9990 +2 99105 006.0274 024.6772 0011464 141.4652 218.6296 14.87770807278305 +1 99105U 12031A 17208.34009421 +.00000904 +00000-0 +21302-4 0 9997 +2 99105 006.0274 017.8673 0011575 155.3356 204.7326 14.87771097278444 +1 99105U 12031A 17209.27857331 +.00000915 +00000-0 +22508-4 0 9992 +2 99105 006.0274 011.0559 0011700 169.0684 190.9699 14.87771598278589 +1 99105U 12031A 17210.01594982 +.00000904 +00000-0 +21336-4 0 9992 +2 99105 006.0272 005.7061 0011817 179.7607 180.2528 14.87771839278690 +1 99105U 12031A 17212.02697476 +.00000902 +00000-0 +21070-4 0 9990 +2 99105 006.0271 351.1098 0012105 208.4229 151.5239 14.87772682278998 +1 99105U 12031A 17213.09952223 +.00000903 +00000-0 +21177-4 0 9991 +2 99105 006.0270 343.3264 0012218 223.4606 136.4559 14.87773114279153 +1 99105U 12031A 17215.37869064 +.00000935 +00000-0 +24567-4 0 9995 +2 99105 006.0276 326.7958 0012409 255.2146 104.6607 14.87774164279492 +1 99105U 12031A 17216.31717247 +.00000938 +00000-0 +24881-4 0 9997 +2 99105 006.0271 319.9832 0012488 267.9804 091.8894 14.87774509279634 +1 99105U 12031A 17217.25565539 +.00000961 +00000-0 +27326-4 0 9999 +2 99105 006.0271 313.1712 0012467 280.9084 078.9642 14.87775079279774 +1 99105U 12031A 17218.26117367 +.00000979 +00000-0 +29220-4 0 9997 +2 99105 006.0269 305.8717 0012373 294.9590 064.9254 14.87775728279929 +1 99105U 12031A 17219.06558940 +.00000953 +00000-0 +26454-4 0 9990 +2 99105 006.0269 300.0330 0012312 306.0737 053.8252 14.87775998280047 +1 99105U 12031A 17220.27221386 +.00000944 +00000-0 +25512-4 0 9990 +2 99105 006.0266 291.2738 0012198 322.7334 037.1949 14.87776562280223 +1 99105U 12031A 17222.14918658 +.00000921 +00000-0 +23043-4 0 9995 +2 99105 006.0263 277.6485 0011991 349.2937 010.6937 14.87777426280504 +1 99105U 12031A 17223.42284666 +.00000929 +00000-0 +23952-4 0 9999 +2 99105 006.0261 268.4032 0011781 007.6000 352.4307 14.87778117280697 +1 99105U 12031A 17224.42836772 +.00000929 +00000-0 +23939-4 0 9996 +2 99105 006.0262 261.1058 0011809 022.0157 338.0479 14.87778488280849 +1 99105U 12031A 17225.36685342 +.00000935 +00000-0 +24668-4 0 9992 +2 99105 006.0259 254.2999 0011575 035.7426 324.3477 14.87779086280988 +1 99105U 12031A 17226.37237155 +.00000928 +00000-0 +23893-4 0 9990 +2 99105 006.0258 246.9994 0011478 050.7724 309.3423 14.87779470281139 +1 99105U 12031A 17227.51195665 .00000943 00000-0 25545-4 0 9999 +2 99105 6.0258 238.7249 0011428 67.9323 292.2019 14.87780058281305 +1 99105U 12031A 17228.51747106 .00000965 00000-0 27865-4 0 9994 +2 99105 6.0258 231.4244 0011389 83.1355 277.0069 14.87780642281452 +1 99105U 12031A 17229.72408554 .00000978 00000-0 29236-4 0 9997 +2 99105 6.0258 222.6634 0011402 101.4516 258.6894 14.87781332281634 +1 99105U 12031A 17231.39993565 +.00001056 +00000-0 +37510-4 0 9996 +2 99105 006.0258 210.5027 0011509 126.6165 233.5023 14.87782659281889 +1 99105U 12031A 17232.33840895 .00001074 00000-0 39409-4 0 9997 +2 99105 6.0258 203.6898 0011600 140.5632 219.5342 14.87783310282024 +1 99105U 12031A 17233.27688076 +.00001084 +00000-0 +40424-4 0 9992 +2 99105 006.0258 196.8763 0011706 154.5140 205.5567 14.87783948282160 +1 99105U 12031A 17234.48348586 +.00001070 +00000-0 +38937-4 0 9996 +2 99105 006.0258 188.1155 0011865 172.0832 187.9485 14.87784652282343 +1 99105U 12031A 17235.48898932 .00001058 00000-0 37655-4 0 9994 +2 99105 6.0258 180.8153 0012009 186.4054 173.5921 14.87785283282499 +1 99105U 12031A 17236.69559276 .00001069 00000-0 38866-4 0 9990 +2 99105 6.0259 172.0548 0012177 203.5680 156.3890 14.87786231282671 +1 99105U 12031A 17237.36592829 +.00001052 +00000-0 +36991-4 0 9995 +2 99105 006.0259 167.1905 0012262 213.0739 146.8622 14.87786683282778 +1 99105U 12031A 17237.96923018 .00001031 00000-0 34777-4 0 9996 +2 99105 6.0260 162.8101 0012341 221.4384 138.4809 14.87787000282865 +1 99105U 12031A 17239.10880156 +.00000997 +00000-0 +31073-4 0 9995 +2 99105 006.0261 154.5382 0012430 237.1792 122.7139 14.87787574283039 +1 99105U 12031A 17240.38244095 .00001035 00000-0 35208-4 0 9991 +2 99105 6.0263 145.2914 0012475 254.7079 105.1671 14.87788700283220 +1 99105U 12031A 17242.52752193 .00001129 00000-0 45213-4 0 9993 +2 99105 6.0265 129.7195 0012450 284.0938 75.7808 14.87790732283540 +1 99105U 12031A 17245.47701693 .00001130 00000-0 45274-4 0 9998 +2 99105 6.0267 108.3159 0012256 325.1105 34.8221 14.87792898283980 +1 99105U 12031A 17246.48252634 .00001172 00000-0 49753-4 0 9998 +2 99105 6.0267 101.0166 0012111 339.2748 20.6890 14.87794156284139 +1 99105U 12031A 17247.01879894 +.00001162 +00000-0 +48778-4 0 9996 +2 99105 006.0267 097.1261 0012003 346.7293 013.2520 14.87794529284219 +1 99105U 12031A 17248.42651124 +.00001318 +00000-0 +65364-4 0 9999 +2 99105 006.0269 086.9072 0011751 006.9723 353.0568 14.87796944284422 +1 99105U 12031A 17249.09684960 +.00001343 +00000-0 +68114-4 0 9991 +2 99105 006.0269 082.0397 0011626 016.7100 343.3411 14.87797807284527 +1 99105U 12031A 17250.43752541 +.00001373 +00000-0 +71326-4 0 9994 +2 99105 006.0270 072.3101 0011460 036.5493 323.5417 14.87799556284726 +1 99105U 12031A 17252.51556239 .00002030 00000-0 14129-3 0 9999 +2 99105 6.0271 57.2273 0011264 67.5357 292.5964 14.87806554285038 +1 99105U 12031A 17253.25292847 +.00001538 +00000-0 +88904-4 0 9998 +2 99105 006.0271 051.8751 0011220 078.7657 281.3733 14.87805080285143 +1 99105U 12031A 17254.05732592 .00001449 00000-0 79381-4 0 9993 +2 99105 6.0271 46.0366 0011196 91.0117 269.1294 14.87805921285265 +1 99105U 12031A 17255.86721468 .00001343 00000-0 68033-4 0 9998 +2 99105 6.0271 32.9008 0011266 118.5152 241.6112 14.87808123285535 +1 99105U 12031A 17258.28038860 +.00001249 +00000-0 +57994-4 0 9993 +2 99105 006.0274 015.3883 0011537 154.7479 205.3214 14.87810918285893 +1 99105U 12031A 17259.28587506 +.00001248 +00000-0 +57896-4 0 9990 +2 99105 006.0275 008.0907 0011678 169.4554 190.5820 14.87811972286040 +1 99105U 12031A 17260.09026327 +.00001203 +00000-0 +53090-4 0 9993 +2 99105 006.0275 002.2518 0011786 181.1281 178.8821 14.87812570286160 +1 99105U 12031A 17261.02871547 .00001191 00000-0 51824-4 0 9998 +2 99105 6.0275 355.4397 0011922 194.5188 165.4598 14.87813530286300 +1 99105U 12031A 17262.03419979 .00001154 00000-0 47784-4 0 9996 +2 99105 6.0275 348.1413 0012062 208.7839 151.1623 14.87814377286452 +1 99105U 12031A 17263.10671708 +.00001119 +00000-0 +44078-4 0 9999 +2 99105 006.0274 340.3586 0012197 223.7876 136.1285 14.87815233286610 +1 99105U 12031A 17264.31330249 +.00001071 +00000-0 +38992-4 0 9990 +2 99105 006.0273 331.6144 0012344 240.3978 119.4921 14.87815933286791 +1 99105U 12031A 17265.18472391 +.00001079 +00000-0 +39777-4 0 9991 +2 99105 006.0272 325.2912 0012379 252.6458 107.2317 14.87816767286924 +1 99105U 12031A 17266.32427577 .00001112 00000-0 43275-4 0 9990 +2 99105 6.0271 317.0200 0012431 268.2624 91.6081 14.87817931287096 +1 99105U 12031A 17267.19569877 +.00001137 +00000-0 +46033-4 0 9998 +2 99105 006.0269 310.6940 0012385 279.9989 079.8742 14.87818758287229 +1 99105U 12031A 17268.20118732 +.00001223 +00000-0 +55106-4 0 9998 +2 99105 006.0268 303.3946 0012348 293.8261 066.0573 14.87820102287372 +1 99105U 12031A 17269.40777497 .00001249 00000-0 57962-4 0 9990 +2 99105 6.0267 294.6361 0012266 310.4491 49.4569 14.87821274287555 +1 99105U 12031A 17270.27920021 +.00001372 +00000-0 +70984-4 0 9996 +2 99105 006.0269 288.3138 0012244 322.9649 036.9635 14.87822748287684 +1 99105U 12031A 17271.41875457 .00001547 00000-0 89638-4 0 9995 +2 99105 6.0268 280.0417 0012102 338.9927 20.9705 14.87825235287855 +1 99105U 12031A 17272.35721106 +.00001460 +00000-0 +80457-4 0 9994 +2 99105 006.0267 273.2287 0011984 352.3039 007.6906 14.87825984287995 +1 99105U 12031A 17273.36269925 +.00001439 +00000-0 +78167-4 0 9993 +2 99105 006.0267 265.9300 0011874 007.2044 352.8255 14.87827333288149 +1 99105U 12031A 17274.36818606 +.00001395 +00000-0 +73567-4 0 9990 +2 99105 006.0266 258.6305 0011701 021.4829 338.5790 14.87828594288291 +1 99105U 12031A 17275.03850966 +.00001387 +00000-0 +72706-4 0 9995 +2 99105 006.0266 253.7645 0011643 031.6184 328.4643 14.87829635288399 +1 99105U 12031A 17276.17805803 .00001297 00000-0 63089-4 0 9995 +2 99105 6.0264 245.4924 0011532 48.5113 311.6005 14.87830870288569 +1 99105U 12031A 17277.31760339 .00001303 00000-0 63748-4 0 9996 +2 99105 6.0263 237.2244 0011341 65.1681 294.9627 14.87832753288739 +1 99105U 12031A 17279.19449424 .00001297 00000-0 63166-4 0 9999 +2 99105 6.0261 223.5969 0011309 93.5810 266.5612 14.87835391289014 +1 99105U 12031A 17280.26699967 .00001284 00000-0 61719-4 0 9996 +2 99105 6.0261 215.8100 0011324 109.8436 250.2914 14.87836665289174 +1 99105U 12031A 17281.13840831 +.00001275 +00000-0 +60715-4 0 9996 +2 99105 006.0261 209.4829 0011402 122.9889 237.1336 14.87837701289305 +1 99105U 12031A 17282.07684670 .00001263 00000-0 59508-4 0 9992 +2 99105 6.0260 202.6695 0011504 136.9141 223.1889 14.87838766289444 +1 99105U 12031A 17283.28340889 .00001222 00000-0 55050-4 0 9990 +2 99105 6.0261 193.9112 0011647 154.8360 205.2337 14.87839722289623 +1 99105U 12031A 17284.42293733 .00001262 00000-0 59256-4 0 9998 +2 99105 6.0261 185.6372 0011813 171.4346 188.5985 14.87841073289795 +1 99105U 12031A 17285.22730892 +.00001300 +00000-0 +63321-4 0 9990 +2 99105 006.0260 179.7942 0011917 183.1518 176.8536 14.87842096289919 +1 99105U 12031A 17286.16574234 .00001293 00000-0 62611-4 0 9992 +2 99105 6.0260 172.9805 0012060 196.5884 163.3850 14.87843084290051 +1 99105U 12031A 17287.50636144 .00001289 00000-0 62086-4 0 9991 +2 99105 6.0260 163.2474 0012248 215.3258 144.6059 14.87844474290258 +1 99105U 12031A 17288.10964218 +.00001241 +00000-0 +56950-4 0 9995 +2 99105 006.0259 158.8777 0012358 223.7062 136.2087 14.87844720290342 +1 99105U 12031A 17289.11510688 .00001208 00000-0 53473-4 0 9995 +2 99105 6.0260 151.5768 0012454 237.6098 122.2825 14.87845575290494 +1 99105U 12031A 17290.12057235 .00001182 00000-0 50677-4 0 9996 +2 99105 6.0260 144.2759 0012518 251.4248 108.4521 14.87846444290647 +1 99105U 12031A 17292.33259965 .00001125 00000-0 44602-4 0 9993 +2 99105 6.0262 128.2145 0012525 281.6771 78.1952 14.87848312290977 +1 99105U 12031A 17293.13697494 +.00001118 +00000-0 +43905-4 0 9992 +2 99105 006.0260 122.3763 0012493 292.9066 066.9744 14.87848999291091 +1 99105U 12031A 17294.34353921 .00001129 00000-0 45060-4 0 9996 +2 99105 6.0262 113.6156 0012393 309.6312 50.2724 14.87849848291279 +1 99105U 12031A 17295.01385283 +.00001137 +00000-0 +45992-4 0 9990 +2 99105 006.0260 108.7519 0012309 318.6941 041.2257 14.87850545291378 +1 99105U 12031A 17296.01932334 .00001158 00000-0 48160-4 0 9999 +2 99105 6.0261 101.4519 0012182 332.7655 27.1835 14.87851517291529 +1 99105U 12031A 17297.42698241 .00001256 00000-0 58606-4 0 9996 +2 99105 6.0262 91.2319 0011975 352.8610 7.1348 14.87853182291737 +1 99105U 12031A 17298.43245230 .00001351 00000-0 68795-4 0 9990 +2 99105 6.0262 83.9335 0011830 7.5015 352.5290 14.87854767291885 +1 99105U 12031A 17299.43792004 +.00001315 +00000-0 +64935-4 0 9999 +2 99105 006.0264 076.6298 0011705 022.0836 337.9797 14.87855714292038 +1 99105U 12031A 17300.64448029 .00001305 00000-0 63860-4 0 9995 +2 99105 6.0267 67.8716 0011534 39.7701 320.3273 14.87857101292211 +1 99105U 12031A 17301.31479101 +.00001235 +00000-0 +56464-4 0 9999 +2 99105 006.0268 063.0054 0011476 050.0025 310.1110 14.87857397292310 +1 99105U 12031A 17302.58837882 .00001208 00000-0 53627-4 0 9991 +2 99105 6.0265 53.7656 0011329 69.2241 290.9101 14.87858748292507 +1 99105U 12031A 17303.12462626 +.00001181 +00000-0 +50775-4 0 9993 +2 99105 006.0270 049.8786 0011267 077.2578 282.8810 14.87859142292580 +1 99105U 12031A 17304.26414810 .00001135 00000-0 45790-4 0 9991 +2 99105 6.0271 41.6095 0011224 94.5379 265.6032 14.87860029292758 +1 99105U 12031A 17305.26960573 +.00001118 +00000-0 +44022-4 0 9992 +2 99105 006.0271 034.3124 0011344 109.5752 250.5602 14.87861002292900 +1 99105U 12031A 17306.34209138 +.00001104 +00000-0 +42528-4 0 9999 +2 99105 006.0272 026.5270 0011423 125.7655 234.3537 14.87861896293062 +1 99105U 12031A 17307.07942392 +.00001108 +00000-0 +42944-4 0 9992 +2 99105 006.0273 021.1749 0011491 136.7168 223.3864 14.87862609293177 +1 99105U 12031A 17308.28596623 .00001131 00000-0 45371-4 0 9991 +2 99105 6.0274 12.4170 0011641 154.6027 205.4674 14.87863778293357 +1 99105U 12031A 17310.02874632 .00001124 00000-0 44624-4 0 9997 +2 99105 6.0275 359.7668 0011864 179.9450 180.0680 14.87865271293618 +1 99105U 12031A 17311.30231483 .00001234 00000-0 56203-4 0 9998 +2 99105 6.0275 350.5223 0012078 198.1509 161.8188 14.87866987293806 +1 99105U 12031A 17312.24073283 +.00001450 +00000-0 +79216-4 0 9993 +2 99105 006.0274 343.7103 0012196 211.4380 148.5019 14.87869159293947 +1 99105U 12031A 17313.71539180 .00001335 00000-0 66984-4 0 9999 +2 99105 6.0276 333.0109 0012252 232.0039 127.8983 14.87869955294169 +1 99105U 12031A 17314.11757181 +.00001257 +00000-0 +58707-4 0 9994 +2 99105 006.0276 330.0930 0012276 237.5459 122.3483 14.87869900294222 +1 99105U 12031A 17315.39114160 +.00001191 +00000-0 +51684-4 0 9993 +2 99105 006.0275 320.8492 0012353 255.1991 104.6768 14.87870802294417 +1 99105U 12031A 17316.32956235 .00001157 00000-0 48018-4 0 9993 +2 99105 6.0275 314.0375 0012369 268.0780 91.7932 14.87871524294554 +1 99105U 12031A 17317.13392390 +.00001157 +00000-0 +47974-4 0 9999 +2 99105 006.0274 308.2006 0012353 279.0611 080.8120 14.87872334294678 +1 99105U 12031A 17318.00531693 .00001066 00000-0 38355-4 0 9999 +2 99105 6.0273 301.8781 0012325 291.2285 68.6528 14.87872573294804 +1 99105U 12031A 17319.34592178 .00001140 00000-0 46198-4 0 9999 +2 99105 6.0271 292.1482 0012236 309.7686 50.1366 14.87874128295000 +1 99105U 12031A 17320.35137584 .00001164 00000-0 48766-4 0 9997 +2 99105 6.0269 284.8490 0012127 323.7620 36.1687 14.87875095295157 +1 99105U 12031A 17321.28980002 +.00001176 +00000-0 +50077-4 0 9992 +2 99105 006.0268 278.0363 0012004 336.9449 023.0142 14.87875909295299 +1 99105U 12031A 17322.02713295 .00001190 00000-0 51532-4 0 9996 +2 99105 6.0267 272.6822 0011886 347.5924 12.3912 14.87876637295409 +1 99105U 12031A 17323.03258784 .00001101 00000-0 42122-4 0 9992 +2 99105 6.0267 265.3850 0011758 1.9411 358.0764 14.87876903295556 +1 99105U 12031A 17324.03804099 .00001082 00000-0 40155-4 0 9991 +2 99105 6.0265 258.0851 0011611 16.4259 343.6246 14.87877715295707 +1 99105U 12031A 17325.37864328 .00001206 00000-0 53368-4 0 9998 +2 99105 6.0264 248.3521 0011457 36.2710 323.8194 14.87879673295901 +1 99105U 12031A 17326.11597367 +.00001188 +00000-0 +51467-4 0 9994 +2 99105 006.0264 242.9999 0011383 047.4643 312.6446 14.87880254296013 +1 99105U 12031A 17326.98736196 .00001197 00000-0 52348-4 0 9993 +2 99105 6.0263 236.6728 0011286 60.4946 299.6307 14.87881237296141 +1 99105U 12031A 17328.05983919 .00001134 00000-0 45650-4 0 9999 +2 99105 6.0264 228.8891 0011225 77.1429 282.9953 14.87881786296309 +1 99105U 12031A 17329.19934229 .00001120 00000-0 44171-4 0 9997 +2 99105 6.0263 220.6148 0011230 94.4879 265.6533 14.87882837296478 +1 99105U 12031A 17330.07072503 +.00001120 +00000-0 +44142-4 0 9992 +2 99105 006.0262 214.2878 0011267 107.7681 252.3678 14.87883725296600 +1 99105U 12031A 17331.07616567 .00001127 00000-0 44912-4 0 9999 +2 99105 6.0258 206.9905 0011346 122.8788 237.2434 14.87884634296756 +1 99105U 12031A 17332.21566161 .00001065 00000-0 38303-4 0 9998 +2 99105 6.0257 198.7155 0011464 139.8330 220.2647 14.87885348296922 +1 99105U 12031A 17333.15406846 +.00001050 +00000-0 +36692-4 0 9998 +2 99105 006.0256 191.9007 0011591 153.7349 206.3368 14.87886059297069 +1 99105U 12031A 17334.42761883 .00001095 00000-0 41458-4 0 9999 +2 99105 6.0256 182.6527 0011868 172.3063 187.7249 14.87887340297255 +1 99105U 12031A 17335.16493700 +.00001066 +00000-0 +38303-4 0 9994 +2 99105 006.0256 177.2986 0012006 182.9790 177.0268 14.87887696297362 +1 99105U 12031A 17336.23739949 .00001012 00000-0 32538-4 0 9996 +2 99105 6.0257 169.5107 0012194 198.2046 161.7646 14.87888118297522 +1 99105U 12031A 17338.11420897 .00000962 00000-0 27257-4 0 9999 +2 99105 6.0257 155.8822 0012426 224.4723 135.4407 14.87889115297808 +1 99105U 12031A 17339.32073120 .00000960 00000-0 27025-4 0 9997 +2 99105 6.0257 147.1197 0012331 241.2357 118.6532 14.87889777297987 +1 99105U 12031A 17340.32616621 .00000994 00000-0 30676-4 0 9991 +2 99105 6.0257 139.8212 0012597 254.7505 105.1230 14.87890750298138 +1 99105U 12031A 17341.26457380 +.00000984 +00000-0 +29528-4 0 9993 +2 99105 006.0258 133.0069 0012626 267.6436 092.2247 14.87891267298277 +1 99105U 12031A 17342.13595327 .00000996 00000-0 30885-4 0 9995 +2 99105 6.0259 126.6795 0012624 279.5872 80.2830 14.87891865298400 +1 99105U 12031A 17343.27545512 .00001172 00000-0 49470-4 0 9999 +2 99105 6.0284 118.4288 0013041 293.3618 66.5162 14.87894034298574 +1 99105U 12031A 17344.14683214 +.00000962 +00000-0 +27279-4 0 9997 +2 99105 006.0261 112.0777 0012527 307.1008 052.7976 14.87892905298707 +1 99105U 12031A 17345.15227713 .00000936 00000-0 24551-4 0 9998 +2 99105 6.0271 104.7925 0012319 320.3929 39.5300 14.87893140298851 +1 99105U 12031A 17346.15771949 .00000977 00000-0 28914-4 0 9996 +2 99105 6.0264 97.4908 0012066 335.3396 24.6156 14.87893834299007 +1 99105U 12031A 17347.43127947 .00000976 00000-0 28849-4 0 9997 +2 99105 6.0266 88.2444 0011879 353.3980 6.5992 14.87894507299197 +1 99105U 12031A 17348.16860355 +.00000979 +00000-0 +29131-4 0 9993 +2 99105 006.0268 082.8919 0011766 004.1506 355.8720 14.87894984299307 +1 99105U 12031A 17349.03998604 .00000974 00000-0 28608-4 0 9992 +2 99105 6.0269 76.5659 0011656 16.8442 343.2073 14.87895488299432 +1 99105U 12031A 17350.31354398 +.00000980 +00000-0 +29244-4 0 9995 +2 99105 006.0270 067.3197 0011481 035.6288 324.4606 14.87896187299623 +1 99105U 12031A 17351.05086596 +.00000952 +00000-0 +26315-4 0 9994 +2 99105 006.0271 061.9670 0011401 046.5670 313.5407 14.87896487299731 +1 99105U 12031A 17352.32442016 .00001034 00000-0 35061-4 0 9998 +2 99105 6.0273 52.7216 0011300 65.7427 294.3882 14.87897645299929 +1 99105U 12031A 17353.46391382 .00001008 00000-0 32304-4 0 9991 +2 99105 6.0275 44.4510 0011269 82.8681 277.2729 14.87898161300097 +1 99105U 12031A 17354.33528945 +.00000995 +00000-0 +30922-4 0 9993 +2 99105 006.0276 038.1274 0011359 096.2160 263.9263 14.87898685300227 +1 99105U 12031A 17355.13963568 .00000985 00000-0 29754-4 0 9996 +2 99105 6.0279 32.2935 0011414 107.5940 252.5436 14.87899079300340 +1 99105U 12031A 17356.21209549 .00000968 00000-0 27950-4 0 9995 +2 99105 6.0281 24.5166 0011424 123.8671 236.2546 14.87899697300501 +1 99105U 12031A 17357.41860932 .00000954 00000-0 26493-4 0 9993 +2 99105 6.0281 15.7596 0011524 141.8869 218.2075 14.87900268300683 +1 99105U 12031A 17358.02186547 +.00000966 +00000-0 +27797-4 0 9995 +2 99105 006.0282 011.3812 0011607 150.5913 209.4869 14.87900665300770 +1 99105U 12031A 17359.16134777 .00000990 00000-0 30310-4 0 9991 +2 99105 6.0282 3.1094 0011736 167.3151 192.7274 14.87901414300940 +1 99105U 12031A 17360.03271613 +.00000992 +00000-0 +30464-4 0 9992 +2 99105 006.0282 356.7847 0011841 179.9055 180.1076 14.87901870301074 +1 99105U 12031A 17361.37328231 +.00001016 +00000-0 +33060-4 0 9992 +2 99105 006.0282 347.0564 0011931 199.2125 160.7554 14.87902808301275 +1 99105U 12031A 17362.37870689 .00000994 00000-0 30666-4 0 9995 +2 99105 6.0278 339.7602 0012195 213.7225 146.2128 14.87903290301426 +1 99105U 12031A 17363.45116021 .00000958 00000-0 26822-4 0 9997 +2 99105 6.0277 331.9755 0012277 228.6467 131.2606 14.87903699301587 +1 99105U 12031A 17364.18847239 +.00000949 +00000-0 +25822-4 0 9999 +2 99105 006.0276 326.6238 0012350 238.8201 121.0716 14.87904048301694 +1 99105U 12031A 17364.99281422 .00000905 00000-0 21197-4 0 9997 +2 99105 6.0274 320.7878 0012349 249.5524 110.3278 14.87904332301816 +1 99105U 12031A 18001.13229887 .00000936 00000-0 24440-4 0 9997 +2 99105 6.0271 312.5167 0012403 265.2433 94.6279 14.87904931301986 +1 99105U 12031A 18002.00367007 +.00000951 +00000-0 +26115-4 0 9999 +2 99105 006.0269 306.1896 0012400 277.1909 082.6810 14.87905439302116 +1 99105U 12031A 18003.27721475 +.00000946 +00000-0 +25567-4 0 9996 +2 99105 006.0268 296.9443 0012286 294.8787 065.0065 14.87906092302305 +1 99105U 12031A 18004.28264589 +.00000939 +00000-0 +24863-4 0 9996 +2 99105 006.0267 289.6447 0012200 308.8187 051.0853 14.87906565302454 +1 99105U 12031A 18005.28807833 .00000933 00000-0 24200-4 0 9991 +2 99105 6.0267 282.3461 0012160 322.9800 36.9490 14.87906979302609 +1 99105U 12031A 18006.36053983 .00000938 00000-0 24711-4 0 9992 +2 99105 6.0265 274.5594 0012019 338.0970 21.8645 14.87907514302764 +1 99105U 12031A 18007.03082854 +.00000932 +00000-0 +24131-4 0 9996 +2 99105 006.0264 269.6932 0011914 347.5973 012.3862 14.87907779302865 +1 99105U 12031A 18009.10872263 .00000932 00000-0 24177-4 0 9997 +2 99105 6.0259 254.6055 0011613 17.6145 342.4386 14.87908634303175 +1 99105U 12031A 18011.38769877 +.00000920 +00000-0 +22901-4 0 9996 +2 99105 006.0256 238.0554 0011341 051.4710 308.6435 14.87909456303514 +1 99105U 12031A 18012.59421250 .00000939 00000-0 24901-4 0 9992 +2 99105 6.0255 229.2938 0011249 69.7200 290.4137 14.87910035303695 +1 99105U 12031A 18014.13586490 .00000956 00000-0 26757-4 0 9992 +2 99105 6.0254 218.0995 0011202 92.9887 267.1524 14.87910891303928 +1 99105U 12031A 18016.34779369 .00000921 00000-0 22936-4 0 9993 +2 99105 6.0254 202.0365 0011330 126.4496 233.6678 14.87911622304251 +1 99105U 12031A 18018.42565932 .00000934 00000-0 24330-4 0 9997 +2 99105 6.0253 186.9521 0011554 157.6023 202.4611 14.87912663304565 +1 99105U 12031A 18020.43649205 .00000950 00000-0 25951-4 0 9995 +2 99105 6.0254 172.3497 0011837 186.7994 173.1974 14.87913778304868 +1 99105U 12031A 18021.03974218 +.00000939 +00000-0 +24787-4 0 9990 +2 99105 006.0254 167.9690 0011918 195.4857 164.4907 14.87913860304959 +1 99105U 12031A 18022.11218601 +.00000949 +00000-0 +25831-4 0 9999 +2 99105 006.0256 160.1807 0012071 210.6853 149.2570 14.87914400305110 +1 99105U 12031A 18024.52518971 .00000914 00000-0 22157-4 0 9999 +2 99105 6.0261 142.6718 0012363 244.3803 115.5048 14.87915281305478 +1 99105U 12031A 18025.12844104 +.00000919 +00000-0 +22677-4 0 9998 +2 99105 006.0260 138.2927 0012372 252.3897 107.4880 14.87915692305562 +1 99105U 12031A 18026.26791690 +.00000938 +00000-0 +24701-4 0 9996 +2 99105 006.0261 130.0200 0012406 267.9921 091.8787 14.87916368305737 +1 99105U 12031A 18027.27333807 .00000936 00000-0 24405-4 0 9999 +2 99105 6.0263 122.7205 0012439 281.7000 78.1733 14.87916858305889 +1 99105U 12031A 18028.27876127 .00000933 00000-0 24192-4 0 9997 +2 99105 6.0266 115.4261 0012366 295.8698 64.0156 14.87917431306030 +1 99105U 12031A 18029.01607170 +.00000944 +00000-0 +25272-4 0 9992 +2 99105 006.0267 110.0727 0012302 306.1037 053.7953 14.87917863306143 +1 99105U 12031A 18030.55772284 +.00000920 +00000-0 +22763-4 0 9998 +2 99105 006.0268 098.8814 0012184 327.5192 032.4188 14.87918423306373 +1 99105U 12031A 18031.16097792 +.00000924 +00000-0 +23182-4 0 9990 +2 99105 006.0270 094.5015 0012115 336.0608 023.8957 14.87918714306464 +1 99105U 12031A 18032.43451678 +.00000920 +00000-0 +22832-4 0 9995 +2 99105 006.0271 085.2559 0011962 354.0748 005.9239 14.87919250306650 +1 99105U 12031A 18034.17725326 .00000925 00000-0 23408-4 0 9997 +2 99105 6.0274 72.6030 0011724 19.3057 340.7515 14.87920067306911 +1 99105U 12031A 18035.11564904 +.00000919 +00000-0 +22709-4 0 9991 +2 99105 006.0275 065.7899 0011607 033.0739 327.0115 14.87920390307057 +1 99105U 12031A 18036.18810004 +.00000935 +00000-0 +24412-4 0 9995 +2 99105 006.0276 058.0054 0011472 049.0388 311.0733 14.87920990307216 +1 99105U 12031A 18037.26054945 +.00000926 +00000-0 +23503-4 0 9999 +2 99105 006.0277 050.2199 0011397 065.1189 295.0124 14.87921358307376 +1 99105U 12031A 18038.06488599 +.00000924 +00000-0 +23266-4 0 9998 +2 99105 006.0279 044.3855 0011345 077.3161 282.8236 14.87921662307498 +1 99105U 12031A 18039.33841454 .00000944 00000-0 25469-4 0 9991 +2 99105 6.0277 35.1376 0011325 96.6557 263.4861 14.87922290307681 +1 99105U 12031A 18040.27680304 +.00000940 +00000-0 +24986-4 0 9990 +2 99105 006.0279 028.3285 0011370 110.7856 249.3491 14.87922645307828 +1 99105U 12031A 18041.34924330 +.00000936 +00000-0 +24514-4 0 9991 +2 99105 006.0277 020.5390 0011436 127.0840 233.0335 14.87923151307981 +1 99105U 12031A 18042.08654610 +.00000926 +00000-0 +23452-4 0 9993 +2 99105 006.0277 015.1866 0011524 138.1020 221.9992 14.87923284308095 +1 99105U 12031A 18043.15898474 .00000931 00000-0 24006-4 0 9992 +2 99105 6.0277 7.4023 0011646 153.9481 206.1235 14.87923778308256 +1 99105U 12031A 18045.37088702 +.00001008 +00000-0 +32112-4 0 9996 +2 99105 006.0279 351.3573 0011972 185.6626 174.3368 14.87925395308580 +1 99105U 12031A 18046.51035037 +.00001026 +00000-0 +34051-4 0 9990 +2 99105 006.0275 343.0869 0012136 201.8398 158.1213 14.87926060308750 +1 99105U 12031A 18047.85089590 .00001114 00000-0 43275-4 0 9997 +2 99105 6.0227 333.4127 0012473 222.4157 137.4467 14.87926452308940 +1 99105U 12031A 18048.58819601 +.00001010 +00000-0 +32287-4 0 9998 +2 99105 006.0274 328.0034 0012362 230.7979 129.1052 14.87927036309061 +1 99105U 12031A 18049.59360652 +.00001010 +00000-0 +32255-4 0 9992 +2 99105 006.0272 320.7060 0012428 244.5874 115.2968 14.87927573309212 +1 99105U 12031A 18050.06279834 +.00001013 +00000-0 +32594-4 0 9996 +2 99105 006.0271 317.3007 0012473 251.0176 108.8600 14.87927872309282 +1 99105U 12031A 18051.26929234 .00001001 00000-0 31328-4 0 9996 +2 99105 6.0270 308.5409 0012496 267.6002 92.2696 14.87928523309460 +1 99105U 12031A 18053.54823075 +.00000953 +00000-0 +26222-4 0 9997 +2 99105 006.0266 291.9957 0012358 298.9992 060.8899 14.87929528309800 +1 99105U 12031A 18054.48661843 +.00000959 +00000-0 +26952-4 0 9990 +2 99105 006.0265 285.1818 0012280 312.0179 047.8904 14.87930099309941 +1 99105U 12031A 18055.15689572 +.00000960 +00000-0 +27050-4 0 9996 +2 99105 006.0264 280.3151 0012227 321.3718 038.5536 14.87930505310043 +1 99105U 12031A 18056.09528447 .00000976 00000-0 28695-4 0 9996 +2 99105 6.0263 273.5012 0012100 334.5330 25.4203 14.87931090310189 +1 99105U 12031A 18057.10070125 .00000977 00000-0 28831-4 0 9997 +2 99105 6.0261 266.2008 0011943 348.8159 11.1705 14.87931628310335 +1 99105U 12031A 18058.10611780 .00000977 00000-0 28916-4 0 9990 +2 99105 6.0261 258.9001 0011778 3.2660 356.7545 14.87932177310487 +1 99105U 12031A 18059.11153376 +.00000985 +00000-0 +29697-4 0 9990 +2 99105 006.0260 251.5997 0011606 017.9408 342.1130 14.87932801310631 +1 99105U 12031A 18060.38505914 +.00000966 +00000-0 +27774-4 0 9993 +2 99105 006.0259 242.3516 0011418 036.8065 323.2847 14.87933409310825 +1 99105U 12031A 18061.45750010 .00000974 00000-0 28552-4 0 9995 +2 99105 6.0258 234.5651 0011300 52.8506 307.2654 14.87934041310983 +1 99105U 12031A 18062.39588440 +.00000968 +00000-0 +27935-4 0 9992 +2 99105 006.0257 227.7509 0011232 066.9971 293.1342 14.87934467311120 +1 99105U 12031A 18063.06615755 .00000975 00000-0 28692-4 0 9996 +2 99105 6.0255 222.8832 0011202 76.9844 283.1535 14.87934920311228 +1 99105U 12031A 18064.07156605 .00000975 00000-0 28730-4 0 9995 +2 99105 6.0256 215.5819 0011170 92.3565 267.7843 14.87935438311370 +1 99105U 12031A 18065.54616194 +.00000981 +00000-0 +29392-4 0 9994 +2 99105 006.0257 204.8757 0011224 114.7718 245.3579 14.87936279311598 +1 99105U 12031A 18066.48453932 +.00000993 +00000-0 +30658-4 0 9997 +2 99105 006.0260 198.0657 0011275 128.8759 231.2376 14.87936916311735 +1 99105U 12031A 18067.48994144 +.00000999 +00000-0 +31179-4 0 9998 +2 99105 006.0260 190.7650 0011377 143.9410 216.1487 14.87937515311889 +1 99105U 12031A 18068.49534395 +.00000992 +00000-0 +30424-4 0 9993 +2 99105 006.0261 183.4714 0011513 158.7117 201.3491 14.87937952312034 +1 99105U 12031A 18069.76885041 .00000996 00000-0 30861-4 0 9990 +2 99105 6.0262 174.2247 0011675 177.2716 182.7477 14.87938692312223 +1 99105U 12031A 18071.04235548 .00001030 00000-0 34412-4 0 9991 +2 99105 6.0260 164.9770 0011896 195.5655 164.4108 14.87939696312412 +1 99105U 12031A 18073.65639337 .00001052 00000-0 36785-4 0 9993 +2 99105 6.0261 145.9959 0012215 232.4094 127.4925 14.87941391312809 +1 99105U 12031A 18074.72882023 .00001064 00000-0 38065-4 0 9990 +2 99105 6.0262 138.2084 0012297 247.3021 112.5807 14.87942171312967 +1 99105U 12031A 18075.13098056 +.00001063 +00000-0 +37884-4 0 9998 +2 99105 006.0263 135.2879 0012314 252.7891 107.0890 14.87942438313027 +1 99105U 12031A 18076.33746238 .00001113 00000-0 43165-4 0 9998 +2 99105 6.0263 126.5266 0012357 269.3915 90.4798 14.87943587313206 +1 99105U 12031A 18077.40989210 .00001121 00000-0 44014-4 0 9999 +2 99105 6.0265 118.7406 0012353 284.1541 75.7216 14.87944491313366 +1 99105U 12031A 18078.08016116 +.00001114 +00000-0 +43321-4 0 9992 +2 99105 006.0266 113.8730 0012325 293.4673 066.4161 14.87944986313467 +1 99105U 12031A 18079.55475501 .00001124 00000-0 44398-4 0 9998 +2 99105 6.0265 103.1672 0012197 313.7827 46.1294 14.87946215313680 +1 99105U 12031A 18080.02394403 +.00001115 +00000-0 +43427-4 0 9995 +2 99105 006.0266 099.7595 0012135 320.3740 039.5503 14.87946558313752 +1 99105U 12031A 18081.49853909 .00001092 00000-0 41016-4 0 9994 +2 99105 6.0268 89.0528 0011967 341.0971 18.8714 14.87947608313975 +1 99105U 12031A 18083.04016394 +.00001111 +00000-0 +43026-4 0 9991 +2 99105 006.0268 077.8746 0011747 003.2499 356.7706 14.87948874314208 +1 99105U 12031A 18083.97854195 .00001096 00000-0 41488-4 0 9994 +2 99105 6.0268 71.0619 0011623 16.8988 343.1527 14.87949511314344 +1 99105U 12031A 18084.98394597 +.00001075 +00000-0 +39251-4 0 9993 +2 99105 006.0269 063.7623 0011477 031.6932 328.3887 14.87950120314497 +1 99105U 12031A 18084.98394601 +.00001073 +00000-0 +39083-4 0 9990 +2 99105 006.0269 063.7623 0011477 031.6932 328.3887 14.87950095314498 +1 99105U 12031A 18087.12880353 .00001049 00000-0 36482-4 0 9996 +2 99105 6.0273 48.1918 0011397 63.6301 296.4998 14.87951516314816 +1 99105U 12031A 18088.60338788 .00001037 00000-0 35231-4 0 9998 +2 99105 6.0275 37.4869 0011286 85.8867 274.2551 14.87952501315033 +1 99105U 12031A 18089.34067844 +.00001014 +00000-0 +32864-4 0 9996 +2 99105 006.0275 032.1344 0011279 097.0581 263.0831 14.87952884315146 +1 99105U 12031A 18090.01094170 .00000999 00000-0 31202-4 0 9990 +2 99105 6.0276 27.2681 0011312 107.2077 252.9290 14.87953217315246 +1 99105U 12031A 18090.94930881 +.00000993 +00000-0 +30589-4 0 9991 +2 99105 006.0276 020.4554 0011378 121.3199 238.8044 14.87953763315385 +1 99105U 12031A 18092.28983033 +.00000996 +00000-0 +30879-4 0 9992 +2 99105 006.0277 010.7231 0011495 141.4193 218.6758 14.87954644315586 +1 99105U 12031A 18093.09414205 .00000998 00000-0 31084-4 0 9993 +2 99105 6.0277 4.8839 0011598 153.3606 206.7119 14.87955136315709 +1 99105U 12031A 18094.50168576 .00000996 00000-0 30803-4 0 9993 +2 99105 6.0277 354.6639 0011787 173.9133 186.1139 14.87955915315917 +1 99105U 12031A 18096.51246085 .00000984 00000-0 29549-4 0 9990 +2 99105 6.0277 340.0650 0012080 202.6646 157.2950 14.87956960316210 +1 99105U 12031A 18097.58487443 .00000979 00000-0 29025-4 0 9996 +2 99105 6.0276 332.2803 0012232 217.7031 142.2240 14.87957559316374 +1 99105U 12031A 18098.05405536 +.00000978 +00000-0 +28911-4 0 9994 +2 99105 006.0276 328.8733 0012281 224.2981 135.6164 14.87957817316444 +1 99105U 12031A 18098.99241792 .00000984 00000-0 29454-4 0 9999 +2 99105 6.0275 322.0604 0012370 237.3109 122.5826 14.87958366316587 +1 99105U 12031A 18100.06483339 .00001000 00000-0 31136-4 0 9996 +2 99105 6.0275 314.2740 0012440 252.1174 107.7598 14.87959009316740 +1 99105U 12031A 18101.67345906 .00001065 00000-0 38068-4 0 9993 +2 99105 6.0271 302.5976 0012508 273.7750 86.0948 14.87960362316989 +1 99105U 12031A 18103.68424353 .00001069 00000-0 38472-4 0 9993 +2 99105 6.0269 287.9986 0012421 301.3310 58.5604 14.87962004317289 +1 99105U 12031A 18104.15342906 +.00001029 +00000-0 +34282-4 0 9990 +2 99105 006.0270 284.5934 0012387 308.1066 051.7946 14.87961763317354 +1 99105U 12031A 18105.02477057 .00001023 00000-0 33607-4 0 9997 +2 99105 6.0269 278.2664 0012292 320.3209 39.6021 14.87962412317487 +1 99105U 12031A 18106.03016540 .00001010 00000-0 32233-4 0 9995 +2 99105 6.0267 270.9666 0012175 334.2792 25.6731 14.87963047317632 +1 99105U 12031A 18107.10258827 .00000996 00000-0 30766-4 0 9991 +2 99105 6.0266 263.1882 0012013 349.2007 10.7863 14.87963644317792 +1 99105U 12031A 18108.44311635 .00001022 00000-0 33563-4 0 9995 +2 99105 6.0265 253.4627 0011789 8.7146 351.3187 14.87964606317994 +1 99105U 12031A 18109.44851025 +.00001014 +00000-0 +32726-4 0 9999 +2 99105 006.0264 246.1618 0011647 023.4211 336.6447 14.87965190318144 +1 99105U 12031A 18110.31985048 +.00000999 +00000-0 +31192-4 0 9997 +2 99105 006.0264 239.8329 0011475 035.9946 324.0955 14.87965572318277 +1 99105U 12031A 18111.45929383 .00001047 00000-0 36257-4 0 9995 +2 99105 6.0263 231.5582 0011325 52.9817 307.1347 14.87966550318449 +1 99105U 12031A 18111.99550222 +.00001049 +00000-0 +36478-4 0 9993 +2 99105 006.0260 227.6672 0011377 061.3914 298.7358 14.87967032318523 +1 99105U 12031A 18113.26899461 +.00001007 +00000-0 +32061-4 0 9994 +2 99105 006.0256 218.4229 0011320 080.6644 279.4765 14.87967598318712 +1 99105U 12031A 18114.54248305 .00000989 00000-0 30182-4 0 9991 +2 99105 6.0254 209.1753 0011307 99.9466 260.1939 14.87968337318902 +1 99105U 12031A 18118.49697963 +.00000974 +00000-0 +28458-4 0 9994 +2 99105 006.0252 180.4579 0011653 158.8440 201.2171 14.87970962319494 +1 99105U 12031A 18119.16723146 .00000954 00000-0 26373-4 0 9994 +2 99105 6.0249 175.5909 0011732 169.0279 191.0106 14.87971142319590 +1 99105U 12031A 18120.17260800 .00000957 00000-0 26681-4 0 9995 +2 99105 6.0250 168.2882 0011865 183.5068 176.4978 14.87971712319749 +1 99105U 12031A 18122.38443724 +.00000956 +00000-0 +26540-4 0 9991 +2 99105 006.0251 152.2268 0012153 214.7917 145.1417 14.87972733320075 +1 99105U 12031A 18123.25576552 +.00000947 +00000-0 +25612-4 0 9993 +2 99105 006.0259 145.9061 0012214 226.8886 133.0221 14.87973156320204 +1 99105U 12031A 18124.59626993 +.00000948 +00000-0 +25634-4 0 9991 +2 99105 006.0259 136.1716 0012316 245.5150 114.3694 14.87973803320405 +1 99105U 12031A 18125.53462444 .00000951 00000-0 25943-4 0 9992 +2 99105 6.0260 129.3576 0012349 258.5007 101.3734 14.87974201320545 +1 99105U 12031A 18126.47297953 +.00001021 +00000-0 +33361-4 0 9999 +2 99105 006.0261 122.5431 0012358 271.3945 088.4768 14.87975094320683 +1 99105U 12031A 18127.14323389 .00001029 00000-0 34228-4 0 9998 +2 99105 6.0262 117.6753 0012319 280.5442 79.3299 14.87975497320783 +1 99105U 12031A 18128.34969301 .00001043 00000-0 35787-4 0 9998 +2 99105 6.0263 108.9143 0012268 297.2382 62.6497 14.87976282320969 +1 99105U 12031A 18129.01994873 +.00001018 +00000-0 +33100-4 0 9994 +2 99105 006.0263 104.0458 0012220 306.4974 053.4029 14.87976473321066 +1 99105U 12031A 18130.42748737 +.00001042 +00000-0 +35650-4 0 9996 +2 99105 006.0263 093.8297 0012088 326.3605 033.5756 14.87977618321277 +1 99105U 12031A 18131.09774380 .00001027 00000-0 34091-4 0 9994 +2 99105 6.0264 88.9621 0011987 335.8267 24.1300 14.87978018321377 +1 99105U 12031A 18132.37123101 .00001061 00000-0 37700-4 0 9993 +2 99105 6.0266 79.7153 0011832 354.0090 5.9897 14.87979205321561 +1 99105U 12031A 18132.97446156 +.00001037 +00000-0 +35211-4 0 9995 +2 99105 006.0267 075.3350 0011726 002.7828 357.2366 14.87979579321659 +1 99105U 12031A 18135.58845910 .00001004 00000-0 31667-4 0 9991 +2 99105 6.0273 56.3623 0011383 40.8039 319.2941 14.87981187322043 +1 99105U 12031A 18136.05763727 +.00000999 +00000-0 +31238-4 0 9995 +2 99105 006.0273 052.9566 0011329 047.8903 312.2188 14.87981446322110 +1 99105U 12031A 18137.26409305 +.00000994 +00000-0 +30682-4 0 9995 +2 99105 006.0274 044.1960 0011229 065.6853 294.4448 14.87982129322291 +1 99105U 12031A 18138.13542076 .00000977 00000-0 28866-4 0 9997 +2 99105 6.0275 37.8699 0011190 78.9092 281.2295 14.87982557322423 +1 99105U 12031A 18139.40889702 .00000958 00000-0 26802-4 0 9995 +2 99105 6.0276 28.6241 0011229 98.5696 261.5705 14.87983209322619 +1 99105U 12031A 18140.28022139 +.00000935 +00000-0 +24403-4 0 9996 +2 99105 006.0277 022.2984 0011288 111.8710 248.2620 14.87983496322745 +1 99105U 12031A 18141.41964329 .00000925 00000-0 23305-4 0 9997 +2 99105 6.0278 14.0259 0011380 129.0320 231.0822 14.87983988322910 +1 99105U 12031A 18145.50814379 +.00000912 +00000-0 +21846-4 0 9994 +2 99105 006.0280 344.3439 0011967 188.9688 171.0227 14.87985777323522 +1 99105U 12031A 18146.51351152 .00000912 00000-0 21803-4 0 9994 +2 99105 6.0279 337.0454 0012101 203.2227 156.7355 14.87986203323670 +1 99105U 12031A 18147.51887977 .00000901 00000-0 20684-4 0 9998 +2 99105 6.0279 329.7467 0012259 217.3210 142.6067 14.87986486323828 +1 99105U 12031A 18148.05507633 +.00000901 +00000-0 +20643-4 0 9991 +2 99105 006.0278 325.8549 0012281 224.8715 135.0421 14.87986741323906 +1 99105U 12031A 18149.26151930 .00000894 00000-0 19896-4 0 9999 +2 99105 6.0278 317.0948 0012378 241.6481 118.2399 14.87987147324084 +1 99105U 12031A 18150.46796471 .00000919 00000-0 22589-4 0 9999 +2 99105 6.0273 308.3431 0012491 257.9279 101.9449 14.87987940324266 +1 99105U 12031A 18151.47333604 +.00000911 +00000-0 +21707-4 0 9999 +2 99105 006.0275 301.0412 0012522 271.7247 088.1447 14.87988303324416 +1 99105U 12031A 18152.54573374 .00000957 00000-0 26624-4 0 9991 +2 99105 6.0274 293.2540 0012499 286.4002 73.4753 14.87989030324576 +1 99105U 12031A 18153.41705804 +.00000975 +00000-0 +28454-4 0 9996 +2 99105 006.0273 286.9283 0012444 298.4363 061.4512 14.87989569324703 +1 99105U 12031A 18155.09268506 +.00000947 +00000-0 +25516-4 0 9990 +2 99105 006.0265 274.7670 0012309 321.3107 038.6140 14.87990082324953 +1 99105U 12031A 18156.63426268 +.00000920 +00000-0 +22710-4 0 9996 +2 99105 006.0262 263.5746 0012102 342.9913 016.9811 14.87990625325188 +1 99105U 12031A 18157.03641364 +.00000919 +00000-0 +22604-4 0 9997 +2 99105 006.0260 260.6569 0012007 348.9436 011.0429 14.87990893325242 +1 99105U 12031A 18158.57799195 +.00000913 +00000-0 +22007-4 0 9994 +2 99105 006.0262 249.4713 0011766 011.0530 348.9857 14.87991710325475 +1 99105U 12031A 18159.51634224 +.00000914 +00000-0 +22108-4 0 9993 +2 99105 006.0260 242.6569 0011635 024.7622 335.3065 14.87992173325616 +1 99105U 12031A 18160.52171671 .00000887 00000-0 19243-4 0 9993 +2 99105 6.0258 235.3554 0011515 39.5860 320.5109 14.87992431325769 +1 99105U 12031A 18161.92923824 +.00000876 +00000-0 +18063-4 0 9990 +2 99105 006.0257 225.1317 0011388 060.5983 299.5282 14.87992945325970 +1 99105U 12031A 18163.53783073 .00000880 00000-0 18589-4 0 9997 +2 99105 6.0255 213.4508 0011292 84.8589 275.2829 14.87993604326214 +1 99105U 12031A 18164.27510071 +.00000870 +00000-0 +17431-4 0 9999 +2 99105 006.0258 208.0960 0011303 095.9149 264.2268 14.87993811326324 +1 99105U 12031A 18165.48154044 +.00000864 +00000-0 +16824-4 0 9995 +2 99105 006.0257 199.3345 0011337 114.0757 246.0558 14.87994188326505 +1 99105U 12031A 18166.48690602 +.00000863 +00000-0 +16719-4 0 9993 +2 99105 006.0255 192.0387 0011426 129.3167 230.7975 14.87994467326651 +1 99105U 12031A 18167.69334140 .00000859 00000-0 16275-4 0 9991 +2 99105 6.0255 183.2760 0011548 147.2261 212.8585 14.87994827326835 +1 99105U 12031A 18168.49763167 +.00000849 +00000-0 +15253-4 0 9999 +2 99105 006.0255 177.4384 0011623 159.0311 201.0295 14.87994951326952 +1 99105U 12031A 18169.30192057 +.00000864 +00000-0 +16754-4 0 9998 +2 99105 006.0255 171.5970 0011733 170.7308 189.3038 14.87995269327076 +1 99105U 12031A 18170.44132903 +.00000883 +00000-0 +18782-4 0 9993 +2 99105 006.0255 163.3214 0011892 187.2125 172.7833 14.87995777327244 +1 99105U 12031A 18171.51371333 .00000881 00000-0 18565-4 0 9991 +2 99105 6.0256 155.5322 0012029 202.4752 157.4849 14.87996130327406 +1 99105U 12031A 18172.38502586 +.00000886 +00000-0 +19083-4 0 9994 +2 99105 006.0257 149.2038 0012133 214.7467 145.1869 14.87996470327533 +1 99105U 12031A 18174.46277340 +.00000904 +00000-0 +20926-4 0 9991 +2 99105 006.0260 134.1136 0012300 243.6365 116.2500 14.87997362327849 +1 99105U 12031A 18175.46813659 .00000939 00000-0 24711-4 0 9995 +2 99105 6.0262 126.8124 0012332 257.4973 102.3776 14.87998044327994 +1 99105U 12031A 18176.00433086 +.00000930 +00000-0 +23726-4 0 9994 +2 99105 006.0263 122.9181 0012336 264.7858 095.0863 14.87998218328070 +1 99105U 12031A 18177.34481815 +.00000985 +00000-0 +29519-4 0 9991 +2 99105 006.0264 113.1844 0012330 283.3137 076.5617 14.87999168328278 +1 99105U 12031A 18178.48423405 .00000965 00000-0 27461-4 0 9990 +2 99105 6.0266 104.9100 0012255 299.0588 60.8314 14.87999629328444 +1 99105U 12031A 18179.55662685 +.00000929 +00000-0 +23672-4 0 9994 +2 99105 006.0268 097.1232 0012155 313.9289 045.9837 14.87999971328600 +1 99105U 12031A 18179.55662685 .00000929 00000-0 23672-4 0 9994 +2 99105 6.0268 97.1232 0012155 313.9289 45.9837 14.87999971328600 +1 99105U 12031A 18181.43331634 +.00000900 +00000-0 +20596-4 0 9992 +2 99105 006.0270 083.4964 0011907 340.3785 019.5885 14.88000706328882 +1 99105U 12031A 18182.70678493 .00000889 00000-0 19460-4 0 9994 +2 99105 6.0273 74.2504 0011701 358.6933 1.3165 14.88001190329072 +1 99105U 12031A 18183.24298260 +.00000869 +00000-0 +17372-4 0 9992 +2 99105 006.0273 070.3579 0011606 006.5079 353.5200 14.88001208329156 +1 99105U 12031A 18184.38240284 +.00000851 +00000-0 +15428-4 0 9992 +2 99105 006.0271 062.0938 0011429 023.3571 336.7076 14.88001524329322 +1 99105U 12031A 18186.52718816 +.00000833 +00000-0 +13582-4 0 9990 +2 99105 006.0276 046.5219 0011253 055.1007 305.0178 14.88001887329647 +1 99105U 12031A 18187.46552949 +.00000859 +00000-0 +16354-4 0 9995 +2 99105 006.0275 039.7091 0011181 069.3545 290.7783 14.88002292329780 +1 99105U 12031A 18188.87303858 .00000863 00000-0 16707-4 0 9991 +2 99105 6.0273 29.4909 0011117 90.5322 269.6080 14.88002758329995 +1 99105U 12031A 18190.68268809 .00000873 00000-0 17840-4 0 9999 +2 99105 6.0274 16.3513 0011201 118.0742 242.0519 14.88003334330263 +1 99105U 12031A 18191.35292742 +.00000880 +00000-0 +18487-4 0 9993 +2 99105 006.0275 011.4867 0011242 128.9066 231.2066 14.88003618330369 +1 99105U 12031A 18193.43066803 +.00000841 +00000-0 +14333-4 0 9994 +2 99105 006.0279 356.4162 0011586 159.2176 200.8425 14.88003994330675 +1 99105U 12031A 18195.84352017 .00000866 00000-0 16997-4 0 9991 +2 99105 6.0276 338.8994 0011886 193.8355 166.1448 14.88004974331034 +1 99105U 12031A 18197.25101646 .00000877 00000-0 18075-4 0 9999 +2 99105 6.0277 328.6791 0012167 214.0689 145.8658 14.88005538331243 +1 99105U 12031A 18198.65851464 .00000880 00000-0 18399-4 0 9997 +2 99105 6.0274 318.4612 0012336 233.6305 126.2685 14.88006012331453 +1 99105U 12031A 18199.12768103 +.00000872 +00000-0 +17555-4 0 9991 +2 99105 006.0270 315.0555 0012351 239.5605 120.3303 14.88006193331523 +1 99105U 12031A 18200.53518149 +.00000892 +00000-0 +19624-4 0 9993 +2 99105 006.0267 304.8369 0012419 258.5860 101.2873 14.88006846331735 +1 99105U 12031A 18201.20542051 +.00000889 +00000-0 +19399-4 0 9992 +2 99105 006.0266 299.9694 0012444 268.0240 091.8463 14.88007085331831 +1 99105U 12031A 18202.41185231 .00000888 00000-0 19267-4 0 9992 +2 99105 6.0263 291.2088 0012415 284.3642 75.5109 14.88007531332018 +1 99105U 12031A 18203.48423766 .00000898 00000-0 20341-4 0 9994 +2 99105 6.0261 283.4223 0012347 299.1985 60.6909 14.88008055332170 +1 99105U 12031A 18205.49496351 +.00000899 +00000-0 +20436-4 0 9993 +2 99105 006.0259 268.8207 0012134 327.0343 032.9030 14.88008754332475 +1 99105U 12031A 18206.29925439 +.00000887 +00000-0 +19209-4 0 9999 +2 99105 006.0256 262.9788 0012000 338.3073 021.6547 14.88008957332594 +1 99105U 12031A 18208.37700882 +.00000849 +00000-0 +15162-4 0 9995 +2 99105 006.0255 247.9029 0011685 009.2401 350.7942 14.88009615332903 +1 99105U 12031A 18209.58344472 +.00000836 +00000-0 +13807-4 0 9998 +2 99105 006.0255 239.1415 0011560 026.8424 333.2302 14.88009929333081 +1 99105U 12031A 18211.46012002 .00000838 00000-0 14017-4 0 9995 +2 99105 6.0255 225.5133 0011379 54.7642 305.3551 14.88010471333368 +1 99105U 12031A 18212.59952704 .00000838 00000-0 14043-4 0 9991 +2 99105 6.0254 217.2369 0011311 71.8325 288.3035 14.88010817333530 +1 99105U 12031A 18213.87297992 +.00000833 +00000-0 +13591-4 0 9996 +2 99105 006.0252 207.9881 0011321 090.8302 269.3124 14.88011095333728 +1 99105U 12031A 18214.14107518 .00000843 00000-0 14640-4 0 9998 +2 99105 6.0252 206.0419 0011349 94.8374 265.3051 14.88011197333760 +1 99105U 12031A 18215.07940621 +.00000840 +00000-0 +14251-4 0 9996 +2 99105 006.0258 199.2271 0011378 109.4140 250.7219 14.88011484333900 +1 99105U 12031A 18216.35285335 .00000853 00000-0 15586-4 0 9998 +2 99105 6.0258 189.9786 0011464 128.2287 231.8875 14.88011960334092 +1 99105U 12031A 18218.29653146 .00000846 00000-0 14908-4 0 9991 +2 99105 6.0261 175.8619 0011689 156.9784 203.0869 14.88012467334383 +1 99105U 12031A 18219.43592690 .00000842 00000-0 14483-4 0 9998 +2 99105 6.0257 167.5836 0011786 173.1429 186.8861 14.88012730334554 +1 99105U 12031A 18220.37425248 +.00000847 +00000-0 +14902-4 0 9998 +2 99105 006.0260 160.7693 0011937 186.6269 173.3702 14.88013039334696 +1 99105U 12031A 18222.51899754 .00000831 00000-0 13216-4 0 9993 +2 99105 6.0262 145.1946 0012221 216.9359 142.9928 14.88013456335012 +1 99105U 12031A 18223.59137214 +.00000838 +00000-0 +13950-4 0 9990 +2 99105 006.0264 137.4140 0012330 231.9650 127.9365 14.88013813335178 +1 99105U 12031A 18224.12755971 +.00000840 +00000-0 +14136-4 0 9996 +2 99105 006.0266 133.5216 0012364 239.4290 120.4619 14.88013933335253 +1 99105U 12031A 18225.40100680 +.00000858 +00000-0 +16021-4 0 9993 +2 99105 006.0264 124.2823 0012389 256.8307 103.0439 14.88014460335445 +1 99105U 12031A 18226.33933642 +.00000868 +00000-0 +17093-4 0 9999 +2 99105 006.0266 117.4672 0012398 269.7577 090.1131 14.88014808335581 +1 99105U 12031A 18228.55111863 .00000892 00000-0 19692-4 0 9992 +2 99105 6.0270 101.4062 0012261 300.0439 59.8474 14.88015682335917 +1 99105U 12031A 18230.02564288 .00000900 00000-0 20551-4 0 9996 +2 99105 6.0271 90.6999 0012160 320.5476 39.3768 14.88016349336137 +1 99105U 12031A 18231.09802499 .00000913 00000-0 21949-4 0 9999 +2 99105 6.0273 82.9124 0012040 335.6914 24.2647 14.88016905336298 +1 99105U 12031A 18232.10338367 .00000923 00000-0 22985-4 0 9992 +2 99105 6.0273 75.6123 0011888 349.9740 10.0152 14.88017409336441 +1 99105U 12031A 18234.44921961 .00000911 00000-0 21806-4 0 9997 +2 99105 6.0275 58.5785 0011534 24.0820 335.9848 14.88018389336792 +1 99105U 12031A 18236.32588492 .00000896 00000-0 20223-4 0 9998 +2 99105 6.0277 44.9517 0011326 52.0081 308.1069 14.88019056337076 +1 99105U 12031A 18238.73873149 .00001015 00000-0 32829-4 0 9990 +2 99105 6.0276 27.4335 0011193 88.4037 271.7374 14.88020908337437 +1 99105U 12031A 18239.47598748 +.00001069 +00000-0 +38540-4 0 9996 +2 99105 006.0276 022.0799 0011216 099.6391 260.5006 14.88021718337545 +1 99105U 12031A 18240.28026575 .00001028 00000-0 34176-4 0 9997 +2 99105 6.0277 16.2402 0011250 111.7732 248.3594 14.88022046337669 +1 99105U 12031A 18241.41965794 .00001043 00000-0 35792-4 0 9990 +2 99105 6.0277 7.9671 0011301 129.0910 231.0224 14.88022843337834 +1 99105U 12031A 18243.43034857 .00000947 00000-0 25564-4 0 9999 +2 99105 6.0275 353.3829 0011572 158.8043 201.2566 14.88023530338138 +1 99105U 12031A 18244.90484987 +.00000914 +00000-0 +22032-4 0 9996 +2 99105 006.0274 342.6765 0011768 180.1899 179.8226 14.88024218338352 +1 99105U 12031A 18246.51339600 .00000923 00000-0 22940-4 0 9994 +2 99105 6.0273 330.9998 0012018 203.1495 156.8092 14.88025134338596 +1 99105U 12031A 18247.65278267 .00000928 00000-0 23465-4 0 9999 +2 99105 6.0272 322.7254 0012185 219.0726 140.8522 14.88025674338766 +1 99105U 12031A 18248.52408054 +.00000944 +00000-0 +25186-4 0 9995 +2 99105 006.0266 316.4048 0012225 231.6082 128.2949 14.88026087338896 +1 99105U 12031A 18250.40072152 .00000932 00000-0 23854-4 0 9998 +2 99105 6.0265 302.7773 0012306 257.3611 102.5141 14.88026886339179 +1 99105U 12031A 18251.07095137 +.00000923 +00000-0 +22979-4 0 9997 +2 99105 006.0264 297.9095 0012323 266.6489 093.2230 14.88027133339277 +1 99105U 12031A 18252.47843615 .00000914 00000-0 21946-4 0 9997 +2 99105 6.0261 287.6889 0012325 285.9235 73.9536 14.88027692339483 +1 99105U 12031A 18253.48378393 .00000909 00000-0 21438-4 0 9995 +2 99105 6.0261 280.3883 0012267 299.8809 60.0101 14.88028099339637 +1 99105U 12031A 18255.56150483 .00000973 00000-0 28277-4 0 9993 +2 99105 6.0258 265.2988 0012062 328.8284 31.1129 14.88029537339945 +1 99105U 12031A 18256.49983114 +.00000973 +00000-0 +28342-4 0 9994 +2 99105 006.0257 258.4841 0011938 342.0812 017.8896 14.88030042340084 +1 99105U 12031A 18257.57220381 .00000998 00000-0 30970-4 0 9991 +2 99105 6.0256 250.6954 0011760 357.3615 2.6451 14.88030816340241 +1 99105U 12031A 18258.64457606 .00001011 00000-0 32369-4 0 9993 +2 99105 6.0255 242.9066 0011627 12.9360 347.1067 14.88031552340407 +1 99105U 12031A 18259.64992432 .00001000 00000-0 31247-4 0 9997 +2 99105 6.0256 235.6055 0011486 27.6899 332.3841 14.88032126340558 +1 99105U 12031A 18260.65527132 .00000985 00000-0 29636-4 0 9991 +2 99105 6.0256 228.3041 0011381 42.6232 317.4779 14.88032686340706 +1 99105U 12031A 18261.45954889 +.00000956 +00000-0 +26526-4 0 9996 +2 99105 006.0259 222.4651 0011444 054.9417 305.1785 14.88032977340829 +1 99105U 12031A 18265.48091541 +.00000961 +00000-0 +27028-4 0 9994 +2 99105 006.0262 193.2700 0011261 115.6918 244.4375 14.88035202341422 diff --git a/nustar_pysolar/planning.py b/nustar_pysolar/planning.py index ebc3d9b..982901e 100644 --- a/nustar_pysolar/planning.py +++ b/nustar_pysolar/planning.py @@ -354,3 +354,135 @@ def sunlight_periods(infile, tstart, tend): return -1 else: return in_range + + +def make_mosaic(orbit, outfile='mosaic.txt', write_output=False, make_regions=False, + reg_pref='testbox'): + ''' + Code to make a mosaic for a 5x5 tiled array on the Sun. + + Input: + + tstart = '2018-05-28T15:37:00' + tend = '2018-05-28T23:10:00' + + positions = make_mosaic(tstart, tend, write_output=True) + + Optional flags: + + write_output = [False] / True + Write the output pointing positions in NuSTAR SOC readable formats in 'outfile' for all of the pointings. + + outfile = ['mosaic.txt'] + Output file if write_output is used. + + make_regions: [False] / True + Make ds9 region files for each tile so that you can see how the FoV moves with each mosaic location. + + reg_pref: 'testbox' + The prefix for the region files. Useful if you want to make this meaningful. + + Output mosaic file has columns of: + "Arrive By Time" RA DEC RA_SUN DEC_SUN + + ''' + import numpy as np + box_pa = get_nustar_roll(orbit[0], 0) + pa = box_pa + 90*u.deg + + print('Step PA', pa) + + base = np.array([-1.45, -0.725, 0, 0.725, 1.45]) + xsteps = np.append(base, np.flip(base, 0)) + xsteps = np.append(xsteps, base) + xsteps = np.append(xsteps, np.flip(base, 0)) + xsteps = np.append(xsteps, base) + + ysteps = np.array(np.zeros(5) + 1.45) + ysteps = np.append(ysteps, np.zeros(5) + 0.725) + ysteps = np.append(ysteps, np.zeros(5)) + ysteps = np.append(ysteps, np.zeros(5)-0.725) + ysteps = np.append(ysteps, np.zeros(5)-1.45) + + + dt = (orbit[1] - orbit[0]) / 25. + + print("Orbit start: {} Orbit end: {}".format(orbit[0].isoformat(), orbit[1].isoformat())) + print("Dwell per position:", dt) + print("") + print("NuSTAR Roll Angle to get DET0 in top right {:.02f} deg".format(box_pa.value)) + print("") + + + if write_output is True: + f = open(outfile, 'w') + + aim_time = orbit[0] + for ind, pair in enumerate(zip(xsteps, ysteps)): + arrive_time = aim_time + aim_time = aim_time + dt + # Rotate to the correct PA angle + + delx = -pair[0]*10/60. + dely = pair[1]*10/60 + + offset = [0., 0.]*u.deg + sun_pos = get_skyfield_position(aim_time, offset, load_path='../data', parallax_correction=True) +# print('Sun time: {} RA (deg): {} Dec (deg): {}'.format(aim_time.isoformat(), sun_pos[0], sun_pos[1])) + + offset = [-delx,dely]*u.deg + sky_pos = get_skyfield_position(aim_time, offset, load_path='../data', parallax_correction=True) +# print('Arrive time: {} RA (deg): {} Dec (deg): {}'.format(arrive_time.isoformat(), sky_pos[0], sky_pos[1])) + + if make_regions: + make_test_region(sky_pos[0], sky_pos[1], box_pa, sun_pos[0], sun_pos[1],box_pa, outname=reg_pref+'{}.reg'.format(ind)) + if write_output: + f.write('{0} {1:.4f} {2:.4f} {3:.4f} {4:.4f}\n'.format(arrive_time.strftime('%Y:%j:%H:%M:%S'), + sky_pos[0].value, sky_pos[1].value, + sun_pos[0].value, sun_pos[1].value)) + + if write_output: + f.close() + + +def make_test_region(boxra, boxdec, boxpa, sunra, sundec, sunpa, outname='testbox.reg'): + ''' + Code to produce a ds9 region file showing the Sun and the FoV on the sky. + + Inputs: + + boxra, boxdec: These are the center of the FoV. Must have astropy units. + boxpa: This is the SKY PA angle. Must have astropy units. + + sunra, sundec: The center of the Sun. + sunpa: The solar north pole angle + + outname: the name of the region file that you want to produce. + Defaults to 'tetsbox.reg' + + The output is a ds9 region file. Go into ds9, enter the Sun ra/dec coordinates and + load an optical stellar background field (under the 'Analysis' tab), then load + the region file to see what's going on. + + ''' + f = open(outname, 'w') + f.write("# Region file format: DS9 version 4.1\n") + f.write('global color=green dashlist=8 3 width=1 font="helvetica 10 normal roman" select=1 highlite=1 dash=0 fixed=0 edit=1 move=1 delete=1 include=1 \ +source=1 \n') + f.write('fk5 \n') + + outstring = 'box({0}, {1}, 720", 720", {2})'.format(boxra.value, boxdec.value, boxpa.value+90 % 360) + f.write(outstring+'\n') + + outstring = 'circle({0}, {1}, 960.5")'.format(sunra.value, sundec.value) + + + f.write(outstring+'\n') + + outstring = 'vector({0}, {1}, 960.5", {2})'.format(sunra.value, sundec.value, sunpa.value+90 % 360 ) + f.write(outstring+'\n') + + f.close + return + +