Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_AS/142/85 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 2024-Apr-19
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+AS/142/85 Recombination coefficients for C II lines (Davey+, 2000)
*================================================================================
*Recombination coefficients for C II
* Davey A.R., Storey P.J., Kisielius R.
* <Astron. Astrophys. Suppl. Ser. 142, 85 (2000)>
* =2000A&AS..142...85D
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Energies of even parity C^+^ bound states relative
to the ionization threshold
integer*4 nr__
parameter (nr__=43) ! Number of records
character*41 ar__ ! Full-size record
character*8 State (nr__) ! Electronic configuration
integer*4 SP1 (nr__) ! 2S+1 of term
integer*4 LP1 (nr__) ! L of term
integer*4 PP1 (nr__) ! Parity (0 for even, 1 for odd) of term
real*8 Ecal (nr__) ! (Ry) Calculated energy
real*8 Eexp (nr__) ! (Ry) ? Experimental energy
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Energies of odd parity C^+^ bound states relative
to the ionization threshold
integer*4 nr__1
parameter (nr__1=29) ! Number of records
character*42 ar__1 ! Full-size record
character*8 State_1 (nr__1) ! Electronic configuration
integer*4 SP1_1 (nr__1) ! 2S+1 of term
integer*4 LP1_1 (nr__1) ! L of term
integer*4 PP1_1 (nr__1) ! Parity (0 for even, 1 for odd) of term
real*8 Ecal_1 (nr__1) ! (Ry) Calculated energy
real*8 Eexp_1 (nr__1) ! (Ry) ? Experimental energy
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Effective recombination coefficients [10^-14^cm^3^/s]
at electron density N_e_=10^4^[cm^-3^],
temperature range T_e_=3500-20000K
integer*4 nr__2
parameter (nr__2=164) ! Number of records
character*99 ar__2 ! Full-size record
integer*4 SP1_2 (nr__2) ! 2S+1 of upper term parent state
integer*4 LP1_2 (nr__2) ! L of upper term parent state
integer*4 PP1_2 (nr__2) ! Parity of upper term parent state (1)
integer*4 n1 (nr__2) ! Principal quantum number of outer electron
character*1 l1 (nr__2) ! Orbital quantum number of outer electron in
* spectroscopic notation
integer*4 S1 (nr__2) ! 2S+1 of upper term
integer*4 L1_1 (nr__2) ! L of upper term
integer*4 P1 (nr__2) ! Parity of upper term (1)
integer*4 SP2 (nr__2) ! 2S+1 of lower term parent state
integer*4 LP2 (nr__2) ! L of lower term parent state
integer*4 PP2 (nr__2) ! Parity of lower term parent state (1)
integer*4 n2 (nr__2) ! Principal quantum number of outer electron
character*1 l2 (nr__2) ! Orbital quantum number of outer electron in
* spectroscopic notation
integer*4 S2 (nr__2) ! 2S+1 of lower term
integer*4 L2_1 (nr__2) ! L of lower term
integer*4 P2 (nr__2) ! Parity of lower term (1)
character*1 Case (nr__2) ! Plasma state case
real*8 Wavel (nr__2) ! (nm) Transition wavelength
real*8 v3500 (nr__2) ! (10-14cm+3/s) Effective recombination coeff. at Te = 3500K
real*8 v5000 (nr__2) ! (10-14cm+3/s) Effective recombination coeff. at Te = 5000K
real*8 v7500 (nr__2) ! (10-14cm+3/s) Effective recombination coeff. at Te = 7500K
real*8 v10000 (nr__2) ! (10-14cm+3/s) Effective recombination coeff. at Te = 10000K
real*8 v12500 (nr__2) ! (10-14cm+3/s) Effective recombination coeff. at Te = 12500K
real*8 v15000 (nr__2) ! (10-14cm+3/s) Effective recombination coeff. at Te = 15000K
real*8 v20000 (nr__2) ! (10-14cm+3/s) Effective recombination coeff. at Te = 20000K
*Note (1): 0 for even, 1 for odd
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Effective recombination coefficients [10^-14^cm^3^/s]
at electron density N_e_=10^4^[cm^-3^],
temperature range T_e_=500-2500K
integer*4 nr__3
parameter (nr__3=231) ! Number of records
character*99 ar__3 ! Full-size record
integer*4 SP1_3 (nr__3) ! 2S+1 of upper term parent state
integer*4 LP1_3 (nr__3) ! L of upper term parent state
integer*4 PP1_3 (nr__3) ! Parity of upper term parent state (1)
integer*4 n1_1 (nr__3) ! Principal quantum number of outer electron
character*1 l1_2 (nr__3) ! Orbital quantum number of outer electron in
* spectroscopic notation
integer*4 S1_1 (nr__3) ! 2S+1 of upper term
integer*4 L1_3 (nr__3) ! L of upper term
integer*4 P1_1 (nr__3) ! Parity of upper term (1)
integer*4 SP2_1 (nr__3) ! 2S+1 of lower term parent state
integer*4 LP2_1 (nr__3) ! L of lower term parent state
integer*4 PP2_1 (nr__3) ! Parity of lower term parent state (1)
integer*4 n2_1 (nr__3) ! Principal quantum number of outer electron
character*1 l2_2 (nr__3) ! Orbital quantum number of outer electron in
* spectroscopic notation
integer*4 S2_1 (nr__3) ! 2S+1 of lower term
integer*4 L2_3 (nr__3) ! L of lower term
integer*4 P2_1 (nr__3) ! Parity of lower term (1)
character*1 Case_1 (nr__3) ! Plasma state case
real*8 Wavel_1 (nr__3) ! (nm) Transition wavelength in nm
real*8 v500 (nr__3) ! (10-14cm+3/s) Effective recombination Ceoff. at Te = 500K
real*8 v750 (nr__3) ! (10-14cm+3/s) Effective recombination Ceoff. at Te = 750K
real*8 v1000 (nr__3) ! (10-14cm+3/s) Effective recombination Ceoff. at Te = 1000K
real*8 v1250 (nr__3) ! (10-14cm+3/s) Effective recombination Ceoff. at Te = 1250K
real*8 v1500 (nr__3) ! (10-14cm+3/s) Effective recombination Ceoff. at Te = 1500K
real*8 v2000 (nr__3) ! (10-14cm+3/s) Effective recombination Ceoff. at Te = 2000K
real*8 v2500 (nr__3) ! (10-14cm+3/s) Effective recombination Ceoff. at Te = 2500K
*Note (1): 0 for even, 1 for odd
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Fitting coefficients and maximum fitting errors (%)
for effective recombination coefficients
Electron density N_e_=10^4^[cm^-3^],
temperature range T_e_=5000-20000K
integer*4 nr__4
parameter (nr__4=143) ! Number of records
character*95 ar__4 ! Full-size record
integer*4 SP1_4 (nr__4) ! 2S+1 of upper term parent state
integer*4 LP1_4 (nr__4) ! L of upper term parent state
integer*4 PP1_4 (nr__4) ! Parity of upper term parent state (1)
integer*4 n1_2 (nr__4) ! Principal quantum number of outer electron
character*1 l1_4 (nr__4) ! Orbital quantum number of outer electron
* in spectroscopic notation
integer*4 S1_2 (nr__4) ! 2S+1 of upper term
integer*4 L1_5 (nr__4) ! L of upper term
integer*4 P1_2 (nr__4) ! Parity of upper term (1)
integer*4 SP2_2 (nr__4) ! 2S+1 of lower term parent state
integer*4 LP2_2 (nr__4) ! L of lower term parent state
integer*4 PP2_2 (nr__4) ! Parity of lower term parent state (1)
integer*4 n2_2 (nr__4) ! Principal quantum number of outer electron
character*1 l2_4 (nr__4) ! Orbital quantum number of outer electron
* in spectroscopic notation
integer*4 S2_2 (nr__4) ! 2S+1 of lower term
integer*4 L2_5 (nr__4) ! L of lower term
integer*4 P2_2 (nr__4) ! Parity of lower term (1)
character*1 Case_2 (nr__4) ! Plasma state case
real*8 Wavel_2 (nr__4) ! (nm) Transition wavelength
real*8 a (nr__4) ! Fitting coefficient a
real*8 b (nr__4) ! Fitting coefficient b
real*8 c (nr__4) ! Fitting coefficient c
real*8 d (nr__4) ! Fitting coefficient d
real*8 f (nr__4) ! Fitting coefficient f
real*8 Err (nr__4) ! (%) Maximum fitting error (in percent)
*Note (1): 0 for even, 1 for odd
C=============================================================================
C Loading file 'table1.dat' ! Energies of even parity C^+^ bound states relative
* to the ionization threshold
C Format for file interpretation
1 format(A8,I3,I3,I3,3X,F9.6,3X,F9.6)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,43
read(1,'(A41)')ar__
read(ar__,1)
+ State(i__),SP1(i__),LP1(i__),PP1(i__),Ecal(i__),Eexp(i__)
if(ar__(33:41) .EQ. '') Eexp(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ State(i__),SP1(i__),LP1(i__),PP1(i__),Ecal(i__),Eexp(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Energies of odd parity C^+^ bound states relative
* to the ionization threshold
C Format for file interpretation
2 format(A8,I3,I3,I3,3X,F9.6,3X,F9.6)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,29
read(1,'(A42)')ar__1
read(ar__1,2)
+ State_1(i__),SP1_1(i__),LP1_1(i__),PP1_1(i__),Ecal_1(i__),
+ Eexp_1(i__)
if(ar__1(33:41) .EQ. '') Eexp_1(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ State_1(i__),SP1_1(i__),LP1_1(i__),PP1_1(i__),Ecal_1(i__),
+ Eexp_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Effective recombination coefficients [10^-14^cm^3^/s]
* at electron density N_e_=10^4^[cm^-3^],
* temperature range T_e_=3500-20000K
C Format for file interpretation
3 format(
+ 2X,I1,I1,I1,1X,I2,A1,2X,I1,I1,I1,4X,I1,I1,I1,1X,I2,A1,2X,I1,
+ I1,I1,3X,A1,F8.1,1X,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,164
read(1,'(A99)')ar__2
read(ar__2,3)
+ SP1_2(i__),LP1_2(i__),PP1_2(i__),n1(i__),l1(i__),S1(i__),
+ L1_1(i__),P1(i__),SP2(i__),LP2(i__),PP2(i__),n2(i__),l2(i__),
+ S2(i__),L2_1(i__),P2(i__),Case(i__),Wavel(i__),v3500(i__),
+ v5000(i__),v7500(i__),v10000(i__),v12500(i__),v15000(i__),
+ v20000(i__)
c ..............Just test output...........
write(6,3)
+ SP1_2(i__),LP1_2(i__),PP1_2(i__),n1(i__),l1(i__),S1(i__),
+ L1_1(i__),P1(i__),SP2(i__),LP2(i__),PP2(i__),n2(i__),l2(i__),
+ S2(i__),L2_1(i__),P2(i__),Case(i__),Wavel(i__),v3500(i__),
+ v5000(i__),v7500(i__),v10000(i__),v12500(i__),v15000(i__),
+ v20000(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Effective recombination coefficients [10^-14^cm^3^/s]
* at electron density N_e_=10^4^[cm^-3^],
* temperature range T_e_=500-2500K
C Format for file interpretation
4 format(
+ 2X,I1,I1,I1,1X,I2,A1,2X,I1,I1,I1,4X,I1,I1,I1,1X,I2,A1,2X,I1,
+ I1,I1,3X,A1,F8.1,1X,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,231
read(1,'(A99)')ar__3
read(ar__3,4)
+ SP1_3(i__),LP1_3(i__),PP1_3(i__),n1_1(i__),l1_2(i__),
+ S1_1(i__),L1_3(i__),P1_1(i__),SP2_1(i__),LP2_1(i__),
+ PP2_1(i__),n2_1(i__),l2_2(i__),S2_1(i__),L2_3(i__),P2_1(i__),
+ Case_1(i__),Wavel_1(i__),v500(i__),v750(i__),v1000(i__),
+ v1250(i__),v1500(i__),v2000(i__),v2500(i__)
c ..............Just test output...........
write(6,4)
+ SP1_3(i__),LP1_3(i__),PP1_3(i__),n1_1(i__),l1_2(i__),
+ S1_1(i__),L1_3(i__),P1_1(i__),SP2_1(i__),LP2_1(i__),
+ PP2_1(i__),n2_1(i__),l2_2(i__),S2_1(i__),L2_3(i__),P2_1(i__),
+ Case_1(i__),Wavel_1(i__),v500(i__),v750(i__),v1000(i__),
+ v1250(i__),v1500(i__),v2000(i__),v2500(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Fitting coefficients and maximum fitting errors (%)
* for effective recombination coefficients
* Electron density N_e_=10^4^[cm^-3^],
* temperature range T_e_=5000-20000K
C Format for file interpretation
5 format(
+ 2X,I1,I1,I1,1X,I2,A1,2X,I1,I1,I1,4X,I1,I1,I1,1X,I2,A1,2X,I1,
+ I1,I1,3X,A1,F8.1,1X,F9.4,F9.4,F9.4,F9.4,F9.4,F7.2)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,143
read(1,'(A95)')ar__4
read(ar__4,5)
+ SP1_4(i__),LP1_4(i__),PP1_4(i__),n1_2(i__),l1_4(i__),
+ S1_2(i__),L1_5(i__),P1_2(i__),SP2_2(i__),LP2_2(i__),
+ PP2_2(i__),n2_2(i__),l2_4(i__),S2_2(i__),L2_5(i__),P2_2(i__),
+ Case_2(i__),Wavel_2(i__),a(i__),b(i__),c(i__),d(i__),f(i__),
+ Err(i__)
c ..............Just test output...........
write(6,5)
+ SP1_4(i__),LP1_4(i__),PP1_4(i__),n1_2(i__),l1_4(i__),
+ S1_2(i__),L1_5(i__),P1_2(i__),SP2_2(i__),LP2_2(i__),
+ PP2_2(i__),n2_2(i__),l2_4(i__),S2_2(i__),L2_5(i__),P2_2(i__),
+ Case_2(i__),Wavel_2(i__),a(i__),b(i__),c(i__),d(i__),f(i__),
+ Err(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end