Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/533/A87 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-20
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/533/A87 Sulphur ions iso-electronic sequence (Liang+, 2011)
*================================================================================
*R-matrix electron-impact excitation data for astrophysically abundant sulphur
*ions.
* Liang G.Y., Badnell N.R., Zhao G., Zhong J.Y., Wang F.L.
* <Astron. Astrophys. 533, A87 (2011)>
* =2011A&A...533A..87L
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'list.dat' ! List of files
integer*4 nr__
parameter (nr__=4) ! Number of records
character*50 ar__ ! Full-size record
character*12 FileName (nr__) ! Name of the file in subdirectory files (1)
character*36 Title (nr__) ! Title of the table
*Note (1): Each file deals with one ion, and includes the 204 fine-structure
* levels and 21809 transitions.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! *Level energies (Ryd) of S^8+^
integer*4 nr__1
parameter (nr__1=92) ! Number of records
character*60 ar__1 ! Full-size record
integer*4 ID (nr__1) ! Identification
character*22 Level (nr__1) ! Level specification
real*8 ENIST (nr__1) ! (Ry) ? NIST level energy (G1)
real*8 EAS (nr__1) ! (Ry) Present AUTOSTRUCTURE calculation level energy (G2)
real*8 EBL03 (nr__1) ! (Ry) ? Bhatia & Landi(2003ApJS..147..409B) level energy
real*8 EMCHF (nr__1) ! (Ry) ? MCHF collection level energy (G3)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! *Level energies (Ryd) of S^9+^
integer*4 nr__2
parameter (nr__2=84) ! Number of records
character*60 ar__2 ! Full-size record
integer*4 ID_1 (nr__2) ! Identification
character*22 Level_1 (nr__2) ! Level specification
real*8 ENIST_1 (nr__2) ! (Ry) ? NIST level energy (G1)
real*8 EAS_1 (nr__2) ! (Ry) Present AUTOSTRUCTURE calculation level energy (G2)
real*8 EBL03_1 (nr__2) ! (Ry) ? Bhatia & Landi(2003ApJS..147..409B) level energy
real*8 EMCHF_1 (nr__2) ! (Ry) ? MCHF collection level energy (G3)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! *Level energies (Ryd) of S^11+^
integer*4 nr__3
parameter (nr__3=204) ! Number of records
character*74 ar__3 ! Full-size record
integer*4 ID_2 (nr__3) ! Identification number
character*17 Level_2 (nr__3) ! Level specification
real*8 ENIST_2 (nr__3) ! (Ry) ? NIST level energy (G1)
real*8 EAS_2 (nr__3) ! (Ry) Present AUTOSTRUCTURE calculation level
* energy (G2)
real*8 ECHIANTI (nr__3) ! (Ry) ? CHIANTI v6 level energy
character*1 n_ECHIANTI (nr__3) ! [x] x for data are found to be observed, after
* checking the original paper (Zhang et al.,
* 1994, Cat. J/A+A/283/319))
real*8 ENSD07 (nr__3) ! (Ry) ? Nataraj et al. (2007, J. Phys. B: At. Mol.
* Opt. Phys., 40, 3153) level energy prediction
* with relativistic coupled-cluster theory
character*1 n_ENSD07 (nr__3) ! [y] y for data from the work of Merkelis et al.
* (1995, Scripta, 51, 233) by using many-body
* perturbation theory (MBPT)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! *Level energies (Ryd) of S^10+^
integer*4 nr__4
parameter (nr__4=254) ! Number of records
character*66 ar__4 ! Full-size record
integer*4 ID_3 (nr__4) ! Identification
character*15 Level_3 (nr__4) ! Level specification
real*8 Enist_3 (nr__4) ! (Ry) ? NIST level energy (G1)
real*8 EAS_3 (nr__4) ! (Ry) Present AUTOSTRUCTURE calculation level energy (G2)
real*8 ELB03 (nr__4) ! (Ry) ? Landi & Bhatia (2003ApJS..149..251L) level energy
real*8 EMCHF_2 (nr__4) ! (Ry) ? MCHF collection level energy (G3)
C=============================================================================
C Loading file 'list.dat' ! List of files
C Format for file interpretation
1 format(A12,2X,A36)
C Effective file loading
open(unit=1,status='old',file=
+'list.dat')
write(6,*) '....Loading file: list.dat'
do i__=1,4
read(1,'(A50)')ar__
read(ar__,1)FileName(i__),Title(i__)
c ..............Just test output...........
write(6,1)FileName(i__),Title(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! *Level energies (Ryd) of S^8+^
C Format for file interpretation
2 format(I3,1X,A22,1X,F9.6,1X,F7.4,1X,F7.4,1X,F7.4)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,92
read(1,'(A60)')ar__1
read(ar__1,2)
+ ID(i__),Level(i__),ENIST(i__),EAS(i__),EBL03(i__),EMCHF(i__)
if(ar__1(28:36) .EQ. '') ENIST(i__) = rNULL__
if(ar__1(46:52) .EQ. '') EBL03(i__) = rNULL__
if(ar__1(54:60) .EQ. '') EMCHF(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ ID(i__),Level(i__),ENIST(i__),EAS(i__),EBL03(i__),EMCHF(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! *Level energies (Ryd) of S^9+^
C Format for file interpretation
3 format(I3,1X,A22,1X,F9.6,1X,F7.4,1X,F7.4,1X,F7.4)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,84
read(1,'(A60)')ar__2
read(ar__2,3)
+ ID_1(i__),Level_1(i__),ENIST_1(i__),EAS_1(i__),EBL03_1(i__),
+ EMCHF_1(i__)
if(ar__2(28:36) .EQ. '') ENIST_1(i__) = rNULL__
if(ar__2(46:52) .EQ. '') EBL03_1(i__) = rNULL__
if(ar__2(54:60) .EQ. '') EMCHF_1(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ ID_1(i__),Level_1(i__),ENIST_1(i__),EAS_1(i__),EBL03_1(i__),
+ EMCHF_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! *Level energies (Ryd) of S^11+^
C Format for file interpretation
4 format(1X,I3,3X,A17,3X,F9.6,2X,F9.6,3X,F9.6,A1,4X,F9.6,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,204
read(1,'(A74)')ar__3
read(ar__3,4)
+ ID_2(i__),Level_2(i__),ENIST_2(i__),EAS_2(i__),ECHIANTI(i__),
+ n_ECHIANTI(i__),ENSD07(i__),n_ENSD07(i__)
if(ar__3(28:36) .EQ. '') ENIST_2(i__) = rNULL__
if(ar__3(51:59) .EQ. '') ECHIANTI(i__) = rNULL__
if(ar__3(65:73) .EQ. '') ENSD07(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ ID_2(i__),Level_2(i__),ENIST_2(i__),EAS_2(i__),ECHIANTI(i__),
+ n_ECHIANTI(i__),ENSD07(i__),n_ENSD07(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! *Level energies (Ryd) of S^10+^
C Format for file interpretation
5 format(1X,I3,3X,A15,3X,F9.6,3X,F7.4,4X,F7.4,4X,F7.4)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,254
read(1,'(A66)')ar__4
read(ar__4,5)
+ ID_3(i__),Level_3(i__),Enist_3(i__),EAS_3(i__),ELB03(i__),
+ EMCHF_2(i__)
if(ar__4(26:34) .EQ. '') Enist_3(i__) = rNULL__
if(ar__4(49:55) .EQ. '') ELB03(i__) = rNULL__
if(ar__4(60:66) .EQ. '') EMCHF_2(i__) = rNULL__
c ..............Just test output...........
write(6,5)
+ ID_3(i__),Level_3(i__),Enist_3(i__),EAS_3(i__),ELB03(i__),
+ EMCHF_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end