Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/135/1276 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-20
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/AJ/135/1276 ATLAS radio observations of ELAIS-S1 (Middelberg+, 2008)
*================================================================================
*Deep Australia Telescope Large Area Survey radio observations of the european
*large area ISO survey S1/Spitzer wide-area infrared extragalactic field.
* Middelberg E., Norris R.P., Cornwell T.J., Voronkov M.A., Siana B.D.,
* Boyle B.J., Ciliegi P., Jackson C.A., Huynh M.T., Berta S., Rubele S.,
* Lonsdale C.J., Ivison R.J., Smail I.
* <Astron. J., 135, 1276-1290 (2008)>
* =2008AJ....135.1276M
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Radio component data
integer*4 nr__
parameter (nr__=1366) ! Number of records
character*144 ar__ ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg ! (deg) Right Ascension J2000
real*8 DEdeg ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*7 CID ! Component number from paper ([MNC2008] CNNNN.N
* in Simbad) (1)
character*19 ATELAIS ! Component designation (JHHMMSS.ss+DDMMSS.s) (2)
integer*4 RAh ! (h) Hour of Right Ascension (J2000)
integer*4 RAm ! (min) Minute of Right Ascension (J2000)
real*4 RAs ! (s) Second of Right Ascension (J2000)
character*1 DE_ ! Sign of the Declination (J2000)
integer*4 DEd ! (deg) Degree of Declination (J2000)
integer*4 DEm ! (arcmin) Arcminute of Declination (J2000)
real*4 DEs ! (arcsec) Arcsecond of Declination (J2000)
real*4 e_RAs ! (arcsec) Error in RAs (3)
real*4 e_DEs ! (arcsec) Error in DEs (3)
real*4 Sp ! (mJy) Peak 20cm flux density
real*4 e_Sp ! (mJy) Error in Sp
real*4 Sint ! (mJy) Integrated 20cm flux density
real*4 e_Sint ! (mJy) Error in Sint
real*4 rms ! (uJy) SExtractor generated RMS at component position
real*4 Bmaj ! (arcsec) FWHM of major axis
real*4 Bmin ! (arcsec) FWHM of minor axis
integer*4 PA ! (deg) Position angle
real*8 Spd ! (mJy) ? Deconvolved peak 20cm flux density (4)
real*4 DBmaj ! (arcsec) ? Deconvolved FWHM of major axis (4)
real*4 DBmin ! (arcsec) ? Deconvolved FWHM of minor axis (4)
integer*4 DPA ! (deg) ? Deconvolved position angle (4)
character*1 L ! [*] Indicates component deemed a sidelobe
*Note (1): In some cases, sources were split up into sub-components,
* resulting in component numbers such as "C5" and "C5.1". However, this
* is no anticipation of the grouping of components to sources, which was
* carried out independently.
*Note (2): In the case of single-component sources, this is identical to the
* source name used in table 5. This is the formal IAU designation and
* should be used in the literature when referring to this component.
*Note (3): These include the formal uncertainties derived from the Gaussian
* fit together with a potential systematic error in the position of the
* calibrator source of 0.1 arcsec.
*Note (4): If a deconvolution was not possible, no value is given.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Radio source data
integer*4 nr__1
parameter (nr__1=1276) ! Number of records
character*491 ar__1 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_1 ! (deg) Right Ascension J2000
real*8 DEdeg_1 ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*7 Src ! Source number from paper (SNNNN.N,
* [MNC2008] SNNNN.N in Simbad)
character*19 ATELAIS_1 ! Source designation (JHHMMSS.ss+DDMMSS.s)
character*48 CID_1 ! Component(s) deemed to belong to source
* ([MNC2008] CNNNN.N in Simbad)
integer*4 RAh_1 ! (h) Hour of Right Ascension (J2000) (1)
integer*4 RAm_1 ! (min) Minute of Right Ascension (J2000) (1)
real*4 RAs_1 ! (s) Second of Right Ascension (J2000) (1)
character*1 DE__1 ! Sign of the Declination (J2000) (1)
integer*4 DEd_1 ! (deg) Degree of Declination (J2000) (1)
integer*4 DEm_1 ! (arcmin) Arcminute of Declination (J2000) (1)
real*4 DEs_1 ! (arcsec) Arcsecond of Declination (J2000) (1)
character*19 SWIRE ! SWIRE source name (JHHMMSS.ss+DDMMSS.s)
real*4 S20cm ! (mJy) Integrated 20cm flux density (2)
real*4 e_S20cm ! (mJy) Error in S20cm (2)
real*8 S3_6 ! (uJy) ? Spitzer/IRAC 3.6um band flux density (3)
real*8 S4_5 ! (uJy) ? Spitzer/IRAC 4.5um band flux density (3)
real*8 S5_8 ! (uJy) ? Spitzer/IRAC 5.8um band flux density (3)
real*8 S8_0 ! (uJy) ? Spitzer/IRAC 8.0um band flux density (3)
real*8 S24 ! (uJy) ? Spitzer/MIPS 24um band flux density (3)
real*4 Bmag ! (mag) ? The B band magnitude of Spitzer counterpart
real*4 Vmag ! (mag) ? The V band magnitude of Spitzer counterpart
real*4 Rmag ! (mag) ? The R band magnitude of Spitzer counterpart
character*2 AGN ! AGN classification code (4)
character*3 M ! [-/-x ] Magliocchetti et al.,
* 1998MNRAS.300..257M, test result (5)
real*4 z ! ? Redshift
character*3 r_z ! Reference for z (6)
character*204 Comm ! Additional comments
character*15 ELAISR ! Gruppioni et al. 1999, Cat. J/MNRAS/305/297
* designation (JHHMMSS+DDMMSSA)
character*30 OName ! Other name(s) from NED
*Note (1): In the case of single-component sources, this is the radio position
* of the component. In the case of multi-component sources with good
* infrared identification, the SWIRE position is used. In the case of
* multi-component sources without infrared identification, the
* coordinates are a flux-weighted mean of the components' coordinates.
*Note (2): In the case of extended or multiple-component sources, the flux
* density has been integrated over the source region, rather than
* taking the sum of its constituent components.
*Note (3): Aperture-corrected flux densities have been used unless the source
* was clearly extended, in which case the flux in a Kron aperture has
* been used.
*Note (4): Code as follows:
* f = based on the far-infrared-radio relation;
* m = based on morphology;
* l = based on classification taken from the literature.
*Note (5): As described in the text, performed for double radio sources.
* A "-" indicates failure, a "x" success of the two parts of the test
* (separation and flux density ratio of the constituents).
*Note (6): References as follows:
* 2df = Colless et al. (2001, Cat. <VII/226>, obsoleted by Cat. <VII/250>);
* 6dF = Jones et al. (2004, Cat. <VII/249>, obsoleted by <VII/259>);
* A01 = Alexander et al. (2001ApJ...554...18A);
* L04 = La Franca et al. (2004, Cat. <J/AJ/127/3075>);
* P06 = Puccetti et al. (2006, Cat. <J/A+A/457/501>);
* S01 = Serjeant et al. (2001MNRAS.322..262S);
* S96 = Shectman et al. (1996, Cat. <VII/203>);
* W03 = Wegner et al. (2003, Cat. <J/AJ/126/2268>).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! A section of the table with SUMSS (Sydney
University Molonglo Sky Survey, Cat. <VIII/81>)
counterparts to 1.4GHz radio sources
integer*4 nr__2
parameter (nr__2=73) ! Number of records
character*105 ar__2 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_2 ! (deg) Right Ascension J2000
real*8 DEdeg_2 ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*6 Src_1 ! Source designation (SNNNN.N)
integer*4 RAh_2 ! (h) SUMMS Hour of Right Ascension (J2000)
integer*4 RAm_2 ! (min) SUMMS Minute of Right Ascension (J2000)
real*4 RAs_2 ! (s) SUMMS Second of Right Ascension (J2000)
character*1 DE__2 ! Sign of the SUMMS Declination (J2000)
integer*4 DEd_2 ! (deg) SUMMS Degree of Declination (J2000)
integer*4 DEm_2 ! (arcmin) SUMMS Arcminute of Declination (J2000)
real*4 DEs_2 ! (arcsec) SUMMS Arcsecond of Declination (J2000)
real*4 S ! (mJy) SUMMS flux density at 843MHz
real*4 Sep ! (arcsec) Separation of SUMMS source to source
real*4 alpha ! Spectral index
character*55 Com ! Comment(s)
C=============================================================================
C Loading file 'table4.dat' ! Radio component data
C Format for file interpretation
1 format(
+ A7,9X,A19,1X,I2,1X,I2,1X,F6.3,1X,A1,I2,1X,I2,1X,F5.2,1X,F4.2,
+ 1X,F4.2,1X,F6.2,1X,F4.2,1X,F6.2,1X,F4.2,1X,F6.2,1X,F5.2,1X,
+ F5.2,1X,I3,1X,F7.2,1X,F5.2,1X,F5.2,1X,I3,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,1366
read(1,'(A144)')ar__
read(ar__,1)
+ CID,ATELAIS,RAh,RAm,RAs,DE_,DEd,DEm,DEs,e_RAs,e_DEs,Sp,e_Sp,
+ Sint,e_Sint,rms,Bmaj,Bmin,PA,Spd,DBmaj,DBmin,DPA,L
if(ar__(120:126) .EQ. '') Spd = rNULL__
if(ar__(128:132) .EQ. '') DBmaj = rNULL__
if(ar__(134:138) .EQ. '') DBmin = rNULL__
if(ar__(140:142) .EQ. '') DPA = iNULL__
RAdeg = rNULL__
DEdeg = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh .GT. -180) RAdeg=RAh*15.
if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
if(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
if(DEd .GE. 0) DEdeg=DEd
if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c ..............Just test output...........
write(6,1)
+ CID,ATELAIS,RAh,RAm,RAs,DE_,DEd,DEm,DEs,e_RAs,e_DEs,Sp,e_Sp,
+ Sint,e_Sint,rms,Bmaj,Bmin,PA,Spd,DBmaj,DBmin,DPA,L
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Radio source data
C Format for file interpretation
2 format(
+ A7,9X,A19,1X,A48,1X,I2,1X,I2,1X,F6.3,1X,A1,I2,1X,I2,1X,F5.2,
+ 8X,A19,2X,F6.2,1X,F5.2,1X,F8.2,1X,F8.2,1X,F8.2,1X,F8.2,1X,
+ F8.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,A2,1X,A3,1X,F4.2,1X,A3,1X,
+ A204,10X,A15,1X,A30)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,1276
read(1,'(A491)')ar__1
read(ar__1,2)
+ Src,ATELAIS_1,CID_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,
+ SWIRE,S20cm,e_S20cm,S3_6,S4_5,S5_8,S8_0,S24,Bmag,Vmag,Rmag,
+ AGN,M,z,r_z,Comm,ELAISR,OName
if(ar__1(153:160) .EQ. '') S3_6 = rNULL__
if(ar__1(162:169) .EQ. '') S4_5 = rNULL__
if(ar__1(171:178) .EQ. '') S5_8 = rNULL__
if(ar__1(180:187) .EQ. '') S8_0 = rNULL__
if(ar__1(189:196) .EQ. '') S24 = rNULL__
if(ar__1(198:202) .EQ. '') Bmag = rNULL__
if(ar__1(204:208) .EQ. '') Vmag = rNULL__
if(ar__1(210:214) .EQ. '') Rmag = rNULL__
if(ar__1(223:226) .EQ. '') z = rNULL__
RAdeg_1 = rNULL__
DEdeg_1 = 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 .GT. -180) RAdeg_1=RAh_1*15.
if(RAm_1 .GT. -180) RAdeg_1=RAdeg_1+RAm_1/4.
if(RAs_1 .GT. -180) RAdeg_1=RAdeg_1+RAs_1/240.
if(DEd_1 .GE. 0) DEdeg_1=DEd_1
if(DEm_1 .GE. 0) DEdeg_1=DEdeg_1+DEm_1/60.
if(DEs_1 .GE. 0) DEdeg_1=DEdeg_1+DEs_1/3600.
if(DE__1.EQ.'-'.AND.DEdeg_1.GE.0) DEdeg_1=-DEdeg_1
c ..............Just test output...........
write(6,2)
+ Src,ATELAIS_1,CID_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,
+ SWIRE,S20cm,e_S20cm,S3_6,S4_5,S5_8,S8_0,S24,Bmag,Vmag,Rmag,
+ AGN,M,z,r_z,Comm,ELAISR,OName
write(6,'(6H Pos: 2F8.4)') RAdeg_1,DEdeg_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! A section of the table with SUMSS (Sydney
* University Molonglo Sky Survey, Cat. <VIII/81>)
* counterparts to 1.4GHz radio sources
C Format for file interpretation
3 format(
+ A6,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F5.2,1X,F5.1,1X,F5.2,
+ 1X,F5.2,1X,A55)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,73
read(1,'(A105)')ar__2
read(ar__2,3)
+ Src_1,RAh_2,RAm_2,RAs_2,DE__2,DEd_2,DEm_2,DEs_2,S,Sep,alpha,
+ Com
RAdeg_2 = rNULL__
DEdeg_2 = 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 .GT. -180) RAdeg_2=RAh_2*15.
if(RAm_2 .GT. -180) RAdeg_2=RAdeg_2+RAm_2/4.
if(RAs_2 .GT. -180) RAdeg_2=RAdeg_2+RAs_2/240.
if(DEd_2 .GE. 0) DEdeg_2=DEd_2
if(DEm_2 .GE. 0) DEdeg_2=DEdeg_2+DEm_2/60.
if(DEs_2 .GE. 0) DEdeg_2=DEdeg_2+DEs_2/3600.
if(DE__2.EQ.'-'.AND.DEdeg_2.GE.0) DEdeg_2=-DEdeg_2
c ..............Just test output...........
write(6,3)
+ Src_1,RAh_2,RAm_2,RAs_2,DE__2,DEd_2,DEm_2,DEs_2,S,Sep,alpha,
+ Com
write(6,'(6H Pos: 2F8.4)') RAdeg_2,DEdeg_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end