High Source Density Bins Description and Links

[Adapted from the Explanatory Supplement.]

The file consists of 80-character ASCII records, with each record containing the information for two bins (Table 1). The flag HSD is hex-encoded by band, as described for the Format of the Point Source Catalog. In this notation each band corresponds to one of the four bits of a hex digit. 12 microns corresponds to bit 0 (Least significant Bit) and 100 microns to bit 3, etc. If a band went through high source density processing, then the appropriate bit in the hex digit is set. Thus, if high source density rules were invoked at 25, 60, and 100 microns, HSD would have the value 1110 (binary) = E (hex). Ecliptic bins start at the ecliptic north pole and step around the sky in bands of constant ecliptic latitude, stepping 1° southward after completing each band. The length of the bin in ecliptic longitude is adjusted for the cosine of the ecliptic latitude to maintain an approximately constant area.

         
Table 1.  Format of the High Source Density Bins File

 Start    Name            Description        Units       Format
 00       BINNUM        SDAS Bin Number      ---           I6
 06       HSD           Bands processed      ---           A2
                        for high source
                        density, hex
                        encoded by band
 08       LAMBDA        Ecliptic longitude   deg           F10.5
                        of bin center
 18       BETA          Ecliptic latitude    deg           F10.5
                        of bin center
 28       LENGTH        Length of bin in     deg           F10.5
                        ecliptic longitude
 38       SPARE         2 spare bytes        ---           2A1
40-79                   Bytes 00-39 are
                        repeated for next bin

FORTRAN programs are given below to convert from bin number to ecliptic coordinates and vice versa. There is a known bug in the computation of bins at ecliptic latitude 60°, causing bin 2842 to be skipped. The following program, while incorrect in this way, will give results consistent with the bin numbers used in the data processing. To avoid problems with roundoff errors within a few arcseconds of bin boundaries, all arithmetic should be calculated in double precision.

Right Click here to download the programs shown below.

C CONVERT AN ECLIPTIC POSITION TO BIN NUMBER
SUBROUTINE P2BIN(LAM,BET,BIN)
IMPLICIT REAL*8 A-H,P-Z
REAL*8 LAM,BET,LSIZE,R2D/57.2957795/
INTEGER*4 BIN,I,J,N,MAXBIN(182)
COMMON /LATBND/MAXBIN
C
I = 90. - BET*R2D + 2.5
N = MAXBIN(I) - MAXBIN(I-1)
LSIZE = 360./FLOAT(N)
J = LAM*R2D/LSIZE +1
BIN = MAXBIN(I-1) + J
RETURN
END

C CONVERT A BIN NUMBER TO A CENTER POSITION AND LONGITUDE LENGTH
SUBROUTINE BIN2P(BIN,CLON,CLAT,LWIDTH)
IMPLICIT REAL*8 A-H,P-Z
INTEGER*4 BIN,I,J,MAXMIN(182),FIRST
REAL*8 CLON,CLAT,LWIDTH,R2D/57.2957795/,NBINS
COMMON /LATBND/MAXBIN
C
DO 100 I=1,182
IF(MAXBIN(I) .GE. BIN) GOTO 150
100 CONTINUE
150 I = I -1
FIRST = MAXBIN(I) + 1
NBINS = MAXBIN(I+1) - FIRST + 1
LWIDTH = 360./NBINS
J = BIN - FIRST
CLON = DFLOAT(J)*LWIDTH + LWIDTH/2.
CLAT = 90. - DFLOAT(I-1)
LWIDTH = LWIDTH/R2D
CLON = CLON/R2D
CLAT = CLAT/R2D
RETURN
END
C
C DATA SUBPROGRAM
BLOCK DATA
INTEGER*4 MAXBIN(182)
DATA MAXBIN
* /0,1,7,19,37,62,93,130,173,223,279,341,409,563,
* 650,743,843,947,1058,1175,1298,1427,1561,1701,1847,1999,
* 2156,2319,2488,2662,2841,3026,3216,3412,3613,3819,4030,
* 4246,4467,4693,4924,5160,5400,5645,5895,6149,6407,6670,
* 6937,7208,7483,7762,8045,8332,8623,8917,9215,9516,9821,
* 10129,10440,10754,11071,11391,11714,12040,12368,12699,
* 13032,13368,13706,14046,14388,14732,15078,15425,15774,
* 16124,16476,16829,17183,17538,17894,18251,18609,18967,
* 19326,19685,20044,20403,20763,21122,21481,21840,22199,
* 22557,22915,23272,23628,23983,24337,24690,25042,25392,
* 25741,26088,26434,26778,27120,27460,27798,28134,28467,
* 28798,29126,29452,29775,30095,30412,30726,31037,31345,
* 31650,31951,32249,32543,32834,33121,33404,33683,33958,
* 34229,34496,34759,35017,35271,35521,35766,36006,36242,
* 38325,38505,38679,38848,39011,39168,39320,39466,39606,
* 39740,39869,39992,40109,40220,40325,40424,40517,40604,
* 40684,40758,40826,40888,40944,40994,41037,41074,41105,
* 41130,41148,41160,41166,41167/
COMMON /LATBND/MAXBIN
END

IRAS High Source Density Bins Description

A service of the HEASARC and of the Astrophysics Science Division at NASA/GSFC
Goddard Space Flight Center, National Aeronautics and Space Administration
HEASARC Director: Dr. Andrew F. Ptak
LAMBDA Director: Dr. Thomas M. Essinger-Hileman
NASA Official: Dr. Thomas M. Essinger-Hileman
Web Curator: Mr. Michael R. Greason