Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/503/293 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/503/293 Stark broadening of He447.1nm line (Gigosos+, 2009)
*================================================================================
*Stark broadening tables for the Helium 447.1 nm line.
*Application to weakly coupled plasmas diagnostics.
* Gigosos M.A., Gonzalez M.A.
* <Astron. Astrophys. 503, 293 (2009)>
* =2009A&A...503..293G
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tables.dat' ! Helium I 447.1 nm Stark profiles
(tables 2-11 of the paper)
integer*4 nr__
parameter (nr__=900) ! Number of records
character*161 ar__ ! Full-size record
real*4 logNe (nr__) ! ([m-3]) Electron density (1)
real*4 Dlambda (nr__) ! (nm) Wavelength detuning from the
* unperturbed wavelength
real*4 IDL0_8_5 (nr__) ! ? Stark profile for mu=0.8 and Te=5000K (2)
real*4 IDL0_8_10 (nr__) ! ? Stark profile for mu=0.8 and Te=10000K (2)
real*4 IDL0_8_20 (nr__) ! Stark profile for mu=0.8 and Te=20000K (2)
real*4 IDL0_8_40 (nr__) ! Stark profile for mu=0.8 and Te=40000K (2)
real*4 IDL2_5 (nr__) ! ? Stark profile for mu=2 and Te=5000K (2)
real*4 IDL2_10 (nr__) ! ? Stark profile for mu=2 and Te=10000K (2)
real*4 IDL2_20 (nr__) ! Stark profile for mu=2 and Te=20000K (2)
real*4 IDL2_40 (nr__) ! Stark profile for mu=2 and Te=40000K (2)
real*4 IDL4_50 (nr__) ! ? Stark profile for mu=4 and Te=5000K (2)
real*4 IDL4_10 (nr__) ! ? Stark profile for mu=4 and Te=10000K (2)
real*4 IDL4_20 (nr__) ! Stark profile for mu=4 and Te=20000K (2)
real*4 IDL4_40 (nr__) ! Stark profile for mu=4 and Te=40000K (2)
*Note (1): Electron density and {mu} (reduced mass of the emitter-ionic
* perturber, in units of proton mass) and T (plasma temperature, same
* for ions and electrons) values:
* 21.00 = data for mu=0.8, 2, 4, Te=5000, 10000, 20000, 40000K
* 21.33 = data for mu=0.8, 2, 4, Te=5000, 10000, 20000, 40000K
* 21.67 = data for mu=0.8, 2, 4, Te=5000, 10000, 20000, 40000K
* 22.00 = data for mu=0.8, 2, 4, Te=5000, 10000, 20000, 40000K
* 22.33 = data for mu=0.8, 2, 4, Te=5000, 10000, 20000, 40000K
* 22.67 = data for mu=0.8, 2, 4, Te=10000, 20000, 40000K
* 23.00 = data for mu=0.8, 2, 4, Te=10000, 20000, 40000K
* 23.33 = data for mu=0.8, 2, 4, Te=10000, 20000, 40000K
* 23.67 = data for mu=0.8, 2, 4, Te=20000, 40000K
* 24.00 = data for mu=0.8, 2, 4, Te=20000, 40000K
*Note (2): Area normalized Stark broadened spectral profiles. Doppler
* effect has not been included in the profiles. mu is the reduced
* mass in units of proton mass.
C=============================================================================
C Loading file 'tables.dat' ! Helium I 447.1 nm Stark profiles
* (tables 2-11 of the paper)
C Format for file interpretation
1 format(
+ F5.2,3X,E9.4,4X,E8.4,4X,E8.4,4X,E8.4,4X,E8.4,4X,E8.4,4X,E8.4,
+ 4X,E8.4,4X,E8.4,4X,E8.4,4X,E8.4,4X,E8.4,4X,E8.4)
C Effective file loading
open(unit=1,status='old',file=
+'tables.dat')
write(6,*) '....Loading file: tables.dat'
do i__=1,900
read(1,'(A161)')ar__
read(ar__,1)
+ logNe(i__),Dlambda(i__),IDL0_8_5(i__),IDL0_8_10(i__),
+ IDL0_8_20(i__),IDL0_8_40(i__),IDL2_5(i__),IDL2_10(i__),
+ IDL2_20(i__),IDL2_40(i__),IDL4_50(i__),IDL4_10(i__),
+ IDL4_20(i__),IDL4_40(i__)
if(ar__(22:29) .EQ. '') IDL0_8_5(i__) = rNULL__
if(ar__(34:41) .EQ. '') IDL0_8_10(i__) = rNULL__
if(ar__(70:77) .EQ. '') IDL2_5(i__) = rNULL__
if(ar__(82:89) .EQ. '') IDL2_10(i__) = rNULL__
if(ar__(118:125) .EQ. '') IDL4_50(i__) = rNULL__
if(ar__(130:137) .EQ. '') IDL4_10(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ logNe(i__),Dlambda(i__),IDL0_8_5(i__),IDL0_8_10(i__),
+ IDL0_8_20(i__),IDL0_8_40(i__),IDL2_5(i__),IDL2_10(i__),
+ IDL2_20(i__),IDL2_40(i__),IDL4_50(i__),IDL4_10(i__),
+ IDL4_20(i__),IDL4_40(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end