Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/164/201 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-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 (nr__) ! Object identifier
real*4 AV (nr__) ! (mag) [0.05/1.44] V band extinction
real*4 E_AV (nr__) ! (mag) [0.02/0.05] Upper uncertainty in AV
real*4 e_AV_1 (nr__) ! (mag) [0.02/0.4] Lower uncertainty in AV
real*4 Rad (nr__) ! (Rsun) [0.84/3.26] Stellar radius
real*4 E_Rad (nr__) ! (Rsun) [0.07/0.4] Upper uncertainty in Rad
real*4 e_Rad_1 (nr__) ! (Rsun) [0.07/0.4] Lower uncertainty in Rad
real*4 Lum (nr__) ! (Lsun) [0.12/1.18] Stellar luminosity
real*4 e_Lum (nr__) ! (Lsun) [0.03/0.3] Uncertainty in Lum
real*4 Lacc (nr__) ! (Lsun) [0.06/0.61] Accretion luminosity
real*4 E_Lacc (nr__) ! (Lsun) [0.00/0.04] Upper uncertainty in Lacc
real*4 e_Lacc_1 (nr__) ! (Lsun) [0.00/0.23] Lower uncertainty in Lacc
real*4 dMass (nr__) ! (10-8Msun/yr) [0.53/17.2] Accretion rate
real*4 E_dMass (nr__) ! (10-8Msun/yr) [0.01/1] Upper uncertainty in dMass
real*4 e_dMass_1 (nr__) ! (10-8Msun/yr) [0.01/7] Lower uncertainty in dMass
real*8 f1E10 (nr__) ! [0.0001/0.24] 1e^10^ accretion column
* stellar surface fraction covered
real*4 E_f1E10 (nr__) ! [0.0004/0.02] Upper uncertainty in f1E10
real*8 e_f1E10_1 (nr__) ! [7e-05/0.03] Lower uncertainty in f1E10
real*8 f1E11 (nr__) ! [9e-05/0.02] 1e^11^ accretion column
* stellar surface fraction covered
real*8 E_f1E11 (nr__) ! [0.0001/0.03] Upper uncertainty in f1E11
real*8 e_f1E11_1 (nr__) ! [3e-05/0.003] Lower uncertainty in f1E11
real*8 f1E12 (nr__) ! [4e-05/0.03] 1e^12^ accretion column
* stellar surface fraction covered
real*8 E_f1E12 (nr__) ! [7e-06/0.002] Upper uncertainty in f1E12
real*8 e_f1E12_1 (nr__) ! [2e-06/0.01] Lower uncertainty in f1E12
integer*4 Tshock (nr__) ! (K) [4800/10542] Shock temperature
integer*4 e_Tshock (nr__) ! (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 (nr__1) ! Object identifier
real*4 vsini (nr__1) ! (km/s) [3.2/18.4] Projected rotational velocity
real*4 e_vsini (nr__1) ! (km/s) [0.8/2] Uncertainty in vsini
real*4 Per (nr__1) ! (d) [4.3/7.1] Period
real*4 Vel (nr__1) ! (km/s) [8.3/23.5] Stellar rotational velocity
real*4 e_Vel (nr__1) ! (km/s) [0.5/2] Uncertainty in Vel
real*4 i (nr__1) ! (deg) [9.3/52.4]? Inclination
real*4 e_i (nr__1) ! (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 (nr__2) ! Object identifier
integer*4 T1E10 (nr__2) ! (K) [3938/4675] 1e^10^ accretion column temperature
integer*4 e_T1E10 (nr__2) ! (K) [1/2] Uncertainty in T1E10
integer*4 T1E11 (nr__2) ! (K) [6310/6880] 1e^11^ accretion column temperature
integer*4 e_T1E11 (nr__2) ! (K) [1/3] Uncertainty in T1E11
integer*4 T1E12 (nr__2) ! (K) [10679/10973] 1e^12^ accretion column temperature
integer*4 e_T1E12 (nr__2) ! (K) [2/3] Uncertainty in T1E12
integer*4 Tshock_1 (nr__2) ! (K) [4800/10542] Weighted-average shock temperature
integer*4 e_Tshock_1 (nr__2) ! (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 (nr__3) ! Object identifier
real*4 Mass_Rad (nr__3) ! [0.08/0.77] Ratio of stellar mass to radius
real*4 dM1E10 (nr__3) ! (10-8Msun/yr) [0.0007/1] 10^10^ accretion column model
* accretion rate
real*4 E_dM1E10 (nr__3) ! (10-8Msun/yr) [0.002/0.3] Upper uncertainty in Mdot1E10
real*4 e_dM1E10_1 (nr__3) ! (10-8Msun/yr) [0.0004/0.4] Lower uncertainty in Mdot1E10
real*4 dM1E11 (nr__3) ! (10-8Msun/yr) [0.001/2.1] 10^11^ accretion column model
* accretion rate
real*4 E_dM1E11 (nr__3) ! (10-8Msun/yr) [0.01/0.7] Upper uncertainty in Mdot1E11
real*4 e_dM1E11_1 (nr__3) ! (10-8Msun/yr) [0.001/0.6] Lower uncertainty in Mdot1E11
real*8 dM1E12 (nr__3) ! (10-8Msun/yr) [0.025/15.1] 10^12^ accretion column model
* accretion rate
real*4 E_dM1E12 (nr__3) ! (10-8Msun/yr) [0.004/2] Upper uncertainty in Mdot1E12
real*4 e_dM1E12_1 (nr__3) ! (10-8Msun/yr) [0.001/5] Lower uncertainty in Mdot1E12
real*4 dMTot (nr__3) ! (10-8Msun/yr) [0.529/17.2] Total accretion rate
real*4 E_dMTot (nr__3) ! (10-8Msun/yr) [0.018/1.0] Upper uncertainty in MdotTot
real*4 e_dMTot_1 (nr__3) ! (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 (nr__4) ! Object identifier
real*4 AV_1 (nr__4) ! (mag) [0.05/1.43] V band extinction
real*4 Rad_1 (nr__4) ! (Rsun) [0.84/3.26] Stellar radius
real*8 dMass_1 (nr__4) ! (10-8Msun/yr) [0.52/17.3] Accretion rate
real*8 f1E10_1 (nr__4) ! [4e-5/0.24] 1e^10^ accretion column stellar
* surface fraction covered
real*8 f1E11_1 (nr__4) ! [4e-5/0.02] 1e^11^ accretion column stellar
* surface fraction covered
real*8 f1E12_1 (nr__4) ! [4e-5/0.03] 1e^12^ accretion column stellar
* surface fraction covered
integer*4 Tshock_2 (nr__4) ! (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(i__),AV(i__),E_AV(i__),e_AV_1(i__),Rad(i__),E_Rad(i__),
+ e_Rad_1(i__),Lum(i__),e_Lum(i__),Lacc(i__),E_Lacc(i__),
+ e_Lacc_1(i__),dMass(i__),E_dMass(i__),e_dMass_1(i__),
+ f1E10(i__),E_f1E10(i__),e_f1E10_1(i__),f1E11(i__),
+ E_f1E11(i__),e_f1E11_1(i__),f1E12(i__),E_f1E12(i__),
+ e_f1E12_1(i__),Tshock(i__),e_Tshock(i__)
c ..............Just test output...........
write(6,1)
+ ID(i__),AV(i__),E_AV(i__),e_AV_1(i__),Rad(i__),E_Rad(i__),
+ e_Rad_1(i__),Lum(i__),e_Lum(i__),Lacc(i__),E_Lacc(i__),
+ e_Lacc_1(i__),dMass(i__),E_dMass(i__),e_dMass_1(i__),
+ f1E10(i__),E_f1E10(i__),e_f1E10_1(i__),f1E11(i__),
+ E_f1E11(i__),e_f1E11_1(i__),f1E12(i__),E_f1E12(i__),
+ e_f1E12_1(i__),Tshock(i__),e_Tshock(i__)
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(i__),vsini(i__),e_vsini(i__),Per(i__),Vel(i__),
+ e_Vel(i__),i(i__),e_i(i__)
if(ar__1(33:36) .EQ. '') i(i__) = rNULL__
if(ar__1(38:40) .EQ. '') e_i(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ ID_1(i__),vsini(i__),e_vsini(i__),Per(i__),Vel(i__),
+ e_Vel(i__),i(i__),e_i(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(i__),T1E10(i__),e_T1E10(i__),T1E11(i__),e_T1E11(i__),
+ T1E12(i__),e_T1E12(i__),Tshock_1(i__),e_Tshock_1(i__)
c ..............Just test output...........
write(6,3)
+ ID_2(i__),T1E10(i__),e_T1E10(i__),T1E11(i__),e_T1E11(i__),
+ T1E12(i__),e_T1E12(i__),Tshock_1(i__),e_Tshock_1(i__)
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(i__),Mass_Rad(i__),dM1E10(i__),E_dM1E10(i__),
+ e_dM1E10_1(i__),dM1E11(i__),E_dM1E11(i__),e_dM1E11_1(i__),
+ dM1E12(i__),E_dM1E12(i__),e_dM1E12_1(i__),dMTot(i__),
+ E_dMTot(i__),e_dMTot_1(i__)
c ..............Just test output...........
write(6,4)
+ ID_3(i__),Mass_Rad(i__),dM1E10(i__),E_dM1E10(i__),
+ e_dM1E10_1(i__),dM1E11(i__),E_dM1E11(i__),e_dM1E11_1(i__),
+ dM1E12(i__),E_dM1E12(i__),e_dM1E12_1(i__),dMTot(i__),
+ E_dMTot(i__),e_dMTot_1(i__)
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(i__),AV_1(i__),Rad_1(i__),dMass_1(i__),f1E10_1(i__),
+ f1E11_1(i__),f1E12_1(i__),Tshock_2(i__)
c ..............Just test output...........
write(6,5)
+ ID_4(i__),AV_1(i__),Rad_1(i__),dMass_1(i__),f1E10_1(i__),
+ f1E11_1(i__),f1E12_1(i__),Tshock_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end