batgrbproduct analysis for trigger=581129

batgrbproduct Analysis for Trigger=581129

Analysis Summary

Notes:

1) No attempt has been made to eliminate the non-significant decimal places.
2) This batgrbproduct analysis is done of the amount of data for this trigger downloaded to date.
3) This batgrbproduct analysis is using the flight position for the mask-weighting calculation.
4) This page will be updated for up to 10 downlinks after the trigger.

======================================================================
                  BAT GRB EVENT DATA PROCESSING REPORT
    Process Script Ver: batgrbproduct v2.46
    Process Start Time (local time): Sun Dec 15 15:17:29 EST 2013
             Host Name: gcn5
   Current Working Dir: /local/data/raid16/prebascript_archive/trigger581129/GCN5_download_cron
======================================================================

  Trigger: 00581129      Segment: 000

  BAT Trigger Time: 408770814.452   [s; MET]
      Trigger Stop: 408770814.516   [s; MET]
               UTC: 2013-12-15T03:26:43.056740    [includes UTCF correction] 
       Where From?: (*SHORT BURST*) TDRSS scaled map message

  BAT 
    RA: 178.341830838294   Dec: 47.9423038496901  [deg; J2000]

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

  Flight Position:   [ source = BAT ]
    RA: 178.341830838294   Dec: 47.9423038496901  [deg; J2000]

  Ground Position:   [ source = BAT pre-slew burst ]
    RA: 178.376472005539   Dec: 47.9193728138365  [deg; J2000]  
        { 11h 53m 30.4s , +47d 55' 09.7" }
         +/- 3.55171098128217 [arcmin] (estimated 90% radius based on SNR)
         +/- 2.34339544460239 [arcmin] (formal 1-sigma fit error)
    SNR: 4.97483780106625
 
  Angular difference between Flight and Ground position is 1.957656 arcmin
 
  Partial Coding Fraction: 0.7734375    [ including projection effects ]

  Duration
     T90: 0.0639999508857727 +/- 
        Measured from: 0
                   to: 0.0639999508857727     [s; relative to TRIGTIME]
     T50: 0.0639999508857727 +/- 
        Measured from: 0
                   to: 0.0639999508857727     [s; relative to TRIGTIME]

    WARNING: battblocks failed.  Used TDRSS fore/background intervals

  Fluence
     Peak Flux  (peak 0.0639999508857727 second)
        Measured from: -2.98023223876953e-08     
                   to: 0.0639999210834503     [s; relative to TRIGTIME]
     Total Fluence        
        Measured from: 0     
                   to: 0.0639999508857727     [s; relative to TRIGTIME]
     
             Band 1     Band 2     Band 3     Band 4
              15-25      25-50     50-100    100-350     keV
 Total     0.001433   0.001514   0.004391   0.005677
           0.001584   0.001647   0.001657   0.001739 [error]
 Peak      0.001433   0.001514   0.004391   0.005677
           0.001584   0.001647   0.001657   0.001739 [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 5.46e-10 4.88e-10 2.96e-10 2.70e-10 1.45e-09 1.65e-09   erg/cm2
              Error 6.14e-10 5.25e-10 4.17e-10 1.88e-09 1.59e-09 1.92e-09   erg/cm2
Peak flux Power-law 2.82e-01 1.41e-01 4.29e-02 1.64e-02 4.77e-01 4.85e-01   ph/cm2/s
              Error 3.20e-01 1.50e-01 5.91e-02 6.88e-02 5.04e-01 5.12e-01   ph/cm2/s
 
Fluence   Cutoff-PL 1.68e-11 8.53e-10 2.08e-09 3.72e-10 3.02e-09 3.44e-09   erg/cm2
              Error 8.96e-11 7.81e-10 1.30e-09 4.63e-10 1.84e-09 1.85e-09   erg/cm2
Peak flux Cutoff-PL 7.28e-03 2.03e-01 3.09e-01 3.22e-02 5.27e-01 5.28e-01   ph/cm2/s
              Error 4.25e-02 1.92e-01 1.90e-01 3.31e-02 3.26e-01 3.26e-01   ph/cm2/s
 
Fluence        Band 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00   erg/cm2
              Error 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00   erg/cm2
Peak flux      Band 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00   ph/cm2/s
              Error 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00   ph/cm2/s

 
=== Time-averaged spectral analysis ===
 
Model Ph_index Ph_90err Epeak  Epeak_90err Norm     Nomr_90err chi2   dof
PL    2.677    0.000    -      -           2.17e-03 2.55e-03   56.290 57
CPL -10.000    0.000      61.9    9.0      2.06e+02 3.78e+02   50.920 56
 
The chi2 of BB fit is  which shows a significant improvement in the fit over a PL.
Check the BB fit.
 
=== 1-s peak spectral analysis ===
 
Model Ph_index Ph_90err Epeak  Epeak_90err Norm     Nomr_90err chi2   dof
PL    2.677    0.000    -      -           2.17e-03 2.55e-03   56.290 57
CPL -10.000    0.000      61.9    9.0      2.06e+02 3.78e+02   50.920 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.064    0.000    0.064
T90     0.064    0.000    0.064
T50     0.064    0.000    0.064
Peak    0.064   -0.000    0.064

Click here for the BAT GCN circular template

Image

Pre-slew 15-350 keV image (Event data; bkg subtracted)

#   RAcent  DECcent POSerr  Theta    Phi     Peak Cts     SNR  Name
  178.3765 +47.9194 0.0387   30.6 -177.2       0.0159     5.0  TRIG_00581129       
------------------------------------------
 
Foreground time interval of the image:
  0.000   0.064 (delta_t =   0.064 [sec])
Background time interval of the image:
-239.452 -47.142 (delta_t = 192.310 [sec])
-46.452  -6.452 (delta_t =  40.000 [sec])

Lightcurves

Notes:

1) All plots contain as much data as has been downloaded to date.
2) The mask-weighted light curves are using the flight position.
3) Multiple plots of different time binning/intervals are shown to cover all scenarios of short/long GRBs, rate/image triggers, and real/false positives.
4) 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.
5) 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.
6) 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.064000

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.064000

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.064000

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

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

