Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/570/A82 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 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 (nr__) ! Object identifier (G1)
real*8 RAdeg (nr__) ! (deg) Right ascension (J2000)
real*8 DEdeg (nr__) ! (deg) Declination (J2000)
real*4 umag (nr__) ! (mag) u-band apparent magnitude
real*4 gmag (nr__) ! (mag) ?=- g-band apparent magnitude
real*4 rmag (nr__) ! (mag) r-band apparent magnitude
real*4 imag (nr__) ! (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 (nr__1) ! Object identifier (G1)
character*2 St (nr__1) ! CTTS (accreting) vs. WTTS (non-accreting)
* classification for members (1)
character*7 SpT (nr__1) ! Spectral type of the object
character*1 r_SpT (nr__1) ! [sp] spectroscopic or photometric type (2)
real*4 Av (nr__1) ! (mag) [0/7.1]?=- Visual extinction
real*8 Lbol (nr__1) ! (Lsun) Bolometric luminosity of the object
real*4 Mass (nr__1) ! (Msun) ?=- Stellar mass
real*4 Rad (nr__1) ! (Rsun) ?=- Stellar radius
real*4 logAge (nr__1) ! ([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 (nr__2) ! Object identifier (G1)
real*4 E_u_g (nr__2) ! (mag) UV excess measured from the u-g color index
real*4 E_u_r (nr__2) ! (mag) ?=- UV excess measured from the u-r
* color index (1)
real*4 logLEu (nr__2) ! ([Lsun]) ?=- Logarithm of the measured u-band excess
* luminosity (1)
real*4 logMa_g (nr__2) ! ([Msun/yr]) ?=- Log of the mass accretion rate derived
* from the (u-g) color excess
character*1 l_logMa_g (nr__2) ! [*] indicates an upper limit at the detection
* threshold (see text)
real*4 logMa_r (nr__2) ! ([Msun/yr]) ?=- Log of the mass accretion rate derived
* from the (u-r) color excess (1)
character*1 l_logMa_r (nr__2) ! [*] * indicates an upper limit at the detection
* threshold (see text)
real*4 var_ (nr__2) ! ([Msun/yr]) ?=- Positive (u-r) accretion variability bar
* log(max/med) (2)
real*4 var__1 (nr__2) ! ([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(i__),RAdeg(i__),DEdeg(i__),umag(i__),gmag(i__),rmag(i__),
+ imag(i__)
if (idig(ar__(38:43)).EQ.0) gmag(i__) = rNULL__
if (idig(ar__(52:57)).EQ.0) imag(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Mon(i__),RAdeg(i__),DEdeg(i__),umag(i__),gmag(i__),rmag(i__),
+ imag(i__)
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(i__),St(i__),SpT(i__),r_SpT(i__),Av(i__),Lbol(i__),
+ Mass(i__),Rad(i__),logAge(i__)
if (idig(ar__1(25:27)).EQ.0) Av(i__) = rNULL__
if (idig(ar__1(37:40)).EQ.0) Mass(i__) = rNULL__
if (idig(ar__1(42:46)).EQ.0) Rad(i__) = rNULL__
if (idig(ar__1(48:51)).EQ.0) logAge(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Mon_1(i__),St(i__),SpT(i__),r_SpT(i__),Av(i__),Lbol(i__),
+ Mass(i__),Rad(i__),logAge(i__)
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(i__),E_u_g(i__),E_u_r(i__),logLEu(i__),logMa_g(i__),
+ l_logMa_g(i__),logMa_r(i__),l_logMa_r(i__),var_(i__),
+ var__1(i__)
if (idig(ar__2(18:22)).EQ.0) E_u_r(i__) = rNULL__
if (idig(ar__2(24:28)).EQ.0) logLEu(i__) = rNULL__
if (idig(ar__2(30:35)).EQ.0) logMa_g(i__) = rNULL__
if (idig(ar__2(38:43)).EQ.0) logMa_r(i__) = rNULL__
if (idig(ar__2(46:49)).EQ.0) var_(i__) = rNULL__
if (idig(ar__2(51:54)).EQ.0) var__1(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ Mon_2(i__),E_u_g(i__),E_u_r(i__),logLEu(i__),logMa_g(i__),
+ l_logMa_g(i__),logMa_r(i__),l_logMa_r(i__),var_(i__),
+ var__1(i__)
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