Conversion of standardized ReadMe file for
file /./ftp/cats/J/AN/347/70071 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-May-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/AN/347/70071 Circumstellar atomic radiation-driven dynamics (Lehtmets+, 2026)
*================================================================================
*Investigating Circumstellar Atomic Radiation-driven Dynamics.
* Lehtmets A., Kama M., Fossati L., Aret A.
* <Astron. Nachr. 347, e70071 (2026)>
* =2026AN....34770071L (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea1.dat' ! Chosen parameters for calculation of beta ratios
and velocity boosts
integer*4 nr__
parameter (nr__=9) ! Number of records
character*23 ar__ ! Full-size record
integer*4 Age ! (Myr) [20] Stellar age
real*4 Mass ! (Msun) [0.75/6.20] Stellar mass
real*4 Radius ! (Rsun) [0.576/2.370] Stellar radius
integer*4 Teff ! (K) [4000/19000] Stellar effective temperature
real*4 logg ! ([cm/s2]) [4.5] Stellar surface gravity in
* logarithmic scale
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea2.dat' ! Beta ratios of all investigated atoms and ions
across 9 different stellar effective temperatures
integer*4 nr__1
parameter (nr__1=76) ! Number of records
character*253 ar__1 ! Full-size record
character*6 Ion ! Ion name (e.g. H_I, He_II)
real*4 B1 ! Beta ratio value at 4000 K
real*4 ea_B1 ! Beta ratio absolute error at 4000 K
real*4 er_B1 ! (%) Beta ratio relative error (%) at 4000 K
real*4 B2 ! Beta ratio value at 4500 K
real*4 ea_B2 ! Beta ratio absolute error at 4500 K
real*4 er_B2 ! (%) Beta ratio relative error (%) at 4500 K
real*4 B3 ! Beta ratio value at 5000 K
real*4 ea_B3 ! Beta ratio absolute error at 5000 K
real*4 er_B3 ! (%) Beta ratio relative error at 5000 K
real*4 B4 ! Beta ratio value at 6000 K
real*4 ea_B4 ! Beta ratio absolute error at 6000 K
real*4 er_B4 ! (%) Beta ratio relative error (%) at 6000 K
real*4 B5 ! Beta ratio value at 8000 K
real*4 ea_B5 ! Beta ratio absolute error at 8000 K
real*4 er_B5 ! (%) Beta ratio relative error (%) at 8000 K
real*4 B6 ! Beta ratio value at 10000 K
real*4 ea_B6 ! Beta ratio absolute error at 10000 K
real*4 er_B6 ! (%) Beta ratio relative error (%) at 10000 K
real*4 B7 ! Beta ratio value at 13000 K
real*4 ea_B7 ! Beta ratio absolute error at 13000 K
real*4 er_B7 ! (%) Beta ratio relative error (%) at 13000 K
real*4 B8 ! Beta ratio value at 17000 K
real*4 ea_B8 ! Beta ratio absolute error at 17000 K
real*4 er_B8 ! (%) Beta ratio relative error (%) at 17000 K
real*4 B9 ! Beta ratio value at 19000 K
real*4 ea_B9 ! Beta ratio absolute error at 19000 K
real*4 er_B9 ! (%) Beta ratio relative error (%) at 19000 K
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea3.dat' ! Ionisation rates of neutral atoms
integer*4 nr__2
parameter (nr__2=28) ! Number of records
character*85 ar__2 ! Full-size record
character*4 Ion_1 ! Ion name (e.g. H_I, He_I)
real*4 I1 ! (10+10s-1) Ionisation rate value at 4000 K
real*4 I2 ! (10+10s-1) Ionisation rate value at 4500 K
real*4 I3 ! (10+10s-1) Ionisation rate value at 5000 K
real*4 I4 ! (10+10s-1) Ionisation rate value at 6000 K
real*4 I5 ! (10+10s-1) Ionisation rate value at 8000 K
real*4 I6 ! (10+10s-1) Ionisation rate value at 10000 K
real*4 I7 ! (10+10s-1) Ionisation rate value at 13000 K
real*4 I8 ! (10+10s-1) Ionisation rate value at 17000 K
real*4 I9 ! (10+10s-1) Ionisation rate value at 19000 K
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea4.dat' ! Velocity boosts of neutral atoms
integer*4 nr__3
parameter (nr__3=28) ! Number of records
character*121 ar__3 ! Full-size record
character*4 Ion_2 ! Ion name (e.g. H_I, He_I)
real*8 Vb1 ! (km/s) ?=- Velocity boosts value at 4000 K (1)
real*8 Vb2 ! (km/s) ?=- Velocity boosts value at 4500 K (1)
real*8 Vb3 ! (km/s) ?=- Velocity boosts value at 5000 K (1)
real*8 Vb4 ! (km/s) ?=- Velocity boosts value at 6000 K (1)
real*8 Vb5 ! (km/s) ?=- Velocity boosts value at 8000 K (1)
real*8 Vb6 ! (km/s) ?=- Velocity boosts value at 10000 K (1)
real*8 Vb7 ! (km/s) ?=- Velocity boosts value at 13000 K (1)
real*8 Vb8 ! (km/s) ?=- Velocity boosts value at 17000 K (1)
real*8 Vb9 ! (km/s) ?=- Velocity boosts value at 19000 K (1)
*Note (1): "---" means that the beta ratio for that species and at that
* stellar effective temperature is under 0.5.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea5.dat' ! Beta ratio values from this work with photospheric
emission and modified stellar EUV compared to
Fernandez et al., 2006ApJ...643...509K results
at 8000K stellar effective temperature
integer*4 nr__4
parameter (nr__4=81) ! Number of records
character*118 ar__4 ! Full-size record
character*6 Ion_3 ! Ion name (e.g. H_I, He_II)
real*4 BPE ! ?=- Beta ratio value from this work with
* photospheric emission (1)
real*4 e_BPE ! ?=- Beta ratio value from this work with
* photospheric emission error (1)
real*4 BMS ! ?=- Beta ratio value from this work with
* modified Stellar EUV (1)
real*4 e_BMS ! ?=- Beta ratio value from this work with
* modified Stellar EUV error (1)
real*8 BF06 ! ?=- Beta ratio value from Fernandez et al.,
* 2006ApJ...643...509K (1)
real*4 e_BF06 ! ?=- Beta ratio value from Fernandez et al.,
* 2006ApJ...643...509K error (1)
*Note (1): '0' means that the value is lower than 10^-30^,
* '---' means that there is no available linelist data.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea6.dat' ! Velocity boost values from this work with
photospheric emission and modified stellar EUV
compared to Fernandez et al., 2006ApJ...643...509K
results at 8000K stellar effective temperature
integer*4 nr__5
parameter (nr__5=19) ! Number of records
character*41 ar__5 ! Full-size record
character*4 Ion_4 ! Ion name (e.g. H_I, He_I)
real*8 VbPE ! (km/s) Velocity boost value from this work with
* photospheric emission
real*8 VbMS ! (km/s) Velocity boost value from this work with
* modified stellar EUV
real*8 VbF06 ! (km/s) Velocity boost value from
* Fernandez et al., 2006ApJ...643...509K
C=============================================================================
C Loading file 'tablea1.dat' ! Chosen parameters for calculation of beta ratios
* and velocity boosts
C Format for file interpretation
1 format(I2,1X,F4.2,1X,F5.3,1X,I5,1X,F3.1)
C Effective file loading
open(unit=1,status='old',file=
+'tablea1.dat')
write(6,*) '....Loading file: tablea1.dat'
do i__=1,9
read(1,'(A23)')ar__
read(ar__,1)Age,Mass,Radius,Teff,logg
c ..............Just test output...........
write(6,1)Age,Mass,Radius,Teff,logg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea2.dat' ! Beta ratios of all investigated atoms and ions
* across 9 different stellar effective temperatures
C Format for file interpretation
2 format(
+ A6,1X,E9.2,1X,E9.2,1X,E8.2,1X,E9.2,1X,E9.2,1X,E8.2,1X,E8.2,1X,
+ E8.2,1X,E8.2,1X,E8.2,1X,E8.2,1X,E8.2,1X,E8.2,1X,E8.2,1X,E8.2,
+ 1X,E8.2,1X,E8.2,1X,E8.2,1X,E8.2,1X,E8.2,1X,E8.2,1X,E8.2,1X,
+ E8.2,1X,E8.2,1X,E8.2,1X,E8.2,1X,E8.2)
C Effective file loading
open(unit=1,status='old',file=
+'tablea2.dat')
write(6,*) '....Loading file: tablea2.dat'
do i__=1,76
read(1,'(A253)')ar__1
read(ar__1,2)
+ Ion,B1,ea_B1,er_B1,B2,ea_B2,er_B2,B3,ea_B3,er_B3,B4,ea_B4,
+ er_B4,B5,ea_B5,er_B5,B6,ea_B6,er_B6,B7,ea_B7,er_B7,B8,ea_B8,
+ er_B8,B9,ea_B9,er_B9
c ..............Just test output...........
write(6,2)
+ Ion,B1,ea_B1,er_B1,B2,ea_B2,er_B2,B3,ea_B3,er_B3,B4,ea_B4,
+ er_B4,B5,ea_B5,er_B5,B6,ea_B6,er_B6,B7,ea_B7,er_B7,B8,ea_B8,
+ er_B8,B9,ea_B9,er_B9
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea3.dat' ! Ionisation rates of neutral atoms
C Format for file interpretation
3 format(
+ A4,1X,E8.2,1X,E8.2,1X,E8.2,1X,E8.2,1X,E8.2,1X,E8.2,1X,E8.2,1X,
+ E8.2,1X,E8.2)
C Effective file loading
open(unit=1,status='old',file=
+'tablea3.dat')
write(6,*) '....Loading file: tablea3.dat'
do i__=1,28
read(1,'(A85)')ar__2
read(ar__2,3)Ion_1,I1,I2,I3,I4,I5,I6,I7,I8,I9
c ..............Just test output...........
write(6,3)Ion_1,I1,I2,I3,I4,I5,I6,I7,I8,I9
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea4.dat' ! Velocity boosts of neutral atoms
C Format for file interpretation
4 format(
+ A4,1X,F10.6,1X,F11.7,F12.7,1X,F15.10,F14.9,1X,F13.8,1X,F13.8,
+ 1X,F12.7,F11.6)
C Effective file loading
open(unit=1,status='old',file=
+'tablea4.dat')
write(6,*) '....Loading file: tablea4.dat'
do i__=1,28
read(1,'(A121)')ar__3
read(ar__3,4)Ion_2,Vb1,Vb2,Vb3,Vb4,Vb5,Vb6,Vb7,Vb8,Vb9
if (idig(ar__3(6:15)).EQ.0) Vb1 = rNULL__
if (idig(ar__3(17:27)).EQ.0) Vb2 = rNULL__
if (idig(ar__3(28:39)).EQ.0) Vb3 = rNULL__
if (idig(ar__3(41:55)).EQ.0) Vb4 = rNULL__
if (idig(ar__3(56:69)).EQ.0) Vb5 = rNULL__
if (idig(ar__3(71:83)).EQ.0) Vb6 = rNULL__
if (idig(ar__3(85:97)).EQ.0) Vb7 = rNULL__
if (idig(ar__3(99:110)).EQ.0) Vb8 = rNULL__
if (idig(ar__3(111:121)).EQ.0) Vb9 = rNULL__
c ..............Just test output...........
write(6,4)Ion_2,Vb1,Vb2,Vb3,Vb4,Vb5,Vb6,Vb7,Vb8,Vb9
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea5.dat' ! Beta ratio values from this work with photospheric
* emission and modified stellar EUV compared to
* Fernandez et al., 2006ApJ...643...509K results
* at 8000K stellar effective temperature
C Format for file interpretation
5 format(
+ A6,1X,E20.16,1X,E20.15,1X,E20.15,1X,E20.15,1X,F15.10,1X,E11.7)
C Effective file loading
open(unit=1,status='old',file=
+'tablea5.dat')
write(6,*) '....Loading file: tablea5.dat'
do i__=1,81
read(1,'(A118)')ar__4
read(ar__4,5)Ion_3,BPE,e_BPE,BMS,e_BMS,BF06,e_BF06
if (idig(ar__4(8:27)).EQ.0) BPE = rNULL__
if (idig(ar__4(29:48)).EQ.0) e_BPE = rNULL__
if (idig(ar__4(50:69)).EQ.0) BMS = rNULL__
if (idig(ar__4(71:90)).EQ.0) e_BMS = rNULL__
if (idig(ar__4(92:106)).EQ.0) BF06 = rNULL__
if (idig(ar__4(108:118)).EQ.0) e_BF06 = rNULL__
c ..............Just test output...........
write(6,5)Ion_3,BPE,e_BPE,BMS,e_BMS,BF06,e_BF06
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea6.dat' ! Velocity boost values from this work with
* photospheric emission and modified stellar EUV
* compared to Fernandez et al., 2006ApJ...643...509K
* results at 8000K stellar effective temperature
C Format for file interpretation
6 format(A4,1X,F14.9,1X,F13.8,1X,F7.4)
C Effective file loading
open(unit=1,status='old',file=
+'tablea6.dat')
write(6,*) '....Loading file: tablea6.dat'
do i__=1,19
read(1,'(A41)')ar__5
read(ar__5,6)Ion_4,VbPE,VbMS,VbF06
c ..............Just test output...........
write(6,6)Ion_4,VbPE,VbMS,VbF06
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