Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/789/2 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/ApJ/789/2 Methane partition function (Rey+, 2014)
*================================================================================
*Theoretical hot methane line lists up to T = 2000 K for astrophysical
*applications.
* Rey M., Nikitin A.V., Tyuterev V.G.
* <Astrophys. J., 789, 2 (2014)>
* =2014ApJ...789....2R (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Methane partition function calculated from
variational rovibrational energy levels
integer*4 nr__
parameter (nr__=156) ! Number of records
character*49 ar__ ! Full-size record
real*4 Temp (nr__) ! (K) Temperature
real*8 Q (nr__) ! Partition function from his work
character*1 f_Q (nr__) ! Flag on Q (1)
real*8 HITRAN (nr__) ! ? HITRAn partition function
real*4 e_HITRAN (nr__) ! (%) ? Percentage error between HITRAN and this work
real*8 Wenger (nr__) ! ? Wenger et al. (2008, JQSRT, 109, 2687)
* partition function
real*4 e_Wenger (nr__) ! (%) ? Percentage error between Wenger and this work
*Note (1): Flag as follows:
* a = 54175.0 extrapolated at J infinity;
* b = 72602.0 extrapolated at J infinity;
* c = 96400.0 extrapolated at J infinity;
* d = 126735.0 extrapolated at J infinity;
* e = 165270.0 extrapolated at J infinity.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'fig2.dat' ! Number of lines and convergence of the
integrated intensity as a function of the
intensity cutoff for different temperatures
integer*4 nr__1
parameter (nr__1=50) ! Number of records
character*30 ar__1 ! Full-size record
integer*4 Temp_1 (nr__1) ! (K) Temperature
integer*4 ICut (nr__1) ! Intensity cutoff; 10^-n^cm/mol
integer*8 Num (nr__1) ! Number of lines
real*4 Sum (nr__1) ! (cm/mol) Integrated intensity
C=============================================================================
C Loading file 'table2.dat' ! Methane partition function calculated from
* variational rovibrational energy levels
C Format for file interpretation
1 format(F6.1,1X,F10.3,1X,A1,1X,F9.2,1X,F3.1,1X,F9.2,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,156
read(1,'(A49)')ar__
read(ar__,1)
+ Temp(i__),Q(i__),f_Q(i__),HITRAN(i__),e_HITRAN(i__),
+ Wenger(i__),e_Wenger(i__)
if(ar__(21:29) .EQ. '') HITRAN(i__) = rNULL__
if(ar__(31:33) .EQ. '') e_HITRAN(i__) = rNULL__
if(ar__(35:43) .EQ. '') Wenger(i__) = rNULL__
if(ar__(45:49) .EQ. '') e_Wenger(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Temp(i__),Q(i__),f_Q(i__),HITRAN(i__),e_HITRAN(i__),
+ Wenger(i__),e_Wenger(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'fig2.dat' ! Number of lines and convergence of the
* integrated intensity as a function of the
* intensity cutoff for different temperatures
C Format for file interpretation
2 format(I4,1X,I2,1X,I11,1X,E10.4)
C Effective file loading
open(unit=1,status='old',file=
+'fig2.dat')
write(6,*) '....Loading file: fig2.dat'
do i__=1,50
read(1,'(A30)')ar__1
read(ar__1,2)Temp_1(i__),ICut(i__),Num(i__),Sum(i__)
c ..............Just test output...........
write(6,2)Temp_1(i__),ICut(i__),Num(i__),Sum(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end