Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/446/3842 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-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__
parameter (rNULL__=--2147483648.) ! NULL real number
parameter (iNULL__=(-2147483647-1)) ! NULL int number
integer idig ! testing NULL number
C=============================================================================
Cat. J/MNRAS/446/3842 Spectral line survey of two LOSs (Armijos-Abendano+, 2015)
*================================================================================
*3-mm spectral line survey of two lines of sight towards two typical cloud
*complexes in the Galactic Centre.
* Armijos-Abendano J., Martin-Pintado J., Requena-Torres M.A., Martin S.,
* Rodriguez-Franco A.
* <Mon. Not. R. Astron. Soc., 446, 3842-3862 (2015)>
* =2015MNRAS.446.3842A (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Line parameters for the LOS+0.693
integer*4 nr__
parameter (nr__=224) ! Number of records
character*126 ar__ ! Full-size record
character*15 Mol ! Molecule
character*1 n_Mol ! [bc] Note on molecule (1)
real*8 Freq ! (MHz) Frequency
character*30 Trans ! Transition
real*4 Area ! (K.km/s) ? Area
real*4 e_Area ! (K.km/s) ? Uncertainty in Area
real*4 b_VLSR ! (km/s) ? Local Standard of Rest velocity interval,
* lower value
real*4 B_VLSR_1 ! (km/s) ? Local Standard of Rest velocity interval,
* upper value
real*4 VLSR ! (km/s) ? Local Standard of Rest velocity
real*4 e_VLSR ! (km/s) ? Uncertainty in VLSR
character*1 n_VLSR ! [a] Note on VLSR (2)
real*4 Deltav1_2 ! (km/s) ? Average value of the FWHM of the line,
* {Delta}v_1/2_
real*4 e_Deltav1_2 ! (km/s) ? Uncertainty in Deltav1/2
character*1 n_Deltav1_2 ! [a] Note on Deltav1/2 (2)
real*4 Ta_ ! (mK) ? Ambient temperature
real*4 e_Ta_ ! (mK) ? Uncertainty in Ta*
character*1 n_Ta_ ! [a] Note on Ta*
character*9 Note ! Note(s) (3)
*Note (1): Note as follows:
* b = Substates EE, AA, EA, AE blended, we show just the most intense
* transition.
* c = Frequency refers to species A.
*Note (2): Note as follows:
* a = Parameter fixed in the Gaussian fit.
*Note (3): Note as follows:
* bl = Blended line;
* m = Multitransition line (frequency refers to the main component of the
* group);
* hf = Hyperfine structure (frequency refers to the main component of the
* group);
* hfa = Hyperfine component, it is possible to resolve this hyperfine component
* since its frequency is sufficiently far from the frequencies of the
* other hyperfine components;
* ot = Transition less affected by opacity;
* cl = This line is contaminated by the emission from an unknown molecular
* species;
* al = Absorption line;
* cd = This transition have been used to derive the column density (although
* several transitions of this molecule have been detected, there is an
* insufficient dynamical range in E_u_ to derive the column density by
* using a rotational diagram).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Line parameters for the LOS-0.11
integer*4 nr__1
parameter (nr__1=150) ! Number of records
character*126 ar__1 ! Full-size record
character*15 Mol_1 ! Molecule
character*1 n_Mol_1 ! [bc] Note on molecule (1)
real*8 Freq_1 ! (MHz) Frequency
character*30 Trans_1 ! Transition
real*4 Area_1 ! (K.km/s) ? Area
real*4 e_Area_1 ! (K.km/s) ? Uncertainty in Area
real*4 b_VLSR_2 ! (km/s) ? Local Standard of Rest velocity interval,
* lower value
real*4 B_VLSR_3 ! (km/s) ? Local Standard of Rest velocity interval,
* upper value
real*4 VLSR_1 ! (km/s) ? Local Standard of Rest velocity
real*4 e_VLSR_1 ! (km/s) ? Uncertainty in VLSR
character*1 n_VLSR_1 ! [a] Note on VLSR (2)
real*4 Deltav1_2_1 ! (km/s) ? Average value of the FWHM of the line,
* {Delta}v_1/2_
real*4 e_Deltav1_2_1 ! (km/s) ? Uncertainty in Deltav1/2
character*1 n_Deltav1_2_1 ! [a] Note on Deltav1/2 (2)
real*4 Ta__1 ! (mK) ? Ambient temperature
real*4 e_Ta__1 ! (mK) ? Uncertainty in Ta*
character*1 n_Ta__1 ! [a] Note on Ta*
character*9 Note_1 ! Note(s) (3)
*Note (1): Note as follows:
* b = Substates EE, AA, EA, AE blended, we show just the most intense
* transition.
* c = Frequency refers to species A.
*Note (2): Note as follows:
* a = Parameter fixed in the Gaussian fit.
*Note (3): Note as follows:
* bl = Blended line;
* m = Multitransition line (frequency refers to the main component of the
* group);
* hf = Hyperfine structure (frequency refers to the main component of the
* group);
* hfa = Hyperfine component, it is possible to resolve this hyperfine component
* since its frequency is sufficiently far from the frequencies of the
* other hyperfine components;
* ot = Transition less affected by opacity;
* cl = This line is contaminated by the emission from an unknown molecular
* species;
* al = Absorption line;
* cd = This transition have been used to derive the column density (although
* several transitions of this molecule have been detected, there is an
* insufficient dynamical range in E_u_ to derive the column density by
* using a rotational diagram).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! T_rot_, column densities and abundances for
both LOSs
integer*4 nr__2
parameter (nr__2=150) ! Number of records
character*99 ar__2 ! Full-size record
character*1 LOS ! [AB] Line of sight identification (1)
character*15 Mol_2 ! Molecule
character*1 n_Mol_2 ! [e] Note on Mol (2)
integer*4 b_VLSR_4 ! (km/s) ? Local Standard of Rest velocity interval,
* lower value
integer*4 B_VLSR_5 ! (km/s) ? Local Standard of Rest velocity interval,
* upper value
character*1 n_B_VLSR ! [g] Note on B_VLSR (3)
character*1 l_VLSR ! [~] Limit flag on VLSR
real*4 VLSR_2 ! (km/s) ? Local Standard of Rest velocity
real*4 e_VLSR_2 ! (km/s) ? Uncertainty in VLSR
character*1 n_VLSR_2 ! [cf] Note on VLSR (4)
real*4 Trot ! (K) ? Rotational temperature (5)
real*4 e_Trot ! (K) ? Uncertainty in Trot
character*1 n_Trot ! [c] Note on Trot (4)
character*2 l_N ! [>~ ] Limit flag on N
real*8 N ! (10+13/cm2) ? Local thermodynamic equilibrium total
* molecular column density
real*8 e_N ! (10+13/cm2) ? Uncertainty in N
character*1 n_N ! [bcd] Note on N (4)
real*4 b_N_NH2 ! (10-9) ? N/NH_2_ column density ratio interval,
* lower value
real*4 B_N_NH2_1 ! (10-9) ? N/NH_2_ column density ratio interval,
* upper value
character*1 n_B_N_NH2 ! [g] Note on B_N/NH2 (3)
character*2 l_N_NH2 ! [>~ ] Limit flag on N/NH2
real*8 N_NH2 ! (10-9) ? N/NH_2_ column density ratio
real*8 e_N_NH2 ! (10-9) ? Uncertainty in N/NH2
*Note (1): LOS as follows:
* A = LOS+0.693;
* B = LOS-0.11.
*Note (2): Note as follows:
* e = The observed transition is contaminated by the emission from an unknown
* molecular species.
*Note (3): Note as follows:
* g = These velocity ranges are chosen for deriving velocity-integrated
* intensities used in the molecular column density estimate. For LOS+0.693
* we have used a velocity range for the ^13^CS(2-1) line as it is affected
* by opacity or self-absorption. For the ^15^N isotopologues of HNC and HCN,
* the velocity ranges are suitable for deriving ^14^N/^15^N ratios (see
* text).
*Note (4): Note as follows:
* b = We have inferred from the ^12^C/^13^C~<15,^14^N/^15^N~<280 and
* ^16^O/^18^O~<186 isotopic ratios given in Table 4 that the column density
* of the most abundant isotopologues of these molecules are biased by
* opacity or self-absorption. Thus here we have derived the column density
* by using either the ^18^O, ^15^N or ^13^C isotopologue for the respective
* velocity component and assuming ^16^O/^18^O=250 or ^14^N/^15^N>600
* (Wilson & Rood, 1994ARA&A..32..191W) and if necessary our value
* ^12^C/^13^C=21. For LOS+0.693, the SiO column density is derived from the
* ^28^SiO isotopologue assuming ^28^Si/^30^Si=18 derived for LOS-0.11;
* c = Only one line of this molecule was detected;
* d = Although several transitions of this molecule are detected, there is an
* insufficient dynamical range for E_u_ to derive the column density from a
* rotational diagram, so we have chosen one transition, usually the less
* affected by opacity;
* f = This velocity is an average of different detected transitions.
*Note (5): T_rot_ derived from rotational diagrams (RDs) or assumed for deriving
* molecular column densities. The assumed T_rot_ for the ^13^C
* isotopologues of CH_3_CCH and CH_3_CN are taken from their other
* isotopologues. The assumed T_rot_ for CH_3_OH and its ^13^C
* isotopologue are taken from Requena-Torres et al.
* (2008ApJ...672..352R). The T_rot_ quoted with uncertainties are
* determined from RDs. When T_rot_ is not listed then T_rot_=10 K is
* assumed, which corresponds to an average value of the low T_rot_
* component derived from other molecules by using RDs.
C=============================================================================
C Loading file 'table1.dat' ! Line parameters for the LOS+0.693
C Format for file interpretation
1 format(
+ A15,A1,1X,F8.1,1X,A30,1X,F6.2,1X,F3.1,1X,F4.1,1X,F5.1,1X,F5.2,
+ 1X,F3.1,A1,1X,F4.1,1X,F5.2,A1,1X,F6.1,1X,F6.2,A1,1X,A9)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,224
read(1,'(A126)')ar__
read(ar__,1)
+ Mol,n_Mol,Freq,Trans,Area,e_Area,b_VLSR,B_VLSR_1,VLSR,e_VLSR,
+ n_VLSR,Deltav1_2,e_Deltav1_2,n_Deltav1_2,Ta_,e_Ta_,n_Ta_,Note
if(ar__(58:63) .EQ. '') Area = rNULL__
if(ar__(65:67) .EQ. '') e_Area = rNULL__
if(ar__(69:72) .EQ. '') b_VLSR = rNULL__
if(ar__(74:78) .EQ. '') B_VLSR_1 = rNULL__
if(ar__(80:84) .EQ. '') VLSR = rNULL__
if(ar__(86:88) .EQ. '') e_VLSR = rNULL__
if(ar__(91:94) .EQ. '') Deltav1_2 = rNULL__
if(ar__(96:100) .EQ. '') e_Deltav1_2 = rNULL__
if(ar__(103:108) .EQ. '') Ta_ = rNULL__
if(ar__(110:115) .EQ. '') e_Ta_ = rNULL__
c ..............Just test output...........
write(6,1)
+ Mol,n_Mol,Freq,Trans,Area,e_Area,b_VLSR,B_VLSR_1,VLSR,e_VLSR,
+ n_VLSR,Deltav1_2,e_Deltav1_2,n_Deltav1_2,Ta_,e_Ta_,n_Ta_,Note
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Line parameters for the LOS-0.11
C Format for file interpretation
2 format(
+ A15,A1,1X,F8.1,1X,A30,1X,F6.2,1X,F3.1,1X,F4.1,1X,F5.1,1X,F5.2,
+ 1X,F3.1,A1,1X,F4.1,1X,F5.2,A1,1X,F6.1,1X,F6.2,A1,1X,A9)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,150
read(1,'(A126)')ar__1
read(ar__1,2)
+ Mol_1,n_Mol_1,Freq_1,Trans_1,Area_1,e_Area_1,b_VLSR_2,
+ B_VLSR_3,VLSR_1,e_VLSR_1,n_VLSR_1,Deltav1_2_1,e_Deltav1_2_1,
+ n_Deltav1_2_1,Ta__1,e_Ta__1,n_Ta__1,Note_1
if(ar__1(58:63) .EQ. '') Area_1 = rNULL__
if(ar__1(65:67) .EQ. '') e_Area_1 = rNULL__
if(ar__1(69:72) .EQ. '') b_VLSR_2 = rNULL__
if(ar__1(74:78) .EQ. '') B_VLSR_3 = rNULL__
if(ar__1(80:84) .EQ. '') VLSR_1 = rNULL__
if(ar__1(86:88) .EQ. '') e_VLSR_1 = rNULL__
if(ar__1(91:94) .EQ. '') Deltav1_2_1 = rNULL__
if(ar__1(96:100) .EQ. '') e_Deltav1_2_1 = rNULL__
if(ar__1(103:108) .EQ. '') Ta__1 = rNULL__
if(ar__1(110:115) .EQ. '') e_Ta__1 = rNULL__
c ..............Just test output...........
write(6,2)
+ Mol_1,n_Mol_1,Freq_1,Trans_1,Area_1,e_Area_1,b_VLSR_2,
+ B_VLSR_3,VLSR_1,e_VLSR_1,n_VLSR_1,Deltav1_2_1,e_Deltav1_2_1,
+ n_Deltav1_2_1,Ta__1,e_Ta__1,n_Ta__1,Note_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! T_rot_, column densities and abundances for
* both LOSs
C Format for file interpretation
3 format(
+ A1,1X,A15,1X,A1,1X,I2,1X,I3,A1,1X,A1,F4.1,1X,F3.1,A1,1X,F4.1,
+ 1X,F5.2,A1,1X,A2,F8.2,1X,F7.2,A1,1X,F4.1,1X,F4.1,A1,1X,A2,
+ F7.2,1X,F7.3)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,150
read(1,'(A99)')ar__2
read(ar__2,3)
+ LOS,Mol_2,n_Mol_2,b_VLSR_4,B_VLSR_5,n_B_VLSR,l_VLSR,VLSR_2,
+ e_VLSR_2,n_VLSR_2,Trot,e_Trot,n_Trot,l_N,N,e_N,n_N,b_N_NH2,
+ B_N_NH2_1,n_B_N_NH2,l_N_NH2,N_NH2,e_N_NH2
if(ar__2(21:22) .EQ. '') b_VLSR_4 = iNULL__
if(ar__2(24:26) .EQ. '') B_VLSR_5 = iNULL__
if(ar__2(30:33) .EQ. '') VLSR_2 = rNULL__
if(ar__2(35:37) .EQ. '') e_VLSR_2 = rNULL__
if(ar__2(40:43) .EQ. '') Trot = rNULL__
if(ar__2(45:49) .EQ. '') e_Trot = rNULL__
if(ar__2(54:61) .EQ. '') N = rNULL__
if(ar__2(63:69) .EQ. '') e_N = rNULL__
if(ar__2(72:75) .EQ. '') b_N_NH2 = rNULL__
if(ar__2(77:80) .EQ. '') B_N_NH2_1 = rNULL__
if(ar__2(85:91) .EQ. '') N_NH2 = rNULL__
if(ar__2(93:99) .EQ. '') e_N_NH2 = rNULL__
c ..............Just test output...........
write(6,3)
+ LOS,Mol_2,n_Mol_2,b_VLSR_4,B_VLSR_5,n_B_VLSR,l_VLSR,VLSR_2,
+ e_VLSR_2,n_VLSR_2,Trot,e_Trot,n_Trot,l_N,N,e_N,n_N,b_N_NH2,
+ B_N_NH2_1,n_B_N_NH2,l_N_NH2,N_NH2,e_N_NH2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end