Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/293/953 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.5, on 2013-May-25
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__
ter (rNULL__=-1.e37) ! NULL real number
parameter (iNULL__=-2147483647) ! NULL int number
C=============================================================================
Cat. J/A+A/293/953 IRON project VI. FeII collision strengths (Zhang+, 1995)
*================================================================================
*Atomic data from the IRON project: VI. Collision strengths and
*rate coefficients for Fe II,
* Zhang H.L., Pradhan A.K.
* <Astron. Astrophys. 293, 953 (1995)>
* =1995A&A...293..953Z
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'config.dat' ! Configuration and observed energy
(log(temp(K))=3.-5, Z = 26, N = 25,
142 fine structure levels)
integer*4 nr__
parameter (nr__=142) ! Number of records
character*41 ar__ ! Full-size record
integer*4 Level (nr__) ! identification
character*11 ConfigI (nr__) ! Initial configuration (with or without one set
* of the intermediate spin and orbital e.g.,
* 5D in the first level means the intermediate
* spin and orbital are S=2 and L=2, and none in
* the 6th level)
character*4 ConfigT (nr__) ! Term configuration
integer*4 J (nr__) ! Total angular momentum
real*8 E (nr__) ! (Ry) Observed energy
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'gamma.dat' ! Effective collision strengths
integer*4 nr__1
parameter (nr__1=10011) ! Number of records
character*229 ar__1 ! Full-size record
integer*4 LevI (nr__1) ! Initial configuration index
integer*4 LevT (nr__1) ! Term configuration index
real*4 Gam1 (nr__1) ! Effective collision strength for T = 1000 K
real*4 Gam3 (nr__1) ! Effective collision strength for T = 3000 K
real*4 Gam5 (nr__1) ! Effective collision strength for T = 5000 K
real*4 Gam7 (nr__1) ! Effective collision strength for T = 7000 K
real*4 Gam10 (nr__1) ! Effective collision strength for T = 10000 K
real*4 Gam12 (nr__1) ! Effective collision strength for T = 12000 K
real*4 Gam15 (nr__1) ! Effective collision strength for T = 15000 K
real*4 Gam17 (nr__1) ! Effective collision strength for T = 17000 K
real*4 Gam20 (nr__1) ! Effective collision strength for T = 20000 K
real*4 Gam25 (nr__1) ! Effective collision strength for T = 25000 K
real*4 Gam30 (nr__1) ! Effective collision strength for T = 30000 K
real*4 Gam35 (nr__1) ! Effective collision strength for T = 35000 K
real*4 Gam40 (nr__1) ! Effective collision strength for T = 40000 K
real*4 Gam45 (nr__1) ! Effective collision strength for T = 45000 K
real*4 Gam50 (nr__1) ! Effective collision strength for T = 50000 K
real*4 Gam60 (nr__1) ! Effective collision strength for T = 60000 K
real*4 Gam70 (nr__1) ! Effective collision strength for T = 70000 K
real*4 Gam80 (nr__1) ! Effective collision strength for T = 80000 K
real*4 Gam90 (nr__1) ! Effective collision strength for T = 90000 K
real*4 Gam100 (nr__1) ! Effective collision strength for T = 100000 K
C=============================================================================
C Loading file 'config.dat' ! Configuration and observed energy
* (log(temp(K))=3.-5, Z = 26, N = 25,
* 142 fine structure levels)
C Format for file interpretation
1 format(2X,I3,2X,A11,2X,A4,5X,I2,3X,F7.5)
C Effective file loading
open(unit=1,file='config.dat', status='old')
write(6,*) '....Loading file: config.dat'
do i__=1,142
read(1,'(A41)')ar__
read(ar__,1)
+ Level(i__),ConfigI(i__),ConfigT(i__),J(i__),E(i__)
c ..............Just test output...........
write(6,1)
+ Level(i__),ConfigI(i__),ConfigT(i__),J(i__),E(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'gamma.dat' ! Effective collision strengths
C Format for file interpretation
2 format(
+ 1X,I3,1X,I3,E11.4,E11.4,E11.4,E11.4,E11.4,E11.4,E11.4,E11.4,
+ 1X,E11.4,E11.4,E11.4,E11.4,E11.4,E11.4,E11.4,E11.4,E11.4,
+ E11.4,E11.4,E11.4)
C Effective file loading
open(unit=1,file='gamma.dat', status='old')
write(6,*) '....Loading file: gamma.dat'
do i__=1,10011
read(1,'(A229)')ar__1
read(ar__1,2)
+ LevI(i__),LevT(i__),Gam1(i__),Gam3(i__),Gam5(i__),Gam7(i__),
+ Gam10(i__),Gam12(i__),Gam15(i__),Gam17(i__),Gam20(i__),
+ Gam25(i__),Gam30(i__),Gam35(i__),Gam40(i__),Gam45(i__),
+ Gam50(i__),Gam60(i__),Gam70(i__),Gam80(i__),Gam90(i__),
+ Gam100(i__)
c ..............Just test output...........
write(6,2)
+ LevI(i__),LevT(i__),Gam1(i__),Gam3(i__),Gam5(i__),Gam7(i__),
+ Gam10(i__),Gam12(i__),Gam15(i__),Gam17(i__),Gam20(i__),
+ Gam25(i__),Gam30(i__),Gam35(i__),Gam40(i__),Gam45(i__),
+ Gam50(i__),Gam60(i__),Gam70(i__),Gam80(i__),Gam90(i__),
+ Gam100(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end