Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/485/2710 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 2026-May-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/MNRAS/485/2710 Radio jets associated with galactic outflows (Jarvis+, 2019)
*================================================================================
*Prevalence of radio jets associated with galactic outflows and feedback from
*quasars.
* Jarvis M.E., Harrison C.M., Thomson A.P., Circosta C., Mainieri V.,
* Alexander D.M., Edge A.C., Lansbury G.B., Molyneux S.J., Mullaney J.R.
* <Mon. Not. R. Astron. Soc., 485, 2710-2730 (2019)>
* =2019MNRAS.485.2710J (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Target list and basic properties
integer*4 nr__
parameter (nr__=10) ! Number of records
character*78 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*10 Name ! Object name (JHHMM+DDMM)
integer*4 RAh ! (h) Right ascension (J2000) (1)
integer*4 RAm ! (min) Right ascension (J2000) (1)
real*4 RAs ! (s) Right ascension (J2000) (1)
character*1 DE_ ! Declination sign (J2000) (1)
integer*4 DEd ! (deg) Declination (J2000) (1)
integer*4 DEm ! (arcmin) Declination (J2000) (1)
real*4 DEs ! (arcsec) Declination (J2000) (1)
real*4 z ! Systemic redshift from the GMOS data
real*4 logLOIII ! ([10-7W]) Total observed [OIII] {lambda}5007
* luminosity from Mullaney et al.
* (2013MNRAS.433..622M) (2)
integer*4 FWHM ! (km/s) Full width at half-maximum (FWHM) of the
* broad component of the [OIII] line fit from
* Mullaney et al. (2013MNRAS.433..622M)
real*4 S1_4 ! (mJy) 1.4GHz flux densities obtained from the
* FIRST survey (Becker et al.
* 1995ApJ...450..559B, Cat. VIII/92)
real*4 e_S1_4 ! (mJy) Error on S1.4 (3)
real*4 logL1_4 ! (W/Hz) Rest-frame radio luminosities using a
* spectral index of {alpha}=-0.7 and assuming
* S_{nu}_{propto}{nu}^{alpha}^ (4)
real*4 Theta ! Radio morphology parameter (5)
real*4 e_Theta ! Error on Theta
*Note (1): Coordinates taken from SDSS (DR7) (Adelman-McCarthy et al.
* 2009ApJS..182..543A, Cat. II/294)
*Note (2): Absorption corrections would increase the values by on average 0.6dex
* (with a maximum increase of 1.4dex)
*Note (3): S1.4 uncertainties are defined as 3x the RMS noise of the radio image
* at the source position
*Note (4): We note that a range of {alpha}=-0.2 to -1.5 introduces a spread of
* +/-0.1dex on the radio luminosity
*Note (5): Sources with Theta>1.06 are classified as extended in the 1.4GHz
* FIRST data (Harrison et al. 2014MNRAS.441.3306H)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Galaxy and AGN parameters derived from SED
fitting for the primary sample
integer*4 nr__1
parameter (nr__1=10) ! Number of records
character*95 ar__1 ! Full-size record
character*10 Name_1 ! Object name (JHHMM+DDMM)
character*1 f_Name ! Flag on name (1)
real*4 logLAGN ! ([10-7W]) Bolometric AGN luminosity
real*4 logM ! ([Msun]) ? Stellar mass
real*4 E_logM ! ([Msun]) ? Upper error on logM
real*4 e_logM_1 ! ([Msun]) ? Lower error on logM
character*1 f_logM ! Flag on logM (2)
real*4 logLIR ! ([10-7W]) Infrared luminosity from star formation in
* the range 8-1000{mu}m
real*4 E_logLIR ! ([10-7W]) Upper error on logLIR
real*4 e_logLIR_1 ! ([10-7W]) Lower error on logLIR
real*4 SFR ! (Msun/yr) Star formation rate
real*4 e_SFR ! (Msun/yr) Error on SFR
real*4 S1_4_1 ! (mJy) 1.4GHz flux predicted from star formation
* following the radio-LIR relation
* (Bell 2003ApJ...586..794B,
* Cat. J/ApJ/586/794)
real*4 e_S1_4_1 ! (mJy) Error on S1.4
real*4 perctSF ! ? Percentage of the FIRST luminosity
* accounted for by star formation in the
* radio excess sources
real*4 e_perctSF ! ? Error on perctSF
real*4 qIR ! qIR ('radio excess') parameter, where
* qIR=<1.8 denotes radio excess
real*4 E_qIR ! Upper error on qIR
real*4 e_qIR_1 ! Lower error on qIR
character*1 Excess ! [YPN] Flag to define if the target is radio
* excess (3)
*Note (1): Flag as follows:
* a = These sources do not have photometric measurements at wavelengths longer
* than 60{mu}m, with J1316+1753 having no photometry above 22{mu}m
*Note (2): Flag as follows:
* b = For this target, the AGN contribution is particularly high in the NIR
* regime and the estimate of the stellar mass is unconstrained, with an
* uncertainty larger than the parameter value itself. We therefore do not
* report a value of logMsun.
*Note (3): Excess as follows:
* Y = radio excess
* P = probably radio excess
* N = not radio excess
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Radio properties of the morphological structures
extracted from resolution matched e-MERLIN and
VLA images
integer*4 nr__2
parameter (nr__2=52) ! Number of records
character*107 ar__2 ! Full-size record
character*10 Name_2 ! Object name (JHHMM+DDMM)
real*4 SFIRST ! (mJy) Flux density from the FIRST survey
real*4 e_SFIRST ! (mJy) Error on SFIRST
character*8 Struc ! Name of structure
character*1 f_Struc ! [ab] Flag on structure (1)
character*22 Interp ! Interpretation of structure or largest linear
* size observed at that resolution in kpc
character*1 f_Interp ! Flag on Interp (2)
character*1 l_S1_5 ! Limit flag on S1.5
real*4 S1_5 ! (mJy) ? Flux density in at 1.5GHz (3)
real*4 e_S1_5 ! (mJy) ? Error on S1.5 (4)
character*1 l_S5_2 ! Limit flag on S5.2
real*4 S5_2 ! (mJy) Flux density in at 5.2GHz (3)
real*4 e_S5_2 ! (mJy) ? Error on S5.2 (4)
character*1 l_S7_2 ! Limit flag on S7.2
real*4 S7_2 ! (mJy) Flux density in at 7.2GHz (3)
real*4 e_S7_2 ! (mJy) ? Error on S7.2 (4)
real*4 alpha ! ? Spectral index (5)
real*4 e_alpha ! ? Error on alpha
*Note (1): Flag as follows:
* a = Although this component is not detected at 5{sigma} significance in any
* individual band, it is detected in all three at 3{sigma} resulting in a
* spectral index of -0.5
* b = Due to the fitting constraints needed to get this component to be well
* fit by a Gaussian in all the three images, the fluxes and spectral
* index for this source are unreliable
*Note (2): Flag as follows:
* c = We determine that HR:B is either a high S/N beam artefact or a variable
* component, which then determines if LR:A is a composite or not
*Note (3): 5{sigma} upper limits are given for non-detections
*Note (4): The errors given are a combination of 1{sigma} random errors and a
* 10 per cent systematic
*Note (5): Spectral index {alpha} defined as S_{nu}_{prop.to}{nu}^{alpha}^ and
* found by fitting a line through all detected points between 1.5 and
* 7.2GHz
C=============================================================================
C Loading file 'table1.dat' ! Target list and basic properties
C Format for file interpretation
1 format(
+ A10,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F6.4,1X,
+ F5.2,1X,I4,1X,F4.1,1X,F3.1,1X,F4.1,1X,F5.3,1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,10
read(1,'(A78)')ar__
read(ar__,1)
+ Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,z,logLOIII,FWHM,S1_4,e_S1_4,
+ logL1_4,Theta,e_Theta
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)
+ Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,z,logLOIII,FWHM,S1_4,e_S1_4,
+ logL1_4,Theta,e_Theta
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Galaxy and AGN parameters derived from SED
* fitting for the primary sample
C Format for file interpretation
2 format(
+ A10,1X,A1,1X,F4.1,1X,F5.2,1X,F4.2,1X,F4.2,1X,A1,1X,F5.2,1X,
+ F4.2,1X,F4.2,1X,F4.1,1X,F4.1,1X,F3.1,1X,F3.1,1X,F4.1,1X,F3.1,
+ 1X,F4.2,1X,F4.2,1X,F4.2,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,10
read(1,'(A95)')ar__1
read(ar__1,2)
+ Name_1,f_Name,logLAGN,logM,E_logM,e_logM_1,f_logM,logLIR,
+ E_logLIR,e_logLIR_1,SFR,e_SFR,S1_4_1,e_S1_4_1,perctSF,
+ e_perctSF,qIR,E_qIR,e_qIR_1,Excess
if(ar__1(19:23) .EQ. '') logM = rNULL__
if(ar__1(25:28) .EQ. '') E_logM = rNULL__
if(ar__1(30:33) .EQ. '') e_logM_1 = rNULL__
if(ar__1(71:74) .EQ. '') perctSF = rNULL__
if(ar__1(76:78) .EQ. '') e_perctSF = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1,f_Name,logLAGN,logM,E_logM,e_logM_1,f_logM,logLIR,
+ E_logLIR,e_logLIR_1,SFR,e_SFR,S1_4_1,e_S1_4_1,perctSF,
+ e_perctSF,qIR,E_qIR,e_qIR_1,Excess
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Radio properties of the morphological structures
* extracted from resolution matched e-MERLIN and
* VLA images
C Format for file interpretation
3 format(
+ A10,1X,F4.1,1X,F3.1,1X,A8,1X,A1,1X,A22,1X,A1,1X,A1,1X,F5.2,1X,
+ F4.1,1X,A1,1X,F5.2,1X,F4.2,1X,A1,1X,F5.2,1X,F4.2,1X,F6.3,1X,
+ F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,52
read(1,'(A107)')ar__2
read(ar__2,3)
+ Name_2,SFIRST,e_SFIRST,Struc,f_Struc,Interp,f_Interp,l_S1_5,
+ S1_5,e_S1_5,l_S5_2,S5_2,e_S5_2,l_S7_2,S7_2,e_S7_2,alpha,
+ e_alpha
if(ar__2(59:63) .EQ. '') S1_5 = rNULL__
if(ar__2(65:68) .EQ. '') e_S1_5 = rNULL__
if(ar__2(78:81) .EQ. '') e_S5_2 = rNULL__
if(ar__2(91:94) .EQ. '') e_S7_2 = rNULL__
if(ar__2(96:101) .EQ. '') alpha = rNULL__
if(ar__2(103:107) .EQ. '') e_alpha = rNULL__
c ..............Just test output...........
write(6,3)
+ Name_2,SFIRST,e_SFIRST,Struc,f_Struc,Interp,f_Interp,l_S1_5,
+ S1_5,e_S1_5,l_S5_2,S5_2,e_S5_2,l_S7_2,S7_2,e_S7_2,alpha,
+ e_alpha
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end