batgrbproduct analysis for trigger=687522

batgrbproduct Analysis for Trigger=687522

Analysis Summary

Notes:

1) No attempt has been made to eliminate the non-significant decimal places.
2) This batgrbproduct analysis is done using the HEASARC data.
3) This batgrbproduct analysis is using the flight position for the mask-weighting calculation.

======================================================================
                  BAT GRB EVENT DATA PROCESSING REPORT
    Process Script Ver: batgrbproduct v2.46
    Process Start Time (local time): Tue Jun 21 08:53:24 EDT 2016
             Host Name: bat1
   Current Working Dir: /local/data/bat1/alien/Swift_3rdBATcatalog/event/batevent_reproc/trigger687522
======================================================================

  Trigger: 00687522      Segment: 000

  BAT Trigger Time: 485861123.968   [s; MET]
      Trigger Stop: 485861124.48   [s; MET]
               UTC: 2016-05-25T09:25:07.649260    [includes UTCF correction] 
       Where From?: TDRSS position message

  BAT 
    RA: 149.442925515661   Dec: 51.1802165201058  [deg; J2000]

     Catalogged Source?: NO
     Point Source?:      YES
     GRB Indicated?:     YES    [ by BAT flight software ]
     Image S/N Ratio:    6.55
     Image Trigger?:     NO
     Rate S/N Ratio:     11.3578166916005  [ if not an image trigger ]
     Image S/N Ratio:    6.55

  Flight Position:   [ source = BAT ]
    RA: 149.442925515661   Dec: 51.1802165201058  [deg; J2000]

  Ground Position:   [ source = BAT pre-slew burst ]
    RA: 149.438236452283   Dec: 51.1967622366146  [deg; J2000]  
        { 09h 57m 45.2s , +51d 11' 48.3" }
         +/- 2.18247287001802 [arcmin] (estimated 90% radius based on SNR)
         +/- 1.15642576875452 [arcmin] (formal 1-sigma fit error)
    SNR: 9.97483525889028
 
  Angular difference between Flight and Ground position is 1.008282 arcmin
 
  Partial Coding Fraction: 0.8125    [ including projection effects ]

  Duration
     T90: 0.291999936103821 +/- 0.0456069787374054
        Measured from: 0.0279999971389771
                   to: 0.319999933242798     [s; relative to TRIGTIME]
     T50: 0.111999988555908 +/- 0.0636553038431823
        Measured from: 0.139999985694885
                   to: 0.251999974250793     [s; relative to TRIGTIME]



  Fluence
     Peak Flux  (peak 1 second)
        Measured from: -0.324000000953674     
                   to: 0.675999999046326     [s; relative to TRIGTIME]
     Total Fluence        
        Measured from: 0.0160000324249268     
                   to: 0.332000017166138     [s; relative to TRIGTIME]
     
             Band 1     Band 2     Band 3     Band 4
              15-25      25-50     50-100    100-350     keV
 Total     0.023671   0.022296   0.015857   0.011736
           0.004328   0.004317   0.003729   0.002999 [error]
 Peak      0.033492   0.037291   0.032721   0.014650
           0.006812   0.007421   0.006300   0.005262 [error]

     [ fluence units of on-axis counts / fully illuminated detector ] 

======================================================================
 
=== Fluence/Peak Flux Summary ===
 
              Model    Band1    Band2    Band3    Band4    Band5    Total      
                       15-25    25-50   50-100  100-350   15-150   15-350   keV
Fluence   Power-law 6.35e-09 9.16e-09 9.83e-09 1.96e-08 3.14e-08 4.49e-08   erg/cm2
              Error 1.79e-09 1.90e-09 3.18e-09 1.24e-08 7.63e-09 1.61e-08   erg/cm2
Peak flux Power-law 2.87e-01 2.95e-01 2.14e-01 2.49e-01 8.93e-01 1.05e+00   ph/cm2/s
              Error 9.26e-02 5.93e-02 5.08e-02 1.22e-01 1.70e-01 1.99e-01   ph/cm2/s
 
Fluence   Cutoff-PL 5.75e-09 1.08e-08 9.42e-09 3.25e-09 2.84e-08 2.92e-08   erg/cm2
              Error 2.02e-09 2.77e-09 4.02e-09 1.17e-08 7.82e-09 1.54e-08   erg/cm2
Peak flux Cutoff-PL 2.88e-01 2.94e-01 2.13e-01 2.55e-01 8.95e-01 1.05e+00   ph/cm2/s
              Error -------- 6.52e-02 5.37e-02 -------- 1.74e-01 2.30e-01   ph/cm2/s
 
Fluence        Band -------- -------- -------- -------- -------- --------   erg/cm2
              Error -------- -------- -------- -------- -------- --------   erg/cm2
Peak flux      Band -------- -------- -------- -------- -------- --------   ph/cm2/s
              Error -------- -------- -------- -------- -------- --------   ph/cm2/s

 
=== Time-averaged spectral analysis ===
 
Model Ph_index Ph_90err Epeak  Epeak_90err Norm     Nomr_90err chi2   dof
PL    1.899    0.350    -      -           1.08e-02 2.66e-03   47.540 57
CPL   0.602    1.839      45.0    0.0      6.12e-02 6.05e-01   45.410 56
 
 
=== 1-s peak spectral analysis ===
 
Model Ph_index Ph_90err Epeak  Epeak_90err Norm     Nomr_90err chi2   dof
PL    1.463    0.327    -      -           7.22e-03 1.38e-03   50.580 57
CPL   1.461    1.519    10000.0    0.0      7.24e-03 -2.93e-03   50.580 56
 
Model      : spectral model: simple power-law (PL) and cutoff power-law (CPL)
Ph_index   : photon index: E^-{Ph_index} (minus sign is not included in the index)
Ph_90err   : 90% error of the photon index
Epeak      : Epeak in keV
Epeak_90err: 90% error of Epeak
Norm       : normalization at 50 keV in PL model and at 1 keV in CPL model
Nomr_90err : 90% error of the normalization
chi2       : Chi-Squared
dof        : degree of freedom

 
=== T90/T50 calculation in the 50-300 keV band (BATSE band) using 64 ms lc ===
 
T90 in the 50-300 keV band: 0.000000 sec.
T50 in the 50-300 keV band: 0.000000 sec.

 
=== Table of the duration information ==
 
        Value   tstart    tstop
T100    0.316    0.016    0.332
T90     0.292    0.028    0.320
T50     0.112    0.140    0.252
Peak    1.000   -0.324    0.676

Click here for the BAT GCN circular template

Estimated redshift from machine learning (Ukwatta et al. 2016) - may take a few minute to calcaulte

Image

Pre-slew 15.0-350.0 keV image (Event data)

#   RAcent    DECcent   POSerr    Theta     Phi              Peak Cts       SNR   Name
  149.4382    51.1968   0.0192    27.8376   172.6543         0.0780268      9.975 TRIG_00687522                 
 
Foreground time interval of the image:
  0.016   0.332 (delta_t =   0.316 [sec])
Background time interval of the image:
-495.767 -480.108 (delta_t =  15.659 [sec])
-239.968   0.016 (delta_t = 239.984 [sec])

Lightcurves

Notes:

1) The mask-weighted light curves are using the flight position.
2) Multiple plots of different time binning/intervals are shown to cover all scenarios of short/long GRBs, rate/image triggers, and real/false positives.
3) For all the mask-weighted lightcurves, the y-axis units are counts/sec/det where a det is 0.4 x 0.4 = 0.16 sq cm.
4) The verticle lines correspond to following: green dotted lines are T50 interval, black dotted lines are T90 interval, blue solid line(s) is a spacecraft slew start time, and orange solid line(s) is a spacecraft slew end time.
5) Time of each bin is in the middle of the bin.

