Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/691/672 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-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/ApJ/691/672 Spectroscopy of YSOs in Serpens molecular cloud (Oliveira+ 2009)
*================================================================================
*Optical characterization of a new young stellar population in the Serpens
*molecular cloud.
* Oliveira I., Merin B., Pontoppidan K.M., van Dishoeck E.F., Overzier R.A.,
* Hernandez J., Sicilia-Aguilar A., Eiroa C., Montesinos B.
* <Astrophys. J., 691, 672-686 (2009)>
* =2009ApJ...691..672O
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Spectral types for the YSO candidates in Serpens
integer*4 nr__
parameter (nr__=78) ! Number of records
character*55 ar__ ! Full-size record
integer*4 Seq (nr__) ! Running number (<[OMP2009] NNN> in Simbad)
integer*4 HMH07 (nr__) ! ? Identification from Harvey et al. 2007,
* Cat. J/ApJ/663/1149 (<[HMH2007] NNN> in Simbad)
character*16 SSTc2d (nr__) ! SSTc2d identification (HHMMSSss+DDMMSSs)
* (<SSTc2d JHHMMSS.s+DDMMSS> in Simbad)
real*4 alpha (nr__) ! {alpha}_2-24um_ slope (1)
character*4 SpT (nr__) ! MK spectral type
real*4 e_SpT (nr__) ! Subtype uncertainty
real*4 AVobs (nr__) ! (mag) Extinction found in this work
real*4 AVCl (nr__) ! (mag) Extinction found by "c2d" over 5' region
real*4 AV (nr__) ! (mag) ? Final extinction
*Note (1): Obtained from a linear fit to the logarithm of the fluxes between
* the 2MASS K (2um) and the MIPS1 (24um) bands (Harvey et al. 2007, Cat.
* J/ApJ/663/1149).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Other stellar parameters for the YSO candidates
integer*4 nr__1
parameter (nr__1=78) ! Number of records
character*111 ar__1 ! Full-size record
integer*4 Seq_1 (nr__1) ! Running number (<[OMP2009] NNN> in Simbad)
character*3 f_Seq (nr__1) ! [a,bc ] Flag on Seq (1)
integer*4 Teff (nr__1) ! (K) Effective temperature
integer*4 E_Teff (nr__1) ! (K) Positive error on Teff
integer*4 e_Teff_1 (nr__1) ! (K) Negative error on Teff
real*4 L_ (nr__1) ! (Lsun) ? Luminosity
real*4 E_L_ (nr__1) ! (Lsun) ? Positive error on L*
real*4 e_L__1 (nr__1) ! (Lsun) ? Negative error on L*
real*4 AgeB (nr__1) ! (Myr) ? Age (2)
real*4 E_AgeB (nr__1) ! (Myr) ? Positive error on AgeB
real*4 e_AgeB_1 (nr__1) ! (Myr) ? Negative error on AgeB
real*4 MassB (nr__1) ! (Msun) ? Mass (2)
real*4 E_MassB (nr__1) ! (Msun) ? Positive error on MassB
real*4 e_MassB_1 (nr__1) ! (Msun) ? Negative error on MassB
real*4 AgeS (nr__1) ! (Myr) ? Age (3)
real*4 E_AgeS (nr__1) ! (Myr) ? Positive error on AgeS
real*4 e_AgeS_1 (nr__1) ! (Myr) ? Negative error on AgeS
real*4 MassS (nr__1) ! (Msun) ? Mass (3)
real*4 E_MassS (nr__1) ! (Msun) ? Positive error on MassS
real*4 e_MassS_1 (nr__1) ! (Msun) ? Negative error on MassS
*Note (1): Flag as follows:
* a = background objects
* b = Post-AGB stars, confirmed with Spitzer IRS spectra
* c = Objects without 2MASS J, H, and K magnitudes
*Note (2): Values derived from the Baraffe et al. models
* (1998, Cat. J/A+A/337/403)
*Note (3): Values derived from the Siess et al. models (2000A&A...358..593S)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! H{alpha} equivalent widths and line widths and mass
accretion rates for YSOs in Serpens
integer*4 nr__2
parameter (nr__2=58) ! Number of records
character*33 ar__2 ! Full-size record
integer*4 Seq_2 (nr__2) ! Running number (<[OMP2009] NNN> in Simbad)
real*4 EWHa (nr__2) ! ? H{alpha} (6563{AA}) equivalent width
* (positive value = emission)
character*1 f_EWHa (nr__2) ! [t] t for classical T Tauri star (CTTS)
* based on EWHa
integer*4 Ha10 (nr__2) ! (km/s) ? H{alpha} full width at 10% of the peak
* intensity
integer*4 e_Ha10 (nr__2) ! (km/s) ? Ha10 uncertainty
character*1 f_Ha10 (nr__2) ! [t] t for CTTS based on Ha10
real*4 logdM_dt (nr__2) ! ([Msun/yr]) ? Mass accretion rate
real*4 e_logdM_dt (nr__2) ! ([Msun/yr]) ? logM uncertainty
C=============================================================================
C Loading file 'table2.dat' ! Spectral types for the YSO candidates in Serpens
C Format for file interpretation
1 format(
+ I3,1X,I3,1X,A16,1X,F5.2,1X,A4,1X,F3.1,1X,F4.1,1X,F5.2,1X,F4.1)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,78
read(1,'(A55)')ar__
read(ar__,1)
+ Seq(i__),HMH07(i__),SSTc2d(i__),alpha(i__),SpT(i__),
+ e_SpT(i__),AVobs(i__),AVCl(i__),AV(i__)
if(ar__(5:7) .EQ. '') HMH07(i__) = iNULL__
if(ar__(52:55) .EQ. '') AV(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Seq(i__),HMH07(i__),SSTc2d(i__),alpha(i__),SpT(i__),
+ e_SpT(i__),AVobs(i__),AVCl(i__),AV(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Other stellar parameters for the YSO candidates
C Format for file interpretation
2 format(
+ I3,A3,1X,I4,1X,I4,1X,I4,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,
+ F6.2,1X,F6.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F5.2,1X,F5.2,
+ 1X,F4.2,1X,F4.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,78
read(1,'(A111)')ar__1
read(ar__1,2)
+ Seq_1(i__),f_Seq(i__),Teff(i__),E_Teff(i__),e_Teff_1(i__),
+ L_(i__),E_L_(i__),e_L__1(i__),AgeB(i__),E_AgeB(i__),
+ e_AgeB_1(i__),MassB(i__),E_MassB(i__),e_MassB_1(i__),
+ AgeS(i__),E_AgeS(i__),e_AgeS_1(i__),MassS(i__),E_MassS(i__),
+ e_MassS_1(i__)
if(ar__1(23:28) .EQ. '') L_(i__) = rNULL__
if(ar__1(30:35) .EQ. '') E_L_(i__) = rNULL__
if(ar__1(37:42) .EQ. '') e_L__1(i__) = rNULL__
if(ar__1(44:49) .EQ. '') AgeB(i__) = rNULL__
if(ar__1(51:56) .EQ. '') E_AgeB(i__) = rNULL__
if(ar__1(58:63) .EQ. '') e_AgeB_1(i__) = rNULL__
if(ar__1(65:68) .EQ. '') MassB(i__) = rNULL__
if(ar__1(70:73) .EQ. '') E_MassB(i__) = rNULL__
if(ar__1(75:78) .EQ. '') e_MassB_1(i__) = rNULL__
if(ar__1(80:84) .EQ. '') AgeS(i__) = rNULL__
if(ar__1(86:90) .EQ. '') E_AgeS(i__) = rNULL__
if(ar__1(92:96) .EQ. '') e_AgeS_1(i__) = rNULL__
if(ar__1(98:101) .EQ. '') MassS(i__) = rNULL__
if(ar__1(103:106) .EQ. '') E_MassS(i__) = rNULL__
if(ar__1(108:111) .EQ. '') e_MassS_1(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Seq_1(i__),f_Seq(i__),Teff(i__),E_Teff(i__),e_Teff_1(i__),
+ L_(i__),E_L_(i__),e_L__1(i__),AgeB(i__),E_AgeB(i__),
+ e_AgeB_1(i__),MassB(i__),E_MassB(i__),e_MassB_1(i__),
+ AgeS(i__),E_AgeS(i__),e_AgeS_1(i__),MassS(i__),E_MassS(i__),
+ e_MassS_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! H{alpha} equivalent widths and line widths and mass
* accretion rates for YSOs in Serpens
C Format for file interpretation
3 format(I3,1X,F5.1,1X,A1,1X,I3,1X,I3,1X,A1,1X,F6.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,58
read(1,'(A33)')ar__2
read(ar__2,3)
+ Seq_2(i__),EWHa(i__),f_EWHa(i__),Ha10(i__),e_Ha10(i__),
+ f_Ha10(i__),logdM_dt(i__),e_logdM_dt(i__)
if(ar__2(5:9) .EQ. '') EWHa(i__) = rNULL__
if(ar__2(13:15) .EQ. '') Ha10(i__) = iNULL__
if(ar__2(17:19) .EQ. '') e_Ha10(i__) = iNULL__
if(ar__2(23:28) .EQ. '') logdM_dt(i__) = rNULL__
if(ar__2(30:33) .EQ. '') e_logdM_dt(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ Seq_2(i__),EWHa(i__),f_EWHa(i__),Ha10(i__),e_Ha10(i__),
+ f_Ha10(i__),logdM_dt(i__),e_logdM_dt(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end