Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/164/201 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-Oct-12
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/AJ/164/201 Accretion in 9 Orion OB1b T-Tauri stars (Pittman+, 2022)
*================================================================================
*Towards a Comprehensive View of Accretion, Inner Disks, and Extinction in
*Classical T Tauri Stars; An ODYSSEUS Study of the Orion OB1b Association.
* Pittman C.V., Espaillat C.C., Robinson C.E., Thanathibodee T., Calvet N.,
* Wendeborn J., Hernandez J., Manara C.F., Walter F., Abraham P., Alcala J.M.,
* Alencar S.H.P., Arulanantham N., Cabrit S., Eisloffel J., Fiorellino E.,
* France K., Gangi M., Grankin K., Herczeg G.J., Kospal A., Mendigutia I.,
* Serna J., Venuti L.
* <Astron. J., 164, 201 (2022)>
* =2022AJ....164..201P
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Best-fit accretion shock model parameters of
Orion OB1 CTTS
integer*4 nr__
parameter (nr__=9) ! Number of records
character*176 ar__ ! Full-size record
character*9 ID ! Object identifier
real*4 AV ! (mag) [0.05/1.44] V band extinction
real*4 E_AV ! (mag) [0.02/0.05] Upper uncertainty in AV
real*4 e_AV_1 ! (mag) [0.02/0.4] Lower uncertainty in AV
real*4 Rad ! (Rsun) [0.84/3.26] Stellar radius
real*4 E_Rad ! (Rsun) [0.07/0.4] Upper uncertainty in Rad
real*4 e_Rad_1 ! (Rsun) [0.07/0.4] Lower uncertainty in Rad
real*4 Lum ! (Lsun) [0.12/1.18] Stellar luminosity
real*4 e_Lum ! (Lsun) [0.03/0.3] Uncertainty in Lum
real*4 Lacc ! (Lsun) [0.06/0.61] Accretion luminosity
real*4 E_Lacc ! (Lsun) [0.00/0.04] Upper uncertainty in Lacc
real*4 e_Lacc_1 ! (Lsun) [0.00/0.23] Lower uncertainty in Lacc
real*4 dMass ! (10-8Msun/yr) [0.53/17.2] Accretion rate
real*4 E_dMass ! (10-8Msun/yr) [0.01/1] Upper uncertainty in dMass
real*4 e_dMass_1 ! (10-8Msun/yr) [0.01/7] Lower uncertainty in dMass
real*8 f1E10 ! [0.0001/0.24] 1e^10^ accretion column
* stellar surface fraction covered
real*4 E_f1E10 ! [0.0004/0.02] Upper uncertainty in f1E10
real*8 e_f1E10_1 ! [7e-05/0.03] Lower uncertainty in f1E10
real*8 f1E11 ! [9e-05/0.02] 1e^11^ accretion column
* stellar surface fraction covered
real*8 E_f1E11 ! [0.0001/0.03] Upper uncertainty in f1E11
real*8 e_f1E11_1 ! [3e-05/0.003] Lower uncertainty in f1E11
real*8 f1E12 ! [4e-05/0.03] 1e^12^ accretion column
* stellar surface fraction covered
real*8 E_f1E12 ! [7e-06/0.002] Upper uncertainty in f1E12
real*8 e_f1E12_1 ! [2e-06/0.01] Lower uncertainty in f1E12
integer*4 Tshock ! (K) [4800/10542] Shock temperature
integer*4 e_Tshock ! (K) [14/1910] Uncertainty in Tshock
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Periods and inclinations of Orion OB1 CTTS
integer*4 nr__1
parameter (nr__1=9) ! Number of records
character*40 ar__1 ! Full-size record
character*9 ID_1 ! Object identifier
real*4 vsini ! (km/s) [3.2/18.4] Projected rotational velocity
real*4 e_vsini ! (km/s) [0.8/2] Uncertainty in vsini
real*4 Per ! (d) [4.3/7.1] Period
real*4 Vel ! (km/s) [8.3/23.5] Stellar rotational velocity
real*4 e_Vel ! (km/s) [0.5/2] Uncertainty in Vel
real*4 i ! (deg) [9.3/52.4]? Inclination
real*4 e_i ! (deg) [2/8]? Uncertainty in i
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! Accretion column temperatures
integer*4 nr__2
parameter (nr__2=9) ! Number of records
character*42 ar__2 ! Full-size record
character*9 ID_2 ! Object identifier
integer*4 T1E10 ! (K) [3938/4675] 1e^10^ accretion column temperature
integer*4 e_T1E10 ! (K) [1/2] Uncertainty in T1E10
integer*4 T1E11 ! (K) [6310/6880] 1e^11^ accretion column temperature
integer*4 e_T1E11 ! (K) [1/3] Uncertainty in T1E11
integer*4 T1E12 ! (K) [10679/10973] 1e^12^ accretion column temperature
integer*4 e_T1E12 ! (K) [2/3] Uncertainty in T1E12
integer*4 Tshock_1 ! (K) [4800/10542] Weighted-average shock temperature
integer*4 e_Tshock_1 ! (K) [14/1910] Uncertainty in Tshock (1)
*Note (1): The error on Tshock comes from the propagated accretion column
* temperature and filling factor uncertainties. Since the filling
* factor uncertainties are not symmetric, the larger of the two
* uncertainties is used in the error propagation.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table8.dat' ! Accretion column mass flux rates
integer*4 nr__3
parameter (nr__3=9) ! Number of records
character*95 ar__3 ! Full-size record
character*9 ID_3 ! Object identifier
real*4 Mass_Rad ! [0.08/0.77] Ratio of stellar mass to radius
real*4 dM1E10 ! (10-8Msun/yr) [0.0007/1] 10^10^ accretion column model
* accretion rate
real*4 E_dM1E10 ! (10-8Msun/yr) [0.002/0.3] Upper uncertainty in Mdot1E10
real*4 e_dM1E10_1 ! (10-8Msun/yr) [0.0004/0.4] Lower uncertainty in Mdot1E10
real*4 dM1E11 ! (10-8Msun/yr) [0.001/2.1] 10^11^ accretion column model
* accretion rate
real*4 E_dM1E11 ! (10-8Msun/yr) [0.01/0.7] Upper uncertainty in Mdot1E11
real*4 e_dM1E11_1 ! (10-8Msun/yr) [0.001/0.6] Lower uncertainty in Mdot1E11
real*8 dM1E12 ! (10-8Msun/yr) [0.025/15.1] 10^12^ accretion column model
* accretion rate
real*4 E_dM1E12 ! (10-8Msun/yr) [0.004/2] Upper uncertainty in Mdot1E12
real*4 e_dM1E12_1 ! (10-8Msun/yr) [0.001/5] Lower uncertainty in Mdot1E12
real*4 dMTot ! (10-8Msun/yr) [0.529/17.2] Total accretion rate
real*4 E_dMTot ! (10-8Msun/yr) [0.018/1.0] Upper uncertainty in MdotTot
real*4 e_dMTot_1 ! (10-8Msun/yr) [0.016/4.0] Lower uncertainty in MdotTot
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'fig1.dat' ! Model parameters for the 9 sources
integer*4 nr__4
parameter (nr__4=9) ! Number of records
character*62 ar__4 ! Full-size record
character*9 ID_4 ! Object identifier
real*4 AV_1 ! (mag) [0.05/1.43] V band extinction
real*4 Rad_1 ! (Rsun) [0.84/3.26] Stellar radius
real*8 dMass_1 ! (10-8Msun/yr) [0.52/17.3] Accretion rate
real*8 f1E10_1 ! [4e-5/0.24] 1e^10^ accretion column stellar
* surface fraction covered
real*8 f1E11_1 ! [4e-5/0.02] 1e^11^ accretion column stellar
* surface fraction covered
real*8 f1E12_1 ! [4e-5/0.03] 1e^12^ accretion column stellar
* surface fraction covered
integer*4 Tshock_2 ! (K) [4802/10650] Weighted-average shock temperature
C=============================================================================
C Loading file 'table3.dat' ! Best-fit accretion shock model parameters of
* Orion OB1 CTTS
C Format for file interpretation
1 format(
+ A9,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,
+ F4.2,1X,F6.4,1X,F6.4,1X,F6.4,1X,F6.3,1X,F5.3,1X,F5.3,1X,F7.5,
+ 1X,F6.4,1X,F7.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,F9.7,1X,F8.6,1X,
+ F9.7,1X,I5,1X,I4)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,9
read(1,'(A176)')ar__
read(ar__,1)
+ ID,AV,E_AV,e_AV_1,Rad,E_Rad,e_Rad_1,Lum,e_Lum,Lacc,E_Lacc,
+ e_Lacc_1,dMass,E_dMass,e_dMass_1,f1E10,E_f1E10,e_f1E10_1,
+ f1E11,E_f1E11,e_f1E11_1,f1E12,E_f1E12,e_f1E12_1,Tshock,
+ e_Tshock
c ..............Just test output...........
write(6,1)
+ ID,AV,E_AV,e_AV_1,Rad,E_Rad,e_Rad_1,Lum,e_Lum,Lacc,E_Lacc,
+ e_Lacc_1,dMass,E_dMass,e_dMass_1,f1E10,E_f1E10,e_f1E10_1,
+ f1E11,E_f1E11,e_f1E11_1,f1E12,E_f1E12,e_f1E12_1,Tshock,
+ e_Tshock
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Periods and inclinations of Orion OB1 CTTS
C Format for file interpretation
2 format(
+ A9,1X,F4.1,1X,F3.1,1X,F3.1,1X,F4.1,1X,F3.1,1X,F4.1,1X,F3.1)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,9
read(1,'(A40)')ar__1
read(ar__1,2)ID_1,vsini,e_vsini,Per,Vel,e_Vel,i,e_i
if(ar__1(33:36) .EQ. '') i = rNULL__
if(ar__1(38:40) .EQ. '') e_i = rNULL__
c ..............Just test output...........
write(6,2)ID_1,vsini,e_vsini,Per,Vel,e_Vel,i,e_i
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7.dat' ! Accretion column temperatures
C Format for file interpretation
3 format(A9,1X,I4,1X,I1,1X,I4,1X,I1,1X,I5,1X,I1,1X,I5,1X,I4)
C Effective file loading
open(unit=1,status='old',file=
+'table7.dat')
write(6,*) '....Loading file: table7.dat'
do i__=1,9
read(1,'(A42)')ar__2
read(ar__2,3)
+ ID_2,T1E10,e_T1E10,T1E11,e_T1E11,T1E12,e_T1E12,Tshock_1,
+ e_Tshock_1
c ..............Just test output...........
write(6,3)
+ ID_2,T1E10,e_T1E10,T1E11,e_T1E11,T1E12,e_T1E12,Tshock_1,
+ e_Tshock_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table8.dat' ! Accretion column mass flux rates
C Format for file interpretation
4 format(
+ A9,1X,F4.2,1X,F6.4,1X,F6.4,1X,F6.4,1X,F6.4,1X,F5.3,1X,F6.4,1X,
+ F7.4,1X,F5.3,1X,F6.4,1X,F6.3,1X,F5.3,1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table8.dat')
write(6,*) '....Loading file: table8.dat'
do i__=1,9
read(1,'(A95)')ar__3
read(ar__3,4)
+ ID_3,Mass_Rad,dM1E10,E_dM1E10,e_dM1E10_1,dM1E11,E_dM1E11,
+ e_dM1E11_1,dM1E12,E_dM1E12,e_dM1E12_1,dMTot,E_dMTot,e_dMTot_1
c ..............Just test output...........
write(6,4)
+ ID_3,Mass_Rad,dM1E10,E_dM1E10,e_dM1E10_1,dM1E11,E_dM1E11,
+ e_dM1E11_1,dM1E12,E_dM1E12,e_dM1E12_1,dMTot,E_dMTot,e_dMTot_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'fig1.dat' ! Model parameters for the 9 sources
C Format for file interpretation
5 format(A9,1X,F4.2,1X,F4.2,1X,F9.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,I5)
C Effective file loading
open(unit=1,status='old',file=
+'fig1.dat')
write(6,*) '....Loading file: fig1.dat'
do i__=1,9
read(1,'(A62)')ar__4
read(ar__4,5)
+ ID_4,AV_1,Rad_1,dMass_1,f1E10_1,f1E11_1,f1E12_1,Tshock_2
c ..............Just test output...........
write(6,5)
+ ID_4,AV_1,Rad_1,dMass_1,f1E10_1,f1E11_1,f1E12_1,Tshock_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end