Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/671/1800 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/671/1800 SCUBA observations of {rho} Oph cloud (Andrews+, 2007)
*================================================================================
*A submillimeter view of circumstellar dust disks in {rho} Ophiuchi.
* Andrews S.M., Williams J.P.
* <Astrophys. J., 671, 1800-1812 (2007)>
* =2007ApJ...671.1800A
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Submillimeter properties of circumstellar disks
in {rho} Ophiuchus
integer*4 nr__
parameter (nr__=147) ! Number of records
character*85 ar__ ! Full-size record
character*19 Name (nr__) ! Object name(s)
character*3 SED (nr__) ! SED classification type (1)
character*2 SpT (nr__) ! MK spectral type from the literature (see
* Section 4.2)
integer*4 F0_35 (nr__) ! (mJy) ? 350{mu}m (850GHz) flux density
integer*4 e_F0_35 (nr__) ! (mJy) ? Uncertainty in F0.35 (2)
character*1 l_F0_45 (nr__) ! Limit flag on F0.45 (3)
integer*4 F0_45 (nr__) ! (mJy) ? 450{mu}m (660GHz) flux density
integer*4 e_F0_45 (nr__) ! (mJy) ? Uncertainty in F0.45 (2)
character*1 l_F0_85 (nr__) ! Limit flag on F0.85 (3)
integer*4 F0_85 (nr__) ! (mJy) ? 850{mu}m (350GHz) flux density
integer*4 e_F0_85 (nr__) ! (mJy) ? Uncertainty in F0.85 (2)
character*1 l_F1_3 (nr__) ! Limit flag on F1.3 (3)
integer*4 F1_3 (nr__) ! (mJy) ? 1.3mm (230GHz) flux density
integer*4 e_F1_3 (nr__) ! (mJy) ? Uncertainty in F1.3 (2)
character*1 l_Mass (nr__) ! Limit flag on Mass (3)
real*4 Mass (nr__) ! (solMass) Disk mass (see section 3.1)
character*1 l_alpha (nr__) ! Limit flag on alpha (3)
real*4 alpha (nr__) ! ? Submillimeter continuum spectral index {alpha}
* (F({nu}){prop.to}{nu}^{alpha}^, see section 3.2)
real*4 e_alpha (nr__) ! ? Uncertainty in Sp+Index
character*3 n_Name (nr__) ! Individual source note(s) (4)
*Note (1): From the spectral index {alpha} defined in the 2-25{mu}m range,
* the SED classification is the following:
* FS = Flat Spectrum, (star + disk + envelope?), -0.3<{alpha}<0.3;
* I = Class I (star + disk + envelope), {alpha}<=-0.3;
* II = Class II (star + disk), 0.3<={alpha}<=1.6;
* III = Class III (star only), {alpha}>1.6
*Note (2): Quoted errors are the 1{sigma} rms noise levels and do not
* include systematic errors in the absolute flux calibration (~25% at
* 350 and 450{mu}m, ~10% at 850{mu}m, and ~20% at 1.3mm).
*Note (3): Upper limits are taken at the 3{sigma} confidence level.
*Note (4): Individual notes as follows:
* 1 = The 1.3mm flux densities are from Andre and Montmerle
* (1994ApJ...420..837A). The 1{sigma} rms uncertainties were
* typically assumed to be ~10mJy
* 2 = 800{mu}m flux densities from Jensen et al. (1996ApJ...458..312J)
* are listed F850um
* 3 = The 1.3mm flux densities are from Nurnberger et al.
* (1998A&A...330..549N)
* 4 = The 1.3mm flux densities are from Motte et al. (1998A&A...336..150M)
* 5 = The 1.3mm flux densities for WL 20 (S) and WSB 60 are from the
* interferometric survey of Andrews and Williams (2007ApJ...659..705A)
* 6 = The 1.3mm flux densities are from Stanke et al. 2006,
* Cat. <J/A+A/447/609>)
C=============================================================================
C Loading file 'table1.dat' ! Submillimeter properties of circumstellar disks
* in {rho} Ophiuchus
C Format for file interpretation
1 format(
+ A19,1X,A3,1X,A2,1X,I4,1X,I3,1X,A1,I4,1X,I3,1X,A1,I4,1X,I2,1X,
+ A1,I3,1X,I2,1X,A1,F6.4,1X,A1,F4.2,1X,F4.2,1X,A3)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,147
read(1,'(A85)')ar__
read(ar__,1)
+ Name(i__),SED(i__),SpT(i__),F0_35(i__),e_F0_35(i__),
+ l_F0_45(i__),F0_45(i__),e_F0_45(i__),l_F0_85(i__),F0_85(i__),
+ e_F0_85(i__),l_F1_3(i__),F1_3(i__),e_F1_3(i__),l_Mass(i__),
+ Mass(i__),l_alpha(i__),alpha(i__),e_alpha(i__),n_Name(i__)
if(ar__(28:31) .EQ. '') F0_35(i__) = iNULL__
if(ar__(33:35) .EQ. '') e_F0_35(i__) = iNULL__
if(ar__(38:41) .EQ. '') F0_45(i__) = iNULL__
if(ar__(43:45) .EQ. '') e_F0_45(i__) = iNULL__
if(ar__(48:51) .EQ. '') F0_85(i__) = iNULL__
if(ar__(53:54) .EQ. '') e_F0_85(i__) = iNULL__
if(ar__(57:59) .EQ. '') F1_3(i__) = iNULL__
if(ar__(61:62) .EQ. '') e_F1_3(i__) = iNULL__
if(ar__(73:76) .EQ. '') alpha(i__) = rNULL__
if(ar__(78:81) .EQ. '') e_alpha(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Name(i__),SED(i__),SpT(i__),F0_35(i__),e_F0_35(i__),
+ l_F0_45(i__),F0_45(i__),e_F0_45(i__),l_F0_85(i__),F0_85(i__),
+ e_F0_85(i__),l_F1_3(i__),F1_3(i__),e_F1_3(i__),l_Mass(i__),
+ Mass(i__),l_alpha(i__),alpha(i__),e_alpha(i__),n_Name(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end