Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/405/2062 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-16
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. J/MNRAS/405/2062 AGN candidates from 2MASS/ROSAT catalogs (Kouzuma+, 2010)
*================================================================================
*Properties of near-infrared-selected AGN candidates with 2MASS/ROSAT catalogues.
* Kouzuma S., Yamaoka H.
* <Mon. Not. R. Astron. Soc., 405, 2062-2068 (2010)>
* =2010MNRAS.405.2062K
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'bsc_k.dat' ! Known AGNs in the Rosat-BSC (Cat. IX/10)
integer*4 nr__
parameter (nr__=2053) ! Number of records
character*134 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
real*8 RAdeg (nr__) ! (deg) Right ascension (J2000) in 2MASS
real*8 DEdeg (nr__) ! (deg) Declination (J2000) in 2MASS
real*4 Jmag (nr__) ! (mag) 2MASS J selected default magnitude
real*4 Hmag (nr__) ! (mag) 2MASS H selected default magnitude
real*4 Ksmag (nr__) ! (mag) 2MASS Ks selected default magnitude
real*8 CRate (nr__) ! (ct/s) Source count rate (1)
real*8 e_CRate (nr__) ! (ct/s) ?=0 Error on CRate (1)
real*4 bgCt (nr__) ! (ct/s/arcmin2) ?=0 Background count rate (1)
integer*4 texp (nr__) ! (s) Exposure time
real*4 HR1 (nr__) ! [-1/1]? Hardness ratio 1 (2)
real*4 e_HR1 (nr__) ! ? Error on HR1 (3)
real*4 HR2 (nr__) ! [-1/1]? Hardness ratio 2 (2)
real*4 e_HR2 (nr__) ! ? Error on HR2 (3)
integer*4 Extent (nr__) ! (arcsec) Source extent, by which the source image
* exceeds the point spread function
real*8 Obs_beg (nr__) ! ("YYMMDD") ? Start time of observation (FSC only)
real*8 Obs_end (nr__) ! ("YYMMDD") ? End time of observation (FSC only)
integer*4 Rel (nr__) ! ? Reliability of source detection
* (FSC data only)
*Note (1): Countrates are vignetting corrected
*Note (2): The hardness ratios are defined by:
* HR1=(B-A)/(B+A), where
* A = count rate in PHA range 11-41 (~0.1-0.4keV)
* B = count rate in PHA range 52-201 (~0.5-2.0keV)
* HR2=(D-C)/(D+C), where
* C = count rate in PHA range 52-90 (~0.5-0.9keV)
* D = count rate in PHA range 91-201 (~0.9-2.0keV)
* (PHA = Pulse Height Amplitude)
*Note (3): The true mean error is probably bigger, and values larger than
* 9.99 are set to 9.99
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'bsc_u.dat' ! Unclassified sources (AGN candidates) in the BSC
integer*4 nr__1
parameter (nr__1=3220) ! Number of records
character*134 ar__1 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
real*8 RAdeg_1 (nr__1) ! (deg) Right ascension (J2000) in 2MASS
real*8 DEdeg_1 (nr__1) ! (deg) Declination (J2000) in 2MASS
real*4 Jmag_1 (nr__1) ! (mag) 2MASS J selected default magnitude
real*4 Hmag_1 (nr__1) ! (mag) 2MASS H selected default magnitude
real*4 Ksmag_1 (nr__1) ! (mag) 2MASS Ks selected default magnitude
real*8 CRate_1 (nr__1) ! (ct/s) Source count rate (1)
real*8 e_CRate_1 (nr__1) ! (ct/s) ?=0 Error on CRate (1)
real*4 bgCt_1 (nr__1) ! (ct/s/arcmin2) ?=0 Background count rate (1)
integer*4 texp_1 (nr__1) ! (s) Exposure time
real*4 HR1_1 (nr__1) ! [-1/1]? Hardness ratio 1 (2)
real*4 e_HR1_1 (nr__1) ! ? Error on HR1 (3)
real*4 HR2_1 (nr__1) ! [-1/1]? Hardness ratio 2 (2)
real*4 e_HR2_1 (nr__1) ! ? Error on HR2 (3)
integer*4 Extent_1 (nr__1) ! (arcsec) Source extent, by which the source image
* exceeds the point spread function
real*8 Obs_beg_1 (nr__1) ! ("YYMMDD") ? Start time of observation (FSC only)
real*8 Obs_end_1 (nr__1) ! ("YYMMDD") ? End time of observation (FSC only)
integer*4 Rel_1 (nr__1) ! ? Reliability of source detection
* (FSC data only)
*Note (1): Countrates are vignetting corrected
*Note (2): The hardness ratios are defined by:
* HR1=(B-A)/(B+A), where
* A = count rate in PHA range 11-41 (~0.1-0.4keV)
* B = count rate in PHA range 52-201 (~0.5-2.0keV)
* HR2=(D-C)/(D+C), where
* C = count rate in PHA range 52-90 (~0.5-0.9keV)
* D = count rate in PHA range 91-201 (~0.9-2.0keV)
* (PHA = Pulse Height Amplitude)
*Note (3): The true mean error is probably bigger, and values larger than
* 9.99 are set to 9.99
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'fsc_k.dat' ! Known AGNs in the Rosat-FSC (Cat. IX/29)
integer*4 nr__2
parameter (nr__2=1008) ! Number of records
character*134 ar__2 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
real*8 RAdeg_2 (nr__2) ! (deg) Right ascension (J2000) in 2MASS
real*8 DEdeg_2 (nr__2) ! (deg) Declination (J2000) in 2MASS
real*4 Jmag_2 (nr__2) ! (mag) 2MASS J selected default magnitude
real*4 Hmag_2 (nr__2) ! (mag) 2MASS H selected default magnitude
real*4 Ksmag_2 (nr__2) ! (mag) 2MASS Ks selected default magnitude
real*8 CRate_2 (nr__2) ! (ct/s) Source count rate (1)
real*8 e_CRate_2 (nr__2) ! (ct/s) ?=0 Error on CRate (1)
real*4 bgCt_2 (nr__2) ! (ct/s/arcmin2) ?=0 Background count rate (1)
integer*4 texp_2 (nr__2) ! (s) Exposure time
real*4 HR1_2 (nr__2) ! [-1/1]? Hardness ratio 1 (2)
real*4 e_HR1_2 (nr__2) ! ? Error on HR1 (3)
real*4 HR2_2 (nr__2) ! [-1/1]? Hardness ratio 2 (2)
real*4 e_HR2_2 (nr__2) ! ? Error on HR2 (3)
integer*4 Extent_2 (nr__2) ! (arcsec) Source extent, by which the source image
* exceeds the point spread function
real*8 Obs_beg_2 (nr__2) ! ("YYMMDD") ? Start time of observation (FSC only)
real*8 Obs_end_2 (nr__2) ! ("YYMMDD") ? End time of observation (FSC only)
integer*4 Rel_2 (nr__2) ! ? Reliability of source detection
* (FSC data only)
*Note (1): Countrates are vignetting corrected
*Note (2): The hardness ratios are defined by:
* HR1=(B-A)/(B+A), where
* A = count rate in PHA range 11-41 (~0.1-0.4keV)
* B = count rate in PHA range 52-201 (~0.5-2.0keV)
* HR2=(D-C)/(D+C), where
* C = count rate in PHA range 52-90 (~0.5-0.9keV)
* D = count rate in PHA range 91-201 (~0.9-2.0keV)
* (PHA = Pulse Height Amplitude)
*Note (3): The true mean error is probably bigger, and values larger than
* 9.99 are set to 9.99
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'fsc_u.dat' ! Unclassified sources (AGN candidates) in the FSC
integer*4 nr__3
parameter (nr__3=9693) ! Number of records
character*134 ar__3 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
real*8 RAdeg_3 (nr__3) ! (deg) Right ascension (J2000) in 2MASS
real*8 DEdeg_3 (nr__3) ! (deg) Declination (J2000) in 2MASS
real*4 Jmag_3 (nr__3) ! (mag) 2MASS J selected default magnitude
real*4 Hmag_3 (nr__3) ! (mag) 2MASS H selected default magnitude
real*4 Ksmag_3 (nr__3) ! (mag) 2MASS Ks selected default magnitude
real*8 CRate_3 (nr__3) ! (ct/s) Source count rate (1)
real*8 e_CRate_3 (nr__3) ! (ct/s) ?=0 Error on CRate (1)
real*4 bgCt_3 (nr__3) ! (ct/s/arcmin2) ?=0 Background count rate (1)
integer*4 texp_3 (nr__3) ! (s) Exposure time
real*4 HR1_3 (nr__3) ! [-1/1]? Hardness ratio 1 (2)
real*4 e_HR1_3 (nr__3) ! ? Error on HR1 (3)
real*4 HR2_3 (nr__3) ! [-1/1]? Hardness ratio 2 (2)
real*4 e_HR2_3 (nr__3) ! ? Error on HR2 (3)
integer*4 Extent_3 (nr__3) ! (arcsec) Source extent, by which the source image
* exceeds the point spread function
real*8 Obs_beg_3 (nr__3) ! ("YYMMDD") ? Start time of observation (FSC only)
real*8 Obs_end_3 (nr__3) ! ("YYMMDD") ? End time of observation (FSC only)
integer*4 Rel_3 (nr__3) ! ? Reliability of source detection
* (FSC data only)
*Note (1): Countrates are vignetting corrected
*Note (2): The hardness ratios are defined by:
* HR1=(B-A)/(B+A), where
* A = count rate in PHA range 11-41 (~0.1-0.4keV)
* B = count rate in PHA range 52-201 (~0.5-2.0keV)
* HR2=(D-C)/(D+C), where
* C = count rate in PHA range 52-90 (~0.5-0.9keV)
* D = count rate in PHA range 91-201 (~0.9-2.0keV)
* (PHA = Pulse Height Amplitude)
*Note (3): The true mean error is probably bigger, and values larger than
* 9.99 are set to 9.99
C=============================================================================
C Loading file 'bsc_k.dat' ! Known AGNs in the Rosat-BSC (Cat. IX/10)
C Format for file interpretation
1 format(
+ F10.6,1X,F10.6,1X,F6.3,1X,F6.3,1X,F6.3,1X,F9.5,1X,F9.5,1X,
+ E11.4,1X,I5,3X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,I3,1X,F9.2,1X,
+ F9.2,2X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'bsc_k.dat')
write(6,*) '....Loading file: bsc_k.dat'
do i__=1,2053
read(1,'(A134)')ar__
read(ar__,1)
+ RAdeg(i__),DEdeg(i__),Jmag(i__),Hmag(i__),Ksmag(i__),
+ CRate(i__),e_CRate(i__),bgCt(i__),texp(i__),HR1(i__),
+ e_HR1(i__),HR2(i__),e_HR2(i__),Extent(i__),Obs_beg(i__),
+ Obs_end(i__),Rel(i__)
if(ar__(84:88) .EQ. '') HR1(i__) = rNULL__
if(ar__(90:94) .EQ. '') e_HR1(i__) = rNULL__
if(ar__(96:100) .EQ. '') HR2(i__) = rNULL__
if(ar__(102:106) .EQ. '') e_HR2(i__) = rNULL__
if(ar__(112:120) .EQ. '') Obs_beg(i__) = rNULL__
if(ar__(122:130) .EQ. '') Obs_end(i__) = rNULL__
if(ar__(133:134) .EQ. '') Rel(i__) = iNULL__
c ..............Just test output...........
write(6,1)
+ RAdeg(i__),DEdeg(i__),Jmag(i__),Hmag(i__),Ksmag(i__),
+ CRate(i__),e_CRate(i__),bgCt(i__),texp(i__),HR1(i__),
+ e_HR1(i__),HR2(i__),e_HR2(i__),Extent(i__),Obs_beg(i__),
+ Obs_end(i__),Rel(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'bsc_u.dat' ! Unclassified sources (AGN candidates) in the BSC
C Format for file interpretation
2 format(
+ F10.6,1X,F10.6,1X,F6.3,1X,F6.3,1X,F6.3,1X,F9.5,1X,F9.5,1X,
+ E11.4,1X,I5,3X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,I3,1X,F9.2,1X,
+ F9.2,2X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'bsc_u.dat')
write(6,*) '....Loading file: bsc_u.dat'
do i__=1,3220
read(1,'(A134)')ar__1
read(ar__1,2)
+ RAdeg_1(i__),DEdeg_1(i__),Jmag_1(i__),Hmag_1(i__),
+ Ksmag_1(i__),CRate_1(i__),e_CRate_1(i__),bgCt_1(i__),
+ texp_1(i__),HR1_1(i__),e_HR1_1(i__),HR2_1(i__),e_HR2_1(i__),
+ Extent_1(i__),Obs_beg_1(i__),Obs_end_1(i__),Rel_1(i__)
if(ar__1(84:88) .EQ. '') HR1_1(i__) = rNULL__
if(ar__1(90:94) .EQ. '') e_HR1_1(i__) = rNULL__
if(ar__1(96:100) .EQ. '') HR2_1(i__) = rNULL__
if(ar__1(102:106) .EQ. '') e_HR2_1(i__) = rNULL__
if(ar__1(112:120) .EQ. '') Obs_beg_1(i__) = rNULL__
if(ar__1(122:130) .EQ. '') Obs_end_1(i__) = rNULL__
if(ar__1(133:134) .EQ. '') Rel_1(i__) = iNULL__
c ..............Just test output...........
write(6,2)
+ RAdeg_1(i__),DEdeg_1(i__),Jmag_1(i__),Hmag_1(i__),
+ Ksmag_1(i__),CRate_1(i__),e_CRate_1(i__),bgCt_1(i__),
+ texp_1(i__),HR1_1(i__),e_HR1_1(i__),HR2_1(i__),e_HR2_1(i__),
+ Extent_1(i__),Obs_beg_1(i__),Obs_end_1(i__),Rel_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'fsc_k.dat' ! Known AGNs in the Rosat-FSC (Cat. IX/29)
C Format for file interpretation
3 format(
+ F10.6,1X,F10.6,1X,F6.3,1X,F6.3,1X,F6.3,1X,F9.5,1X,F9.5,1X,
+ E11.4,1X,I5,3X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,I3,1X,F9.2,1X,
+ F9.2,2X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'fsc_k.dat')
write(6,*) '....Loading file: fsc_k.dat'
do i__=1,1008
read(1,'(A134)')ar__2
read(ar__2,3)
+ RAdeg_2(i__),DEdeg_2(i__),Jmag_2(i__),Hmag_2(i__),
+ Ksmag_2(i__),CRate_2(i__),e_CRate_2(i__),bgCt_2(i__),
+ texp_2(i__),HR1_2(i__),e_HR1_2(i__),HR2_2(i__),e_HR2_2(i__),
+ Extent_2(i__),Obs_beg_2(i__),Obs_end_2(i__),Rel_2(i__)
if(ar__2(84:88) .EQ. '') HR1_2(i__) = rNULL__
if(ar__2(90:94) .EQ. '') e_HR1_2(i__) = rNULL__
if(ar__2(96:100) .EQ. '') HR2_2(i__) = rNULL__
if(ar__2(102:106) .EQ. '') e_HR2_2(i__) = rNULL__
if(ar__2(112:120) .EQ. '') Obs_beg_2(i__) = rNULL__
if(ar__2(122:130) .EQ. '') Obs_end_2(i__) = rNULL__
if(ar__2(133:134) .EQ. '') Rel_2(i__) = iNULL__
c ..............Just test output...........
write(6,3)
+ RAdeg_2(i__),DEdeg_2(i__),Jmag_2(i__),Hmag_2(i__),
+ Ksmag_2(i__),CRate_2(i__),e_CRate_2(i__),bgCt_2(i__),
+ texp_2(i__),HR1_2(i__),e_HR1_2(i__),HR2_2(i__),e_HR2_2(i__),
+ Extent_2(i__),Obs_beg_2(i__),Obs_end_2(i__),Rel_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'fsc_u.dat' ! Unclassified sources (AGN candidates) in the FSC
C Format for file interpretation
4 format(
+ F10.6,1X,F10.6,1X,F6.3,1X,F6.3,1X,F6.3,1X,F9.5,1X,F9.5,1X,
+ E11.4,1X,I5,3X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,I3,1X,F9.2,1X,
+ F9.2,2X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'fsc_u.dat')
write(6,*) '....Loading file: fsc_u.dat'
do i__=1,9693
read(1,'(A134)')ar__3
read(ar__3,4)
+ RAdeg_3(i__),DEdeg_3(i__),Jmag_3(i__),Hmag_3(i__),
+ Ksmag_3(i__),CRate_3(i__),e_CRate_3(i__),bgCt_3(i__),
+ texp_3(i__),HR1_3(i__),e_HR1_3(i__),HR2_3(i__),e_HR2_3(i__),
+ Extent_3(i__),Obs_beg_3(i__),Obs_end_3(i__),Rel_3(i__)
if(ar__3(84:88) .EQ. '') HR1_3(i__) = rNULL__
if(ar__3(90:94) .EQ. '') e_HR1_3(i__) = rNULL__
if(ar__3(96:100) .EQ. '') HR2_3(i__) = rNULL__
if(ar__3(102:106) .EQ. '') e_HR2_3(i__) = rNULL__
if(ar__3(112:120) .EQ. '') Obs_beg_3(i__) = rNULL__
if(ar__3(122:130) .EQ. '') Obs_end_3(i__) = rNULL__
if(ar__3(133:134) .EQ. '') Rel_3(i__) = iNULL__
c ..............Just test output...........
write(6,4)
+ RAdeg_3(i__),DEdeg_3(i__),Jmag_3(i__),Hmag_3(i__),
+ Ksmag_3(i__),CRate_3(i__),e_CRate_3(i__),bgCt_3(i__),
+ texp_3(i__),HR1_3(i__),e_HR1_3(i__),HR2_3(i__),e_HR2_3(i__),
+ Extent_3(i__),Obs_beg_3(i__),Obs_end_3(i__),Rel_3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end