Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/792/118 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-Jun-13
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/792/118 Protonated oxirane characterization (Puzzarini+, 2014)
*================================================================================
*Accurate spectroscopic characterization of protonated oxirane:
*a potential prebiotic species in Titan's atmosphere.
* Puzzarini C., Ali A., Biczysko M., Barone V.
* <Astrophys. J., 792, 118 (2014)>
* =2014ApJ...792..118P (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Harmonic wavenumbers
integer*4 nr__
parameter (nr__=18) ! Number of records
character*9 ar__ ! Full-size record
real*8 Wave ! (cm-1) Harmonic wavenumber
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Anharmonicity Constants X Matrix (Xmn = Xnm)
integer*4 nr__1
parameter (nr__1=18) ! Number of records
character*254 ar__1 ! Full-size record
integer*4 Index ! Normal mode label
real*4 X1 ! (cm-1) X1n (n=1-18) Anharmonicity constant
real*4 X2 ! (cm-1) ? X2n (n=2-18) Anharmonicity constant
real*4 X3 ! (cm-1) ? X3n (n=3-18) Anharmonicity constant
real*4 X4 ! (cm-1) ? X4n (n=4-18) Anharmonicity constant
real*4 X5 ! (cm-1) ? X5n (n=5-18) Anharmonicity constant
real*4 X6 ! (cm-1) ? X6n (n=6-18) Anharmonicity constant
real*4 X7 ! (cm-1) ? X7n (n=7-18) Anharmonicity constant
real*4 X8 ! (cm-1) ? X8n (n=8-18) Anharmonicity constant
real*4 X9 ! (cm-1) ? X9n (n=9-18) Anharmonicity constant
real*4 X10 ! (cm-1) ? X10n (n=10-18) Anharmonicity constant
real*4 X11 ! (cm-1) ? X11n (n=11-18) Anharmonicity constant
real*4 X12 ! (cm-1) ? X12n (n=12-18) Anharmonicity constant
real*4 X13 ! (cm-1) ? X13n (n=13-18) Anharmonicity constant
real*4 X14 ! (cm-1) ? X14n (n=14-18) Anharmonicity constant
real*4 X15 ! (cm-1) ? X15n (n=15-18) Anharmonicity constant
real*4 X16 ! (cm-1) ? X16n (n=16-18) Anharmonicity constant
real*4 X17 ! (cm-1) ? X17n (n=17-18) Anharmonicity constant
real*4 X18 ! (cm-1) ? X1818 Anharmonicity constant
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Cubic, semi-diagonal quartic force constants
integer*4 nr__2
parameter (nr__2=3425) ! Number of records
character*27 ar__2 ! Full-size record
integer*4 I1 ! Normal mode numbering index 1
integer*4 I2 ! Normal mode numbering index 2
integer*4 I3 ! Normal mode numbering index 3
integer*4 I4 ! ? Normal mode numbering index 4 (1)
real*8 FC ! (cm-1) Cubic or semidiagonal quartic force constant
*Note (1): Blank for cubic force constants.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Coriolis couplings
integer*4 nr__3
parameter (nr__3=918) ! Number of records
character*17 ar__3 ! Full-size record
character*1 Index_1 ! Inertial axis numbering index
integer*4 I1_1 ! Normal mode numbering index 1
integer*4 I2_1 ! Normal mode numbering index 2
real*8 CC ! Coriolis coupling constant
C=============================================================================
C Loading file 'table1.dat' ! Harmonic wavenumbers
C Format for file interpretation
1 format(F9.4)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,18
read(1,'(A9)')ar__
read(ar__,1)Wave
c ..............Just test output...........
write(6,1)Wave
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Anharmonicity Constants X Matrix (Xmn = Xnm)
C Format for file interpretation
2 format(
+ I2,1X,E13.6,1X,E13.6,1X,E13.6,1X,E13.6,1X,E13.6,1X,E13.6,1X,
+ E13.6,1X,E13.6,1X,E13.6,1X,E13.6,1X,E13.6,1X,E13.6,1X,E13.6,
+ 1X,E13.6,1X,E13.6,1X,E13.6,1X,E13.6,1X,E13.6)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,18
read(1,'(A254)')ar__1
read(ar__1,2)
+ Index,X1,X2,X3,X4,X5,X6,X7,X8,X9,X10,X11,X12,X13,X14,X15,X16,
+ X17,X18
if(ar__1(18:30) .EQ. '') X2 = rNULL__
if(ar__1(32:44) .EQ. '') X3 = rNULL__
if(ar__1(46:58) .EQ. '') X4 = rNULL__
if(ar__1(60:72) .EQ. '') X5 = rNULL__
if(ar__1(74:86) .EQ. '') X6 = rNULL__
if(ar__1(88:100) .EQ. '') X7 = rNULL__
if(ar__1(102:114) .EQ. '') X8 = rNULL__
if(ar__1(116:128) .EQ. '') X9 = rNULL__
if(ar__1(130:142) .EQ. '') X10 = rNULL__
if(ar__1(144:156) .EQ. '') X11 = rNULL__
if(ar__1(158:170) .EQ. '') X12 = rNULL__
if(ar__1(172:184) .EQ. '') X13 = rNULL__
if(ar__1(186:198) .EQ. '') X14 = rNULL__
if(ar__1(200:212) .EQ. '') X15 = rNULL__
if(ar__1(214:226) .EQ. '') X16 = rNULL__
if(ar__1(228:240) .EQ. '') X17 = rNULL__
if(ar__1(242:254) .EQ. '') X18 = rNULL__
c ..............Just test output...........
write(6,2)
+ Index,X1,X2,X3,X4,X5,X6,X7,X8,X9,X10,X11,X12,X13,X14,X15,X16,
+ X17,X18
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Cubic, semi-diagonal quartic force constants
C Format for file interpretation
3 format(I2,1X,I2,1X,I2,1X,I2,1X,F15.10)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,3425
read(1,'(A27)')ar__2
read(ar__2,3)I1,I2,I3,I4,FC
if(ar__2(10:11) .EQ. '') I4 = iNULL__
c ..............Just test output...........
write(6,3)I1,I2,I3,I4,FC
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Coriolis couplings
C Format for file interpretation
4 format(A1,1X,I2,1X,I2,1X,F9.6)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,918
read(1,'(A17)')ar__3
read(ar__3,4)Index_1,I1_1,I2_1,CC
c ..............Just test output...........
write(6,4)Index_1,I1_1,I2_1,CC
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end