Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/117/427 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.5, on 2013-May-19
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__
ter (rNULL__=-1.e37) ! NULL real number
parameter (iNULL__=-2147483647) ! NULL int number
C=============================================================================
Cat. J/ApJS/117/427 Sgr B2 spectral survey (Nummelin+, 1998)
*================================================================================
*A three-position spectral line survey of Sagittarius B2 between 218 and 263 GHz.
*I. The observational data.
* Nummelin A., Bergman P., Hjalmarson A., Friberg P., Irvine W.M.,
* Millar T.J., Ohishi M., Saito S.
* <Astrophys. J. Suppl. Ser. 117, 427 (1998)>
* =1998ApJS..117..427N (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Line identifications toward Sagittarius B2(N)
integer*4 nr__
parameter (nr__=1726) ! Number of records
character*95 ar__ ! Full-size record
integer*4 Freq (nr__) ! (MHz) Observed frequency
real*4 Tmb (nr__) ! (K) ? Peak main-beam brightness temperature
integer*4 DeltaV (nr__) ! (km/s) ? Velocity width
character*7 Ref (nr__) ! Orion references (1)
character*36 Iden (nr__) ! Identification
character*26 Rem (nr__) ! Remarks
*Note (1): B86: Blake et al., 1986ApJS...60..357B
* S85: Sutton et al., 1985ApJS...58..341S
* G91: Greaves & White, 1991A&AS...91..237G
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Line identifications toward Sagittarius B2(M)
integer*4 nr__1
parameter (nr__1=664) ! Number of records
character*95 ar__1 ! Full-size record
integer*4 Freq_1 (nr__1) ! (MHz) Observed frequency
real*4 Tmb_1 (nr__1) ! (K) ? Peak main-beam brightness temperature
integer*4 DeltaV_1 (nr__1) ! (km/s) ? Velocity width
character*7 Ref_1 (nr__1) ! Orion references (1)
character*36 Iden_1 (nr__1) ! Identification
character*26 Rem_1 (nr__1) ! Remarks
*Note (1): B86: Blake et al., 1986ApJS...60..357B
* S85: Sutton et al., 1985ApJS...58..341S
* G91: Greaves & White, 1991A&AS...91..237G
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Line identifications toward Sagittarius B2(NW)
integer*4 nr__2
parameter (nr__2=112) ! Number of records
character*95 ar__2 ! Full-size record
integer*4 Freq_2 (nr__2) ! (MHz) Observed frequency
real*4 Tmb_2 (nr__2) ! (K) ? Peak main-beam brightness temperature
integer*4 DeltaV_2 (nr__2) ! (km/s) ? Velocity width
character*7 Ref_2 (nr__2) ! Orion references (1)
character*36 Iden_2 (nr__2) ! Identification
character*26 Rem_2 (nr__2) ! Remarks
*Note (1): B86: Blake et al., 1986ApJS...60..357B
* S85: Sutton et al., 1985ApJS...58..341S
* G91: Greaves & White, 1991A&AS...91..237G
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tab6-45.dat' ! Lines listed species by species
integer*4 nr__3
parameter (nr__3=1861) ! Number of records
character*171 ar__3 ! Full-size record
character*16 El (nr__3) ! Element (1)
real*8 RFreq (nr__3) ! (MHz) Rest frequency
character*15 State (nr__3) ! State (2)
character*27 Trans (nr__3) ! Transition (3)
character*6 n_Trans (nr__3) ! Note on transition (4)
integer*4 Eu (nr__3) ! (K) Upper state energy (5)
real*4 Aul (nr__3) ! (s-1) A-coefficient (5)
real*4 N (nr__3) ! (K.km/s) ? Observed integrated intensity at Sgr B2(N) (6)
character*1 n_N (nr__3) ! [uda] Note (7)
real*4 M (nr__3) ! (K.km/s) ? Observed integrated intensity at Sgr B2(M)(6)
character*1 n_M (nr__3) ! [udba] Note (7)
real*4 NW (nr__3) ! (K.km/s) ? Observed integrated intensity at Sgr B2(NW)(6)
character*1 n_NW (nr__3) ! [uda] Note (7)
character*60 Rem_3 (nr__3) ! Remarks (8)
*Note (1): SO_2_ nu_2_=1 : SO_2_ in the {nu}_2_ = 1 Vibrational State
* NH_2_CN nu=1 : NH_2_CN in the v = 1 Vibrational State
* HC_3_N 1nu_5_ : HC_3_N in the 1{nu}_5_ Vibrational State
* HC_3_N 2nu_7_ : HC_3_N in the 2{nu}_7_ Vibrational State
* HC_3_N 1nu_6_ : HC_3_N in the 1{nu}_6_ Vibrational State
* HC_3_N 1nu_7_ : HC_3_N in the 1{nu}_7_ Vibrational State
* CH_3_OH vt=1 : CH_3_OH in the v_t_ = 1 Torsional State
* CH_3_OH vt=2 : CH_3_OH in the v_t_ = 2 Torsional State
* ^13^CH_3_OH vt=1 : ^13^CH_3_OH in the v_t_ = 1 Torsional State
* CH_3_CN nu_8_=1 : CH_3_CN in the {nu}_8_ = 1 Vibrational State
* CH_3_CHO vt=1 : CH_3_CHO in the v_t_ = 1 Torsional State
* CH_3_CHO vt=2 : CH_3_CHO in the v_t_ = 2 Torsional State
* C_2_H_3_CN 1nu11 : C_2_H_3_CN in the 1{nu}_11_ Vibrational State
* C_2_H_3_CN 1nu15 : C_2_H_3_CN in the 1{nu}_15_ Vibrational State
* C_2_H_3_CN 2nu11 : C_2_H_3_CN in the 2{nu}_11_ Vibrational State
* C_2_H_5_OH tran : C_2_H_5_OH in the tran Substates
* C_2_H_5_OH gau : C_2_H_5_OH in the gauche Substates
* CH_3_OCH_3_ For most transitions, the splitting of CH_3_OCH_3_
* into its four torsional substates AA,EE, EA and AE has not been
* resolved, and for brevity we only list J_Ka,Kc_ in these cases
*Note (2): For HCN and HCO+, species
* For CH_3_OH, CH_3_SH, CH_3_NH_2_, CH_3_CHO, CH_3_OCHO: symmetry
* Symmetry : A, E; a = asymmetric state, s = symmetric state
*Note (3): Type of transitions listed:
* J: SO+, HCN, HN^13^C, HCO, N_2_O, HCS+, OCS, NCNH+,
* _1_{SIGMA} diatomic species (CO, CS, SiO)
* J_F_: NO, NS
* N_J_: SO, C_2_H
* J_l_: HC_3_N
* J_K_: CH_3_CN, CH_3_C_2_H
* J_Ka,Kc_: HDO, SO_2_, H_2_CO, HNCO, HOCO+, H_2_CS+, CH2NH, c-C_3_H_2_,
* CH_2_CO, NH_2_CN, HCOOH, NH_2_CHO, c-c_2_H_4_O, C_2_H_3_CN,
* CH_3_OCHO, C_2_H_5_OH, CH_3_OCH_3_, C_2_H_5_CN
* J^p^_K_: CH_3_SH
* J^p^_Ka_: CH_3_OH, CH_3_NH_2_
* J^p^_Ka,Kc_: CH_3_CHO
* N_Ka,Kc_: CH_2_CN
*Note (4): F1 corresponds to J = (N + 3/2) -> (N + 1/2)
* F2 to J = (N + 1/2) -> (N - 1/2)
* a: a-Type Transitions
* b: b-type transitions
* c: c-type transitions
* a/b: a/b-Type transitions
* x: x-Type Transitions, See Plummer, Herbst, & De Lucia,
* 1987ApJ...318..873P
* For C_2_H_5_OH gau: + -> +; + -> -; - -> -; - -> +
*Note (5): ^33^SO: The E_u_ and A_ul_ values listed for are derived using the
* ^34^SO and SO values
* C_2_H_3_CN: In case of K-doublets, the Aul values listed refer to the
* individual transitions
* CH_3_OCH_3_: The rest frequencies and Aul values listed always refer to the
* subcomponent with highest statistical weight, EE, unless
* otherwise stated.
*Note (6): Absorption lines have been indicated in the tables by negative
* integrated intensities
*Note (7): u: up
* d: down
* a: line intensity incorrect due to emission in reference beam
* b: hyperfine-splitting partially resolved
*Note (8): a: ambiguous
* b: blend
* o: overlap
* gnd: ground vibrational state of C_2_H_3_CN
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table46.dat' ! Unidentified lines
integer*4 nr__4
parameter (nr__4=370) ! Number of records
character*80 ar__4 ! Full-size record
integer*4 Freq_3 (nr__4) ! (MHz) Observed frequency (1)
real*4 N_1 (nr__4) ! (K.km/s) ? Observed integrated intensity at Sgr B2(N)
real*4 M_1 (nr__4) ! (K.km/s) ? Observed integrated intensity at Sgr B2(M)
real*4 NW_1 (nr__4) ! (K.km/s) ? Observed integrated intensity at Sgr B2(NW) (2)
character*58 Rem_4 (nr__4) ! Remarks (3)
*Note (1): Assuming V_LSR_=+62km/s for all three positions.
* In case the U-line is common to N and M, the observed frequency
* listed is valid for N, and the frequency for M is noted in the
* Remarks column, if different. For the U-line that is common to M and
* NW (U247563) the listed observed frequency is for M.
*Note (2): Absorption lines have been indicated in the tables by negative
* integrated intensities
*Note (3): S85: Detected by Sutton et al., 1985ApJS...58..341S
* B86: Detected by Blake et al., 1986ApJS...60..357B
* G91: Detected by Greaves & White, 1991A&AS...91..237G
C=============================================================================
C Loading file 'table3.dat' ! Line identifications toward Sagittarius B2(N)
C Format for file interpretation
1 format(1X,I6,1X,F5.2,7X,I2,1X,A7,2X,A36,1X,A26)
C Effective file loading
open(unit=1,file='table3.dat', status='old')
write(6,*) '....Loading file: table3.dat'
do i__=1,1726
read(1,'(A95)')ar__
read(ar__,1)
+ Freq(i__),Tmb(i__),DeltaV(i__),Ref(i__),Iden(i__),Rem(i__)
if(ar__(9:13) .EQ. '') Tmb(i__) = rNULL__
if(ar__(21:22) .EQ. '') DeltaV(i__) = iNULL__
c ..............Just test output...........
write(6,1)
+ Freq(i__),Tmb(i__),DeltaV(i__),Ref(i__),Iden(i__),Rem(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Line identifications toward Sagittarius B2(M)
C Format for file interpretation
2 format(1X,I6,1X,F5.2,7X,I2,1X,A7,2X,A36,1X,A26)
C Effective file loading
open(unit=1,file='table4.dat', status='old')
write(6,*) '....Loading file: table4.dat'
do i__=1,664
read(1,'(A95)')ar__1
read(ar__1,2)
+ Freq_1(i__),Tmb_1(i__),DeltaV_1(i__),Ref_1(i__),Iden_1(i__),
+ Rem_1(i__)
if(ar__1(9:13) .EQ. '') Tmb_1(i__) = rNULL__
if(ar__1(21:22) .EQ. '') DeltaV_1(i__) = iNULL__
c ..............Just test output...........
write(6,2)
+ Freq_1(i__),Tmb_1(i__),DeltaV_1(i__),Ref_1(i__),Iden_1(i__),
+ Rem_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Line identifications toward Sagittarius B2(NW)
C Format for file interpretation
3 format(1X,I6,1X,F5.2,7X,I2,1X,A7,2X,A36,1X,A26)
C Effective file loading
open(unit=1,file='table5.dat', status='old')
write(6,*) '....Loading file: table5.dat'
do i__=1,112
read(1,'(A95)')ar__2
read(ar__2,3)
+ Freq_2(i__),Tmb_2(i__),DeltaV_2(i__),Ref_2(i__),Iden_2(i__),
+ Rem_2(i__)
if(ar__2(9:13) .EQ. '') Tmb_2(i__) = rNULL__
if(ar__2(21:22) .EQ. '') DeltaV_2(i__) = iNULL__
c ..............Just test output...........
write(6,3)
+ Freq_2(i__),Tmb_2(i__),DeltaV_2(i__),Ref_2(i__),Iden_2(i__),
+ Rem_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tab6-45.dat' ! Lines listed species by species
C Format for file interpretation
4 format(
+ A16,1X,F8.1,1X,A15,1X,A27,A6,I4,1X,E7.2,1X,F6.1,A1,1X,F6.1,A1,
+ 1X,F5.1,A1,1X,A60)
C Effective file loading
open(unit=1,file='tab6-45.dat', status='old')
write(6,*) '....Loading file: tab6-45.dat'
do i__=1,1861
read(1,'(A171)')ar__3
read(ar__3,4)
+ El(i__),RFreq(i__),State(i__),Trans(i__),n_Trans(i__),Eu(i__),
+ Aul(i__),N(i__),n_N(i__),M(i__),n_M(i__),NW(i__),n_NW(i__),
+ Rem_3(i__)
if(ar__3(89:94) .EQ. '') N(i__) = rNULL__
if(ar__3(97:102) .EQ. '') M(i__) = rNULL__
if(ar__3(105:109) .EQ. '') NW(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ El(i__),RFreq(i__),State(i__),Trans(i__),n_Trans(i__),Eu(i__),
+ Aul(i__),N(i__),n_N(i__),M(i__),n_M(i__),NW(i__),n_NW(i__),
+ Rem_3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table46.dat' ! Unidentified lines
C Format for file interpretation
5 format(I6,1X,F4.1,1X,F4.1,1X,F4.1,1X,A58)
C Effective file loading
open(unit=1,file='table46.dat', status='old')
write(6,*) '....Loading file: table46.dat'
do i__=1,370
read(1,'(A80)')ar__4
read(ar__4,5)
+ Freq_3(i__),N_1(i__),M_1(i__),NW_1(i__),Rem_4(i__)
if(ar__4(8:11) .EQ. '') N_1(i__) = rNULL__
if(ar__4(13:16) .EQ. '') M_1(i__) = rNULL__
if(ar__4(18:21) .EQ. '') NW_1(i__) = rNULL__
c ..............Just test output...........
write(6,5)
+ Freq_3(i__),N_1(i__),M_1(i__),NW_1(i__),Rem_4(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end