Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/570/A82 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 2026-Sep-15
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/A+A/570/A82 Mapping accretion variability in NGC 2264 (Venuti+, 2014)
*================================================================================
*Mapping accretion and its variability in the young open cluster NGC 2264:
*a study based on u-band photometry.
* Venuti L., Bouvier J., Flaccomio E., Alencar S.H.P., Irwin J.,
* Stauffer J.R., Cody A.M., Teixeira P.S., Sousa A.P., Micela G.,
* Cuillandre J.-C., Peres G.
* <Astron. Astrophys. 570, A82 (2014)>
* =2014A&A...570A..82V
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Photometry for NGC 2264 members monitored at CFHT
integer*4 nr__
parameter (nr__=757) ! Number of records
character*57 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 Mon ! Object identifier (G1)
real*8 RAdeg ! (deg) Right ascension (J2000)
real*8 DEdeg ! (deg) Declination (J2000)
real*4 umag ! (mag) u-band apparent magnitude
real*4 gmag ! (mag) ?=- g-band apparent magnitude
real*4 rmag ! (mag) r-band apparent magnitude
real*4 imag ! (mag) ?=- i-band apparent magnitude
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Spectral type, extinction and stellar parameters
for NGC 2264 members monitored at CFHT
integer*4 nr__1
parameter (nr__1=757) ! Number of records
character*51 ar__1 ! Full-size record
integer*4 Mon_1 ! Object identifier (G1)
character*2 St ! CTTS (accreting) vs. WTTS (non-accreting)
* classification for members (1)
character*7 SpT ! Spectral type of the object
character*1 r_SpT ! [sp] spectroscopic or photometric type (2)
real*4 Av ! (mag) [0/7.1]?=- Visual extinction
real*8 Lbol ! (Lsun) Bolometric luminosity of the object
real*4 Mass ! (Msun) ?=- Stellar mass
real*4 Rad ! (Rsun) ?=- Stellar radius
real*4 logAge ! ([yr]) ?=- Logarithm of the age of the object
*Note (1): Classification code as follows:
* c = CTTS (classical T Tau star)
* w = WTTS (weak-line T Tau star)
* cc = CTTS candidate.
* c* = the object is a member that was previously classified as WTTS and
* is re-classified as CTTS in this study
*Note (2): spectral type estimate from:
* s = spectroscopic SpType estimate (retrieved from literature)
* p = photometric SpType estimate (from CFHT optical colors)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! UV excess, u-band excess luminosity, accretion
rate and accretion rate variability estimates
for NGC 2264 accreting members
integer*4 nr__2
parameter (nr__2=237) ! Number of records
character*54 ar__2 ! Full-size record
integer*4 Mon_2 ! Object identifier (G1)
real*4 E_u_g ! (mag) UV excess measured from the u-g color index
real*4 E_u_r ! (mag) ?=- UV excess measured from the u-r
* color index (1)
real*4 logLEu ! ([Lsun]) ?=- Logarithm of the measured u-band excess
* luminosity (1)
real*4 logMa_g ! ([Msun/yr]) ?=- Log of the mass accretion rate derived
* from the (u-g) color excess
character*1 l_logMa_g ! [*] indicates an upper limit at the detection
* threshold (see text)
real*4 logMa_r ! ([Msun/yr]) ?=- Log of the mass accretion rate derived
* from the (u-r) color excess (1)
character*1 l_logMa_r ! [*] * indicates an upper limit at the detection
* threshold (see text)
real*4 var_ ! ([Msun/yr]) ?=- Positive (u-r) accretion variability bar
* log(max/med) (2)
real*4 var__1 ! ([Msun/yr]) ?=- Negative (u-r) accretion variability bar
* log(med/min) (2)
*Note (1): Value measured at the median luminosity state of the system.
*Note (2): Variability range (from maximum to minimum value measured) detected
* around the median mass accretion rate level during the 2 week-long u,r
* monitoring. If no accretion activity is detected above the detection limit
* for the object at the minimum/median luminosity level, but some significant
* accretion is detected at the brightest state, only the upper variability bar
* is reported; if the object is a non-detection in accretion, no variability
* measurement is reported.
C=============================================================================
C Loading file 'table2.dat' ! Photometry for NGC 2264 members monitored at CFHT
C Format for file interpretation
1 format(4X,I6,1X,F9.5,1X,F8.5,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,757
read(1,'(A57)')ar__
read(ar__,1)Mon,RAdeg,DEdeg,umag,gmag,rmag,imag
if (idig(ar__(38:43)).EQ.0) gmag = rNULL__
if (idig(ar__(52:57)).EQ.0) imag = rNULL__
c ..............Just test output...........
write(6,1)Mon,RAdeg,DEdeg,umag,gmag,rmag,imag
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Spectral type, extinction and stellar parameters
* for NGC 2264 members monitored at CFHT
C Format for file interpretation
2 format(
+ 4X,I6,1X,A2,1X,A7,1X,A1,1X,F3.1,1X,F7.4,1X,F4.2,1X,F5.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,757
read(1,'(A51)')ar__1
read(ar__1,2)Mon_1,St,SpT,r_SpT,Av,Lbol,Mass,Rad,logAge
if (idig(ar__1(25:27)).EQ.0) Av = rNULL__
if (idig(ar__1(37:40)).EQ.0) Mass = rNULL__
if (idig(ar__1(42:46)).EQ.0) Rad = rNULL__
if (idig(ar__1(48:51)).EQ.0) logAge = rNULL__
c ..............Just test output...........
write(6,2)Mon_1,St,SpT,r_SpT,Av,Lbol,Mass,Rad,logAge
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! UV excess, u-band excess luminosity, accretion
* rate and accretion rate variability estimates
* for NGC 2264 accreting members
C Format for file interpretation
3 format(
+ 4X,I6,1X,F5.2,1X,F5.2,1X,F5.2,1X,F6.2,A1,1X,F6.2,A1,1X,F4.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,237
read(1,'(A54)')ar__2
read(ar__2,3)
+ Mon_2,E_u_g,E_u_r,logLEu,logMa_g,l_logMa_g,logMa_r,l_logMa_r,
+ var_,var__1
if (idig(ar__2(18:22)).EQ.0) E_u_r = rNULL__
if (idig(ar__2(24:28)).EQ.0) logLEu = rNULL__
if (idig(ar__2(30:35)).EQ.0) logMa_g = rNULL__
if (idig(ar__2(38:43)).EQ.0) logMa_r = rNULL__
if (idig(ar__2(46:49)).EQ.0) var_ = rNULL__
if (idig(ar__2(51:54)).EQ.0) var__1 = rNULL__
c ..............Just test output...........
write(6,3)
+ Mon_2,E_u_g,E_u_r,logLEu,logMa_g,l_logMa_g,logMa_r,l_logMa_r,
+ var_,var__1
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