Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/446/3842 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-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 (nr__) ! Molecule
character*1 n_Mol (nr__) ! [bc] Note on molecule (1)
real*8 Freq (nr__) ! (MHz) Frequency
character*30 Trans (nr__) ! Transition
real*4 Area (nr__) ! (K.km/s) ? Area
real*4 e_Area (nr__) ! (K.km/s) ? Uncertainty in Area
real*4 b_VLSR (nr__) ! (km/s) ? Local Standard of Rest velocity interval,
* lower value
real*4 B_VLSR_1 (nr__) ! (km/s) ? Local Standard of Rest velocity interval,
* upper value
real*4 VLSR (nr__) ! (km/s) ? Local Standard of Rest velocity
real*4 e_VLSR (nr__) ! (km/s) ? Uncertainty in VLSR
character*1 n_VLSR (nr__) ! [a] Note on VLSR (2)
real*4 Deltav1_2 (nr__) ! (km/s) ? Average value of the FWHM of the line,
* {Delta}v_1/2_
real*4 e_Deltav1_2(nr__) ! (km/s) ? Uncertainty in Deltav1/2
character*1 n_Deltav1_2(nr__) ! [a] Note on Deltav1/2 (2)
real*4 Ta_ (nr__) ! (mK) ? Ambient temperature
real*4 e_Ta_ (nr__) ! (mK) ? Uncertainty in Ta*
character*1 n_Ta_ (nr__) ! [a] Note on Ta*
character*9 Note (nr__) ! 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 (nr__1) ! Molecule
character*1 n_Mol_1 (nr__1) ! [bc] Note on molecule (1)
real*8 Freq_1 (nr__1) ! (MHz) Frequency
character*30 Trans_1 (nr__1) ! Transition
real*4 Area_1 (nr__1) ! (K.km/s) ? Area
real*4 e_Area_1 (nr__1) ! (K.km/s) ? Uncertainty in Area
real*4 b_VLSR_2 (nr__1) ! (km/s) ? Local Standard of Rest velocity interval,
* lower value
real*4 B_VLSR_3 (nr__1) ! (km/s) ? Local Standard of Rest velocity interval,
* upper value
real*4 VLSR_1 (nr__1) ! (km/s) ? Local Standard of Rest velocity
real*4 e_VLSR_1 (nr__1) ! (km/s) ? Uncertainty in VLSR
character*1 n_VLSR_1 (nr__1) ! [a] Note on VLSR (2)
real*4 Deltav1_2_1(nr__1) ! (km/s) ? Average value of the FWHM of the line,
* {Delta}v_1/2_
real*4 e_Deltav1_2_1(nr__1) ! (km/s) ? Uncertainty in Deltav1/2
character*1 n_Deltav1_2_1(nr__1) ! [a] Note on Deltav1/2 (2)
real*4 Ta__1 (nr__1) ! (mK) ? Ambient temperature
real*4 e_Ta__1 (nr__1) ! (mK) ? Uncertainty in Ta*
character*1 n_Ta__1 (nr__1) ! [a] Note on Ta*
character*9 Note_1 (nr__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 (nr__2) ! [AB] Line of sight identification (1)
character*15 Mol_2 (nr__2) ! Molecule
character*1 n_Mol_2 (nr__2) ! [e] Note on Mol (2)
integer*4 b_VLSR_4 (nr__2) ! (km/s) ? Local Standard of Rest velocity interval,
* lower value
integer*4 B_VLSR_5 (nr__2) ! (km/s) ? Local Standard of Rest velocity interval,
* upper value
character*1 n_B_VLSR (nr__2) ! [g] Note on B_VLSR (3)
character*1 l_VLSR (nr__2) ! [~] Limit flag on VLSR
real*4 VLSR_2 (nr__2) ! (km/s) ? Local Standard of Rest velocity
real*4 e_VLSR_2 (nr__2) ! (km/s) ? Uncertainty in VLSR
character*1 n_VLSR_2 (nr__2) ! [cf] Note on VLSR (4)
real*4 Trot (nr__2) ! (K) ? Rotational temperature (5)
real*4 e_Trot (nr__2) ! (K) ? Uncertainty in Trot
character*1 n_Trot (nr__2) ! [c] Note on Trot (4)
character*2 l_N (nr__2) ! [>~ ] Limit flag on N
real*8 N (nr__2) ! (10+13/cm2) ? Local thermodynamic equilibrium total
* molecular column density
real*8 e_N (nr__2) ! (10+13/cm2) ? Uncertainty in N
character*1 n_N (nr__2) ! [bcd] Note on N (4)
real*4 b_N_NH2 (nr__2) ! (10-9) ? N/NH_2_ column density ratio interval,
* lower value
real*4 B_N_NH2_1 (nr__2) ! (10-9) ? N/NH_2_ column density ratio interval,
* upper value
character*1 n_B_N_NH2 (nr__2) ! [g] Note on B_N/NH2 (3)
character*2 l_N_NH2 (nr__2) ! [>~ ] Limit flag on N/NH2
real*8 N_NH2 (nr__2) ! (10-9) ? N/NH_2_ column density ratio
real*8 e_N_NH2 (nr__2) ! (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(i__),n_Mol(i__),Freq(i__),Trans(i__),Area(i__),
+ e_Area(i__),b_VLSR(i__),B_VLSR_1(i__),VLSR(i__),e_VLSR(i__),
+ n_VLSR(i__),Deltav1_2(i__),e_Deltav1_2(i__),n_Deltav1_2(i__),
+ Ta_(i__),e_Ta_(i__),n_Ta_(i__),Note(i__)
if(ar__(58:63) .EQ. '') Area(i__) = rNULL__
if(ar__(65:67) .EQ. '') e_Area(i__) = rNULL__
if(ar__(69:72) .EQ. '') b_VLSR(i__) = rNULL__
if(ar__(74:78) .EQ. '') B_VLSR_1(i__) = rNULL__
if(ar__(80:84) .EQ. '') VLSR(i__) = rNULL__
if(ar__(86:88) .EQ. '') e_VLSR(i__) = rNULL__
if(ar__(91:94) .EQ. '') Deltav1_2(i__) = rNULL__
if(ar__(96:100) .EQ. '') e_Deltav1_2(i__) = rNULL__
if(ar__(103:108) .EQ. '') Ta_(i__) = rNULL__
if(ar__(110:115) .EQ. '') e_Ta_(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Mol(i__),n_Mol(i__),Freq(i__),Trans(i__),Area(i__),
+ e_Area(i__),b_VLSR(i__),B_VLSR_1(i__),VLSR(i__),e_VLSR(i__),
+ n_VLSR(i__),Deltav1_2(i__),e_Deltav1_2(i__),n_Deltav1_2(i__),
+ Ta_(i__),e_Ta_(i__),n_Ta_(i__),Note(i__)
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(i__),n_Mol_1(i__),Freq_1(i__),Trans_1(i__),Area_1(i__),
+ e_Area_1(i__),b_VLSR_2(i__),B_VLSR_3(i__),VLSR_1(i__),
+ e_VLSR_1(i__),n_VLSR_1(i__),Deltav1_2_1(i__),
+ e_Deltav1_2_1(i__),n_Deltav1_2_1(i__),Ta__1(i__),e_Ta__1(i__),
+ n_Ta__1(i__),Note_1(i__)
if(ar__1(58:63) .EQ. '') Area_1(i__) = rNULL__
if(ar__1(65:67) .EQ. '') e_Area_1(i__) = rNULL__
if(ar__1(69:72) .EQ. '') b_VLSR_2(i__) = rNULL__
if(ar__1(74:78) .EQ. '') B_VLSR_3(i__) = rNULL__
if(ar__1(80:84) .EQ. '') VLSR_1(i__) = rNULL__
if(ar__1(86:88) .EQ. '') e_VLSR_1(i__) = rNULL__
if(ar__1(91:94) .EQ. '') Deltav1_2_1(i__) = rNULL__
if(ar__1(96:100) .EQ. '') e_Deltav1_2_1(i__) = rNULL__
if(ar__1(103:108) .EQ. '') Ta__1(i__) = rNULL__
if(ar__1(110:115) .EQ. '') e_Ta__1(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Mol_1(i__),n_Mol_1(i__),Freq_1(i__),Trans_1(i__),Area_1(i__),
+ e_Area_1(i__),b_VLSR_2(i__),B_VLSR_3(i__),VLSR_1(i__),
+ e_VLSR_1(i__),n_VLSR_1(i__),Deltav1_2_1(i__),
+ e_Deltav1_2_1(i__),n_Deltav1_2_1(i__),Ta__1(i__),e_Ta__1(i__),
+ n_Ta__1(i__),Note_1(i__)
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(i__),Mol_2(i__),n_Mol_2(i__),b_VLSR_4(i__),B_VLSR_5(i__),
+ n_B_VLSR(i__),l_VLSR(i__),VLSR_2(i__),e_VLSR_2(i__),
+ n_VLSR_2(i__),Trot(i__),e_Trot(i__),n_Trot(i__),l_N(i__),
+ N(i__),e_N(i__),n_N(i__),b_N_NH2(i__),B_N_NH2_1(i__),
+ n_B_N_NH2(i__),l_N_NH2(i__),N_NH2(i__),e_N_NH2(i__)
if(ar__2(21:22) .EQ. '') b_VLSR_4(i__) = iNULL__
if(ar__2(24:26) .EQ. '') B_VLSR_5(i__) = iNULL__
if(ar__2(30:33) .EQ. '') VLSR_2(i__) = rNULL__
if(ar__2(35:37) .EQ. '') e_VLSR_2(i__) = rNULL__
if(ar__2(40:43) .EQ. '') Trot(i__) = rNULL__
if(ar__2(45:49) .EQ. '') e_Trot(i__) = rNULL__
if(ar__2(54:61) .EQ. '') N(i__) = rNULL__
if(ar__2(63:69) .EQ. '') e_N(i__) = rNULL__
if(ar__2(72:75) .EQ. '') b_N_NH2(i__) = rNULL__
if(ar__2(77:80) .EQ. '') B_N_NH2_1(i__) = rNULL__
if(ar__2(85:91) .EQ. '') N_NH2(i__) = rNULL__
if(ar__2(93:99) .EQ. '') e_N_NH2(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ LOS(i__),Mol_2(i__),n_Mol_2(i__),b_VLSR_4(i__),B_VLSR_5(i__),
+ n_B_VLSR(i__),l_VLSR(i__),VLSR_2(i__),e_VLSR_2(i__),
+ n_VLSR_2(i__),Trot(i__),e_Trot(i__),n_Trot(i__),l_N(i__),
+ N(i__),e_N(i__),n_N(i__),b_N_NH2(i__),B_N_NH2_1(i__),
+ n_B_N_NH2(i__),l_N_NH2(i__),N_NH2(i__),e_N_NH2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end