Commits

monoid committed 1bac55f

Display C/2009 P1 (Garradd) instead of C/2009 R1 (McNaught).

Comments (0)

Files changed (1)

     (new StarMap.Telrad(0.901, 0.451)).draw(ctx, this['proj']);
     (new StarMap.Telrad(0.301, 2.151)).draw(ctx, this['proj']);
 
+    /*
+Ephemerides from http://www.minorplanetcenter.net/iau/MPEph/MPEph.html
+
+C/2009 P1 (Garradd)
+
+Perturbed ephemeris below is based on elements from MPEC 2012-A52.
+
+   CK09P010
+Date       UT      R.A. (J2000) Decl.    Delta     r     El.    Ph.   m1     Sky Motion
+            h m s                                                            "/min    P.A.
+2012 01 17 000000   17.4533  +32.839     1.749   1.586   64.0  33.9   7.2    1.17    350.9
+2012 01 22 000000   17.4175  +35.301     1.683   1.602   68.1  34.8   7.2    1.34    348.6
+2012 01 27 000000   17.3649  +38.107     1.615   1.621   72.6  35.4   7.1    1.54    346.0
+2012 02 01 000000   17.2890  +41.290     1.548   1.642   77.3  35.8   7.1    1.77    343.2
+2012 02 06 000000   17.1802  +44.877     1.482   1.666   82.2  35.9   7.1    2.02    340.1
+2012 02 11 000000   17.0241  +48.874     1.421   1.692   87.3  35.6   7.0    2.29    336.3
+2012 02 16 000000   16.7981  +53.248     1.367   1.720   92.4  35.0   7.0    2.57    331.7
+2012 02 21 000000   16.4653  +57.889     1.321   1.750   97.4  34.1   7.0    2.85    325.6
+2012 02 26 000000   15.9670  +62.547     1.288   1.782  102.1  32.9   7.1    3.09    317.2
+2012 03 02 000000   15.2189  +66.752     1.269   1.815  106.2  31.6   7.1    3.26    305.0
+2012 03 07 000000   14.1497  +69.747     1.267   1.851  109.4  30.4   7.2    3.35    288.1
+2012 03 12 000000   12.8391  +70.676     1.283   1.888  111.4  29.4   7.3    3.33    267.6
+2012 03 17 000000   11.5897  +69.266     1.317   1.926  112.0  28.6   7.4    3.20    247.8
+2012 03 22 000000   10.6405  +66.145     1.368   1.965  111.4  28.2   7.6    3.00    232.4
+2012 03 27 000000   09.9973  +62.191     1.436   2.005  109.7  27.9   7.8    2.75    221.2
+2012 04 01 000000   09.5739  +58.016     1.517   2.047  107.0  27.8   8.0    2.48    213.1
+2012 04 06 000000   09.2944  +53.951     1.611   2.089  103.8  27.7   8.2    2.21    206.8
+2012 04 11 000000   09.1087  +50.147     1.714   2.132  100.1  27.6   8.5    1.96    201.6
+2012 04 16 000000   08.9858  +46.657     1.825   2.176   96.2  27.3   8.7    1.74    196.9
+2012 04 21 000000   08.9067  +43.487     1.943   2.221   92.2  26.9   8.9    1.54    192.6
+2012 04 26 000000   08.8588  +40.618     2.065   2.266   88.0  26.4   9.1    1.38    188.4
+2012 05 01 000000   08.8340  +38.022     2.190   2.311   83.9  25.7   9.3    1.24    184.2
+2012 05 06 000000   08.8263  +35.669     2.317   2.358   79.8  24.9   9.5    1.12    180.2
+2012 05 11 000000   08.8317  +33.526     2.446   2.404   75.7  24.0   9.8    1.03    176.2
+2012 05 16 000000   08.8474  +31.567     2.574   2.451   71.6  23.0   9.9    0.95    172.2
+2012 05 21 000000   08.8713  +29.766     2.703   2.498   67.6  22.0  10.1    0.88    168
+    */
+
     // Draw path of C/2009 R1 (McNaught)
