Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/712/925 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 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/ApJ/712/925 Transition circumstellar disks in Ophiuchus (Cieza+, 2010)
*================================================================================
*The nature of transition circumstellar disks.
*I. The Ophiuchus molecular cloud.
* Cieza L.A., Schreiber M.R., Romero G.A., Mora M.D., Merin B., Swift J.J.,
* Orellana M., Williams J.P., Harvey P.M., Evans N.J.
* <Astrophys. J., 712, 925-941 (2010)>
* =2010ApJ...712..925C
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Transition disk sample
integer*4 nr__
parameter (nr__=34) ! Number of records
character*136 ar__ ! Full-size record
integer*4 Seq (nr__) ! [1/34] Running sequence number
character*16 SSTc2d (nr__) ! Object Name from "Cores to Disks" catalog
* (JHHMMSS.s+DDMMSS)
character*14 AltName (nr__) ! Alternative name
real*4 R1mag (nr__) ! (mag) First R mag from USNO-B1 catalog
real*4 R2mag (nr__) ! (mag) Second R mag from USNO-B1 catalog
real*4 Jmag (nr__) ! (mJy) 2MASS J band magnitude (1)
real*4 Hmag (nr__) ! (mJy) 2MASS H band magnitude (1)
real*4 Ksmag (nr__) ! (mJy) 2MASS Ks band magnitude (1)
real*4 F3_6 (nr__) ! (mJy) Spitzer/IRAC 3.6 micron band flux density (1)
real*4 F4_5 (nr__) ! (mJy) Spitzer/IRAC 4.5 micron band flux density (1)
real*4 F5_8 (nr__) ! (mJy) Spitzer/IRAC 5.8 micron band flux density (1)
real*4 F8_0 (nr__) ! (mJy) Spitzer/IRAC 8.0 micron band flux density (1)
real*4 F24 (nr__) ! (mJy) Spitzer/MIPS 24 micron band flux density (1)
character*1 l_F70 (nr__) ! The 5{sigma} limit flag on F70
real*4 F70 (nr__) ! (mJy) Spitzer/MIPS 70 micron band flux density (2)
*Note (1): Detections are >7{sigma} (i.e., the photometric uncertainties
* are <15%)
*Note (2): The 70 micron fluxes are either >5{sigma} detections or
* 5{sigma} upper limits.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Observed properties
integer*4 nr__1
parameter (nr__1=34) ! Number of records
character*88 ar__1 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 Seq_1 (nr__1) ! [1/34] Running sequence number
real*8 RAdeg (nr__1) ! (deg) Right Ascension in decimal degrees (J2000)
real*8 DEdeg (nr__1) ! (deg) Declination in decimal degrees (J2000)
character*7 Tel (nr__1) ! Telescope used to obtain spectrum (Clay, CFHT
* or Du Pont)
character*2 SpT (nr__1) ! MK spectral type
character*4 LiI (nr__1) ! (0.1nm) ? LiI equivalent width; in Angstroms (3)
character*3 CaII (nr__1) ! [Yes/No ]? Presence of CaII emission
integer*4 Halpha (nr__1) ! (km/s) ? H{alpha} velocity dispersion at 10% peak
* intensity
character*2 f_Halpha (nr__1) ! [-1/No] H{alpha} in emission (4)
real*4 mm (nr__1) ! (mm) The (sub)millimeter wavelength
character*1 l_Fmm (nr__1) ! ? Upper limit flag for Fmm (3{sigma})
real*4 Fmm (nr__1) ! (mJy) The (sub)millimeter flux (5)
real*4 e_Fmm (nr__1) ! (mJy) ? The 1{sigma} uncertainty for Fmm
real*4 Sep (nr__1) ! (arcsec) ? Companion separation (G1)
real*4 PA (nr__1) ! (deg) ? Position angle of companion
character*1 l_dK (nr__1) ! ? Upper limit flag for dk (5{sigma})
real*4 dK (nr__1) ! (mag) ? Flux ratio of companion in the Ks band
*Note (3): The absence of data implies that the signal to noise in this
* region of the spectrum is too low to measure the width or establish
* the presence of the LiI line.
*Note (4): A "-1" implies that H{alpha} is seen in absorption.
*Note (5): The 1.3mm data for source #14 and 17 comes from Andrews &
* Williams (2007, Cat. J/ApJ/671/1800), while the 1.3mm and 850um data
* for sources #12, 13 and 27 comes from Cieza et al. (2008, Cat.
* J/ApJ/686/L115).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Derived properties
integer*4 nr__2
parameter (nr__2=34) ! Number of records
character*84 ar__2 ! Full-size record
integer*4 Seq_2 (nr__2) ! [1/34] Running sequence number
character*1 l_accRate (nr__2) ! ? Upper limit flag for AccRate
real*4 accRate (nr__2) ! ([Msun/yr]) ? Accretion Rate (log scale)
character*1 l_Mdisk (nr__2) ! ? Upper limit flag for DiskMass
real*4 Mdisk (nr__2) ! (Mjup) ? Disk Mass (6)
real*4 rProj (nr__2) ! (AU) ? Projected separation (G1)
real*4 lam_to (nr__2) ! (um) {lambda}_turnoff_ parameter (7)
real*4 alphaE (nr__2) ! {alpha}_excess_ parameter (7)
character*1 l_log_Ld_L_(nr__2) ! ? Upper limit flag for log(Ld/L*)
character*6 log_Ld_L_ (nr__2) ! ([-]) Fractional disk luminosity (log scale)
real*4 AJ (nr__2) ! (mag) Extinction in the J-band
character*35 OType (nr__2) ! Object type
*Note (6): The disk mass upper limits for targets #12 and 27 come from
* Cieza et al. (2008, Cat. J/ApJ/686/L115).
*Note (7): the SED shape is quantified by the 2 parameters {lambda}_turnoff_
* (longest wavelength at which the observed flux is dominated by the
* stellar photosphere) and {alpha}_excess_, slope of the IR excess.
C=============================================================================
C Loading file 'table1.dat' ! Transition disk sample
C Format for file interpretation
1 format(
+ I2,8X,A16,1X,A14,1X,F5.2,1X,F5.2,1X,E8.2,1X,E8.2,1X,E8.2,1X,
+ E8.2,1X,E8.2,1X,E8.2,1X,E8.2,1X,E8.2,1X,A1,1X,E8.2)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,34
read(1,'(A136)')ar__
read(ar__,1)
+ Seq(i__),SSTc2d(i__),AltName(i__),R1mag(i__),R2mag(i__),
+ Jmag(i__),Hmag(i__),Ksmag(i__),F3_6(i__),F4_5(i__),F5_8(i__),
+ F8_0(i__),F24(i__),l_F70(i__),F70(i__)
c ..............Just test output...........
write(6,1)
+ Seq(i__),SSTc2d(i__),AltName(i__),R1mag(i__),R2mag(i__),
+ Jmag(i__),Hmag(i__),Ksmag(i__),F3_6(i__),F4_5(i__),F5_8(i__),
+ F8_0(i__),F24(i__),l_F70(i__),F70(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Observed properties
C Format for file interpretation
2 format(
+ I2,1X,F9.5,1X,F9.5,1X,A7,1X,A2,1X,A4,1X,A3,1X,I3,A2,1X,F4.2,
+ 1X,A1,1X,F6.2,1X,F5.2,1X,F5.3,1X,F6.2,1X,A1,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,34
read(1,'(A88)')ar__1
read(ar__1,2)
+ Seq_1(i__),RAdeg(i__),DEdeg(i__),Tel(i__),SpT(i__),LiI(i__),
+ CaII(i__),Halpha(i__),f_Halpha(i__),mm(i__),l_Fmm(i__),
+ Fmm(i__),e_Fmm(i__),Sep(i__),PA(i__),l_dK(i__),dK(i__)
if(ar__1(44:46) .EQ. '') Halpha(i__) = iNULL__
if(ar__1(64:68) .EQ. '') e_Fmm(i__) = rNULL__
if(ar__1(70:74) .EQ. '') Sep(i__) = rNULL__
if(ar__1(76:81) .EQ. '') PA(i__) = rNULL__
if(ar__1(85:88) .EQ. '') dK(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Seq_1(i__),RAdeg(i__),DEdeg(i__),Tel(i__),SpT(i__),LiI(i__),
+ CaII(i__),Halpha(i__),f_Halpha(i__),mm(i__),l_Fmm(i__),
+ Fmm(i__),e_Fmm(i__),Sep(i__),PA(i__),l_dK(i__),dK(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Derived properties
C Format for file interpretation
3 format(
+ I2,1X,A1,1X,F5.1,1X,A1,1X,F4.1,1X,F6.2,1X,F4.2,1X,F5.2,1X,A1,
+ 1X,A6,1X,F3.1,1X,A35)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,34
read(1,'(A84)')ar__2
read(ar__2,3)
+ Seq_2(i__),l_accRate(i__),accRate(i__),l_Mdisk(i__),
+ Mdisk(i__),rProj(i__),lam_to(i__),alphaE(i__),
+ l_log_Ld_L_(i__),log_Ld_L_(i__),AJ(i__),OType(i__)
if(ar__2(6:10) .EQ. '') accRate(i__) = rNULL__
if(ar__2(14:17) .EQ. '') Mdisk(i__) = rNULL__
if(ar__2(19:24) .EQ. '') rProj(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ Seq_2(i__),l_accRate(i__),accRate(i__),l_Mdisk(i__),
+ Mdisk(i__),rProj(i__),lam_to(i__),alphaE(i__),
+ l_log_Ld_L_(i__),log_Ld_L_(i__),AJ(i__),OType(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end