Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/257/21 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-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 (nr__) ! Galaxy name
character*3 f_Name (nr__) ! Flag(s) on Name (1)
real*8 RAdeg (nr__) ! (deg) [181/194] NED right ascension of the galaxy
* optical center (J2000)
real*8 DEdeg (nr__) ! (deg) [2/19.4] NED declination of the galaxy optical
* center (J2000)
integer*4 Vel (nr__) ! (km/s) [-228/2538] NED optical heliocentric recession
* velocity, {nu}_opt_
real*4 Inc (nr__) ! (deg) [13.6/90] Optical r-band inclination, i (2)
real*4 PA (nr__) ! (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 (nr__1) ! Galaxy name
character*1 n_Name (nr__1) ! Flag on Name (3)
real*4 S_N (nr__1) ! [1.6/232.6]? Signal-to-noise ratio of
* the integrated spectrum
real*4 Diamb (nr__1) ! (arcsec) [7.1/10.2] Diameter of the circularized
* synthesized beam, {theta}_b_
integer*4 rms (nr__1) ! (mJy) [21/680] Integrated spectrum rms
* in 10.6km/s channels
integer*4 Vlsr (nr__1) ! (km/s) [-222/2501]? Local Standard of Rest
* recession velocity of the line center (4)
integer*4 Delv (nr__1) ! (km/s) [92/621]? Velocity width of the line
* at zero intensity, {Delta}{nu}_line_
character*1 l_Sco (nr__1) ! 3{sigma} upper limit flag on Sco
real*8 Sco (nr__1) ! (Jy.km/s) [7.55/8369] Velocity-integrated CO flux
integer*4 e_Sco (nr__1) ! (Jy.km/s) [6/45]? Sco uncertainty (5)
character*1 l_logLCO (nr__1) ! 3{sigma} upper limit flag on logLCO
real*4 logLCO (nr__1) ! ([K.km/s.pc2]) [6/9.2] CO line luminosity given
* by Equation 5 (6)
real*4 e_logLCO (nr__1) ! ([K.km/s.pc2]) [0/0.27]? logLCO uncertainty
character*1 l_TCO (nr__1) ! 3{sigma} upper limit flag on TCO
real*4 TCO (nr__1) ! (K.km/s) [0.04/11.44] Mean velocity-integrated CO
* intensity
real*4 e_TCO (nr__1) ! (K.km/s) [0.02/0.11]? TCO uncertainty
character*1 l_logMmol (nr__1) ! 3{sigma} upper limit flag on logMmol
real*4 logMmol (nr__1) ! ([Msun]) [6.84/9.88] logarithm of molecular gas mass
* given by Equation 7 (7)
real*4 e_logMmol (nr__1) ! ([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 (nr__2) ! Galaxy name
real*4 R50 (nr__2) ! (kpc) [0.49/6.2] Radius containing 50% of the global
* CO(2-1) flux
real*4 R90 (nr__2) ! (kpc) [1.3/11.2] Radius containing 90% of the global
* CO(2-1) flux
real*4 Riso (nr__2) ! (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(i__),f_Name(i__),RAdeg(i__),DEdeg(i__),Vel(i__),Inc(i__),
+ PA(i__)
c ..............Just test output...........
write(6,1)
+ Name(i__),f_Name(i__),RAdeg(i__),DEdeg(i__),Vel(i__),Inc(i__),
+ PA(i__)
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(i__),n_Name(i__),S_N(i__),Diamb(i__),rms(i__),
+ Vlsr(i__),Delv(i__),l_Sco(i__),Sco(i__),e_Sco(i__),
+ l_logLCO(i__),logLCO(i__),e_logLCO(i__),l_TCO(i__),TCO(i__),
+ e_TCO(i__),l_logMmol(i__),logMmol(i__),e_logMmol(i__)
if(ar__1(11:15) .EQ. '') S_N(i__) = rNULL__
if(ar__1(26:29) .EQ. '') Vlsr(i__) = iNULL__
if(ar__1(31:33) .EQ. '') Delv(i__) = iNULL__
if(ar__1(45:46) .EQ. '') e_Sco(i__) = iNULL__
if(ar__1(55:58) .EQ. '') e_logLCO(i__) = rNULL__
if(ar__1(68:71) .EQ. '') e_TCO(i__) = rNULL__
if(ar__1(80:83) .EQ. '') e_logMmol(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1(i__),n_Name(i__),S_N(i__),Diamb(i__),rms(i__),
+ Vlsr(i__),Delv(i__),l_Sco(i__),Sco(i__),e_Sco(i__),
+ l_logLCO(i__),logLCO(i__),e_logLCO(i__),l_TCO(i__),TCO(i__),
+ e_TCO(i__),l_logMmol(i__),logMmol(i__),e_logMmol(i__)
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(i__),R50(i__),R90(i__),Riso(i__)
if(ar__2(22:27) .EQ. '') Riso(i__) = rNULL__
c ..............Just test output...........
write(6,3)Name_2(i__),R50(i__),R90(i__),Riso(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end