Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/674/A201 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-Aug-11
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/674/A201 He^+^(1s)-He(1s^2^) electron-capture processes (Wang+, 2023)
*================================================================================
*Elastic and electron capture processes in slow He+-He collision.
* Wang Y., Wang F., Lin X., Liu L., Wu Y., Wang J.
* <Astron. Astrophys. 674, A201 (2023)>
* =2023A&A...674A.201W (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablee1.dat' ! QMOCC calculations for total, spin-resolved
state-selective charge transfer cross-sections
for electron capture to the singlet and triplet
states of He(1snl, n=1-3) atom in He+-He
collision system.
Resluts of 34 molecular states.
integer*4 nr__
parameter (nr__=26) ! Number of records
character*105 ar__ ! Full-size record
real*4 Energy ! (eV/u) Collision energy
real*4 TotalCT ! (10-16cm2) Total charge transfer cross-section
* for 34 states
real*4 v1s2_1S ! (10-16cm2) Electron capture to He(1s2 1S)
real*8 v1s2s_1S ! (10-16cm2) Electron capture to He(1s2s 1S)
real*8 v1s2s_3S ! (10-16cm2) Electron capture to He(1s2s 3S)
real*8 v1s2p_3P ! (10-16cm2) Electron capture to He(1s2p 3P)
real*8 v1s2p_1P ! (10-16cm2) Electron capture to He(1s2p 1P)
real*8 v1s3s_3S ! (10-16cm2) Electron capture to He(1s3s 3S)
real*4 v1s3s_1S ! (10-16cm2) Electron capture to He(1s3s 1S)
real*4 v1s3p_3P ! (10-16cm2) Electron capture to He(1s3p 3P)
real*4 v1s3d_3D ! (10-16cm2) Electron capture to He(1s3d 3D)
real*4 v1s3s_1D ! (10-16cm2) Electron capture to He(1s3d 1D)
real*4 v1s3p_1P ! (10-16cm2) Electron capture to He(1s3p 1P)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablee2.dat' ! QMOCC calculations for resonant charge transfer
(or ground state transfer), elastic, and
transfer cross-sections in He+-He collision
system. Resluts of 2 molecular states
integer*4 nr__1
parameter (nr__1=750) ! Number of records
character*60 ar__1 ! Full-size record
real*8 Energy_1 ! (eV/u) Collision energy
real*4 ResonantCT ! (10-16cm2) Resonant charge transfer cross-section
real*4 Elastic ! (10-16cm2) Elastic cross-section
real*4 Total ! (10-16cm2) Elastic integer total cross-sections
real*4 ct_el ! (10-16cm2) Sum of ResonantCT and elastic cross-sections
real*4 Diffusion ! (10-16cm2) Diffusion cross-sections
real*8 Viscous ! (10-16cm2) Viscous cross-sections
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablee3.dat' ! nl-resolved charge transfer reaction rate
coefficients (RC) for He+(1s)-He(1s2) collisions
integer*4 nr__2
parameter (nr__2=64) ! Number of records
character*81 ar__2 ! Full-size record
integer*4 Temp ! (K) Temperature
real*4 TotalRC_2s ! (cm3/s) Total rate coefficient in case of 2-states
real*4 TotalRC_34 ! (cm3/s) Total rate coefficient in case of 34-states
real*4 v1s2 ! (cm3/s) Rate coefficient for transfer to He(1s2)
real*4 v1s2s ! (cm3/s) Rate coefficient for transfer to He(1s2s)
real*4 v1s2p ! (cm3/s) Rate coefficient for transfer to He(1s2p)
real*4 v1s3s ! (cm3/s) Rate coefficient for transfer to He(1s3s)
real*4 v1s3p ! (cm3/s) Rate coefficient for transfer to He(1s3p)
real*4 v1s3d ! (cm3/s) Rate coefficient for transfer to He(1s3d)
C=============================================================================
C Loading file 'tablee1.dat' ! QMOCC calculations for total, spin-resolved
* state-selective charge transfer cross-sections
* for electron capture to the singlet and triplet
* states of He(1snl, n=1-3) atom in He+-He
* collision system.
* Resluts of 34 molecular states.
C Format for file interpretation
1 format(
+ F6.1,1X,F5.2,1X,F5.2,1X,F7.5,1X,F7.5,1X,F7.5,1X,F8.6,1X,F8.6,
+ 1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3)
C Effective file loading
open(unit=1,status='old',file=
+'tablee1.dat')
write(6,*) '....Loading file: tablee1.dat'
do i__=1,26
read(1,'(A105)')ar__
read(ar__,1)
+ Energy,TotalCT,v1s2_1S,v1s2s_1S,v1s2s_3S,v1s2p_3P,v1s2p_1P,
+ v1s3s_3S,v1s3s_1S,v1s3p_3P,v1s3d_3D,v1s3s_1D,v1s3p_1P
c ..............Just test output...........
write(6,1)
+ Energy,TotalCT,v1s2_1S,v1s2s_1S,v1s2s_3S,v1s2p_3P,v1s2p_1P,
+ v1s3s_3S,v1s3s_1S,v1s3p_3P,v1s3d_3D,v1s3s_1D,v1s3p_1P
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablee2.dat' ! QMOCC calculations for resonant charge transfer
* (or ground state transfer), elastic, and
* transfer cross-sections in He+-He collision
* system. Resluts of 2 molecular states
C Format for file interpretation
2 format(F14.9,1X,F6.3,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F10.6)
C Effective file loading
open(unit=1,status='old',file=
+'tablee2.dat')
write(6,*) '....Loading file: tablee2.dat'
do i__=1,750
read(1,'(A60)')ar__1
read(ar__1,2)
+ Energy_1,ResonantCT,Elastic,Total,ct_el,Diffusion,Viscous
c ..............Just test output...........
write(6,2)
+ Energy_1,ResonantCT,Elastic,Total,ct_el,Diffusion,Viscous
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablee3.dat' ! nl-resolved charge transfer reaction rate
* coefficients (RC) for He+(1s)-He(1s2) collisions
C Format for file interpretation
3 format(
+ I9,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,
+ E8.3)
C Effective file loading
open(unit=1,status='old',file=
+'tablee3.dat')
write(6,*) '....Loading file: tablee3.dat'
do i__=1,64
read(1,'(A81)')ar__2
read(ar__2,3)
+ Temp,TotalRC_2s,TotalRC_34,v1s2,v1s2s,v1s2p,v1s3s,v1s3p,v1s3d
c ..............Just test output...........
write(6,3)
+ Temp,TotalRC_2s,TotalRC_34,v1s2,v1s2s,v1s2p,v1s3s,v1s3p,v1s3d
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end