Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/609/A129 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 2024-Apr-20
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/A+A/609/A129 HCN, HNC and DNC spectra of 27 sources (Colzi+, 2018)
*================================================================================
*Nitrogen and hydrogen fractionation in high-mass star forming cores from
*observations of HCN and HNC.
* Colzi L., Fontani F., Caselli P., Ceccarelli C., Hily-Blant P.,
* Bizzocchi L.
* <Astron. Astrophys. 609, A129 (2018)>
* =2018A&A...609A.129C (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'sources.dat' ! *Total column densities (beam-averaged), of
H^15^NC, HN^13^C, DNC (table 3) and HC^15^N and
H^13^CN(1-0) transitions (table 4)
integer*4 nr__
parameter (nr__=27) ! Number of records
character*142 ar__ ! Full-size record
character*5 Type (nr__) ! Type (HMPO, HMSC or UCHII)
character*13 Source (nr__) ! Source name
character*1 n_Source (nr__) ! [w] Note on Source (1)
character*2 l_N_H15NC (nr__) ! [<= ] Limit flag on N(H15NC)
integer*4 N_H15NC (nr__) ! (10+10cm-2) Total H^15^NC column density
* (beam-averaged)
integer*4 e_N_H15NC (nr__) ! (10+10cm-2) ? rms uncertainty on N(H15NC) (2)
character*1 n_N_H15NC (nr__) ! [ut] Note on N(H15NC) (1)
integer*4 DN1 (nr__) ! (10+10cm-2) ?=- Error on N(H15NC) without considering
* the calibration error
real*4 N_HN13C (nr__) ! (10+11cm-2) Total HN^13^C column density
* (beam-averaged)
real*4 e_N_HN13C (nr__) ! (10+11cm-2) rms uncertainty on N(HN13C) (2)
real*4 DN2 (nr__) ! (10+11cm-2) Error on N(HN13C) without considering the
* calibration error
real*4 N_DNC (nr__) ! (10+11cm-2) Total DNC column density (beam-averaged)
real*4 e_N_DNC (nr__) ! (10+11cm-2) rms uncertainty on N(DNC) (2)
character*2 l_HNC_H15NC(nr__) ! [>= ] Limit flag on HNC/H15NC
integer*4 HNC_H15NC (nr__) ! HNC/H^15^NC isotopic ratio (3)
integer*4 e_HNC_H15NC(nr__) ! ? rms uncertainty on HNC/H15NC (2)
integer*4 HNC_DNC (nr__) ! HNC/DNC isotopic ratio (4)
integer*4 e_HNC_DNC (nr__) ! rms uncertainty on HNC/DNC (2)
real*4 DNC_HNC (nr__) ! DNC/HNC isotopic ratio
real*4 e_DNC_HNC (nr__) ! rms uncertainty on DNC/HNC (2)
character*2 l_N_HC15N (nr__) ! [<= ] Limit flag on N(HC15N)
integer*4 N_HC15N (nr__) ! (10+10cm-2) Total N(HC15N) column density
* (beam-averaged)
integer*4 e_N_HC15N (nr__) ! (10+10cm-2) ? rms uncertainty on N(HC15N)
character*1 n_N_HC15N (nr__) ! [ut] Note on N(HC15N) (1)
integer*4 DN3 (nr__) ! (10+10cm-2) ?=- Error on N(HC15N) without considering
* the calibration error
integer*4 N_H13CN (nr__) ! (10+11cm-2) Total H^13^CN column density
* (beam-averaged)
integer*4 e_N_H13CN (nr__) ! (10+11cm-2) rms uncertainty on N(H13CN) (2)
real*4 DN4 (nr__) ! (10+11cm-2) Error on N(H13CN) without considering the
* calibration error
character*2 l_HCN_HC15N(nr__) ! [>= ] Limit flag on HCN/HC15N
integer*4 HCN_HC15N (nr__) ! HCN/HC15N isotopic ratio (3)
integer*4 e_HCN_HC15N(nr__) ! ? rms uncertainty on HCN/HC15N (2)
integer*4 Tk (nr__) ! (K) Kinetic temperature of the clump (5)
character*1 n_Tk (nr__) ! [*] Note on Tk (6)
character*12 FName (nr__) ! Name used in the spectra files (7)
*Note (1): Flag as follows:
* w = "warm" HMSC
* u = upper limit
* t = tentative detection
*Note (2): Uncertainties in the column densities and in the isotope ratios
* have been computed as explained in Sect. 3.1 and 3.2.
*Note (3): it has been multiplied by ^12^C/^13^C as described in Section (3.1).
*Note (4): it has been multiplied by ^12^C/^13^C and by the correction of the
* different beams 3.09.
*Note (5): kinetic temperatures of the clumps derived from Fontani et al.
* (2015ApJ...808L..46F): for the sources without a derivation of Tk, the mean
* value for that evolutionary stage was taken (for the HMSCs the average was
* done without the "warm" ones, i.e those with a Tk<20K).
*Note (6): * for average value for the specific evolutionary stage
*Note (7): The spectra files are named :
* FName-dnc.spt in dnc subdirectory,
* FName-h13cn.spt in h13cn subdirectory,
* FName-h15cn.spt in 15cnn subdirectory,
* FName-hc15n.spt in hc15n subdirectory,
* FName-hn13c.spt in hn13c subdirectory.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'dnc/*' ! DNC spectra
integer*4 nr__1
parameter (nr__1=27) ! Number of records
character*89 ar__1 ! Full-size record
real*4 Vlsr (nr__1) ! (km/s) LSR velocity
real*4 Tmb (nr__1) ! (K) Beam brightness temperature
real*4 Tmbf1 (nr__1) ! ? Tmb computed in the fit
real*4 Tmbf2 (nr__1) ! ? Tmb computed in the fit
real*4 Vlsr2 (nr__1) ! ? LSR velocity for the possible second
* component
real*4 Tmbs1 (nr__1) ! ? Beam brightness temperature for the possible
* second component
real*4 Tmbs2 (nr__1) ! ? Tmb computed in the fit for the possible
* second component
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'h13cn/*' ! H^13^CN spectra
integer*4 nr__2
parameter (nr__2=27) ! Number of records
character*89 ar__2 ! Full-size record
real*4 Vlsr_1 (nr__2) ! (km/s) LSR velocity
real*4 Tmb_1 (nr__2) ! (K) Beam brightness temperature
real*4 Tmbf1_1 (nr__2) ! ? Tmb computed in the fit
real*4 Tmbf2_1 (nr__2) ! ? Tmb computed in the fit
real*4 Vlsr2_1 (nr__2) ! ? LSR velocity for the possible second
* component
real*4 Tmbs1_1 (nr__2) ! ? Beam brightness temperature for the possible
* second component
real*4 Tmbs2_1 (nr__2) ! ? Tmb computed in the fit for the possible
* second component
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'h15nc/*' ! H^15^NC spectra
integer*4 nr__3
parameter (nr__3=27) ! Number of records
character*89 ar__3 ! Full-size record
real*4 Vlsr_2 (nr__3) ! (km/s) LSR velocity
real*4 Tmb_2 (nr__3) ! (K) Beam brightness temperature
real*4 Tmbf1_2 (nr__3) ! ? Tmb computed in the fit
real*4 Tmbf2_2 (nr__3) ! ? Tmb computed in the fit
real*4 Vlsr2_2 (nr__3) ! ? LSR velocity for the possible second
* component
real*4 Tmbs1_2 (nr__3) ! ? Beam brightness temperature for the possible
* second component
real*4 Tmbs2_2 (nr__3) ! ? Tmb computed in the fit for the possible
* second component
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hc15n/*' ! HC^15^N spectra
integer*4 nr__4
parameter (nr__4=27) ! Number of records
character*89 ar__4 ! Full-size record
real*4 Vlsr_3 (nr__4) ! (km/s) LSR velocity
real*4 Tmb_3 (nr__4) ! (K) Beam brightness temperature
real*4 Tmbf1_3 (nr__4) ! ? Tmb computed in the fit
real*4 Tmbf2_3 (nr__4) ! ? Tmb computed in the fit
real*4 Vlsr2_3 (nr__4) ! ? LSR velocity for the possible second
* component
real*4 Tmbs1_3 (nr__4) ! ? Beam brightness temperature for the possible
* second component
real*4 Tmbs2_3 (nr__4) ! ? Tmb computed in the fit for the possible
* second component
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hn13c/*' ! HN^13^C spectra
integer*4 nr__5
parameter (nr__5=27) ! Number of records
character*89 ar__5 ! Full-size record
real*4 Vlsr_4 (nr__5) ! (km/s) LSR velocity
real*4 Tmb_4 (nr__5) ! (K) Beam brightness temperature
real*4 Tmbf1_4 (nr__5) ! ? Tmb computed in the fit
real*4 Tmbf2_4 (nr__5) ! ? Tmb computed in the fit
real*4 Vlsr2_4 (nr__5) ! ? LSR velocity for the possible second
* component
real*4 Tmbs1_4 (nr__5) ! ? Beam brightness temperature for the possible
* second component
real*4 Tmbs2_4 (nr__5) ! ? Tmb computed in the fit for the possible
* second component
C=============================================================================
C Loading file 'sources.dat' ! *Total column densities (beam-averaged), of
* H^15^NC, HN^13^C, DNC (table 3) and HC^15^N and
* H^13^CN(1-0) transitions (table 4)
C Format for file interpretation
1 format(
+ A5,1X,A13,A1,A2,I3,1X,I2,A1,1X,I2,1X,F5.1,1X,F4.1,1X,F3.1,1X,
+ F4.1,1X,F3.1,1X,A2,I3,1X,I3,1X,I4,1X,I3,1X,E6.2,1X,E6.2,1X,A2,
+ I4,1X,I3,A1,1X,I2,1X,I3,1X,I2,1X,F3.1,1X,A2,I3,1X,I3,1X,I2,A1,
+ 1X,A12)
C Effective file loading
open(unit=1,status='old',file=
+'sources.dat')
write(6,*) '....Loading file: sources.dat'
do i__=1,27
read(1,'(A142)')ar__
read(ar__,1)
+ Type(i__),Source(i__),n_Source(i__),l_N_H15NC(i__),
+ N_H15NC(i__),e_N_H15NC(i__),n_N_H15NC(i__),DN1(i__),
+ N_HN13C(i__),e_N_HN13C(i__),DN2(i__),N_DNC(i__),e_N_DNC(i__),
+ l_HNC_H15NC(i__),HNC_H15NC(i__),e_HNC_H15NC(i__),HNC_DNC(i__),
+ e_HNC_DNC(i__),DNC_HNC(i__),e_DNC_HNC(i__),l_N_HC15N(i__),
+ N_HC15N(i__),e_N_HC15N(i__),n_N_HC15N(i__),DN3(i__),
+ N_H13CN(i__),e_N_H13CN(i__),DN4(i__),l_HCN_HC15N(i__),
+ HCN_HC15N(i__),e_HCN_HC15N(i__),Tk(i__),n_Tk(i__),FName(i__)
if(ar__(27:28) .EQ. '') e_N_H15NC(i__) = iNULL__
if (DN1(i__) .EQ. 45) DN1(i__) = iNULL__
if(ar__(64:66) .EQ. '') e_HNC_H15NC(i__) = iNULL__
if(ar__(98:100) .EQ. '') e_N_HC15N(i__) = iNULL__
if (DN3(i__) .EQ. 45) DN3(i__) = iNULL__
if(ar__(123:125) .EQ. '') e_HCN_HC15N(i__) = iNULL__
c ..............Just test output...........
write(6,1)
+ Type(i__),Source(i__),n_Source(i__),l_N_H15NC(i__),
+ N_H15NC(i__),e_N_H15NC(i__),n_N_H15NC(i__),DN1(i__),
+ N_HN13C(i__),e_N_HN13C(i__),DN2(i__),N_DNC(i__),e_N_DNC(i__),
+ l_HNC_H15NC(i__),HNC_H15NC(i__),e_HNC_H15NC(i__),HNC_DNC(i__),
+ e_HNC_DNC(i__),DNC_HNC(i__),e_DNC_HNC(i__),l_N_HC15N(i__),
+ N_HC15N(i__),e_N_HC15N(i__),n_N_HC15N(i__),DN3(i__),
+ N_H13CN(i__),e_N_H13CN(i__),DN4(i__),l_HCN_HC15N(i__),
+ HCN_HC15N(i__),e_HCN_HC15N(i__),Tk(i__),n_Tk(i__),FName(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'dnc/*' ! DNC spectra
C Format for file interpretation
2 format(
+ E17.13,1X,E11.5,2X,E10.5,2X,E10.5,2X,E10.5,2X,E10.5,2X,E10.5)
C Effective file loading
open(unit=1,status='old',file=
+'dnc/*')
write(6,*) '....Loading file: dnc/*'
do i__=1,27
read(1,'(A89)')ar__1
read(ar__1,2)
+ Vlsr(i__),Tmb(i__),Tmbf1(i__),Tmbf2(i__),Vlsr2(i__),
+ Tmbs1(i__),Tmbs2(i__)
if(ar__1(32:41) .EQ. '') Tmbf1(i__) = rNULL__
if(ar__1(44:53) .EQ. '') Tmbf2(i__) = rNULL__
if(ar__1(56:65) .EQ. '') Vlsr2(i__) = rNULL__
if(ar__1(68:77) .EQ. '') Tmbs1(i__) = rNULL__
if(ar__1(80:89) .EQ. '') Tmbs2(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Vlsr(i__),Tmb(i__),Tmbf1(i__),Tmbf2(i__),Vlsr2(i__),
+ Tmbs1(i__),Tmbs2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'h13cn/*' ! H^13^CN spectra
C Format for file interpretation
3 format(
+ E17.13,1X,E11.5,2X,E10.5,2X,E10.5,2X,E10.5,2X,E10.5,2X,E10.5)
C Effective file loading
open(unit=1,status='old',file=
+'h13cn/*')
write(6,*) '....Loading file: h13cn/*'
do i__=1,27
read(1,'(A89)')ar__2
read(ar__2,3)
+ Vlsr_1(i__),Tmb_1(i__),Tmbf1_1(i__),Tmbf2_1(i__),Vlsr2_1(i__),
+ Tmbs1_1(i__),Tmbs2_1(i__)
if(ar__2(32:41) .EQ. '') Tmbf1_1(i__) = rNULL__
if(ar__2(44:53) .EQ. '') Tmbf2_1(i__) = rNULL__
if(ar__2(56:65) .EQ. '') Vlsr2_1(i__) = rNULL__
if(ar__2(68:77) .EQ. '') Tmbs1_1(i__) = rNULL__
if(ar__2(80:89) .EQ. '') Tmbs2_1(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ Vlsr_1(i__),Tmb_1(i__),Tmbf1_1(i__),Tmbf2_1(i__),Vlsr2_1(i__),
+ Tmbs1_1(i__),Tmbs2_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'h15nc/*' ! H^15^NC spectra
C Format for file interpretation
4 format(
+ E17.13,1X,E11.5,2X,E10.5,2X,E10.5,2X,E10.5,2X,E10.5,2X,E10.5)
C Effective file loading
open(unit=1,status='old',file=
+'h15nc/*')
write(6,*) '....Loading file: h15nc/*'
do i__=1,27
read(1,'(A89)')ar__3
read(ar__3,4)
+ Vlsr_2(i__),Tmb_2(i__),Tmbf1_2(i__),Tmbf2_2(i__),Vlsr2_2(i__),
+ Tmbs1_2(i__),Tmbs2_2(i__)
if(ar__3(32:41) .EQ. '') Tmbf1_2(i__) = rNULL__
if(ar__3(44:53) .EQ. '') Tmbf2_2(i__) = rNULL__
if(ar__3(56:65) .EQ. '') Vlsr2_2(i__) = rNULL__
if(ar__3(68:77) .EQ. '') Tmbs1_2(i__) = rNULL__
if(ar__3(80:89) .EQ. '') Tmbs2_2(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ Vlsr_2(i__),Tmb_2(i__),Tmbf1_2(i__),Tmbf2_2(i__),Vlsr2_2(i__),
+ Tmbs1_2(i__),Tmbs2_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'hc15n/*' ! HC^15^N spectra
C Format for file interpretation
5 format(
+ E17.13,1X,E11.5,2X,E10.5,2X,E10.5,2X,E10.5,2X,E10.5,2X,E10.5)
C Effective file loading
open(unit=1,status='old',file=
+'hc15n/*')
write(6,*) '....Loading file: hc15n/*'
do i__=1,27
read(1,'(A89)')ar__4
read(ar__4,5)
+ Vlsr_3(i__),Tmb_3(i__),Tmbf1_3(i__),Tmbf2_3(i__),Vlsr2_3(i__),
+ Tmbs1_3(i__),Tmbs2_3(i__)
if(ar__4(32:41) .EQ. '') Tmbf1_3(i__) = rNULL__
if(ar__4(44:53) .EQ. '') Tmbf2_3(i__) = rNULL__
if(ar__4(56:65) .EQ. '') Vlsr2_3(i__) = rNULL__
if(ar__4(68:77) .EQ. '') Tmbs1_3(i__) = rNULL__
if(ar__4(80:89) .EQ. '') Tmbs2_3(i__) = rNULL__
c ..............Just test output...........
write(6,5)
+ Vlsr_3(i__),Tmb_3(i__),Tmbf1_3(i__),Tmbf2_3(i__),Vlsr2_3(i__),
+ Tmbs1_3(i__),Tmbs2_3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'hn13c/*' ! HN^13^C spectra
C Format for file interpretation
6 format(
+ E17.13,1X,E11.5,2X,E10.5,2X,E10.5,2X,E10.5,2X,E10.5,2X,E10.5)
C Effective file loading
open(unit=1,status='old',file=
+'hn13c/*')
write(6,*) '....Loading file: hn13c/*'
do i__=1,27
read(1,'(A89)')ar__5
read(ar__5,6)
+ Vlsr_4(i__),Tmb_4(i__),Tmbf1_4(i__),Tmbf2_4(i__),Vlsr2_4(i__),
+ Tmbs1_4(i__),Tmbs2_4(i__)
if(ar__5(32:41) .EQ. '') Tmbf1_4(i__) = rNULL__
if(ar__5(44:53) .EQ. '') Tmbf2_4(i__) = rNULL__
if(ar__5(56:65) .EQ. '') Vlsr2_4(i__) = rNULL__
if(ar__5(68:77) .EQ. '') Tmbs1_4(i__) = rNULL__
if(ar__5(80:89) .EQ. '') Tmbs2_4(i__) = rNULL__
c ..............Just test output...........
write(6,6)
+ Vlsr_4(i__),Tmb_4(i__),Tmbf1_4(i__),Tmbf2_4(i__),Vlsr2_4(i__),
+ Tmbs1_4(i__),Tmbs2_4(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end