Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/602/A115 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/A+A/602/A115 X-ray survey of NGC7000/IC5070 (Damiani+, 2017)
*================================================================================
*X-ray survey of the North-America and Pelican star-forming complex
*(NGC7000/IC5070).
* Damiani F., Pillitteri I., Prisinzano L.
* <Astron. Astrophys. 602, A115 (2017)>
* =2017A&A...602A.115D (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! X-ray source properties
integer*4 nr__
parameter (nr__=721) ! Number of records
character*107 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 X_rayId (nr__) ! [1/721] Unique X-ray source ID
character*16 CXOId (nr__) ! IAU designation of Chandra X-ray detections
* (JHHMMSS.s+DDMMSS)
character*16 XMMId (nr__) ! IAU designation of XMM X-ray detections
* (JHHMMSS.s+DDMMSS)
real*8 RAdeg (nr__) ! (deg) Right ascension (J2000)
real*8 DEdeg (nr__) ! (deg) Declination (J2000)
real*4 e_pos (nr__) ! (arcsec) Position error of X-ray source
real*8 CRate (nr__) ! (ct/ks) X-ray count-rate
real*4 e_CRate (nr__) ! (ct/ks) Error on CRate
real*4 ACIS_HR1 (nr__) ! ? Hardness ratio 1 for Chandra ACIS sources
real*4 ACIS_HR2 (nr__) ! ? Hardness ratio 2 for Chandra ACIS sources
real*4 XMM_HR1 (nr__) ! ? Hardness ratio 1 for XMM EPIC sources
real*4 XMM_HR2 (nr__) ! ? Hardness ratio 2 for XMM EPIC sources
real*4 logLX (nr__) ! ([10-7W]) ? log of X-ray luminosity
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! X-ray source identifications
integer*4 nr__1
parameter (nr__1=721) ! Number of records
character*46 ar__1 ! Full-size record
integer*4 X_rayId_1 (nr__1) ! [1/721] Unique X-ray source ID
real*4 imag (nr__1) ! (mag) ? Magnitude in i band from IPHAS
real*4 r_i (nr__1) ! (mag) ? Color index r-i from IPHAS
real*4 r_Ha (nr__1) ! (mag) ? Color index r-Halpha from IPHAS
real*4 Jmag (nr__1) ! (mag) ? Magnitude in J band from UKIDSS/2MASS
real*4 Hmag (nr__1) ! (mag) ? Magnitude in H band from UKIDSS/2MASS
real*4 Kmag (nr__1) ! (mag) ? Magnitude in K band from UKIDSS/2MASS
character*1 IR_excess (nr__1) ! [Y] Star with NIR excess
character*1 Ha_emiss (nr__1) ! [Y] Star with Halpha emission
integer*4 v_RGS2011_ (nr__1) ! ? Sequential number in Rebull et al.
* (2011, Cat. J/ApJS/193/25)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! NIR-excess and Halpha-emission stars
integer*4 nr__2
parameter (nr__2=273) ! Number of records
character*65 ar__2 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 Seq (nr__2) ! [722/994] Unique X-ray undetected member ID
real*8 RAdeg_1 (nr__2) ! (deg) Right ascension (J2000)
real*8 DEdeg_1 (nr__2) ! (deg) Declination (J2000)
real*4 imag_1 (nr__2) ! (mag) ? Magnitude in i band from IPHAS
real*4 r_i_1 (nr__2) ! (mag) ? Color index r-i from IPHAS
real*4 r_Ha_1 (nr__2) ! (mag) ? Color index r-Halpha from IPHAS
real*4 Jmag_1 (nr__2) ! (mag) ? Magnitude in J band from UKIDSS/2MASS
real*4 Hmag_1 (nr__2) ! (mag) ? Magnitude in H band from UKIDSS/2MASS
real*4 Kmag_1 (nr__2) ! (mag) ? Magnitude in K band from UKIDSS/2MASS
character*1 IR_excess_1(nr__2) ! [Y] Star with NIR excess
character*1 Ha_emiss_1 (nr__2) ! [Y] Star with Halpha emission
integer*4 v_RGS2011__1(nr__2) ! ? Sequential number in Rebull et al.
* (2011, Cat. J/ApJS/193/25)
C=============================================================================
C Loading file 'table3.dat' ! X-ray source properties
C Format for file interpretation
1 format(
+ I3,1X,A16,1X,A16,1X,F9.5,1X,F8.5,1X,F4.2,1X,F8.3,1X,F6.3,1X,
+ F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,721
read(1,'(A107)')ar__
read(ar__,1)
+ X_rayId(i__),CXOId(i__),XMMId(i__),RAdeg(i__),DEdeg(i__),
+ e_pos(i__),CRate(i__),e_CRate(i__),ACIS_HR1(i__),
+ ACIS_HR2(i__),XMM_HR1(i__),XMM_HR2(i__),logLX(i__)
if(ar__(79:83) .EQ. '') ACIS_HR1(i__) = rNULL__
if(ar__(85:89) .EQ. '') ACIS_HR2(i__) = rNULL__
if(ar__(91:95) .EQ. '') XMM_HR1(i__) = rNULL__
if(ar__(97:101) .EQ. '') XMM_HR2(i__) = rNULL__
if(ar__(103:107) .EQ. '') logLX(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ X_rayId(i__),CXOId(i__),XMMId(i__),RAdeg(i__),DEdeg(i__),
+ e_pos(i__),CRate(i__),e_CRate(i__),ACIS_HR1(i__),
+ ACIS_HR2(i__),XMM_HR1(i__),XMM_HR2(i__),logLX(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! X-ray source identifications
C Format for file interpretation
2 format(
+ I3,1X,F5.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,A1,1X,
+ A1,1X,I4)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,721
read(1,'(A46)')ar__1
read(ar__1,2)
+ X_rayId_1(i__),imag(i__),r_i(i__),r_Ha(i__),Jmag(i__),
+ Hmag(i__),Kmag(i__),IR_excess(i__),Ha_emiss(i__),
+ v_RGS2011_(i__)
if(ar__1(5:9) .EQ. '') imag(i__) = rNULL__
if(ar__1(11:14) .EQ. '') r_i(i__) = rNULL__
if(ar__1(16:19) .EQ. '') r_Ha(i__) = rNULL__
if(ar__1(21:25) .EQ. '') Jmag(i__) = rNULL__
if(ar__1(27:31) .EQ. '') Hmag(i__) = rNULL__
if(ar__1(33:37) .EQ. '') Kmag(i__) = rNULL__
if(ar__1(43:46) .EQ. '') v_RGS2011_(i__) = iNULL__
c ..............Just test output...........
write(6,2)
+ X_rayId_1(i__),imag(i__),r_i(i__),r_Ha(i__),Jmag(i__),
+ Hmag(i__),Kmag(i__),IR_excess(i__),Ha_emiss(i__),
+ v_RGS2011_(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! NIR-excess and Halpha-emission stars
C Format for file interpretation
3 format(
+ I3,1X,F9.5,1X,F8.5,1X,F5.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F5.2,1X,
+ F5.2,1X,A1,1X,A1,1X,I4)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,273
read(1,'(A65)')ar__2
read(ar__2,3)
+ Seq(i__),RAdeg_1(i__),DEdeg_1(i__),imag_1(i__),r_i_1(i__),
+ r_Ha_1(i__),Jmag_1(i__),Hmag_1(i__),Kmag_1(i__),
+ IR_excess_1(i__),Ha_emiss_1(i__),v_RGS2011__1(i__)
if(ar__2(24:28) .EQ. '') imag_1(i__) = rNULL__
if(ar__2(30:33) .EQ. '') r_i_1(i__) = rNULL__
if(ar__2(35:38) .EQ. '') r_Ha_1(i__) = rNULL__
if(ar__2(40:44) .EQ. '') Jmag_1(i__) = rNULL__
if(ar__2(46:50) .EQ. '') Hmag_1(i__) = rNULL__
if(ar__2(52:56) .EQ. '') Kmag_1(i__) = rNULL__
if(ar__2(62:65) .EQ. '') v_RGS2011__1(i__) = iNULL__
c ..............Just test output...........
write(6,3)
+ Seq(i__),RAdeg_1(i__),DEdeg_1(i__),imag_1(i__),r_i_1(i__),
+ r_Ha_1(i__),Jmag_1(i__),Hmag_1(i__),Kmag_1(i__),
+ IR_excess_1(i__),Ha_emiss_1(i__),v_RGS2011__1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end