The mask-weighted light curves based on cleaned DPHs

Spectra

Notes:

1) The fitting includes the systematic errors.
2) For long bursts, a spectral fit of 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.
4) For short bursts (T90<1sec), the specrtum is also fit with Blackbody, OTTB, and Double Blackbody.

Time averaged spectrum fit using the pre-slew DRM

Power-law model

Time interval is from    0.000 sec. to    0.064 sec.
 
Spectral model in power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 2.67743     ( -2.7557   -2.7557 )
Norm@50keV  : 2.16729E-03    ( -0.00214664   0.00295702 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          56.29 using 59 PHA bins.
# Reduced chi-squared =         0.9876 for     57 degrees of freedom 
# Null hypothesis probability =   5.016058e-01
 
Photon flux (15-150 keV)    in 0.064 sec: 0.477201 ( -0.477201  0.531092 ) ph/cm2/s
Energy fluence (15-150 keV)             : 1.44605e-09 ( -1.44605e-09  1.73437e-09 ) ergs/cm2

Cutoff power-law model

Time interval is from    0.000 sec. to    0.064 sec.
Spectral model in the cutoff power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: -9.99996     (     )
Epeak [keV] : 61.9128      ( -9.56592 8.49554 )
Norm@50keV  : 206.267    ( -206.278   549.644 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          50.92 using 59 PHA bins.
# Reduced chi-squared =         0.9093 for     56 degrees of freedom 
# Null hypothesis probability =   6.669382e-01
 
Photon flux (15-150 keV)    in 0.064 sec: 0.527324 ( -0.324074  0.3274 ) ph/cm2/s
Energy fluence (15-150 keV)             : 3.02132e-09 ( -1.84086e-09  1.84588e-09 ) ergs/cm2

Band function

Time interval is from    0.000 sec. to    0.064 sec.
 
Spectral model in the Band function:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
alpha       : 5.00000   ( -4.03567   -5 )
beta        : -10.0000    ( 10   10 )
Epeak [keV] : 56.3915      ( -16.456 14.3245 )
Norm@50keV   : 5.02673    ( -4.98873   31.6571 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          53.33 using 59 PHA bins.
# Reduced chi-squared =         0.9695 for     55 degrees of freedom 
# Null hypothesis probability =   5.388556e-01
 
Photon flux (15-150 keV)    in 0.064 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]    : 12.4509  (-12.393  6.22579)
R^2/D10^2   : 0.151172  (-0.151369   14.3648)
(R is the radius in km and D10 is the distance to the source in units of 10 kpc)
------------------------------------------------------------
#Fit statistic : Chi-Squared =          54.38 using 59 PHA bins.
# Reduced chi-squared =         0.9540 for     57 degrees of freedom 
# Null hypothesis probability =   5.741260e-01
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  3.179e-10   2.532e-10
25- 50  1.044e-09   7.472e-10
50-150  9.968e-10   1.069e-09
15-150  2.359e-09   1.900e-09

Thermal bremsstrahlung (OTTB)

 
 
Spectral model: thermal bremsstrahlung
------------------------------------------------------------
Parameters  : value         Lower 90%    Upper 90%
kT [keV]    : 37.4317       (-36.7483      -36.7483)
Norm        : 3.97307       (-3.99825     1672.22)
------------------------------------------------------------
#Fit statistic : Chi-Squared =          55.54 using 59 PHA bins.
# Reduced chi-squared =         0.9743 for     57 degrees of freedom 
# Null hypothesis probability =   5.301828e-01
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  5.703e-10   4.326e-10
25- 50  7.285e-10   6.045e-10
50-150  5.233e-10   6.314e-10
15-150  1.822e-09   1.489e-09

Double BB

 
 
Spectral model: bbodyrad<1> + bbodyrad<2>
------------------------------------------------------------
Parameters  : value         Lower 90%    Upper 90%
kT1 [keV]   : 3.19768        (-3.19266      15.1846)
R1^2/D10^2  : 9.50609        (-9.65033      -9.65033)
kT2 [keV]   : 12.6811        (-12.6815      -12.6815)
R2^2/D10^2  : 0.130786    (-0.130776   9.74741)
------------------------------------------------------------
#Fit statistic : Chi-Squared =          54.27 using 59 PHA bins.
# Reduced chi-squared =         0.9867 for     55 degrees of freedom 
# Null hypothesis probability =   5.025667e-01
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  4.527e-10   3.312e-10
25- 50  9.871e-10   8.173e-10
50-150  9.585e-10   8.250e-10
15-150  2.398e-09   2.239e-09

Peak spectrum fit

Power-law model

Time interval is from    0.000 sec. to    0.064 sec.
 
Spectral model in power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: 2.67743     ( -2.7557   -2.7557 )
Norm@50keV  : 2.16729E-03    ( -0.00214664   0.00295702 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          56.29 using 59 PHA bins.
# Reduced chi-squared =         0.9876 for     57 degrees of freedom 
# Null hypothesis probability =   5.016058e-01
 
Photon flux (15-150 keV)    in 0.064 sec: 0.477201 ( -0.477201  0.531092 ) ph/cm2/s
Energy fluence (15-150 keV)             : 1.44605e-09 ( -1.44605e-09  1.73437e-09 ) ergs/cm2

Cutoff power-law model

Time interval is from    0.000 sec. to    0.064 sec.
Spectral model in the cutoff power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: -9.99996     (     )
Epeak [keV] : 61.9128      ( -9.56592 8.49554 )
Norm@50keV  : 206.267    ( -206.278   549.644 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          50.92 using 59 PHA bins.
# Reduced chi-squared =         0.9093 for     56 degrees of freedom 
# Null hypothesis probability =   6.669382e-01
 
Photon flux (15-150 keV)    in 0.064 sec: 0.527324 ( -0.324074  0.3274 ) ph/cm2/s
Energy fluence (15-150 keV)             : 3.02132e-09 ( -1.84086e-09  1.84588e-09 ) ergs/cm2

Band function

Time interval is from    0.000 sec. to    0.064 sec.
 
Spectral model in the Band function:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
alpha       : 5.00000   ( -4.03567   -5 )
beta        : -10.0000    ( 10   10 )
Epeak [keV] : 56.3915      ( -16.456 14.3245 )
Norm@50keV   : 5.02673    ( -4.98873   31.6571 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          53.33 using 59 PHA bins.
# Reduced chi-squared =         0.9695 for     55 degrees of freedom 
# Null hypothesis probability =   5.388556e-01
 
Photon flux (15-150 keV)    in 0.064 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]    : 12.4509  (-12.393  6.22579)
R^2/D10^2   : 0.151172  (-0.151369   14.3648)
(R is the radius in km and D10 is the distance to the source in units of 10 kpc)
------------------------------------------------------------
#Fit statistic : Chi-Squared =          54.38 using 59 PHA bins.
# Reduced chi-squared =         0.9540 for     57 degrees of freedom 
# Null hypothesis probability =   5.741260e-01
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  3.179e-10   2.564e-10
25- 50  1.044e-09   8.464e-10
50-150  9.968e-10   1.034e-09
15-150  2.359e-09   1.958e-09

Thermal bremsstrahlung (OTTB)

 
 
Spectral model: thermal bremsstrahlung
------------------------------------------------------------
Parameters  : value         Lower 90%    Upper 90%
kT [keV]    : 37.4317       (-36.7483      -36.7483)
Norm        : 3.97307       (-3.99825     1672.22)
------------------------------------------------------------
#Fit statistic : Chi-Squared =          55.54 using 59 PHA bins.
# Reduced chi-squared =         0.9743 for     57 degrees of freedom 
# Null hypothesis probability =   5.301828e-01
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  5.703e-10   4.582e-10
25- 50  7.285e-10   6.493e-10
50-150  5.233e-10   6.867e-10
15-150  1.822e-09   1.700e-09

Double BB

 
 
Spectral model: bbodyrad<1> + bbodyrad<2>
------------------------------------------------------------
Parameters  : value         Lower 90%    Upper 90%
kT1 [keV]   : 3.19768        (-3.19266      15.1846)
R1^2/D10^2  : 9.50609        (-9.65033      -9.65033)
kT2 [keV]   : 12.6811        (-12.6815      -12.6815)
R2^2/D10^2  : 0.130786    (-0.130776   9.74741)
------------------------------------------------------------
#Fit statistic : Chi-Squared =          54.27 using 59 PHA bins.
# Reduced chi-squared =         0.9867 for     55 degrees of freedom 
# Null hypothesis probability =   5.025667e-01
 
Energy  Fluence     90% Error
[keV]   [erg/cm2]   [erg/cm2]
15- 25  4.527e-10   2.568e-10
25- 50  9.871e-10   7.962e-10
50-150  9.585e-10   9.117e-10
15-150  2.398e-09   2.042e-09

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: 2.67743     ( -2.7557   -2.7557 )
Norm@50keV  : 2.16729E-03    ( -0.00214664   0.00295702 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          56.29 using 59 PHA bins.
# Reduced chi-squared =         0.9876 for     57 degrees of freedom 
# Null hypothesis probability =   5.016058e-01
 
Photon flux (15-150 keV)    in 0.064 sec: 0.477201 ( -0.477201  0.531092 ) ph/cm2/s
Energy fluence (15-150 keV)             : 1.44605e-09 ( -1.44605e-09  1.73437e-09 ) ergs/cm2

Cutoff power-law model

Spectral model in the cutoff power-law:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
Photon index: -9.99996     (     )
Epeak [keV] : 61.9128      ( -9.56592 8.49554 )
Norm@50keV  : 206.267    ( -206.278   549.644 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          50.92 using 59 PHA bins.
# Reduced chi-squared =         0.9093 for     56 degrees of freedom 
# Null hypothesis probability =   6.669382e-01
 
Photon flux (15-150 keV)    in 0.064 sec: 0.527324 ( -0.324074  0.3274 ) ph/cm2/s
Energy fluence (15-150 keV)             : 3.02132e-09 ( -1.84086e-09  1.84588e-09 ) ergs/cm2

Band function

 
Spectral model in the Band function:
------------------------------------------------------------
Parameters  : value       lower 90%   higher 90%
alpha       : 5.00000   ( -4.03567   -5 )
beta        : -10.0000    ( 10   10 )
Epeak [keV] : 56.3915      ( -16.456 14.3245 )
Norm@50keV   : 5.02673    ( -4.98873   31.6571 )
------------------------------------------------------------
#Fit statistic : Chi-Squared =          53.33 using 59 PHA bins.
# Reduced chi-squared =         0.9695 for     55 degrees of freedom 
# Null hypothesis probability =   5.388556e-01
 
Photon flux (15-150 keV)    in 0.064 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)  = 4.87561e-10
S(50-100 keV) = 2.96193e-10

T90 vs. Hardness ratio plot

T90                                   = 0.0639999508857727 sec.
Hardness ratio (energy fluence ratio) = 0.607499

Color vs. Color plot

Count Ratio  (25-50 keV) / (15-25 keV) = 1.05652
Count Ratio (50-100 keV) / (15-25 keV) = 3.0642

Mask shadow pattern

IMX = 0, IMY = 0

TIME vs. PHA plot around the trigger time

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

TIME vs. DetID plot around the trigger time

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

100 us light curve (15-350 keV)

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

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:
-239.452000 -47.142000
-46.452000 -6.452000
 
# RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
#  [deg]    [deg]    ['] [deg]  [deg]                   [']
214.8265  58.9856    3.6  54.3 -171.7   2.0284   3.2 ------	UNKNOWN                       
211.8137  55.1151    3.1  51.5 -167.7   1.9578   3.7 ------	UNKNOWN                       
159.9985 -19.4279    2.4  41.1  -25.3   1.3107   4.8 ------	UNKNOWN                       
200.2340  20.5895    2.9  38.6 -119.7   1.6968   4.0 ------	UNKNOWN                       
198.3246   6.6754    3.0  41.1  -98.1   3.8731   3.9 ------	UNKNOWN

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

Time interval of the image:
0.000000 0.064000
 
# RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
#  [deg]    [deg]    ['] [deg]  [deg]                   [']
165.4808 -17.2938    2.6  39.5  -33.6   0.0191   4.5 ------	UNKNOWN

Post-slew image of 15-350 keV band

Time interval of the image:
64.448600 962.578000
 
# RAcent  DECcent POSerr Theta    Phi  PeakCts   SNR AngSep     Name
#  [deg]    [deg]    ['] [deg]  [deg]                   [']
 55.1871  77.1716    3.5  50.0  146.3   8.4462   3.3 ------	UNKNOWN                       
167.4777   6.6449    3.0  42.3   -3.6   2.8593   3.9 ------	UNKNOWN                       
178.8206  38.4895    2.6   9.5  -22.1   2.0925   4.4 ------	UNKNOWN                       
217.4376  58.1239    3.4  25.2 -148.2   1.7467   3.4 ------	UNKNOWN

Plot creation:
Sun Dec 29 23:53:18 EST 2013