Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/792/118 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 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 (nr__) ! (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 (nr__1) ! Normal mode label
real*4 X1 (nr__1) ! (cm-1) X1n (n=1-18) Anharmonicity constant
real*4 X2 (nr__1) ! (cm-1) ? X2n (n=2-18) Anharmonicity constant
real*4 X3 (nr__1) ! (cm-1) ? X3n (n=3-18) Anharmonicity constant
real*4 X4 (nr__1) ! (cm-1) ? X4n (n=4-18) Anharmonicity constant
real*4 X5 (nr__1) ! (cm-1) ? X5n (n=5-18) Anharmonicity constant
real*4 X6 (nr__1) ! (cm-1) ? X6n (n=6-18) Anharmonicity constant
real*4 X7 (nr__1) ! (cm-1) ? X7n (n=7-18) Anharmonicity constant
real*4 X8 (nr__1) ! (cm-1) ? X8n (n=8-18) Anharmonicity constant
real*4 X9 (nr__1) ! (cm-1) ? X9n (n=9-18) Anharmonicity constant
real*4 X10 (nr__1) ! (cm-1) ? X10n (n=10-18) Anharmonicity constant
real*4 X11 (nr__1) ! (cm-1) ? X11n (n=11-18) Anharmonicity constant
real*4 X12 (nr__1) ! (cm-1) ? X12n (n=12-18) Anharmonicity constant
real*4 X13 (nr__1) ! (cm-1) ? X13n (n=13-18) Anharmonicity constant
real*4 X14 (nr__1) ! (cm-1) ? X14n (n=14-18) Anharmonicity constant
real*4 X15 (nr__1) ! (cm-1) ? X15n (n=15-18) Anharmonicity constant
real*4 X16 (nr__1) ! (cm-1) ? X16n (n=16-18) Anharmonicity constant
real*4 X17 (nr__1) ! (cm-1) ? X17n (n=17-18) Anharmonicity constant
real*4 X18 (nr__1) ! (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 (nr__2) ! Normal mode numbering index 1
integer*4 I2 (nr__2) ! Normal mode numbering index 2
integer*4 I3 (nr__2) ! Normal mode numbering index 3
integer*4 I4 (nr__2) ! ? Normal mode numbering index 4 (1)
real*8 FC (nr__2) ! (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 (nr__3) ! Inertial axis numbering index
integer*4 I1_1 (nr__3) ! Normal mode numbering index 1
integer*4 I2_1 (nr__3) ! Normal mode numbering index 2
real*8 CC (nr__3) ! 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(i__)
c ..............Just test output...........
write(6,1)Wave(i__)
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(i__),X1(i__),X2(i__),X3(i__),X4(i__),X5(i__),X6(i__),
+ X7(i__),X8(i__),X9(i__),X10(i__),X11(i__),X12(i__),X13(i__),
+ X14(i__),X15(i__),X16(i__),X17(i__),X18(i__)
if(ar__1(18:30) .EQ. '') X2(i__) = rNULL__
if(ar__1(32:44) .EQ. '') X3(i__) = rNULL__
if(ar__1(46:58) .EQ. '') X4(i__) = rNULL__
if(ar__1(60:72) .EQ. '') X5(i__) = rNULL__
if(ar__1(74:86) .EQ. '') X6(i__) = rNULL__
if(ar__1(88:100) .EQ. '') X7(i__) = rNULL__
if(ar__1(102:114) .EQ. '') X8(i__) = rNULL__
if(ar__1(116:128) .EQ. '') X9(i__) = rNULL__
if(ar__1(130:142) .EQ. '') X10(i__) = rNULL__
if(ar__1(144:156) .EQ. '') X11(i__) = rNULL__
if(ar__1(158:170) .EQ. '') X12(i__) = rNULL__
if(ar__1(172:184) .EQ. '') X13(i__) = rNULL__
if(ar__1(186:198) .EQ. '') X14(i__) = rNULL__
if(ar__1(200:212) .EQ. '') X15(i__) = rNULL__
if(ar__1(214:226) .EQ. '') X16(i__) = rNULL__
if(ar__1(228:240) .EQ. '') X17(i__) = rNULL__
if(ar__1(242:254) .EQ. '') X18(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Index(i__),X1(i__),X2(i__),X3(i__),X4(i__),X5(i__),X6(i__),
+ X7(i__),X8(i__),X9(i__),X10(i__),X11(i__),X12(i__),X13(i__),
+ X14(i__),X15(i__),X16(i__),X17(i__),X18(i__)
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(i__),I2(i__),I3(i__),I4(i__),FC(i__)
if(ar__2(10:11) .EQ. '') I4(i__) = iNULL__
c ..............Just test output...........
write(6,3)I1(i__),I2(i__),I3(i__),I4(i__),FC(i__)
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(i__),I1_1(i__),I2_1(i__),CC(i__)
c ..............Just test output...........
write(6,4)Index_1(i__),I1_1(i__),I2_1(i__),CC(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end