Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/405/2062 into FORTRAN code for reading data files line by line.
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. 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 ! (deg) Right ascension (J2000) in 2MASS
real*8 DEdeg ! (deg) Declination (J2000) in 2MASS
real*4 Jmag ! (mag) 2MASS J selected default magnitude
real*4 Hmag ! (mag) 2MASS H selected default magnitude
real*4 Ksmag ! (mag) 2MASS Ks selected default magnitude
real*8 CRate ! (ct/s) Source count rate (1)
real*8 e_CRate ! (ct/s) ?=0 Error on CRate (1)
real*4 bgCt ! (ct/s/arcmin2) ?=0 Background count rate (1)
integer*4 texp ! (s) Exposure time
real*4 HR1 ! [-1/1]? Hardness ratio 1 (2)
real*4 e_HR1 ! ? Error on HR1 (3)
real*4 HR2 ! [-1/1]? Hardness ratio 2 (2)
real*4 e_HR2 ! ? Error on HR2 (3)
integer*4 Extent ! (arcsec) Source extent, by which the source image
* exceeds the point spread function
real*8 Obs_beg ! ("YYMMDD") ? Start time of observation (FSC only)
real*8 Obs_end ! ("YYMMDD") ? End time of observation (FSC only)
integer*4 Rel ! ? 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 ! (deg) Right ascension (J2000) in 2MASS
real*8 DEdeg_1 ! (deg) Declination (J2000) in 2MASS
real*4 Jmag_1 ! (mag) 2MASS J selected default magnitude
real*4 Hmag_1 ! (mag) 2MASS H selected default magnitude
real*4 Ksmag_1 ! (mag) 2MASS Ks selected default magnitude
real*8 CRate_1 ! (ct/s) Source count rate (1)
real*8 e_CRate_1 ! (ct/s) ?=0 Error on CRate (1)
real*4 bgCt_1 ! (ct/s/arcmin2) ?=0 Background count rate (1)
integer*4 texp_1 ! (s) Exposure time
real*4 HR1_1 ! [-1/1]? Hardness ratio 1 (2)
real*4 e_HR1_1 ! ? Error on HR1 (3)
real*4 HR2_1 ! [-1/1]? Hardness ratio 2 (2)
real*4 e_HR2_1 ! ? Error on HR2 (3)
integer*4 Extent_1 ! (arcsec) Source extent, by which the source image
* exceeds the point spread function
real*8 Obs_beg_1 ! ("YYMMDD") ? Start time of observation (FSC only)
real*8 Obs_end_1 ! ("YYMMDD") ? End time of observation (FSC only)
integer*4 Rel_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 ! (deg) Right ascension (J2000) in 2MASS
real*8 DEdeg_2 ! (deg) Declination (J2000) in 2MASS
real*4 Jmag_2 ! (mag) 2MASS J selected default magnitude
real*4 Hmag_2 ! (mag) 2MASS H selected default magnitude
real*4 Ksmag_2 ! (mag) 2MASS Ks selected default magnitude
real*8 CRate_2 ! (ct/s) Source count rate (1)
real*8 e_CRate_2 ! (ct/s) ?=0 Error on CRate (1)
real*4 bgCt_2 ! (ct/s/arcmin2) ?=0 Background count rate (1)
integer*4 texp_2 ! (s) Exposure time
real*4 HR1_2 ! [-1/1]? Hardness ratio 1 (2)
real*4 e_HR1_2 ! ? Error on HR1 (3)
real*4 HR2_2 ! [-1/1]? Hardness ratio 2 (2)
real*4 e_HR2_2 ! ? Error on HR2 (3)
integer*4 Extent_2 ! (arcsec) Source extent, by which the source image
* exceeds the point spread function
real*8 Obs_beg_2 ! ("YYMMDD") ? Start time of observation (FSC only)
real*8 Obs_end_2 ! ("YYMMDD") ? End time of observation (FSC only)
integer*4 Rel_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 ! (deg) Right ascension (J2000) in 2MASS
real*8 DEdeg_3 ! (deg) Declination (J2000) in 2MASS
real*4 Jmag_3 ! (mag) 2MASS J selected default magnitude
real*4 Hmag_3 ! (mag) 2MASS H selected default magnitude
real*4 Ksmag_3 ! (mag) 2MASS Ks selected default magnitude
real*8 CRate_3 ! (ct/s) Source count rate (1)
real*8 e_CRate_3 ! (ct/s) ?=0 Error on CRate (1)
real*4 bgCt_3 ! (ct/s/arcmin2) ?=0 Background count rate (1)
integer*4 texp_3 ! (s) Exposure time
real*4 HR1_3 ! [-1/1]? Hardness ratio 1 (2)
real*4 e_HR1_3 ! ? Error on HR1 (3)
real*4 HR2_3 ! [-1/1]? Hardness ratio 2 (2)
real*4 e_HR2_3 ! ? Error on HR2 (3)
integer*4 Extent_3 ! (arcsec) Source extent, by which the source image
* exceeds the point spread function
real*8 Obs_beg_3 ! ("YYMMDD") ? Start time of observation (FSC only)
real*8 Obs_end_3 ! ("YYMMDD") ? End time of observation (FSC only)
integer*4 Rel_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,DEdeg,Jmag,Hmag,Ksmag,CRate,e_CRate,bgCt,texp,HR1,e_HR1,
+ HR2,e_HR2,Extent,Obs_beg,Obs_end,Rel
if(ar__(84:88) .EQ. '') HR1 = rNULL__
if(ar__(90:94) .EQ. '') e_HR1 = rNULL__
if(ar__(96:100) .EQ. '') HR2 = rNULL__
if(ar__(102:106) .EQ. '') e_HR2 = rNULL__
if(ar__(112:120) .EQ. '') Obs_beg = rNULL__
if(ar__(122:130) .EQ. '') Obs_end = rNULL__
if(ar__(133:134) .EQ. '') Rel = iNULL__
c ..............Just test output...........
write(6,1)
+ RAdeg,DEdeg,Jmag,Hmag,Ksmag,CRate,e_CRate,bgCt,texp,HR1,e_HR1,
+ HR2,e_HR2,Extent,Obs_beg,Obs_end,Rel
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,DEdeg_1,Jmag_1,Hmag_1,Ksmag_1,CRate_1,e_CRate_1,
+ bgCt_1,texp_1,HR1_1,e_HR1_1,HR2_1,e_HR2_1,Extent_1,Obs_beg_1,
+ Obs_end_1,Rel_1
if(ar__1(84:88) .EQ. '') HR1_1 = rNULL__
if(ar__1(90:94) .EQ. '') e_HR1_1 = rNULL__
if(ar__1(96:100) .EQ. '') HR2_1 = rNULL__
if(ar__1(102:106) .EQ. '') e_HR2_1 = rNULL__
if(ar__1(112:120) .EQ. '') Obs_beg_1 = rNULL__
if(ar__1(122:130) .EQ. '') Obs_end_1 = rNULL__
if(ar__1(133:134) .EQ. '') Rel_1 = iNULL__
c ..............Just test output...........
write(6,2)
+ RAdeg_1,DEdeg_1,Jmag_1,Hmag_1,Ksmag_1,CRate_1,e_CRate_1,
+ bgCt_1,texp_1,HR1_1,e_HR1_1,HR2_1,e_HR2_1,Extent_1,Obs_beg_1,
+ Obs_end_1,Rel_1
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,DEdeg_2,Jmag_2,Hmag_2,Ksmag_2,CRate_2,e_CRate_2,
+ bgCt_2,texp_2,HR1_2,e_HR1_2,HR2_2,e_HR2_2,Extent_2,Obs_beg_2,
+ Obs_end_2,Rel_2
if(ar__2(84:88) .EQ. '') HR1_2 = rNULL__
if(ar__2(90:94) .EQ. '') e_HR1_2 = rNULL__
if(ar__2(96:100) .EQ. '') HR2_2 = rNULL__
if(ar__2(102:106) .EQ. '') e_HR2_2 = rNULL__
if(ar__2(112:120) .EQ. '') Obs_beg_2 = rNULL__
if(ar__2(122:130) .EQ. '') Obs_end_2 = rNULL__
if(ar__2(133:134) .EQ. '') Rel_2 = iNULL__
c ..............Just test output...........
write(6,3)
+ RAdeg_2,DEdeg_2,Jmag_2,Hmag_2,Ksmag_2,CRate_2,e_CRate_2,
+ bgCt_2,texp_2,HR1_2,e_HR1_2,HR2_2,e_HR2_2,Extent_2,Obs_beg_2,
+ Obs_end_2,Rel_2
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,DEdeg_3,Jmag_3,Hmag_3,Ksmag_3,CRate_3,e_CRate_3,
+ bgCt_3,texp_3,HR1_3,e_HR1_3,HR2_3,e_HR2_3,Extent_3,Obs_beg_3,
+ Obs_end_3,Rel_3
if(ar__3(84:88) .EQ. '') HR1_3 = rNULL__
if(ar__3(90:94) .EQ. '') e_HR1_3 = rNULL__
if(ar__3(96:100) .EQ. '') HR2_3 = rNULL__
if(ar__3(102:106) .EQ. '') e_HR2_3 = rNULL__
if(ar__3(112:120) .EQ. '') Obs_beg_3 = rNULL__
if(ar__3(122:130) .EQ. '') Obs_end_3 = rNULL__
if(ar__3(133:134) .EQ. '') Rel_3 = iNULL__
c ..............Just test output...........
write(6,4)
+ RAdeg_3,DEdeg_3,Jmag_3,Hmag_3,Ksmag_3,CRate_3,e_CRate_3,
+ bgCt_3,texp_3,HR1_3,e_HR1_3,HR2_3,e_HR2_3,Extent_3,Obs_beg_3,
+ Obs_end_3,Rel_3
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end