-    var C2009R1_DA = ["2010 05 04", "2010 05 09", "2010 05 14",
-                      "2010 05 19", "2010 05 24", "2010 05 29",
-                      "2010 06 03", "2010 06 08", "2010 06 13",
-                      "2010 06 18", "2010 06 23", "2010 06 28",
-                      "2010 07 03", "2010 07 08", "2010 07 13",
-                      "2010 07 18", "2010 07 23", "2010 07 28"];
-    var C2009R1_RA = [23+35.99/60.0, 23+48.09/60.0,  0+ 2.04/60.0,
-                       0+18.56/60.0,  0+38.69/60.0,  1+ 4.06/60.0,
-                       1+37.04/60.0,  2+20.73/60.0,  3+17.79/60.0,
-                       4+26.49/60.0,  5+37.00/60.0,  6+36.84/60.0,
-                       7+20.69/60.0,  7+50.73/60.0,  8+11.73/60.0,
-                       8+27.57/60.0,  8+40.53/60.0,  8+51.84/60.0];
-    var C2009R1_DE = [ 9+56.1/60.0, 13+11.8/60.0, 16+53.6/60.0,
-                      21+ 5.1/60.0, 25+48.9/60.0, 31+ 3.2/60.0,
-                      36+37.2/60.0, 42+ 1.6/60.0, 46+19.9/60.0,
-                      48+14.0/60.0, 46+44.8/60.0, 42+ 0.7/60.0,
-                      35+ 8.7/60.0, 27+30.3/60.0, 20+ 4.7/60.0,
-                      13+16.4/60.0, 07+ 8.3/60.0,  1+35.8/60.0];
-    var C2009R1 = new StarMap.Path(C2009R1_DA, C2009R1_RA, C2009R1_DE, {
+    var C2009P1_DA = ["2012 01 17", "2012 01 22", "2012 01 27",
+                      "2012 02 01", "2012 02 06", "2012 02 11",
+                      "2012 02 16", "2012 02 21", "2012 02 26",
+                      "2012 03 02", "2012 03 07", "2012 03 12",
+                      "2012 03 17", "2012 03 22", "2012 03 27",
+                      "2012 04 01", "2012 04 06", "2012 04 11"];
+    var C2009P1_RA = [
+        17.4533, 17.4175, 17.3649,
+        17.2890, 17.1802, 17.0241,
+        16.7981, 16.4653, 15.9670,
+        15.2189, 14.1497, 12.8391,
+        11.5897, 10.6405,  9.9973,
+         9.5739,  9.2944,  9.1087
+    ];
+    var C2009P1_DE = [
+            +32.839, +35.301, +38.107,
+            +41.290, +44.877, +48.874,
+            +53.248, +57.889, +62.547,
+            +66.752, +69.747, +70.676,
+            +69.266, +66.145, +62.191,
+            +58.016, +53.951, +50.147
+    ];
+    var C2009P1 = new StarMap.Path(C2009P1_DA, C2009P1_RA, C2009P1_DE, {
         'fillStyle': 'green',
         'strokeStyle': 'green',
         'lineWidth': 1
         'fillStyle': '#8F8',
         'strokeStyle': '#8F8'
     });
-    C2009R1.draw(ctx, this['proj']);
+    C2009P1.draw(ctx, this['proj']);
     
+    /*
+    http://www.minorplanetcenter.net/mpec/K12/K12A52.html
+
+    C/2009 P1 (Garradd)
+Epoch 2011 Dec. 25.0 TT = JDT 2455920.5
+T 2011 Dec. 23.67758 TT                                 MPC
+q   1.5505367            (2000.0)            P               Q
+z  -0.0006782      Peri.   90.74770     -0.16661319     -0.82691127
+ +/-0.0000002      Node   325.99770     -0.58719595     +0.52078702
+e   1.0010516      Incl.  106.17749     +0.79211171     +0.21212879
+From 4943 observations 2009 Aug. 13-2011 Dec. 4, mean residual 0".4.
+    */
     // Draw current position of the comet
-    var C2009R1_param = {
-        title: 'C/2009 R1 (McNaught)',
-        t0: 2455379.6792-2400000.5,
-        q: 0.405011,
-        z: -0.000808,
-        e: 1.000327,
-        peri: 130.7013*DEG2RAD,
-        node: 322.6220*DEG2RAD,
-        incl: 77.0319*DEG2RAD
+    var C2009P1_param = {
+        title: 'C/2009 P1 (Garradd)',
+        t0: 2455920.5-2400000.5,
+        q: 1.5505367,
+        z: -0.0006782,
+        e: 1.0010516,
+        peri: 90.74770*DEG2RAD,
+        node: 325.99770*DEG2RAD,
+        incl: 106.17749*DEG2RAD
     };
 
-    (new StarMap.Object(C2009R1_param, 'rgba(220, 255, 220, 0.7)', mjd, jct)).draw(ctx, this['proj']);
+    (new StarMap.Object(C2009P1_param, 'rgba(220, 255, 220, 0.7)', mjd, jct)).draw(ctx, this['proj']);
 
 
     var lutetia = {