1 s binning

From T0-310 s to T0+310 s

From T100_start - 20 s to T100_end + 30 s

Full time interval

64 ms binning

Full time interval

From T100_start to T100_end

From T0-5 sec to T0+5 sec

16 ms binning

From T100_start-5 s to T100_start+5 s

Blue+Cyan dotted lines: Time interval of the scaled map
From 0.000000 to 0.512000

8 ms binning

From T100_start-3 s to T100_start+3 s

Blue+Cyan dotted lines: Time interval of the scaled map
From 0.000000 to 0.512000

2 ms binning

From T100_start-1 s to T100_start+1 s

Blue+Cyan dotted lines: Time interval of the scaled map
From 0.000000 to 0.512000

SN=5 or 10 sec. binning (whichever binning is satisfied first)

SN=5 or 10 sec. binning (T < 200 sec)

Quad-rate summed light curves (from T0-300s to T0+1000s)

Spectral Evolution

Spectra

Notes:

1) The fitting includes the systematic errors.
2) When the burst includes telescope slew time periods, the fitting uses an average response file made from multiple 5-s response files through out the slew time plus single time preiod for the non-slew times, and weighted by the total counts in the corresponding time period using addrmf. An average response file is needed becuase a spectral fit using the pre-slew DRM will introduce some errors in both a spectral shape and a normalization if the PHA file contains a lot of the slew/post-slew time interval.
3) For fits to more complicated models (e.g. a power-law over a cutoff power-law), the BAT team has decided to require a chi-square improvement of more than 6 for each extra dof.

