Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/415/32 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-19
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/415/32 COLD GASS survey (Saintonge+, 2011)
*================================================================================
*COLD GASS, an IRAM legacy survey of molecular gas in massive galaxies.
*I. Relations between H_2_, H I, stellar content and structural properties.
* Saintonge A., Kauffmann G., Kramer C., Tacconi L.J., Buchbender C.,
* Catinella B., Fabello S., Gracia-Carpio J., Wang J., Cortese L., Fu J.,
* Genzel R., Giovanelli R., Guo Q., Haynes M.P., Heckman T.M., Krumholz M.R.,
* Lemonias J., Li C., Moran S., Rodriguez-Fernandez N., Schiminovich D.,
* Schuster K., Sievers A.
* <Mon. Not. R. Astron. Soc., 415, 32-60 (2011)>
* =2011MNRAS.415...32S
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'catalog.dat' ! *COLD GASS Second Data Release - 2 February 2012
integer*4 nr__
parameter (nr__=308) ! Number of records
character*143 ar__ ! Full-size record
integer*4 GASS ! [1977/57017] GASS Identification number
character*19 SDSS ! SDSS designation (JHHMMSS.ss+DDMMSS.s)
real*8 z ! SDSS redshift
real*4 logM_ ! ([Msun]) Stellar mass
real*4 logmu_ ! ([Msun/kpc2]) ?=- Stellar mass surface density
real*4 D25 ! (arcsec) ?=- g-band optical diameter
real*4 C ! ?=- Concentration index (C=R_90_/R_50_)
real*4 NUV_r ! (mag) ?=- NUV-r color index, corrected for
* Galactic extinction
real*4 rmag ! (mag) ?=- r-band model magnitude
real*4 sigma ! (mK) rms per 20km/s-wide channel in the IRAM
* spectrum
real*4 S_N ! S/N of the CO(1-0) line (1)
real*4 SCO ! (Jy.km/s) CO(1-0) line flux, observed
real*4 SCOc ! (Jy.km/s) CO(1-0) line flux, corrected
real*4 foff ! ?=0 Flux ratio between offset and central
* pointings
real*4 logMH2 ! ([Msun]) Total H2 gas mass, measured value, or
* upper limit if flag CO=2 (2)
real*4 logMH2_M_ ! ([-]) H2 to stellar mass ratio,
* or upper limit if FlagCO=2
integer*4 CO ! [1/2] CO emission line flag,
* 1: detection, 2: non-detection
real*4 logMHI ! ([Msun]) Total neutral H mass (measured value,
* or upper limit if HI=4)
integer*4 HI ! [0/4] HI emission line flag (3)
real*4 Weight ! Weights (4)
*Note (1): use with caution the MH2 values when S/N<5
*Note (2): computed with alpha_CO_=3.2M_{sun}_/(K.(km/s).pc^2^), which does NOT
* include the contribution of He.
*Note (3): Flags as follows:
* 0 = no public data
* 1 = detection,, values of MHI from ALFALFA (Giovanelli et al. 2007,
* Cat. J/AJ/133/2569; Saintonge et al. 2008, Cat. J/AJ/135/588;
* Kent et al. 2008, Cat. J/AJ/136/713; Martin et al. 2009,
* Cat. J/ApJS/183/214; Stierwalt et al. 2009, Cat. J/AJ/138/338;
* Haynes et al. 2011AJ....142..170H)
* 2 = detection, values of MHI from Cornell HI archive
* (Springob et al., 2005, Cat. VIII/77)
* 3 = detection, values of MHI from GASS DR1 (Catinella et al., 2010,
* Cat. J/MNRAS/403/683)
* 4 = non-detection, values of MHI from GASS DR1 (Catinella et al.,
* 2010, Cat. J/MNRAS/403/683)
* The ALFALFA & HI archive HI masses have been recomputed in the same
* system as the GASS HI masses.
*Note (4): Weights that should be used when measuring sample-wide averages or
* fitting relations to the sample, to correct for the flat stellar mass
* distribution of the COLD GASS sample compared to the complete volume-limited
* parent sample (see paper for a more complete description).
C=============================================================================
C Loading file 'catalog.dat' ! *COLD GASS Second Data Release - 2 February 2012
C Format for file interpretation
1 format(
+ I5,3X,A19,2X,F7.5,2X,F5.2,3X,F4.2,2X,F4.1,2X,F4.2,2X,F4.2,2X,
+ F5.2,2X,F5.3,3X,F4.1,2X,F5.2,2X,F5.2,2X,F4.2,3X,F4.2,2X,F5.2,
+ 2X,I1,2X,F5.2,2X,I1,2X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'catalog.dat')
write(6,*) '....Loading file: catalog.dat'
do i__=1,308
read(1,'(A143)')ar__
read(ar__,1)
+ GASS,SDSS,z,logM_,logmu_,D25,C,NUV_r,rmag,sigma,S_N,SCO,SCOc,
+ foff,logMH2,logMH2_M_,CO,logMHI,HI,Weight
if (idig(ar__(47:50)).EQ.0) logmu_ = rNULL__
if (idig(ar__(53:56)).EQ.0) D25 = rNULL__
if (idig(ar__(59:62)).EQ.0) C = rNULL__
if (idig(ar__(65:68)).EQ.0) NUV_r = rNULL__
if (idig(ar__(71:75)).EQ.0) rmag = rNULL__
c ..............Just test output...........
write(6,1)
+ GASS,SDSS,z,logM_,logmu_,D25,C,NUV_r,rmag,sigma,S_N,SCO,SCOc,
+ foff,logMH2,logMH2_M_,CO,logMHI,HI,Weight
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end
C Locate position of first digit in string; or return 0
integer function idig(c)
character*(*) c
character*1 c1
integer lc,i
lc=len(c)
idig=0
do i=1,lc
if(c(i:i).ne.' ') go to 1
end do
1 if(i.gt.lc) return
c1=c(i:i)
if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1
if(i.gt.lc) return
c1=c(i:i)
if(c1.ge.'0'.and.c1.le.'9') idig=i
return
end