Conversion of standardized ReadMe file for
file /./ftp/cats/IX/22 into FORTRAN code for loading all data files into arrays.
Note that special values are assigned to unknown or unspecified
numbers (also called NULL numbers);
when necessary, the coordinate components making up the right ascension
and declination are converted into floating-point numbers
representing these angles in degrees.
program load_ReadMe
C=============================================================================
C F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2024-Apr-19
C=============================================================================
* This code was generated from the ReadMe file documenting a catalogue
* according to the "Standard for Documentation of Astronomical Catalogues"
* currently in use by the Astronomical Data Centers (CDS, ADC, A&A)
* (see full documentation at URL http://vizier.u-strasbg.fr/doc/catstd.htx)
* Please report problems or questions to
C=============================================================================
implicit none
* Unspecified or NULL values, generally corresponding to blank columns,
* are assigned one of the following special values:
* rNULL__ for unknown or NULL floating-point values
* iNULL__ for unknown or NULL integer values
real*4 rNULL__
integer*4 iNULL__
parameter (rNULL__=--2147483648.) ! NULL real number
parameter (iNULL__=(-2147483647-1)) ! NULL int number
integer idig ! testing NULL number
C=============================================================================
Cat. IX/22 Einstein extended deep survey (Primini+, 1991)
*================================================================================
*The CfA Eisntein observatory extended deep X-ray survey
* Primini F.A., Murray S.S., Huchra J., Schild R., Burg R., Giacconi R.
* <Astrophys. J. 374, 440 (1991)>
* =1991ApJ...374..440P
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'targets.dat' ! EDS Target Information
integer*4 nr__
parameter (nr__=43) ! Number of records
character*40 ar__ ! Full-size record
C B1950 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg (nr__) ! (deg) Right Ascension B1950
real*8 DEdeg (nr__) ! (deg) Declination B1950
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 Seqno (nr__) ! Sequence number (G1)
character*3 detector (nr__) ! Detector (IPC or HRI)
integer*4 RAh (nr__) ! (h) Right Ascension (B1950) (hour)
integer*4 RAm (nr__) ! (min) Right Ascension (B1950) (min)
integer*4 RAs (nr__) ! (s) Right Ascension (B1950) (sec)
character*1 DE_ (nr__) ! Declination sign (B1950)
integer*4 DEd (nr__) ! (deg) Declination (B1950) (deg)
integer*4 DEm (nr__) ! (arcmin) Declination (B1950) (arcmin)
integer*4 DEs (nr__) ! (arcsec) Declination (B1950) (arcsec)
real*8 exptime (nr__) ! (s) Exposure time (2)
*Note (2): This column gives the net time, in seconds, for the processed
* image (uncorrected for deadtime effects).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hfield.dat' ! EDS Field Parameters for HRI Observations
integer*4 nr__1
parameter (nr__1=34) ! Number of records
character*110 ar__1 ! Full-size record
C B1950 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_1 (nr__1) ! (deg) Right Ascension B1950
real*8 DEdeg_1 (nr__1) ! (deg) Declination B1950
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 Hseqno (nr__1) ! HRI Sequence Number (G1)
real*8 ontime (nr__1) ! (s) Exposure Time
real*4 globk (nr__1) ! (ct) Background in 12"x12" detect cell (in 9' circle)
real*4 detbk (nr__1) ! (ct) Background in 12"x12" detect cell (in 5' circle)
real*8 th12 (nr__1) ! (ct) Threshold counts for source detection (12"x12"
* cell)
real*8 th24 (nr__1) ! (ct) Threshold counts in 24"x24" cell
real*8 th36 (nr__1) ! (ct) Threshold counts in 36"x36" cell
real*8 th48 (nr__1) ! (ct) Threshold counts in 48"x48" cell
real*8 th72 (nr__1) ! (ct) Threshold counts in 72"x72" cell
real*8 th120 (nr__1) ! (ct) Threshold counts in 120"x120" cell
integer*4 RAh_1 (nr__1) ! (h) Right Ascension (B1950) (hour) (1)
integer*4 RAm_1 (nr__1) ! (min) Right Ascension (B1950) (min)
integer*4 RAs_1 (nr__1) ! (s) Right Ascension (B1950) (sec)
character*1 DE__1 (nr__1) ! Declination sign (B1950)
integer*4 DEd_1 (nr__1) ! (deg) Declination (B1950) (deg)
integer*4 DEm_1 (nr__1) ! (arcmin) Declination (B1950) (arcmin)
integer*4 DEs_1 (nr__1) ! (arcsec) Declination (B1950) (arcsec)
real*8 nroll (nr__1) ! (deg) Nominal Roll Angle (G2)
*Note (1): These columns give the nominal right ascension and declination
* (B1950) of the field. The nominal position is the center of the image
* which differs slightly from the center of the field of view due to
* aspect and boresight offsets.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hfluxes.dat' ! EDS Secondary HRI Source Parameters
integer*4 nr__2
parameter (nr__2=210) ! Number of records
character*332 ar__2 ! Full-size record
integer*4 Hseqno_1 (nr__2) ! HRI sequence number (G1)
integer*4 Hsrcno (nr__2) ! HRI source number (1)
integer*4 cts0to3 (nr__2) ! (ct) Counts within 0-3" of source position
integer*4 cts3to9 (nr__2) ! (ct) Counts within 3-9" of source position
integer*4 cts9to18 (nr__2) ! (ct) Counts within 9-18" of source position
real*4 USProb (nr__2) ! Probability that the source is unresolved
integer*4 ctR (10,nr__2) ! (ct) Counts in annuli containing ~10% of
* total src cts
real*4 radR (10,nr__2) ! (arcsec) Outer radii of the annular regions
real*4 denR (10,nr__2) ! (ct/arcsec2) Photon densities in the annular regions
real*4 CR_S_PS (nr__2) ! (ct/ks) Count rate for assumed soft thermal point
* source (2)
real*4 FL_S_PS (nr__2) ! (mW/m2) Corresponding flux
real*4 CR_H_PS (nr__2) ! (ct/ks) C-rate for assumed hard thermal pt src (3)
real*4 FL_H_PS (nr__2) ! (mW/m2) Corresponding flux
real*4 xtnt (nr__2) ! (arcsec) Best estimate of source extent
real*4 CR_S_XS (nr__2) ! (ct/ks) C-rate for soft thermal source
* within estimated extent (2)
real*4 FL_S_XS (nr__2) ! (mW/m2) Corresponding flux
real*4 CR_H_XS (nr__2) ! (ct/ks) Count rate for hard thermal source
* within estimated extent (3)
real*4 FL_H_XS (nr__2) ! (mW/m2) Corresponding flux
real*4 ulxtnt (nr__2) ! (arcsec) Upper limit on source extent
real*4 CR_S_UXS (nr__2) ! (ct/ks) Count rate for soft thermal source
* within upper limit extent (2)
real*4 FL_S_UXS (nr__2) ! (mW/m2) Corresponding flux
real*4 CR_H_UXS (nr__2) ! (ct/ks) Count rate for hard thermal source
* within upper limit extent (3)
real*4 FL_H_UXS (nr__2) ! (mW/m2) Corresponding flux
real*4 llxtnt (nr__2) ! (arcsec) Lower limit on source extent
real*4 CR_S_LXS (nr__2) ! (ct/ks) Count rate for soft thermal source
* within lower limit extent (2)
real*4 FL_S_LXS (nr__2) ! (mW/m2) Corresponding flux
real*4 CR_H_LXS (nr__2) ! (ct/ks) Count rate for hard thermal source
* within lower limit extent (3)
real*8 FL_H_LXS (nr__2) ! (mW/m2) Corresponding flux
*Note (1): The 'Hsrcno' column gives the HRI source number within the field.
*Note (2): Assumed Thermal Bremmstrahlung spectrum at kT=0.3 keV
*Note (3): Assumed Thermal Bremmstrahlung spectrum at kT=3.0 keV
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hsrc.dat' ! EDS HRI Source Positions and Count Rates
integer*4 nr__3
parameter (nr__3=210) ! Number of records
character*174 ar__3 ! Full-size record
C B1950 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_2 (nr__3) ! (deg) Right Ascension B1950
real*8 DEdeg_2 (nr__3) ! (deg) Declination B1950
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 Hseqno_2 (nr__3) ! HRI Sequence Number (G1)
integer*4 Hsrcno_1 (nr__3) ! HRI Source Number (1)
integer*4 csize (nr__3) ! (arcsec) Detection cell size
real*8 y (nr__3) ! (pix) Position in y
real*8 z (nr__3) ! (pix) Position in z
real*4 poserr (nr__3) ! (arcsec) Position error
integer*4 RAh_2 (nr__3) ! (h) Right Ascension (B1950) (hour)
integer*4 RAm_2 (nr__3) ! (min) Right Ascension (B1950) (min)
integer*4 RAs_2 (nr__3) ! (s) Right Ascension (B1950) (sec)
character*1 DE__2 (nr__3) ! Declination sign (B1950)
integer*4 DEd_2 (nr__3) ! (deg) Declination (B1950) (deg)
integer*4 DEm_2 (nr__3) ! (arcmin) Declination (B1950) (arcmin)
integer*4 DEs_2 (nr__3) ! (arcsec) Declination (B1950) (arcsec)
real*4 netcr (nr__3) ! (ct/ks) Corrected Count Rate per 1000 sec
real*4 e_netcr (nr__3) ! (ct/ks) Error on Count Rate
real*4 bpcts (nr__3) ! (ct) Counts in detect cell centered at best
* position (2)
real*4 eqcts (nr__3) ! (ct) Counts in circle of equivalent area (2)
real*4 bp18cts (nr__3) ! (ct) Counts in standard 18" circle at best
* position (3)
real*8 bp18bg (nr__3) ! (ct) Background in 18" circle (3)
real*4 bp18scat (nr__3) ! Fraction of src cnts scattered out of
* 18" circle (3)
real*4 vcorr (nr__3) ! Telescope vignetting correction at
* source position (4)
real*4 dtfac (nr__3) ! Dead time correction factor (4)
real*4 qec (nr__3) ! Detector quantum efficiency at source
* position (4)
real*8 exptime_1 (nr__3) ! (s) Exposure time at source position (5)
integer*4 xccts (nr__3) ! (ct) Counts in extended circle (6)
real*8 xcbg (nr__3) ! (ct) Background in extended circle (6)
real*8 xcscat (nr__3) ! (ct) Fraction of src cts scattered out of
* extended circle (6)
real*8 xcrad (nr__3) ! (arcsec) Radius of extended circle (6)
real*8 sizcor (nr__3) ! Source xtent correction factor for cts/ks (7)
*Note (1): The 'Hsrcno' column gives the HRI source number within the field.
*Note (2): The number of counts in a detect cell centered at the best
* position for each HRI source is given in the 'bpcts' columns. In the
* 'eqcts' column, the number of counts in a circle of equivalent area is
* given.
*Note (3): Source parameters are also calculated for each source at its
* best position in a standard circle with a radius of 18". The 'bp18cts'
* column lists source counts in the 18" circle, and the 'bp18bg' lists
* the expected background value in this circle determined from the
* background in a 5' region about the target center, scaled by the
* circle size and the exposure at the source location. The 'bp18scat'
* column lists the fraction of total source counts scattered out of the
* 18" circle.
*Note (4): The 'vcorr' and 'dtfac' columns list corrections made to the
* observation time; 'vcorr' gives the telescope vignetting correction at
* the source position and 'dtfac' gives the dead time correction factor
* at the time of the observation. The 'qec' column gives the detector
* quantum efficiency at the source position.
*Note (5): This column gives the total exposure time, in seconds, at the
* source position. The exposure time is corrected for intervals lost due
* to unacceptable aspect solution, detector dead time, and other event
* screening. It may differ from the net time in the processed image if
* the source is not in the field of view for the entire observation.
*Note (6): These columns give information about extended sources. The
* number of counts in an extended circle (that circle which
* circumscribes all detect cells above threshold) is given in 'xccts'.
* The background in this extended circle is given by 'xcbg'. 'xcscat'
* gives the fraction of source counts scattered out of the extended
* circle; and 'xcrad' gives the radius of the extended circle in
* arcseconds.
*Note (7): This column give the size correction factor, 'sizcor', a
* parameter developed to adjust intensities for extended sources. Since
* Rev 1 did not attempt to provide estimates of count rates or fluxes of
* extended sources, we have developed a new parameter, sizcor, which is
* conceptually similar to the classical beam broadening correction used
* by radio astronomers when performing drift scans with single-dish
* telescopes. Sizcor serves to correct all intensity measurements
* (optimized for point sources) of resolved sources for the X-ray
* brightness which is outside of the detection box. 'sizcor' can be used
* with any intensity measurement (counts, count rate, or flux) to obtain
* an estimate for extended sources. For a more complete description of
* sizcor, see the table <IX/23/source.dat>.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hris_sep.dat' ! EDS Separations Between HRI Srcs in a Single Field
integer*4 nr__4
parameter (nr__4=380) ! Number of records
character*80 ar__4 ! Full-size record
integer*4 Hseqno_3 (nr__4) ! HRI sequence number (G1)
integer*4 H1srcno (nr__4) ! HRI source number (1)
integer*4 H1cell (nr__4) ! (arcsec) Detection cell size
real*8 H1y (nr__4) ! (pix) Source y position
real*8 H1z (nr__4) ! (pix) Source z position
real*4 H1poserr (nr__4) ! (arcsec) Position uncertainty
integer*4 H2srcno (nr__4) ! # of another source in same sequence
* detected with different cell size (1)
integer*4 H2cell (nr__4) ! (arcsec) Detection cell size
real*8 H2y (nr__4) ! (pix) Second source y position
real*8 H2z (nr__4) ! (pix) Second source z position
real*4 H2poserr (nr__4) ! (arcsec) Second source position uncertainty
real*8 separation (nr__4) ! (arcsec) Separation between the two sources (2)
real*8 nsep (nr__4) ! Normalized separation (2)
*Note (1): These columns give the HRI source numbers, within a given field,
* of two sources which were detected with different cell sizes.
*Note (2): The 'separation' column gives the pure separation between the
* two sources in arcseconds. The nsep column gives the normalized
* separation, defined to be the ratio of separation to total positional
* uncertainty. Total positional uncertainty is defined as the sum of
* individual positional uncertainties, with a systematic error of 3.5"
* added in quadrature to each individual uncertainty.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'ifield.dat' ! EDS IPC Field Parameters
integer*4 nr__5
parameter (nr__5=9) ! Number of records
character*153 ar__5 ! Full-size record
C B1950 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_3 (nr__5) ! (deg) Right Ascension B1950
real*8 DEdeg_3 (nr__5) ! (deg) Declination B1950
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 Iseqno (nr__5) ! IPC Sequence Number (G1)
real*8 ontime_1 (nr__5) ! (s) Uncorrected exposure time (1)
real*8 ltime (nr__5) ! (s) Corrected livetime (2)
real*4 sfbk (nr__5) ! (ct/pix2) Field average background for soft band (3)
real*4 hfbk (nr__5) ! (ct/pix2) Field average background for hard band (3)
real*4 bfbk (nr__5) ! (ct/pix2) Field average background for broad band (3)
integer*4 RAh_3 (nr__5) ! (h) Right Ascension (B1950) (hour)
integer*4 RAm_3 (nr__5) ! (min) Right Ascension (B1950) (min)
integer*4 RAs_3 (nr__5) ! (s) Right Ascension (B1950) (sec)
character*1 DE__3 (nr__5) ! Declination sign (B1950)
integer*4 DEd_3 (nr__5) ! (deg) Declination (B1950) (deg)
integer*4 DEm_3 (nr__5) ! (arcmin) Declination (B1950) (arcmin)
integer*4 DEs_3 (nr__5) ! (arcsec) Declination (B1950) (arcsec)
real*8 nroll_1 (nr__5) ! (deg) Nominal Roll Angle (G2)
real*4 nH (nr__5) ! (cm-2) Neutral Hydrogen Column Density (4)
real*4 sr2f (nr__5) ! (mJ/m2/ct) Flux conversion factor for soft band (5)
real*4 hr2f (nr__5) ! (mJ/m2/ct) Flux conversion factor for hard band (5)
real*4 br2f (nr__5) ! (mJ/m2/ct) Flux conversion factor for broad band (5)
real*4 alpha (nr__5) ! Backgrnd map normaliz. factor for DSMAP (6)
real*4 sbeta (nr__5) ! BG map normaliz. factor for soft band BEMAP (6)
real*4 hbeta (nr__5) ! BG map normaliz. factor for hard band BEMAP (6)
real*4 bbeta (nr__5) ! BG map norm. factor for broad band BEMAP (6)
real*4 bal (nr__5) ! Time-averaged BAL for the sequence (7)
*Note (1): The 'ontime' column gives the total uncorrected exposure time,
* in seconds.
*Note (2): The 'ltime' column gives the corrected livetime (in seconds);
* this has been corrected for intervals lost because of Earth
* occultation, passage through the South Atlantic Anomaly, unacceptable
* aspect solution, and detector dead time. It is the net time, in
* seconds, during which the detectors gathered data.
*Note (3): These columns give the field-averaged background, in
* counts/square pixel, for the soft (0.16-0.8 keV), hard (0.8-3.5 keV),
* and broad (0.16-3.5 keV) bands. The background has been obtained from
* an annulus between 8 and 15 arcminutes about the target center.
*Note (4): This column gives the column density of neutral hydrogen in
* atoms per cm^2^, integrated along the line of sight through our Galaxy
* in the direction of the field center of the observation. Specified in
* exponential notation, 'nH' values are obtained from sky surveys of
* neutral hydrogen using 21 cm radio observations. The data are from the
* Bell Lab survey of Stark et al., private communication.
*Note (5): These columns give the conversion factor from count rate
* (cts/sec) to fluxes (ergs/cm**2/sec) for the soft, hard, and broad
* bands.
*Note (6): These quantities describe the relative contributions of the two
* template background maps in the background map generated for the
* image. The DSMAP template is generated by summing several deep survey
* exposures. In constructing the image background map, the DSMAP
* template is normalized by the ratio of image livetime to DSMAP
* livetime. This normalization factor, 'alpha', is the same for all
* energy bands. The BEMAP template is generated by summing several
* exposures of the bright earth. In constructing the image background
* map, the BEMAP template is normalized so that the total counts in the
* image (with sources subtracted) and the total counts in the sum of the
* normalized DSMAP and BEMAP are the same. The BEMAP normalization
* factors, 'sbeta', 'hbeta', and 'bbeta', thus represent the ratio of
* BEMAP counts in the image to counts in the BEMAP template. They differ
* in the three energy bands. For more information on the construction of
* the image background map from the DSMAP and BEMAP templates, see
* Primini, et al. 1991ApJ...374..440P.
*Note (7): This column contains the time-averaged gain of the instrument
* for the observation.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'isrc.dat' ! EDS IPC Source Parameters
integer*4 nr__6
parameter (nr__6=178) ! Number of records
character*271 ar__6 ! Full-size record
C B1950 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_4 (nr__6) ! (deg) Right Ascension B1950
real*8 DEdeg_4 (nr__6) ! (deg) Declination B1950
C ---------------------------------- ! (position vector(s) in degrees)
character*4 Iseqno_1 (nr__6) ! IPC sequence number (G1)
integer*4 Isrcno (nr__6) ! IPC source number in sequence (1)
integer*4 reco (nr__6) ! Rib and Edge Code (2)
character*2 pcod (nr__6) ! Detection method and band (3)
integer*4 poscode (nr__6) ! Position code (1 = within central 32'x32') (4)
integer*4 poserr_1 (nr__6) ! (arcsec) Position uncertainty
integer*4 RAh_4 (nr__6) ! (h) Right Ascension (B1950) (hour)
integer*4 RAm_4 (nr__6) ! (min) Right Ascension (B1950) (min)
integer*4 RAs_4 (nr__6) ! (s) Right Ascension (B1950) (sec)
character*1 DE__4 (nr__6) ! Declination sign (B1950)
integer*4 DEd_4 (nr__6) ! (deg) Declination (B1950) (deg)
integer*4 DEm_4 (nr__6) ! (arcmin) Declination (B1950) (arcmin)
integer*4 DEs_4 (nr__6) ! (arcsec) Declination (B1950) (arcsec)
real*4 theta (nr__6) ! (arcmin) Off-axis angle (5)
real*4 vcorr_1 (nr__6) ! Rev. 1B vignetting correction (6)
integer*4 Scell (nr__6) ! (ct) Soft counts in detect cell at source posn (7)
real*4 e_Scell (nr__6) ! (ct) Error in detect cell counts in soft band (7)
integer*4 Hcell (nr__6) ! (ct) Hard counts in detect cell at source posn (7)
real*4 e_Hcell (nr__6) ! (ct) Error in detect cell counts in hard band (7)
integer*4 Bcell (nr__6) ! (ct) Broad counts in detect cell at source posn (7)
real*4 e_Bcell (nr__6) ! (ct) Error in detect cell counts in broad band (7)
real*4 Scellbk (nr__6) ! (ct) Background in detect cell in soft band (8)
integer*4 e_Scellbk (nr__6) ! (ct) Error in detect cell bg in soft band (8)
real*4 Hcellbk (nr__6) ! (ct) Background in detect cell in hard band (8)
real*4 e_Hcellbk (nr__6) ! (ct) Error in detect cell bg in hard band (8)
real*4 Bcellbk (nr__6) ! (ct) Background in detect cell in broad band (8)
integer*4 e_Bcellbk (nr__6) ! (ct) Error in detect cell bg in broad band (8)
real*4 Sccr (nr__6) ! (ct/s) Soft band corrected count rate (9)
real*4 e_Sccr (nr__6) ! (ct/s) Error in soft band corrected count rate (9)
real*4 Hccr (nr__6) ! (ct/s) Hard band corrected count rate (9)
real*4 e_Hccr (nr__6) ! (ct/s) Error in hard band corrected count rate (9)
real*4 Bccr (nr__6) ! (ct/s) Broad band corrected count rate (9)
real*4 e_Bccr (nr__6) ! (ct/s) Error in broad band corrected count rate (9)
real*4 Sflux (nr__6) ! (mW/m2) Soft band flux (10)
real*4 e_Sflux (nr__6) ! (mW/m2) Error in soft band flux (10)
real*4 Hflux (nr__6) ! (mW/m2) Hard band flux (10)
real*4 e_Hflux (nr__6) ! (mW/m2) Error in hard band flux (10)
real*4 Bflux (nr__6) ! (mW/m2) Broad band flux (10)
real*4 e_Bflux (nr__6) ! (mW/m2) Error in broad band flux (10)
real*4 Scorrsnr (nr__6) ! Soft band signal-to-noise ratio
real*4 Hcorrsnr (nr__6) ! Hard band signal-to-noise ratio
real*4 Bcorrsnr (nr__6) ! Broad band signal-to-noise ratio
real*4 f (nr__6) ! Deep survey count rate correction (11)
*Note (1): The 'srcno' column gives the IPC source number within the field.
*Note (2): This column gives the ribs and edges code ('reco'), a numeric
* flag with nonzero values indicating that the detection cell (and/or
* frame defined for "local detect") falls near or on a rib shadow or
* detector edge. The code is formed by concatenating a pair of two-digit
* numbers: the first (thousands and hundreds columns) gives the number
* (up to 16) of frame subcells potentially shadowed by the ribs or
* masked-field edge; the second (tens and units columns) gives the
* number of subcells (up to 09) within the detect box which are
* potentially shadowed.
*Note (3): The 'pcod' column gives the detection method and band which
* describes the origin of the position for each IPC source detected
* during standard Rev 1B processing. There are six combinations which
* describe the detection method and band for each IPC source detected;
* these combinations are given by a two-digit code, with the first digit
* representing the detection method (L for L detect, or M for M detect),
* and the second digit representing the spectral band used in detection
* (S for soft, H for hard, or B for broad). 'pcod' describes from which
* method/band combination the position is given; it is the first
* available one from the following prioritized order: MH, MB, MS, LH,
* LB, LS.
*Note (4): The 'poscode' column gives the position code for each source.
* The position code will have bit 0 (the least significant bit) set if
* the source is within the central 32'x32' of the field. It will have
* bit 1 set if 'reco' > 0 (see above for description of 'reco').
*Note (5): This is the distance of the source from the pointing position of
* the observation.
*Note (6): This column gives the Rev 1B vignetting correction and is used
* in determining the corrected count rates.
*Note (7): The 'Scell' and 'e_Scell' columns give the number of source
* counts in the detect cell and the corresponding error at the source
* position in the soft band (0.16-0.8 keV); 'Hcell' and 'e_Hcell' give
* these quantities for the hard band (0.8-3.5 keV); 'Bcell' and
* 'e_Bcell' give these quantities for the broad band (0.16-3.5 keV).
*Note (8): The 'Scellbk' and 'e_Scellbk' columns give the number of
* background counts and the corresponding error in the detect cell in
* the soft band. 'Hcellbk' and 'e_Hcellbk' give these quantities for the
* hard band; 'Bcellbk' and 'e_Bcellbk' give these quantities for the
* broad band. The background counts are determined either from counts in
* a region surrounding the detect cell (LDETECT) or from counts in the
* corresponding detect cell of the background map (MDETECT). In the
* latter case, the error in the background counts is assumed to be zero.
*Note (9): The 'Sccr' and 'e_Sccr' columns give the soft band corrected
* count rate (in counts/sec) and the corresponding error. 'Hccr' and
* 'e_Hccr' give these quantities for the hard band; 'Bccr' and 'e_Bccr'
* give these quantities for the broad band. In general, the corrected
* count rate is determined as follows:
*
* ccr = (cell-cellbk) * vcorr * CMS * CPRF / livetime
*
* vcorr is described above. CPRF corrects for source photons which fall
* outside the detect cell. Detect cell sizes are chosen so that
* CPRF=1.13 for all energy bands. CMS corrects for large angle
* scattering intrinsic in the HEAO-2 mirror. This factor is
* energy-dependent and has a value of 1.09 for the soft band and 1.18
* for the hard and broad bands. For deep survey sources, additional
* minor count rate corrections may be necessary. These are described in
* the "f" column below.
*Note (10): The 'Sflux' and 'e_Sflux' columns give the soft band flux (in
* ergs/cm**2/s or mW/m2) and the corresponding error. 'Hflux' and
* 'e_Hflux' give these quantities for the hard band; 'Bflux' and
* 'e_Bflux' give these quantities for the broad band.
*Note (11): This column gives the deep survey count rate correction factor;
* this number will be 0 if the source was excluded from the DSMAP.
* Otherwise, it will be equal to the ratio of the livetime of the
* sequence number to the livetime of the DSMAP. This correction is
* necessary because the DSMAP templates used to construct background
* maps are generated from the deep survey observations themselves. Thus
* calculation of count rates for deep survey sources may be biased by
* inclusion of source counts in the background map. Correction for this
* bias is described in detail in Primini et al. 1991.
* (=1991ApJ...374..440P)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'ispec.dat' ! EDS Spectral Parameters for IPC Sources
integer*4 nr__7
parameter (nr__7=178) ! Number of records
character*375 ar__7 ! Full-size record
character*4 Iseqno_2 (nr__7) ! IPC sequence number (G1)
integer*4 Isrcno_1 (nr__7) ! IPC source no. within sequence (1)
real*4 Barea (nr__7) ! (pix2) pixels Area of background region (2)
real*4 Sarea (nr__7) ! (pix2) pixels Area of source region (2)
real*8 hr (nr__7) ! Hardness ratio (3)
real*8 E_hr (nr__7) ! Positive hardness ratio error (3)
real*8 e_hr_1 (nr__7) ! Negative hardness ratio error (real) (3)
real*4 phbk (16,nr__7) ! Pulse height spectrum of background region (4)
real*4 pibk (16,nr__7) ! Pi spectrum of background region (5)
real*4 ph (16,nr__7) ! Pulse height spectrum of source (6)
real*4 pi (16,nr__7) ! Pi spectrum of source (7)
*Note (1): This column gives the IPC source number within the field, for
* each IPC source detected in the Einstein Deep Survey.
*Note (2): These columns give the area, in square pixels, for the
* background region and the source region used to accumulate counts for
* the IP/RHA spectra. The background region is nominally a 5'-6'annulus
* around the source, with any other sources excluded. The source region
* is nominally a 3' circle around the source, with any other sources
* excluded. An IPC pixel is 8"x8".
*Note (3): These columns give the hardness ratio and its plus and minus
* error, from Rev1B processing. The hardness ratio for a source is
* defined as the ratio of the difference between the number of counts in
* the hard and soft bands to the sum of these quantities:
* counts(hard) - counts(soft)
* hr = ---------------------------
* counts(hard) + counts(soft)
* Note that the error values are assumed symmetric and are thus only
* approximate.
*Note (4): These columns give the pulse height spectrum of the background
* region for each of the 16 pha channels.
*Note (5): These columns give the pi spectrum of the background region for
* each of the 16 pi channels.
*Note (6): These columns give the pulse height spectrum of the source for
* each of the 16 pha channels.
*Note (7): These columns give the pi spectrum of the source for each of the
* 16 pi channels.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'icomp.dat' ! EDS IPC Source Component Parameters
integer*4 nr__8
parameter (nr__8=453) ! Number of records
character*98 ar__8 ! Full-size record
integer*4 Iseqno_3 (nr__8) ! IPC sequence number (G1)
integer*4 Isrcno_2 (nr__8) ! IPC source number in sequence (1)
character*2 comp (nr__8) ! [LMBHS] Detection method and band (2)
real*8 ypix (nr__8) ! (pix) Source y position
real*4 e_ypix (nr__8) ! (pix) Error in y position
real*8 zpix (nr__8) ! (pix) Source z position
real*4 e_zpix (nr__8) ! (pix) Error in z position
integer*4 cellcts (nr__8) ! (ct) Counts in detect cell centered at y,z (3)
real*4 e_cellcts (nr__8) ! (ct) Error in detect cell counts (3)
real*8 cellbk (nr__8) ! (ct) Background in detect cell (3)
real*4 e_cellbk (nr__8) ! (ct) Error in background (3)
real*8 netccts (nr__8) ! (ct) Net source counts in detect cell
real*4 e_netccts (nr__8) ! (ct) Error in net counts
real*8 thr (nr__8) ! (ct) Threshold counts in detect cell for source
* detection.
real*4 snr (nr__8) ! Count rate signal-to-noise ratio
*Note (1): The 'Isrcno' column gives the IPC source number within the field.
*Note (2): This column gives a two-digit code which identifies the
* detection method and band for each IPC source. The two detection
* method codes are L and M, for LDETECT and MDETECT. The band will be
* identified by an S, H, or B (for soft, hard, or broad). Thus, a
* typical entry in this column might be MH for an MDETECT hard band
* detection.
*Note (3): The 'cellcts' and 'e_cellcts' columns give the counts and
* corresponding error in the detect cell centered at the given source
* position, y,z. Likewise, the 'cellbk' and 'e_cellbk' give the
* background counts and corresponding error in the source detect cell.
* The background counts are determined either from counts in a region
* surrounding the detect cell (LDETECT) or from counts in the
* corresponding detect cell of the background map (MDETECT). In the
* latter case, the error in the background counts is assumed to be zero.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'ipchri.dat' ! EDS Separations Between IPC and HRI Sources
integer*4 nr__9
parameter (nr__9=5861) ! Number of records
character*89 ar__9 ! Full-size record
C B1950 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_5 (nr__9) ! (deg) Right Ascension B1950
real*8 DEdeg_5 (nr__9) ! (deg) Declination B1950
C ---------------------------------- ! (position vector(s) in degrees)
character*4 Iseqno_4 (nr__9) ! IPC sequence number (G1)
integer*4 Isrcno_3 (nr__9) ! IPC source number (1)
integer*4 RAh_5 (nr__9) ! (h) Right Ascension (B1950) (hour)
integer*4 RAm_5 (nr__9) ! (min) Right Ascension (B1950) (min)
real*4 RAs_5 (nr__9) ! (s) Right Ascension (B1950) (sec)
character*1 DE__5 (nr__9) ! Declination sign (B1950)
integer*4 DEd_5 (nr__9) ! (deg) Declination (B1950) (deg)
integer*4 DEm_5 (nr__9) ! (arcmin) Declination (B1950) (arcmin)
real*4 DEs_5 (nr__9) ! (arcsec) Declination (B1950) (arcsec)
integer*4 ipcperr (nr__9) ! (arcsec) IPC source position uncertainty (2)
real*4 ipchsnr (nr__9) ! IPC Hard band signal-to-noise ratio
integer*4 Hseqno_4 (nr__9) ! ?=0 HRI sequence number (G1)
integer*4 Hsrcno_2 (nr__9) ! ?=0 HRI source number
integer*4 RA_HRI_h (nr__9) ! (h) HRI Right Ascension (B1950) (hour)
integer*4 RA_HRI_m (nr__9) ! (min) HRI Right Ascension (B1950) (min)
real*4 RA_HRI_s (nr__9) ! (s) HRI Right Ascension (B1950) (sec)
character*1 DE_HRI__ (nr__9) ! HRI Declination sign (B1950)
integer*4 DE_HRI_d (nr__9) ! (deg) HRI Declination (B1950) (deg)
integer*4 DE_HRI_m (nr__9) ! (arcmin) HRI Declination (B1950) (arcmin)
integer*4 DE_HRI_s (nr__9) ! (arcsec) HRI Declination (B1950) (arcsec)
real*4 hriperr (nr__9) ! (arcsec) HRI source position uncertainty (3)
real*4 theta_1 (nr__9) ! (arcsec) Separation between IPC and HRI sources (4)
real*4 nsep_1 (nr__9) ! Normalized separation (4)
*Note (1): These columns give the IPC sequence number and source number
* within that field for Einstein Deep Survey IPC sources.
*Note (2): This column gives the positional uncertainty, in arcseconds, of
* the given IPC source.
*Note (3): This column gives the positional uncertainty, in arcseconds, of
* the given HRI source.
*Note (4): 'theta' gives the separation between IPC and HRI sources
* detected in corresponding fields for each target in the Einstein Deep
* Survey. The 'nsep' column gives the normalized separation between IPC
* and HRI sources, defined to be the ratio of separation to sum of IPC
* and HRI positional uncertainties. A systematic uncertainty of 3.5" is
* added in quadrature to the HRI uncertainty before calculating nsep.
* The IPC position error also includes a systematic uncertainty, which
* depends on source position, and ranges from 25" to 45".
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'ihsample.dat' ! EDS Cross-Reference of HRI and IPC Sources
integer*4 nr__10
parameter (nr__10=45) ! Number of records
character*15 ar__10 ! Full-size record
integer*4 Iseqno_5 (nr__10) ! IPC sequence number (1)
integer*4 Isrcno_4 (nr__10) ! IPC source number (1)
integer*4 Hseqno_5 (nr__10) ! ?=0 HRI sequence number (2)
integer*4 Hsrcno_3 (nr__10) ! ?=0 HRI source number of matching HRI source (2)
*Note (1): These columns give the IPC sequence number and source number
* within that field for Einstein Deep Survey IPC sources greater than
* 4.5 sigma which also have a coincident HRI source detection.
*Note (2): These columns give the HRI sequence number and source number
* within that field for Einstein Deep Survey HRI sources which also have
* a coincident IPC source detection greater than 4.5 sigma. Coincidences
* are defined as an IPC-HRI pair with a normalized separation of two or
* less.
* The normalized separation is described as:
*
* IPC-HRI separation
* ------------------------------------
* sum of IPC and HRI positional errors
*
* This threshold was determined empirically by examining the
* distribution of normalized separations for all IPC-HRI pairs.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'sample.dat' ! EDS IPC, HRI + Optical Data for 4.5Sigma IPC Sample
integer*4 nr__11
parameter (nr__11=45) ! Number of records
character*136 ar__11 ! Full-size record
C 1950 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_6 (nr__11) ! (deg) Right Ascension 1950
real*8 DEdeg_6 (nr__11) ! (deg) Declination 1950
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 Iseqno_6 (nr__11) ! IPC sequence number (G1)
integer*4 Isrcno_5 (nr__11) ! IPC source number (1)
integer*4 RAh_6 (nr__11) ! (h) Right Ascension (1950) (hour)
integer*4 RAm_6 (nr__11) ! (min) Right Ascension (1950) (min)
integer*4 RAs_6 (nr__11) ! (s) Right Ascension (1950) (sec)
character*1 DE__6 (nr__11) ! Declination Sign (1950)
integer*4 DEd_6 (nr__11) ! (deg) Declination (1950) (deg)
integer*4 DEm_6 (nr__11) ! (arcmin) Declination (1950) (arcmin)
integer*4 DEs_6 (nr__11) ! (arcsec) Declination (1950) (arcsec)
integer*4 Iposerr (nr__11) ! (arcsec) IPC position uncertainty
real*8 Iflux (nr__11) ! (10-17W/m2) IPC source flux (*1.0E14) (2)
real*4 e_Iflux (nr__11) ! (10-17W/m2) IPC flux error (*1.0E14)
real*4 Isnr (nr__11) ! IPC count rate signal-to-noise ratio
integer*4 Hseqno_6 (nr__11) ! ?=0 HRI sequence number (3)
integer*4 Hsrcno_4 (nr__11) ! ?=0 HRI src # of HRI src (if any) that
* matches IPC src (3)
integer*4 Hcell_1 (nr__11) ! (arcsec) HRI detection cell size
character*27 Hpos (nr__11) ! HRI B1950 source position (sexag.)
real*4 e_Hpos (nr__11) ! (arcsec) HRI position uncertainty
real*4 Hrate (nr__11) ! (ct/ks) HRI source count rate per 1000 s (4)
real*4 e_Hrate (nr__11) ! (ct/ks) HRI rate error
character*1 l_Vmag (nr__11) ! [> ] uncertainty in Vmag
real*4 Vmag (nr__11) ! (mag) ? Optical visual magnitude
character*1 l_fxfv (nr__11) ! [> ] uncertainty in fxfv
real*4 fxfv (nr__11) ! ? Log of ratio of broad x-ray flux/optical
* flux (5)
real*4 z_1 (nr__11) ! ? OPT redshift
character*6 ID (nr__11) ! Identification (6)
*Note (1): This column gives the IPC source number within that field for
* Einstein Deep Survey IPC sources contained in the 4.5 sigma signal-to-
* noise-limited sample.
*Note (2): This column give the hard band IPC source flux in units of
* 10^-14^ ergs/cm^2^/s.
*Note (3): These columns give the HRI sequence number and source number
* within that field for Einstein Deep Survey HRI sources contained in
* the 4.5 sigma signal-to-noise-limited sample. Note that there may be
* more than one HRI detection for each IPC source.
*Note (4): These columns give the HRI count rate per 1000 sec and the
* corresponding error for each HRI source detected in the 4.5 sigma
* sample. The count rates have been corrected for vignetting,
* scattering, dead time, and detector efficiency.
*Note (5): This column gives the logarithm of the ratio of x-ray flux in
* the IPC broad band (0.16-3.5 keV) to the optical flux.
*Note (6): This column gives the optical identification for the IPC source.
* Current identifications include: Galaxy, Star, and QSO.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'optinfo.dat' ! EDS Optical Data for Srcs in the 4.5 Sigma Sample
integer*4 nr__12
parameter (nr__12=25) ! Number of records
character*45 ar__12 ! Full-size record
integer*4 Iseqno_7 (nr__12) ! IPC sequence number (G1)
integer*4 Isrcno_6 (nr__12) ! IPC source number (1)
character*1 l_Vmag_1 (nr__12) ! [> ] limit to the visual magnitude
real*4 Vmag_1 (nr__12) ! (mag) ? Visual magnitude
character*1 l_fxfv_1 (nr__12) ! [> ] '>' when fxfv is a lower limit
real*4 fxfv_1 (nr__12) ! ? Log of ratio of broad x-ray flux/optical flux (2)
real*4 z_2 (nr__12) ! ? Redshift
character*6 ID_1 (nr__12) ! Galaxy Identification (3)
*Note (1): The source number within that field for Einstein Deep Survey IPC
* sources for which there is optical information.
*Note (2): This column gives the logarithm of the ratio of x-ray flux in
* the broad band (0.16-3.5 keV) to the optical flux.
*Note (3): This column gives the optical identification for the IPC source.
* Current identifications include: Galaxy, Star, and QSO.
C=============================================================================
C Loading file 'targets.dat' ! EDS Target Information
C Format for file interpretation
1 format(
+ 1X,I5,1X,A3,1X,I2,1X,I2,1X,I2,1X,A1,I2,1X,I2,1X,I2,1X,F10.3)
C Effective file loading
open(unit=1,status='old',file=
+'targets.dat')
write(6,*) '....Loading file: targets.dat'
do i__=1,43
read(1,'(A40)')ar__
read(ar__,1)
+ Seqno(i__),detector(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
+ DEd(i__),DEm(i__),DEs(i__),exptime(i__)
RAdeg(i__) = rNULL__
DEdeg(i__) = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c ..............Just test output...........
write(6,1)
+ Seqno(i__),detector(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
+ DEd(i__),DEm(i__),DEs(i__),exptime(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'hfield.dat' ! EDS Field Parameters for HRI Observations
C Format for file interpretation
2 format(
+ 1X,I5,1X,F10.3,1X,F6.3,1X,F6.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,
+ F8.3,1X,F8.3,1X,F8.3,1X,I2,1X,I2,1X,I2,1X,A1,I2,1X,I2,1X,I2,
+ 1X,F8.3)
C Effective file loading
open(unit=1,status='old',file=
+'hfield.dat')
write(6,*) '....Loading file: hfield.dat'
do i__=1,34
read(1,'(A110)')ar__1
read(ar__1,2)
+ Hseqno(i__),ontime(i__),globk(i__),detbk(i__),th12(i__),
+ th24(i__),th36(i__),th48(i__),th72(i__),th120(i__),RAh_1(i__),
+ RAm_1(i__),RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),
+ DEs_1(i__),nroll(i__)
RAdeg_1(i__) = rNULL__
DEdeg_1(i__) = rNULL__
c Derive coordinates RAdeg_1 and DEdeg_1 from input data
c (RAdeg_1 and DEdeg_1 are set to rNULL__ when unknown)
if(RAh_1(i__) .GT. -180) RAdeg_1(i__)=RAh_1(i__)*15.
if(RAm_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAm_1(i__)/4.
if(RAs_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAs_1(i__)/240.
if(DEd_1(i__) .GE. 0) DEdeg_1(i__)=DEd_1(i__)
if(DEm_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEm_1(i__)/60.
if(DEs_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEs_1(i__)/3600.
if(DE__1(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__)
c ..............Just test output...........
write(6,2)
+ Hseqno(i__),ontime(i__),globk(i__),detbk(i__),th12(i__),
+ th24(i__),th36(i__),th48(i__),th72(i__),th120(i__),RAh_1(i__),
+ RAm_1(i__),RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),
+ DEs_1(i__),nroll(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg_1(i__),DEdeg_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'hfluxes.dat' ! EDS Secondary HRI Source Parameters
C Format for file interpretation
3 format(
+ 1X,I5,1X,I2,1X,I2,1X,I3,1X,I3,1X,F5.3,10I3,10F5.1,10F6.3,1X,
+ F6.3,1X,E10.5,1X,E10.5,1X,E10.5,1X,F5.3,1X,F5.2,1X,E10.5,1X,
+ E10.5,1X,E10.5,1X,F2.0,1X,F5.2,1X,E10.5,1X,E10.5,1X,E10.5,1X,
+ F2.0,1X,F5.2,1X,E10.5,1X,E10.5,1X,F7.3)
C Effective file loading
open(unit=1,status='old',file=
+'hfluxes.dat')
write(6,*) '....Loading file: hfluxes.dat'
do i__=1,210
read(1,'(A332)')ar__2
read(ar__2,3)
+ Hseqno_1(i__),Hsrcno(i__),cts0to3(i__),cts3to9(i__),
+ cts9to18(i__),USProb(i__),(ctR(j__,i__),j__=1,10),(radR(j__,
+ i__),j__=1,10),(denR(j__,i__),j__=1,10),CR_S_PS(i__),
+ FL_S_PS(i__),CR_H_PS(i__),FL_H_PS(i__),xtnt(i__),CR_S_XS(i__),
+ FL_S_XS(i__),CR_H_XS(i__),FL_H_XS(i__),ulxtnt(i__),
+ CR_S_UXS(i__),FL_S_UXS(i__),CR_H_UXS(i__),FL_H_UXS(i__),
+ llxtnt(i__),CR_S_LXS(i__),FL_S_LXS(i__),CR_H_LXS(i__),
+ FL_H_LXS(i__)
c ..............Just test output...........
write(6,3)
+ Hseqno_1(i__),Hsrcno(i__),cts0to3(i__),cts3to9(i__),
+ cts9to18(i__),USProb(i__),(ctR(j__,i__),j__=1,10),(radR(j__,
+ i__),j__=1,10),(denR(j__,i__),j__=1,10),CR_S_PS(i__),
+ FL_S_PS(i__),CR_H_PS(i__),FL_H_PS(i__),xtnt(i__),CR_S_XS(i__),
+ FL_S_XS(i__),CR_H_XS(i__),FL_H_XS(i__),ulxtnt(i__),
+ CR_S_UXS(i__),FL_S_UXS(i__),CR_H_UXS(i__),FL_H_UXS(i__),
+ llxtnt(i__),CR_S_LXS(i__),FL_S_LXS(i__),CR_H_LXS(i__),
+ FL_H_LXS(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'hsrc.dat' ! EDS HRI Source Positions and Count Rates
C Format for file interpretation
4 format(
+ 1X,I5,1X,I2,1X,I3,1X,F8.3,1X,F8.3,1X,F4.1,1X,I2,1X,I2,1X,I2,
+ 1X,A1,I2,1X,I2,1X,I2,1X,F6.3,1X,F5.3,1X,F5.0,1X,F5.0,1X,F4.0,
+ 1X,F7.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F9.2,1X,I4,1X,F7.2,
+ 1X,F10.7,1X,F11.7,1X,F10.7)
C Effective file loading
open(unit=1,status='old',file=
+'hsrc.dat')
write(6,*) '....Loading file: hsrc.dat'
do i__=1,210
read(1,'(A174)')ar__3
read(ar__3,4)
+ Hseqno_2(i__),Hsrcno_1(i__),csize(i__),y(i__),z(i__),
+ poserr(i__),RAh_2(i__),RAm_2(i__),RAs_2(i__),DE__2(i__),
+ DEd_2(i__),DEm_2(i__),DEs_2(i__),netcr(i__),e_netcr(i__),
+ bpcts(i__),eqcts(i__),bp18cts(i__),bp18bg(i__),bp18scat(i__),
+ vcorr(i__),dtfac(i__),qec(i__),exptime_1(i__),xccts(i__),
+ xcbg(i__),xcscat(i__),xcrad(i__),sizcor(i__)
RAdeg_2(i__) = rNULL__
DEdeg_2(i__) = rNULL__
c Derive coordinates RAdeg_2 and DEdeg_2 from input data
c (RAdeg_2 and DEdeg_2 are set to rNULL__ when unknown)
if(RAh_2(i__) .GT. -180) RAdeg_2(i__)=RAh_2(i__)*15.
if(RAm_2(i__) .GT. -180) RAdeg_2(i__)=RAdeg_2(i__)+RAm_2(i__)/4.
if(RAs_2(i__) .GT. -180) RAdeg_2(i__)=RAdeg_2(i__)+RAs_2(i__)/240.
if(DEd_2(i__) .GE. 0) DEdeg_2(i__)=DEd_2(i__)
if(DEm_2(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEm_2(i__)/60.
if(DEs_2(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEs_2(i__)/3600.
if(DE__2(i__).EQ.'-'.AND.DEdeg_2(i__).GE.0) DEdeg_2(i__)=-DEdeg_2(i__)
c ..............Just test output...........
write(6,4)
+ Hseqno_2(i__),Hsrcno_1(i__),csize(i__),y(i__),z(i__),
+ poserr(i__),RAh_2(i__),RAm_2(i__),RAs_2(i__),DE__2(i__),
+ DEd_2(i__),DEm_2(i__),DEs_2(i__),netcr(i__),e_netcr(i__),
+ bpcts(i__),eqcts(i__),bp18cts(i__),bp18bg(i__),bp18scat(i__),
+ vcorr(i__),dtfac(i__),qec(i__),exptime_1(i__),xccts(i__),
+ xcbg(i__),xcscat(i__),xcrad(i__),sizcor(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg_2(i__),DEdeg_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'hris_sep.dat' ! EDS Separations Between HRI Srcs in a Single Field
C Format for file interpretation
5 format(
+ I5,1X,I1,1X,I2,1X,F8.3,1X,F8.3,1X,F3.1,1X,I2,1X,I3,1X,F8.3,1X,
+ F8.3,1X,F4.1,1X,F8.3,1X,F7.3)
C Effective file loading
open(unit=1,status='old',file=
+'hris_sep.dat')
write(6,*) '....Loading file: hris_sep.dat'
do i__=1,380
read(1,'(A80)')ar__4
read(ar__4,5)
+ Hseqno_3(i__),H1srcno(i__),H1cell(i__),H1y(i__),H1z(i__),
+ H1poserr(i__),H2srcno(i__),H2cell(i__),H2y(i__),H2z(i__),
+ H2poserr(i__),separation(i__),nsep(i__)
c ..............Just test output...........
write(6,5)
+ Hseqno_3(i__),H1srcno(i__),H1cell(i__),H1y(i__),H1z(i__),
+ H1poserr(i__),H2srcno(i__),H2cell(i__),H2y(i__),H2z(i__),
+ H2poserr(i__),separation(i__),nsep(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'ifield.dat' ! EDS IPC Field Parameters
C Format for file interpretation
6 format(
+ I5,1X,F8.2,1X,F8.2,1X,F5.3,1X,F5.3,1X,F5.3,1X,I2,1X,I2,1X,I2,
+ 1X,A1,I2,1X,I2,1X,I2,1X,F7.2,1X,E8.4,1X,E9.4,1X,E9.4,1X,E9.4,
+ 1X,F5.3,1X,E10.4,1X,E10.4,1X,E10.4,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'ifield.dat')
write(6,*) '....Loading file: ifield.dat'
do i__=1,9
read(1,'(A153)')ar__5
read(ar__5,6)
+ Iseqno(i__),ontime_1(i__),ltime(i__),sfbk(i__),hfbk(i__),
+ bfbk(i__),RAh_3(i__),RAm_3(i__),RAs_3(i__),DE__3(i__),
+ DEd_3(i__),DEm_3(i__),DEs_3(i__),nroll_1(i__),nH(i__),
+ sr2f(i__),hr2f(i__),br2f(i__),alpha(i__),sbeta(i__),
+ hbeta(i__),bbeta(i__),bal(i__)
RAdeg_3(i__) = rNULL__
DEdeg_3(i__) = rNULL__
c Derive coordinates RAdeg_3 and DEdeg_3 from input data
c (RAdeg_3 and DEdeg_3 are set to rNULL__ when unknown)
if(RAh_3(i__) .GT. -180) RAdeg_3(i__)=RAh_3(i__)*15.
if(RAm_3(i__) .GT. -180) RAdeg_3(i__)=RAdeg_3(i__)+RAm_3(i__)/4.
if(RAs_3(i__) .GT. -180) RAdeg_3(i__)=RAdeg_3(i__)+RAs_3(i__)/240.
if(DEd_3(i__) .GE. 0) DEdeg_3(i__)=DEd_3(i__)
if(DEm_3(i__) .GE. 0) DEdeg_3(i__)=DEdeg_3(i__)+DEm_3(i__)/60.
if(DEs_3(i__) .GE. 0) DEdeg_3(i__)=DEdeg_3(i__)+DEs_3(i__)/3600.
if(DE__3(i__).EQ.'-'.AND.DEdeg_3(i__).GE.0) DEdeg_3(i__)=-DEdeg_3(i__)
c ..............Just test output...........
write(6,6)
+ Iseqno(i__),ontime_1(i__),ltime(i__),sfbk(i__),hfbk(i__),
+ bfbk(i__),RAh_3(i__),RAm_3(i__),RAs_3(i__),DE__3(i__),
+ DEd_3(i__),DEm_3(i__),DEs_3(i__),nroll_1(i__),nH(i__),
+ sr2f(i__),hr2f(i__),br2f(i__),alpha(i__),sbeta(i__),
+ hbeta(i__),bbeta(i__),bal(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg_3(i__),DEdeg_3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'isrc.dat' ! EDS IPC Source Parameters
C Format for file interpretation
7 format(
+ 1X,A4,1X,I2,1X,I4,1X,A2,1X,I1,1X,I2,1X,I2,1X,I2,1X,I2,1X,A1,
+ I2,1X,I2,1X,I2,1X,F6.3,1X,F5.3,1X,I3,1X,F4.1,1X,I3,1X,F4.1,1X,
+ I4,1X,F4.1,1X,F5.1,1X,I2,1X,F6.2,1X,F4.2,1X,F6.2,1X,I2,1X,
+ E11.5,1X,E10.5,1X,E11.5,1X,E10.5,1X,E10.5,1X,E10.5,1X,E11.5,
+ 1X,E10.5,1X,E11.5,1X,E10.5,1X,E10.5,1X,E10.5,1X,F5.2,1X,F5.2,
+ 1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'isrc.dat')
write(6,*) '....Loading file: isrc.dat'
do i__=1,178
read(1,'(A271)')ar__6
read(ar__6,7)
+ Iseqno_1(i__),Isrcno(i__),reco(i__),pcod(i__),poscode(i__),
+ poserr_1(i__),RAh_4(i__),RAm_4(i__),RAs_4(i__),DE__4(i__),
+ DEd_4(i__),DEm_4(i__),DEs_4(i__),theta(i__),vcorr_1(i__),
+ Scell(i__),e_Scell(i__),Hcell(i__),e_Hcell(i__),Bcell(i__),
+ e_Bcell(i__),Scellbk(i__),e_Scellbk(i__),Hcellbk(i__),
+ e_Hcellbk(i__),Bcellbk(i__),e_Bcellbk(i__),Sccr(i__),
+ e_Sccr(i__),Hccr(i__),e_Hccr(i__),Bccr(i__),e_Bccr(i__),
+ Sflux(i__),e_Sflux(i__),Hflux(i__),e_Hflux(i__),Bflux(i__),
+ e_Bflux(i__),Scorrsnr(i__),Hcorrsnr(i__),Bcorrsnr(i__),f(i__)
RAdeg_4(i__) = rNULL__
DEdeg_4(i__) = rNULL__
c Derive coordinates RAdeg_4 and DEdeg_4 from input data
c (RAdeg_4 and DEdeg_4 are set to rNULL__ when unknown)
if(RAh_4(i__) .GT. -180) RAdeg_4(i__)=RAh_4(i__)*15.
if(RAm_4(i__) .GT. -180) RAdeg_4(i__)=RAdeg_4(i__)+RAm_4(i__)/4.
if(RAs_4(i__) .GT. -180) RAdeg_4(i__)=RAdeg_4(i__)+RAs_4(i__)/240.
if(DEd_4(i__) .GE. 0) DEdeg_4(i__)=DEd_4(i__)
if(DEm_4(i__) .GE. 0) DEdeg_4(i__)=DEdeg_4(i__)+DEm_4(i__)/60.
if(DEs_4(i__) .GE. 0) DEdeg_4(i__)=DEdeg_4(i__)+DEs_4(i__)/3600.
if(DE__4(i__).EQ.'-'.AND.DEdeg_4(i__).GE.0) DEdeg_4(i__)=-DEdeg_4(i__)
c ..............Just test output...........
write(6,7)
+ Iseqno_1(i__),Isrcno(i__),reco(i__),pcod(i__),poscode(i__),
+ poserr_1(i__),RAh_4(i__),RAm_4(i__),RAs_4(i__),DE__4(i__),
+ DEd_4(i__),DEm_4(i__),DEs_4(i__),theta(i__),vcorr_1(i__),
+ Scell(i__),e_Scell(i__),Hcell(i__),e_Hcell(i__),Bcell(i__),
+ e_Bcell(i__),Scellbk(i__),e_Scellbk(i__),Hcellbk(i__),
+ e_Hcellbk(i__),Bcellbk(i__),e_Bcellbk(i__),Sccr(i__),
+ e_Sccr(i__),Hccr(i__),e_Hccr(i__),Bccr(i__),e_Bccr(i__),
+ Sflux(i__),e_Sflux(i__),Hflux(i__),e_Hflux(i__),Bflux(i__),
+ e_Bflux(i__),Scorrsnr(i__),Hcorrsnr(i__),Bcorrsnr(i__),f(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg_4(i__),DEdeg_4(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'ispec.dat' ! EDS Spectral Parameters for IPC Sources
C Format for file interpretation
8 format(
+ 1X,A4,1X,I2,1X,F5.0,1X,F5.0,1X,F10.7,1X,F11.7,1X,F11.7,16F5.0,
+ 16F5.0,16F5.0,16F5.0)
C Effective file loading
open(unit=1,status='old',file=
+'ispec.dat')
write(6,*) '....Loading file: ispec.dat'
do i__=1,178
read(1,'(A375)')ar__7
read(ar__7,8)
+ Iseqno_2(i__),Isrcno_1(i__),Barea(i__),Sarea(i__),hr(i__),
+ E_hr(i__),e_hr_1(i__),(phbk(j__,i__),j__=1,16),(pibk(j__,i__),
+ j__=1,16),(ph(j__,i__),j__=1,16),(pi(j__,i__),j__=1,16)
c ..............Just test output...........
write(6,8)
+ Iseqno_2(i__),Isrcno_1(i__),Barea(i__),Sarea(i__),hr(i__),
+ E_hr(i__),e_hr_1(i__),(phbk(j__,i__),j__=1,16),(pibk(j__,i__),
+ j__=1,16),(ph(j__,i__),j__=1,16),(pi(j__,i__),j__=1,16)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'icomp.dat' ! EDS IPC Source Component Parameters
C Format for file interpretation
9 format(
+ I5,1X,I2,1X,A2,1X,F7.3,1X,F6.3,1X,F7.3,1X,F6.3,1X,I4,1X,F6.3,
+ 1X,F7.3,1X,F6.3,1X,F8.3,1X,F6.3,1X,F7.3,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'icomp.dat')
write(6,*) '....Loading file: icomp.dat'
do i__=1,453
read(1,'(A98)')ar__8
read(ar__8,9)
+ Iseqno_3(i__),Isrcno_2(i__),comp(i__),ypix(i__),e_ypix(i__),
+ zpix(i__),e_zpix(i__),cellcts(i__),e_cellcts(i__),cellbk(i__),
+ e_cellbk(i__),netccts(i__),e_netccts(i__),thr(i__),snr(i__)
c ..............Just test output...........
write(6,9)
+ Iseqno_3(i__),Isrcno_2(i__),comp(i__),ypix(i__),e_ypix(i__),
+ zpix(i__),e_zpix(i__),cellcts(i__),e_cellcts(i__),cellbk(i__),
+ e_cellbk(i__),netccts(i__),e_netccts(i__),thr(i__),snr(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'ipchri.dat' ! EDS Separations Between IPC and HRI Sources
C Format for file interpretation
10 format(
+ 1X,A4,1X,I2,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,I2,
+ 1X,F5.2,1X,I5,1X,I2,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,
+ 1X,F4.1,1X,F6.1,1X,F5.1)
C Effective file loading
open(unit=1,status='old',file=
+'ipchri.dat')
write(6,*) '....Loading file: ipchri.dat'
do i__=1,5861
read(1,'(A89)')ar__9
read(ar__9,10)
+ Iseqno_4(i__),Isrcno_3(i__),RAh_5(i__),RAm_5(i__),RAs_5(i__),
+ DE__5(i__),DEd_5(i__),DEm_5(i__),DEs_5(i__),ipcperr(i__),
+ ipchsnr(i__),Hseqno_4(i__),Hsrcno_2(i__),RA_HRI_h(i__),
+ RA_HRI_m(i__),RA_HRI_s(i__),DE_HRI__(i__),DE_HRI_d(i__),
+ DE_HRI_m(i__),DE_HRI_s(i__),hriperr(i__),theta_1(i__),
+ nsep_1(i__)
RAdeg_5(i__) = rNULL__
DEdeg_5(i__) = rNULL__
c Derive coordinates RAdeg_5 and DEdeg_5 from input data
c (RAdeg_5 and DEdeg_5 are set to rNULL__ when unknown)
if(RAh_5(i__) .GT. -180) RAdeg_5(i__)=RAh_5(i__)*15.
if(RAm_5(i__) .GT. -180) RAdeg_5(i__)=RAdeg_5(i__)+RAm_5(i__)/4.
if(RAs_5(i__) .GT. -180) RAdeg_5(i__)=RAdeg_5(i__)+RAs_5(i__)/240.
if(DEd_5(i__) .GE. 0) DEdeg_5(i__)=DEd_5(i__)
if(DEm_5(i__) .GE. 0) DEdeg_5(i__)=DEdeg_5(i__)+DEm_5(i__)/60.
if(DEs_5(i__) .GE. 0) DEdeg_5(i__)=DEdeg_5(i__)+DEs_5(i__)/3600.
if(DE__5(i__).EQ.'-'.AND.DEdeg_5(i__).GE.0) DEdeg_5(i__)=-DEdeg_5(i__)
c ..............Just test output...........
write(6,10)
+ Iseqno_4(i__),Isrcno_3(i__),RAh_5(i__),RAm_5(i__),RAs_5(i__),
+ DE__5(i__),DEd_5(i__),DEm_5(i__),DEs_5(i__),ipcperr(i__),
+ ipchsnr(i__),Hseqno_4(i__),Hsrcno_2(i__),RA_HRI_h(i__),
+ RA_HRI_m(i__),RA_HRI_s(i__),DE_HRI__(i__),DE_HRI_d(i__),
+ DE_HRI_m(i__),DE_HRI_s(i__),hriperr(i__),theta_1(i__),
+ nsep_1(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg_5(i__),DEdeg_5(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'ihsample.dat' ! EDS Cross-Reference of HRI and IPC Sources
C Format for file interpretation
11 format(I4,1X,I2,1X,I5,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'ihsample.dat')
write(6,*) '....Loading file: ihsample.dat'
do i__=1,45
read(1,'(A15)')ar__10
read(ar__10,11)
+ Iseqno_5(i__),Isrcno_4(i__),Hseqno_5(i__),Hsrcno_3(i__)
c ..............Just test output...........
write(6,11)
+ Iseqno_5(i__),Isrcno_4(i__),Hseqno_5(i__),Hsrcno_3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'sample.dat' ! EDS IPC, HRI + Optical Data for 4.5Sigma IPC Sample
C Format for file interpretation
12 format(
+ I4,1X,I2,1X,I2,1X,I2,1X,I2,1X,A1,I2,1X,I2,1X,I2,1X,I2,1X,F7.4,
+ 1X,F6.4,1X,F6.3,1X,I5,1X,I1,1X,I3,1X,A27,1X,F5.2,1X,F6.3,1X,
+ F5.3,1X,A1,1X,F4.1,1X,A1,1X,F4.1,1X,F4.2,1X,A6)
C Effective file loading
open(unit=1,status='old',file=
+'sample.dat')
write(6,*) '....Loading file: sample.dat'
do i__=1,45
read(1,'(A136)')ar__11
read(ar__11,12)
+ Iseqno_6(i__),Isrcno_5(i__),RAh_6(i__),RAm_6(i__),RAs_6(i__),
+ DE__6(i__),DEd_6(i__),DEm_6(i__),DEs_6(i__),Iposerr(i__),
+ Iflux(i__),e_Iflux(i__),Isnr(i__),Hseqno_6(i__),Hsrcno_4(i__),
+ Hcell_1(i__),Hpos(i__),e_Hpos(i__),Hrate(i__),e_Hrate(i__),
+ l_Vmag(i__),Vmag(i__),l_fxfv(i__),fxfv(i__),z_1(i__),ID(i__)
if(ar__11(114:117) .EQ. '') Vmag(i__) = rNULL__
if(ar__11(121:124) .EQ. '') fxfv(i__) = rNULL__
if(ar__11(126:129) .EQ. '') z_1(i__) = rNULL__
RAdeg_6(i__) = rNULL__
DEdeg_6(i__) = rNULL__
c Derive coordinates RAdeg_6 and DEdeg_6 from input data
c (RAdeg_6 and DEdeg_6 are set to rNULL__ when unknown)
if(RAh_6(i__) .GT. -180) RAdeg_6(i__)=RAh_6(i__)*15.
if(RAm_6(i__) .GT. -180) RAdeg_6(i__)=RAdeg_6(i__)+RAm_6(i__)/4.
if(RAs_6(i__) .GT. -180) RAdeg_6(i__)=RAdeg_6(i__)+RAs_6(i__)/240.
if(DEd_6(i__) .GE. 0) DEdeg_6(i__)=DEd_6(i__)
if(DEm_6(i__) .GE. 0) DEdeg_6(i__)=DEdeg_6(i__)+DEm_6(i__)/60.
if(DEs_6(i__) .GE. 0) DEdeg_6(i__)=DEdeg_6(i__)+DEs_6(i__)/3600.
if(DE__6(i__).EQ.'-'.AND.DEdeg_6(i__).GE.0) DEdeg_6(i__)=-DEdeg_6(i__)
c ..............Just test output...........
write(6,12)
+ Iseqno_6(i__),Isrcno_5(i__),RAh_6(i__),RAm_6(i__),RAs_6(i__),
+ DE__6(i__),DEd_6(i__),DEm_6(i__),DEs_6(i__),Iposerr(i__),
+ Iflux(i__),e_Iflux(i__),Isnr(i__),Hseqno_6(i__),Hsrcno_4(i__),
+ Hcell_1(i__),Hpos(i__),e_Hpos(i__),Hrate(i__),e_Hrate(i__),
+ l_Vmag(i__),Vmag(i__),l_fxfv(i__),fxfv(i__),z_1(i__),ID(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg_6(i__),DEdeg_6(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'optinfo.dat' ! EDS Optical Data for Srcs in the 4.5 Sigma Sample
C Format for file interpretation
13 format(3X,I4,5X,I2,1X,A1,1X,F4.1,1X,A1,1X,F4.1,3X,F4.2,3X,A6)
C Effective file loading
open(unit=1,status='old',file=
+'optinfo.dat')
write(6,*) '....Loading file: optinfo.dat'
do i__=1,25
read(1,'(A45)')ar__12
read(ar__12,13)
+ Iseqno_7(i__),Isrcno_6(i__),l_Vmag_1(i__),Vmag_1(i__),
+ l_fxfv_1(i__),fxfv_1(i__),z_2(i__),ID_1(i__)
if(ar__12(18:21) .EQ. '') Vmag_1(i__) = rNULL__
if(ar__12(25:28) .EQ. '') fxfv_1(i__) = rNULL__
if(ar__12(32:35) .EQ. '') z_2(i__) = rNULL__
c ..............Just test output...........
write(6,13)
+ Iseqno_7(i__),Isrcno_6(i__),l_Vmag_1(i__),Vmag_1(i__),
+ l_fxfv_1(i__),fxfv_1(i__),z_2(i__),ID_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end