Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/257/21 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-Sep-17
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/ApJS/257/21 VERTICO: the Virgo EnviRonment Traced In CO survey (Brown+, 2021)
*================================================================================
*VERTICO: the Virgo EnviRonment Traced In CO survey.
* Brown T., Wilson C.D., Zabel N., Davis T.A., Boselli A., Chung A.,
* Ellison S.L., Lagos C.D.P., Stevens A.R.H., Cortese L., Bahe Y.M.,
* Bisaria D., Bolatto A.D., Cashmore C.R., Catinella B., Chown R., Diemer B.,
* Elahi P.J., Hani M.H., Jimenez-Donaire M.J., Lee B., Leidig K., Mok A.,
* Olsen K.P., Parker L.C., Roberts I.D., Smith R., Spekkens K., Thorp M.,
* Tonnesen S., Vienneau E., Villanueva V., Vogel S.N., Wadsley J., Welker C.,
* Yoon H.
* <Astrophys. J. Suppl. Ser., 257, 21 (2021)>
* =2021ApJS..257...21B
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! The VERTICO target sample
integer*4 nr__
parameter (nr__=51) ! Number of records
character*44 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*7 Name ! Galaxy name
character*3 f_Name ! Flag(s) on Name (1)
real*8 RAdeg ! (deg) [181/194] NED right ascension of the galaxy
* optical center (J2000)
real*8 DEdeg ! (deg) [2/19.4] NED declination of the galaxy optical
* center (J2000)
integer*4 Vel ! (km/s) [-228/2538] NED optical heliocentric recession
* velocity, {nu}_opt_
real*4 Inc ! (deg) [13.6/90] Optical r-band inclination, i (2)
real*4 PA ! (deg) [1.8/356] Optical r-band position angle of the
* kinematically redshifted half of the galaxy,
* calculated east of north (2)
*Note (1): Flag as follows:
* a = Data are from 7m and Total Power arrays. Other observations are 7m array
* only.
* b = Data are archival ALMA ACA (7m+total power) ^12^CO(2-1) observations at
* comparable sensitivity to our Cycle 7 observations. The sources of
* these data are the PHANGS-ALMA program (Leroy+ 2021ApJS..257...43L ;
* 14 galaxies) and one regular PI program (Cramer+ 2019ApJ...870...63C ;
* NGC 4402).
*
* See Section 3 and the Acknowledgments section for further details.
*Note (2): Inclination and PA calculated using fits to SDSS photometry described
* in Section 4.2.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Global CO(2-1) line properties for the
51 VERTICO galaxies fits
integer*4 nr__1
parameter (nr__1=51) ! Number of records
character*83 ar__1 ! Full-size record
character*7 Name_1 ! Galaxy name
character*1 n_Name ! Flag on Name (3)
real*4 S_N ! [1.6/232.6]? Signal-to-noise ratio of
* the integrated spectrum
real*4 Diamb ! (arcsec) [7.1/10.2] Diameter of the circularized
* synthesized beam, {theta}_b_
integer*4 rms ! (mJy) [21/680] Integrated spectrum rms
* in 10.6km/s channels
integer*4 Vlsr ! (km/s) [-222/2501]? Local Standard of Rest
* recession velocity of the line center (4)
integer*4 Delv ! (km/s) [92/621]? Velocity width of the line
* at zero intensity, {Delta}{nu}_line_
character*1 l_Sco ! 3{sigma} upper limit flag on Sco
real*8 Sco ! (Jy.km/s) [7.55/8369] Velocity-integrated CO flux
integer*4 e_Sco ! (Jy.km/s) [6/45]? Sco uncertainty (5)
character*1 l_logLCO ! 3{sigma} upper limit flag on logLCO
real*4 logLCO ! ([K.km/s.pc2]) [6/9.2] CO line luminosity given
* by Equation 5 (6)
real*4 e_logLCO ! ([K.km/s.pc2]) [0/0.27]? logLCO uncertainty
character*1 l_TCO ! 3{sigma} upper limit flag on TCO
real*4 TCO ! (K.km/s) [0.04/11.44] Mean velocity-integrated CO
* intensity
real*4 e_TCO ! (K.km/s) [0.02/0.11]? TCO uncertainty
character*1 l_logMmol ! 3{sigma} upper limit flag on logMmol
real*4 logMmol ! ([Msun]) [6.84/9.88] logarithm of molecular gas mass
* given by Equation 7 (7)
real*4 e_logMmol ! ([Msun]) [0/0.27]? logMmol uncertainty
*Note (3): Flag as follows:
* a = Values are 3{sigma} upper limits calculated over an on-sky circle
* with radius=30" (projected radius ~2.4kpc) and a line width of 100km/s.
* The positions and systemic velocities used are listed in Table 1.
*Note (4): Local standard of rest recession velocity of the line center,
* calculated as the midpoint of {Delta}{nu}_line_
*Note (5): We do not include the typical Band 6 calibration uncertainty of
* 5%-10% (0.02-0.04dex) in the flux measurement uncertainties. One may
* add this in quadrature to account for this.
*Note (6): Following Solomon & Vanden Bout (2005ARA&A..43..677S), we calculate
* the CO line luminosity, L_CO_, in K.km/s.pc^2^, expressed as the
* product of the velocity-integrated source brightness temperature and
* the source area, Equation (5):
* L_CO_=3.25x10^7^S_CO_{nu}_obs_^-2^D_L_^2^,
* where S_CO_ is the velocity-integrated flux in Jy.km/s, {nu}_obs_ is
* the observed frequency in GHz, and D_L_ is the luminosity distance to
* the source in Mpc. See Section 4.3.
*Note (7): We convert CO luminosity to molecular gas mass, including the
* contribution from heavy elements, in units of solar mass using the
* relation, Equation (7):
* M_mol_={alpha}_CO_/R_21_L_CO_,
* where {alpha}_CO_=4.35M_{sun}/pc^2^(K.km/s)^-1^, the molecular gas
* mass-to-CO (1-0) luminosity ratio calculated for the Milky Way disk
* by Bolatto+ (2013ARA&A..51..207B), and R_21_=0.8.
* See Section 4.4.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Molecular gas disk radii estimates fits
integer*4 nr__2
parameter (nr__2=49) ! Number of records
character*27 ar__2 ! Full-size record
character*7 Name_2 ! Galaxy name
real*4 R50 ! (kpc) [0.49/6.2] Radius containing 50% of the global
* CO(2-1) flux
real*4 R90 ! (kpc) [1.3/11.2] Radius containing 90% of the global
* CO(2-1) flux
real*4 Riso ! (kpc) [-0.08/9.8]? Inclination-corrected molecular
* gas isodensity radius at
* {Sigma}_mol_(r)=5M_{sun}_/pc^2^, at
* the distance of Virgo 1"=80pc
C=============================================================================
C Loading file 'table1.dat' ! The VERTICO target sample
C Format for file interpretation
1 format(A7,1X,A3,1X,F7.3,1X,F7.4,1X,I5,1X,F4.1,1X,F5.1)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,51
read(1,'(A44)')ar__
read(ar__,1)Name,f_Name,RAdeg,DEdeg,Vel,Inc,PA
c ..............Just test output...........
write(6,1)Name,f_Name,RAdeg,DEdeg,Vel,Inc,PA
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Global CO(2-1) line properties for the
* 51 VERTICO galaxies fits
C Format for file interpretation
2 format(
+ A7,1X,A1,1X,F5.1,1X,F4.1,1X,I3,1X,I4,1X,I3,1X,A1,1X,F7.2,1X,
+ I2,1X,A1,1X,F4.2,1X,F4.2,1X,A1,1X,F5.2,1X,F4.2,1X,A1,1X,F4.2,
+ 1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,51
read(1,'(A83)')ar__1
read(ar__1,2)
+ Name_1,n_Name,S_N,Diamb,rms,Vlsr,Delv,l_Sco,Sco,e_Sco,
+ l_logLCO,logLCO,e_logLCO,l_TCO,TCO,e_TCO,l_logMmol,logMmol,
+ e_logMmol
if(ar__1(11:15) .EQ. '') S_N = rNULL__
if(ar__1(26:29) .EQ. '') Vlsr = iNULL__
if(ar__1(31:33) .EQ. '') Delv = iNULL__
if(ar__1(45:46) .EQ. '') e_Sco = iNULL__
if(ar__1(55:58) .EQ. '') e_logLCO = rNULL__
if(ar__1(68:71) .EQ. '') e_TCO = rNULL__
if(ar__1(80:83) .EQ. '') e_logMmol = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1,n_Name,S_N,Diamb,rms,Vlsr,Delv,l_Sco,Sco,e_Sco,
+ l_logLCO,logLCO,e_logLCO,l_TCO,TCO,e_TCO,l_logMmol,logMmol,
+ e_logMmol
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Molecular gas disk radii estimates fits
C Format for file interpretation
3 format(A7,1X,F5.3,1X,F6.3,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,49
read(1,'(A27)')ar__2
read(ar__2,3)Name_2,R50,R90,Riso
if(ar__2(22:27) .EQ. '') Riso = rNULL__
c ..............Just test output...........
write(6,3)Name_2,R50,R90,Riso
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end