Time averaged spectrum fit using the average DRM

Power-law model

Time interval is from    0.016 sec. to    0.332 sec.
 
Spectral model in power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 1.89868     ( -0.335711   0.364993 )
Norm@50keV  : 1.08048E-02    ( -0.00269763   0.00261349 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          47.54 using 59 PHA bins.
# Reduced chi-squared =         0.8341 for     57 degrees of freedom 
# Null hypothesis probability =   8.097246e-01
 
Photon flux (15-150 keV)    in 0.316 sec: 1.54982 ( -0.32097  0.32105 ) ph/cm2/s
Energy fluence (15-150 keV)             : 3.13962e-08 ( -7.52484e-09  7.73625e-09 ) ergs/cm2

Cutoff power-law model

Time interval is from    0.016 sec. to    0.332 sec.
Spectral model in the cutoff power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 0.601590     ( -2.19439   1.48417 )
Epeak [keV] : 45.0010      ( -44.999 -44.999 )
Norm@50keV  : 6.12397E-02    ( -0.0612829   1.14849 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          45.41 using 59 PHA bins.
# Reduced chi-squared =         0.8109 for     56 degrees of freedom 
# Null hypothesis probability =   8.432209e-01
 
Photon flux (15-150 keV)    in 0.316 sec: 1.51022 ( -0.32336  0.32409 ) ph/cm2/s
Energy fluence (15-150 keV)             : 2.84457e-08 ( -7.3416e-09  8.30511e-09 ) ergs/cm2

Band function

Time interval is from    0.016 sec. to    0.332 sec.
 
Spectral model in the Band function:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
alpha       : -0.585791   ( -0.720236   -0.720236 )
beta        : -3.91332    ( 3.94099   13.4556 )
Epeak [keV] : 44.8903      ( -44.8907 409.13 )
Norm@50keV   : 6.26023E-02    ( -0.0625976   991.385 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          45.39 using 59 PHA bins.
# Reduced chi-squared =         0.8253 for     55 degrees of freedom 
# Null hypothesis probability =   8.188457e-01
 
Photon flux (15-150 keV)    in 0.316 sec:  (    ) ph/cm2/s
Energy fluence (15-150 keV)             : 0 ( 0  0 ) ergs/cm2

Single BB

 
Spectral model blackbody:
------------------------------------------------------------
Parameters  : value    Lower 90%  Upper 90%
kT [keV]    : 10.6785  (-2.19855  2.75507)
R^2/D10^2   : 0.601472  (-0.352364   0.831772)
(R is the radius in km and D10 is the distance to the source in units of 10 kpc)
------------------------------------------------------------
#Fit statistic : Chi-Squared =          48.11 using 59 PHA bins.
# Reduced chi-squared =         0.8440 for     57 degrees of freedom 
# Null hypothesis probability =   7.930744e-01
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  4.499e-09   1.611e-09
25- 50  1.224e-08   4.040e-09
50-150  7.701e-09   3.860e-09
15-150  2.444e-08   7.916e-09

Thermal bremsstrahlung (OTTB)

 
 
Spectral model: thermal bremsstrahlung
------------------------------------------------------------
Parameters  : value         Lower 90%    Upper 90%
kT [keV]    : 74.6709       (-32.7104      92.6849)
Norm        : 7.60641       (-1.98202     2.76623)
------------------------------------------------------------
#Fit statistic : Chi-Squared =          45.92 using 59 PHA bins.
# Reduced chi-squared =         0.8056 for     57 degrees of freedom 
# Null hypothesis probability =   8.531877e-01
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  6.205e-09   1.856e-09
25- 50  1.004e-08   3.541e-09
50-150  1.352e-08   9.316e-09
15-150  2.976e-08   1.317e-08

Double BB

 
 
Spectral model: bbodyrad<1> + bbodyrad<2>
------------------------------------------------------------
Parameters  : value         Lower 90%    Upper 90%
kT1 [keV]   : 5.65593        (-4.09159      -5.64171)
R1^2/D10^2  : 3.44299        (-3.33631      1648.73)
kT2 [keV]   : 15.6092        (-15.6263      -15.6272)
R2^2/D10^2  : 0.102169    (-0.102955   0.606767)
------------------------------------------------------------
#Fit statistic : Chi-Squared =          45.19 using 59 PHA bins.
# Reduced chi-squared =         0.8216 for     55 degrees of freedom 
# Null hypothesis probability =   8.246061e-01
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  5.862e-09   3.222e-09
25- 50  1.043e-08   5.356e-09
50-150  1.162e-08   7.131e-09
15-150  2.791e-08   1.534e-08

Peak spectrum fit

Power-law model

Time interval is from   -0.324 sec. to    0.676 sec.
 
Spectral model in power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 1.46340     ( -0.32656   0.328343 )
Norm@50keV  : 7.22167E-03    ( -0.00138615   0.0013718 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          50.58 using 59 PHA bins.
# Reduced chi-squared =         0.8874 for     57 degrees of freedom 
# Null hypothesis probability =   7.130228e-01
 
Photon flux (15-150 keV)    in 1 sec: 0.892987 ( -0.169495  0.169813 ) ph/cm2/s
Energy fluence (15-150 keV)             : 6.89477e-08 ( -1.44284e-08  1.46145e-08 ) ergs/cm2

Cutoff power-law model

Time interval is from   -0.324 sec. to    0.676 sec.
Spectral model in the cutoff power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 1.46083     ( -1.28125   1.75589 )
Epeak [keV] : 1.00000E+04      ( -10000 -10000 )
Norm@50keV  : 7.24496E-03    ( -0.00139504   -0.00724506 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          50.58 using 59 PHA bins.
# Reduced chi-squared =         0.9033 for     56 degrees of freedom 
# Null hypothesis probability =   6.793377e-01
 
Photon flux (15-150 keV)    in 1 sec: 0.894736 ( -0.176433  0.171934 ) ph/cm2/s
Energy fluence (15-150 keV)             : 6.92261e-08 ( -1.60128e-08  1.47973e-08 ) ergs/cm2

Band function

Time interval is from   -0.324 sec. to    0.676 sec.
 
Spectral model in the Band function:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
alpha       : 4.70776   ( -4.70776   -4.70776 )
beta        : -1.50358    ( -0.371102   0.34545 )
Epeak [keV] : 21.7613      ( -21.7613 25.3316 )
Norm@50keV   : 1026.29    ( -190.519   4.5553e+06 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          50.38 using 59 PHA bins.
# Reduced chi-squared =         0.9160 for     55 degrees of freedom 
# Null hypothesis probability =   6.514728e-01
 
Photon flux (15-150 keV)    in 1 sec:  (    ) ph/cm2/s
Energy fluence (15-150 keV)             : 0 ( 0  0 ) ergs/cm2

Single BB

 
Spectral model blackbody:
------------------------------------------------------------
Parameters  : value    Lower 90%  Upper 90%
kT [keV]    : 15.1666  (-3.69517  5.11703)
R^2/D10^2   : 9.70691E-02  (-0.0614547   0.161601)
(R is the radius in km and D10 is the distance to the source in units of 10 kpc)
------------------------------------------------------------
#Fit statistic : Chi-Squared =          59.97 using 59 PHA bins.
# Reduced chi-squared =          1.052 for     57 degrees of freedom 
# Null hypothesis probability =   3.684179e-01
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  4.666e-09   1.880e-09
25- 50  1.873e-08   6.296e-09
50-150  2.946e-08   1.215e-08
15-150  5.286e-08   1.672e-08

Thermal bremsstrahlung (OTTB)

 
 
Spectral model: thermal bremsstrahlung
------------------------------------------------------------
Parameters  : value         Lower 90%    Upper 90%
kT [keV]    : 199.363       (-92.3823      -199.363)
Norm        : 3.67091       (-0.684491     0.684478)
------------------------------------------------------------
#Fit statistic : Chi-Squared =          51.05 using 59 PHA bins.
# Reduced chi-squared =         0.8956 for     57 degrees of freedom 
# Null hypothesis probability =   6.967888e-01
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  9.141e-09   5.320e-09
25- 50  1.758e-08   1.013e-08
50-150  3.856e-08   2.178e-08
15-150  6.529e-08   3.676e-08

Double BB

 
 
Spectral model: bbodyrad<1> + bbodyrad<2>
------------------------------------------------------------
Parameters  : value         Lower 90%    Upper 90%
kT1 [keV]   : 7.07591        (-3.47572      5.32026)
R1^2/D10^2  : 0.955902        (-0.763996      1.66683)
kT2 [keV]   : 30.5475        (-13.1838      -30.9724)
R2^2/D10^2  : 7.03811E-03    (-0.00665414   0.0319733)
------------------------------------------------------------
#Fit statistic : Chi-Squared =          50.77 using 59 PHA bins.
# Reduced chi-squared =         0.9231 for     55 degrees of freedom 
# Null hypothesis probability =   6.368313e-01
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  8.790e-09   5.145e-09
25- 50  1.718e-08   8.955e-09
50-150  4.319e-08   2.463e-08
15-150  6.916e-08   3.624e-08

Pre-slew spectrum (Pre-slew PHA with pre-slew DRM)

Power-law model

 
Spectral model in power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 1.89868     ( -0.335711   0.364993 )
Norm@50keV  : 1.08048E-02    ( -0.00269763   0.00261349 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          47.54 using 59 PHA bins.
# Reduced chi-squared =         0.8341 for     57 degrees of freedom 
# Null hypothesis probability =   8.097246e-01
 
Photon flux (15-150 keV)    in 0.316 sec: 1.54972 ( -0.320896  0.321024 ) ph/cm2/s
Energy fluence (15-150 keV)             : 3.13914e-08 ( -7.51767e-09  7.74845e-09 ) ergs/cm2

Cutoff power-law model

Spectral model in the cutoff power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 0.601779     ( -2.19441   1.48414 )
Epeak [keV] : 45.0017      ( -44.9991 -44.9991 )
Norm@50keV  : 6.12238E-02    ( -0.0612828   1.14846 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          45.41 using 59 PHA bins.
# Reduced chi-squared =         0.8109 for     56 degrees of freedom 
# Null hypothesis probability =   8.432209e-01
 
Photon flux (15-150 keV)    in 0.316 sec: 1.50998 ( -0.32325  0.324 ) ph/cm2/s
Energy fluence (15-150 keV)             : 2.84527e-08 ( -7.34239e-09  8.31734e-09 ) ergs/cm2

Band function

 
Spectral model in the Band function:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
alpha       : -0.585791   ( -0.720236   -0.720236 )
beta        : -3.91332    ( 3.94099   13.4556 )
Epeak [keV] : 44.8903      ( -44.8907 409.13 )
Norm@50keV   : 6.26023E-02    ( -0.0625976   991.385 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          45.39 using 59 PHA bins.
# Reduced chi-squared =         0.8253 for     55 degrees of freedom 
# Null hypothesis probability =   8.188457e-01
 
Photon flux (15-150 keV)    in 0.316 sec:  (    ) ph/cm2/s
Energy fluence (15-150 keV)             : 0 ( 0  0 ) ergs/cm2

S(25-50 keV) vs. S(50-100 keV) plot

S(25-50 keV)  = 9.16e-09
S(50-100 keV) = 9.83e-09

T90 vs. Hardness ratio plot

T90                                   = 0.291999936103821 sec.
Hardness ratio (energy fluence ratio) = 1.07314

Color vs. Color plot

Count Ratio  (25-50 keV) / (15-25 keV) = 0.941912
Count Ratio (50-100 keV) / (15-25 keV) = 0.669891

Mask shadow pattern

IMX = -5.237457822660342E-01, IMY = -6.751795674553224E-02

TIME vs. PHA plot around the trigger time

Blue+Cyan dotted lines: Time interval of the scaled map
From 0.000000 to 0.512000

TIME vs. DetID plot around the trigger time

Blue+Cyan dotted lines: Time interval of the scaled map
From 0.000000 to 0.512000

100 us light curve (15-350 keV)

Blue+Cyan dotted lines: Time interval of the scaled map
From 0.000000 to 0.512000

Spacecraft aspect plot

Bright source in the pre-burst/pre-slew/post-slew images

Pre-burst image of 15-350 keV band

Time interval of the image:
-495.767400 -480.108000
-239.968000 0.016000
 
# RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
#  [deg]    [deg]    ['] [deg]  [deg]                   [']
278.3160  47.1942    3.0  50.4   36.7   4.0124   3.8 ------	UNKNOWN                       
205.8270  80.7764    2.5   9.8   47.0   0.9179   4.5 ------	UNKNOWN                       
147.1494  23.3067    3.0  55.7  175.5   4.4386   3.9 ------	UNKNOWN                       
117.4420  62.8873    2.1  18.5 -138.4   1.1217   5.4 ------	UNKNOWN                       
 69.7350  61.7442    2.7  28.0  -86.5   1.2111   4.3 ------	UNKNOWN

Pre-slew background subtracted image of 15-350 keV band

Time interval of the image:
0.016000 0.332000
 
# RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
#  [deg]    [deg]    ['] [deg]  [deg]                   [']
202.7531  55.1971    2.3  29.7  104.3   0.0543   5.1 ------	UNKNOWN                       
149.4383  51.1967    1.2  27.8  172.7   0.0780  10.0 ------	UNKNOWN                       
113.9782  30.9359    2.1  50.2 -146.8   0.1035   5.6 ------	UNKNOWN

Post-slew image of 15-350 keV band

Time interval of the image:
67.532620 962.058500
 
# RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
#  [deg]    [deg]    ['] [deg]  [deg]                   [']
190.8693  25.4207    3.6  40.6  102.4   7.3309   3.3 ------	UNKNOWN                       
287.9650  69.8051    2.3  55.3    5.1   9.8009   5.0 ------	UNKNOWN                       
121.0864  -1.4623    5.0  57.9 -156.9   9.2229   2.3 ------	UNKNOWN                       
 13.0447  67.6314    4.0  56.8  -29.3  10.4589   2.9 ------	UNKNOWN                       
 77.1301  53.0980    4.4  42.5  -68.9   6.7599   2.7 ------	UNKNOWN

Plot creation:
Tue Jun 21 10:20:39 EDT 2016