Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/753/27 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 2026-Oct-15
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/ApJ/753/27 New point sources in the WMAP 7yr data (Scodeller+, 2012)
*================================================================================
*Detection of new point sources in WMAP 7 year data using internal templates and
*needlets.
* Scodeller S., Hansen F.K., Marinucci D.
* <Astrophys. J., 753, 27 (2012)>
* =2012ApJ...753...27S
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! New point sources in WMAP
integer*4 nr__
parameter (nr__=1116) ! Number of records
character*122 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 Seq (nr__) ! Running sequence number
real*8 RAdeg (nr__) ! (deg) Right Ascension in decimal degrees (J2000) (1)
real*8 DEdeg (nr__) ! (deg) Declination in decimal degrees (J2000) (1)
real*4 FK (nr__) ! (Jy) [0.1/24]? WMAP K band (23GHz) flux density (2)
real*4 e_FK (nr__) ! (Jy) [0.04/0.3]? Uncertainty in FK
character*1 f_FK (nr__) ! [i] very weak FK (3)
real*4 FKa (nr__) ! (Jy) [0.1/22]? WMAP Ka band (30GHz) flux density (2)
real*4 e_FKa (nr__) ! (Jy) ? Uncertainty in FKa
character*1 f_FKa (nr__) ! [i] very weak FKa (3)
real*4 FQ (nr__) ! (Jy) [0.1/20]? WMAP Q band (41GHz) flux density (2)
real*4 e_FQ (nr__) ! (Jy) ? Uncertainty in FQ
character*1 f_FQ (nr__) ! [i] very weak FQ (3)
real*4 FV (nr__) ! (Jy) [0.1/18]? WMAP V band (61GHz) flux density (2)
real*4 e_FV (nr__) ! (Jy) ? Uncertainty in FV
character*1 f_FV (nr__) ! [i] very weak FV (3)
real*4 FW (nr__) ! (Jy) [0.2/15]? WMAP W band (94GHz) flux density (2)
real*4 e_FW (nr__) ! (Jy) ? Uncertainty in FW
character*1 f_FW (nr__) ! [i] very weak FW (3)
character*5 Fflag (nr__) ! The 5{sigma} detection flag (4)
character*3 Cflag (nr__) ! Catalog flag (5)
character*19 Cat_ID (nr__) ! Source identification (6)
real*4 Off (nr__) ! (arcmin) [0.1/15]? Offset from the Cat-ID position
character*2 f_Off (nr__) ! [MA] Multiple of weak source (7)
*Note (1): Given by the weighted average of the positions in all frequencies
* where the source is found.
*Note (2): For converting the source temperature and errors from Kelvin to
* Jansky we use the effective beam area and relative error as presented
* in table 4 of Jarosik et al. (2011ApJS..192...14J).
*Note (3):
* i = A very weak source where we are unable to correct the flux for
* Eddington bias due to the low signal-to-noise level (see section 4.3).
* We use the internal template detection to estimate the unbiased flux in
* the given channel as explained in section 4.3.
*Note (4): One flag for each band that indicates whether we detected the
* source at 5{sigma} in the individual channel ("C"), or in the template
* with a 3{sigma} detection in the channel ("F") or not at all (".").
*Note (5): Flags indicating if the found source is identified in a given
* catalog in at least one frequency:
* W = WMAP-catalogs (Gold et al. 2011, Cat. J/ApJS/192/15);
* N = NEWPS 5yr 5s-catalog (Massardi et al. 2009, Cat. J/MNRAS/392/733);
* P = ERCSC-catalogs (at 30,44,70,100 GHz only; Planck Collaboration 2011,
* Cat. VIII/88).
*Note (6): Identification with the GB6 (Gregory et al. 1996, Cat. VIII/40),
* PMN (Griffith et al. 1993-1996, Cat. VIII/38);
* 1Jy (Kuehr et al. 1981, Cat. VIII/5) when possible or with the lower
* frequency catalogs NVSS (Condon et al. 1998, Cat. VIII/65) and SUMSS
* (Mauch et al. 2003, Cat. VIII/81) catalog counterparts.
*Note (7): A flag of "M" is used if there are multiple identifications possible
* (where we took the brightest one) and a flag "A" if the listed
* counterpart has a flux density below 100 mJy.
C=============================================================================
C Loading file 'table6.dat' ! New point sources in WMAP
C Format for file interpretation
1 format(
+ I4,1X,F7.3,1X,F7.3,1X,F5.2,1X,F4.2,1X,A1,1X,F5.2,1X,F4.2,1X,
+ A1,1X,F5.2,1X,F4.2,1X,A1,1X,F5.2,1X,F4.2,1X,A1,1X,F5.2,1X,
+ F4.2,1X,A1,1X,A5,1X,A3,1X,A19,1X,F4.1,A2)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,1116
read(1,'(A122)')ar__
read(ar__,1)
+ Seq(i__),RAdeg(i__),DEdeg(i__),FK(i__),e_FK(i__),f_FK(i__),
+ FKa(i__),e_FKa(i__),f_FKa(i__),FQ(i__),e_FQ(i__),f_FQ(i__),
+ FV(i__),e_FV(i__),f_FV(i__),FW(i__),e_FW(i__),f_FW(i__),
+ Fflag(i__),Cflag(i__),Cat_ID(i__),Off(i__),f_Off(i__)
if(ar__(22:26) .EQ. '') FK(i__) = rNULL__
if(ar__(28:31) .EQ. '') e_FK(i__) = rNULL__
if(ar__(35:39) .EQ. '') FKa(i__) = rNULL__
if(ar__(41:44) .EQ. '') e_FKa(i__) = rNULL__
if(ar__(48:52) .EQ. '') FQ(i__) = rNULL__
if(ar__(54:57) .EQ. '') e_FQ(i__) = rNULL__
if(ar__(61:65) .EQ. '') FV(i__) = rNULL__
if(ar__(67:70) .EQ. '') e_FV(i__) = rNULL__
if(ar__(74:78) .EQ. '') FW(i__) = rNULL__
if(ar__(80:83) .EQ. '') e_FW(i__) = rNULL__
if(ar__(117:120) .EQ. '') Off(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Seq(i__),RAdeg(i__),DEdeg(i__),FK(i__),e_FK(i__),f_FK(i__),
+ FKa(i__),e_FKa(i__),f_FKa(i__),FQ(i__),e_FQ(i__),f_FQ(i__),
+ FV(i__),e_FV(i__),f_FV(i__),FW(i__),e_FW(i__),f_FW(i__),
+ Fflag(i__),Cflag(i__),Cat_ID(i__),Off(i__),f_Off(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end