Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/609/A129 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 2024-Apr-17
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 ! Type (HMPO, HMSC or UCHII)
character*13 Source ! Source name
character*1 n_Source ! [w] Note on Source (1)
character*2 l_N_H15NC ! [<= ] Limit flag on N(H15NC)
integer*4 N_H15NC ! (10+10cm-2) Total H^15^NC column density
* (beam-averaged)
integer*4 e_N_H15NC ! (10+10cm-2) ? rms uncertainty on N(H15NC) (2)
character*1 n_N_H15NC ! [ut] Note on N(H15NC) (1)
integer*4 DN1 ! (10+10cm-2) ?=- Error on N(H15NC) without considering
* the calibration error
real*4 N_HN13C ! (10+11cm-2) Total HN^13^C column density
* (beam-averaged)
real*4 e_N_HN13C ! (10+11cm-2) rms uncertainty on N(HN13C) (2)
real*4 DN2 ! (10+11cm-2) Error on N(HN13C) without considering the
* calibration error
real*4 N_DNC ! (10+11cm-2) Total DNC column density (beam-averaged)
real*4 e_N_DNC ! (10+11cm-2) rms uncertainty on N(DNC) (2)
character*2 l_HNC_H15NC ! [>= ] Limit flag on HNC/H15NC
integer*4 HNC_H15NC ! HNC/H^15^NC isotopic ratio (3)
integer*4 e_HNC_H15NC ! ? rms uncertainty on HNC/H15NC (2)
integer*4 HNC_DNC ! HNC/DNC isotopic ratio (4)
integer*4 e_HNC_DNC ! rms uncertainty on HNC/DNC (2)
real*4 DNC_HNC ! DNC/HNC isotopic ratio
real*4 e_DNC_HNC ! rms uncertainty on DNC/HNC (2)
character*2 l_N_HC15N ! [<= ] Limit flag on N(HC15N)
integer*4 N_HC15N ! (10+10cm-2) Total N(HC15N) column density
* (beam-averaged)
integer*4 e_N_HC15N ! (10+10cm-2) ? rms uncertainty on N(HC15N)
character*1 n_N_HC15N ! [ut] Note on N(HC15N) (1)
integer*4 DN3 ! (10+10cm-2) ?=- Error on N(HC15N) without considering
* the calibration error
integer*4 N_H13CN ! (10+11cm-2) Total H^13^CN column density
* (beam-averaged)
integer*4 e_N_H13CN ! (10+11cm-2) rms uncertainty on N(H13CN) (2)
real*4 DN4 ! (10+11cm-2) Error on N(H13CN) without considering the
* calibration error
character*2 l_HCN_HC15N ! [>= ] Limit flag on HCN/HC15N
integer*4 HCN_HC15N ! HCN/HC15N isotopic ratio (3)
integer*4 e_HCN_HC15N ! ? rms uncertainty on HCN/HC15N (2)
integer*4 Tk ! (K) Kinetic temperature of the clump (5)
character*1 n_Tk ! [*] Note on Tk (6)
character*12 FName ! 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 ! (km/s) LSR velocity
real*4 Tmb ! (K) Beam brightness temperature
real*4 Tmbf1 ! ? Tmb computed in the fit
real*4 Tmbf2 ! ? Tmb computed in the fit
real*4 Vlsr2 ! ? LSR velocity for the possible second
* component
real*4 Tmbs1 ! ? Beam brightness temperature for the possible
* second component
real*4 Tmbs2 ! ? 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 ! (km/s) LSR velocity
real*4 Tmb_1 ! (K) Beam brightness temperature
real*4 Tmbf1_1 ! ? Tmb computed in the fit
real*4 Tmbf2_1 ! ? Tmb computed in the fit
real*4 Vlsr2_1 ! ? LSR velocity for the possible second
* component
real*4 Tmbs1_1 ! ? Beam brightness temperature for the possible
* second component
real*4 Tmbs2_1 ! ? 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 ! (km/s) LSR velocity
real*4 Tmb_2 ! (K) Beam brightness temperature
real*4 Tmbf1_2 ! ? Tmb computed in the fit
real*4 Tmbf2_2 ! ? Tmb computed in the fit
real*4 Vlsr2_2 ! ? LSR velocity for the possible second
* component
real*4 Tmbs1_2 ! ? Beam brightness temperature for the possible
* second component
real*4 Tmbs2_2 ! ? 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 ! (km/s) LSR velocity
real*4 Tmb_3 ! (K) Beam brightness temperature
real*4 Tmbf1_3 ! ? Tmb computed in the fit
real*4 Tmbf2_3 ! ? Tmb computed in the fit
real*4 Vlsr2_3 ! ? LSR velocity for the possible second
* component
real*4 Tmbs1_3 ! ? Beam brightness temperature for the possible
* second component
real*4 Tmbs2_3 ! ? 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 ! (km/s) LSR velocity
real*4 Tmb_4 ! (K) Beam brightness temperature
real*4 Tmbf1_4 ! ? Tmb computed in the fit
real*4 Tmbf2_4 ! ? Tmb computed in the fit
real*4 Vlsr2_4 ! ? LSR velocity for the possible second
* component
real*4 Tmbs1_4 ! ? Beam brightness temperature for the possible
* second component
real*4 Tmbs2_4 ! ? 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,Source,n_Source,l_N_H15NC,N_H15NC,e_N_H15NC,n_N_H15NC,
+ DN1,N_HN13C,e_N_HN13C,DN2,N_DNC,e_N_DNC,l_HNC_H15NC,HNC_H15NC,
+ e_HNC_H15NC,HNC_DNC,e_HNC_DNC,DNC_HNC,e_DNC_HNC,l_N_HC15N,
+ N_HC15N,e_N_HC15N,n_N_HC15N,DN3,N_H13CN,e_N_H13CN,DN4,
+ l_HCN_HC15N,HCN_HC15N,e_HCN_HC15N,Tk,n_Tk,FName
if(ar__(27:28) .EQ. '') e_N_H15NC = iNULL__
if (DN1 .EQ. 45) DN1 = iNULL__
if(ar__(64:66) .EQ. '') e_HNC_H15NC = iNULL__
if(ar__(98:100) .EQ. '') e_N_HC15N = iNULL__
if (DN3 .EQ. 45) DN3 = iNULL__
if(ar__(123:125) .EQ. '') e_HCN_HC15N = iNULL__
c ..............Just test output...........
write(6,1)
+ Type,Source,n_Source,l_N_H15NC,N_H15NC,e_N_H15NC,n_N_H15NC,
+ DN1,N_HN13C,e_N_HN13C,DN2,N_DNC,e_N_DNC,l_HNC_H15NC,HNC_H15NC,
+ e_HNC_H15NC,HNC_DNC,e_HNC_DNC,DNC_HNC,e_DNC_HNC,l_N_HC15N,
+ N_HC15N,e_N_HC15N,n_N_HC15N,DN3,N_H13CN,e_N_H13CN,DN4,
+ l_HCN_HC15N,HCN_HC15N,e_HCN_HC15N,Tk,n_Tk,FName
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,Tmb,Tmbf1,Tmbf2,Vlsr2,Tmbs1,Tmbs2
if(ar__1(32:41) .EQ. '') Tmbf1 = rNULL__
if(ar__1(44:53) .EQ. '') Tmbf2 = rNULL__
if(ar__1(56:65) .EQ. '') Vlsr2 = rNULL__
if(ar__1(68:77) .EQ. '') Tmbs1 = rNULL__
if(ar__1(80:89) .EQ. '') Tmbs2 = rNULL__
c ..............Just test output...........
write(6,2)Vlsr,Tmb,Tmbf1,Tmbf2,Vlsr2,Tmbs1,Tmbs2
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,Tmb_1,Tmbf1_1,Tmbf2_1,Vlsr2_1,Tmbs1_1,Tmbs2_1
if(ar__2(32:41) .EQ. '') Tmbf1_1 = rNULL__
if(ar__2(44:53) .EQ. '') Tmbf2_1 = rNULL__
if(ar__2(56:65) .EQ. '') Vlsr2_1 = rNULL__
if(ar__2(68:77) .EQ. '') Tmbs1_1 = rNULL__
if(ar__2(80:89) .EQ. '') Tmbs2_1 = rNULL__
c ..............Just test output...........
write(6,3)
+ Vlsr_1,Tmb_1,Tmbf1_1,Tmbf2_1,Vlsr2_1,Tmbs1_1,Tmbs2_1
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,Tmb_2,Tmbf1_2,Tmbf2_2,Vlsr2_2,Tmbs1_2,Tmbs2_2
if(ar__3(32:41) .EQ. '') Tmbf1_2 = rNULL__
if(ar__3(44:53) .EQ. '') Tmbf2_2 = rNULL__
if(ar__3(56:65) .EQ. '') Vlsr2_2 = rNULL__
if(ar__3(68:77) .EQ. '') Tmbs1_2 = rNULL__
if(ar__3(80:89) .EQ. '') Tmbs2_2 = rNULL__
c ..............Just test output...........
write(6,4)
+ Vlsr_2,Tmb_2,Tmbf1_2,Tmbf2_2,Vlsr2_2,Tmbs1_2,Tmbs2_2
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,Tmb_3,Tmbf1_3,Tmbf2_3,Vlsr2_3,Tmbs1_3,Tmbs2_3
if(ar__4(32:41) .EQ. '') Tmbf1_3 = rNULL__
if(ar__4(44:53) .EQ. '') Tmbf2_3 = rNULL__
if(ar__4(56:65) .EQ. '') Vlsr2_3 = rNULL__
if(ar__4(68:77) .EQ. '') Tmbs1_3 = rNULL__
if(ar__4(80:89) .EQ. '') Tmbs2_3 = rNULL__
c ..............Just test output...........
write(6,5)
+ Vlsr_3,Tmb_3,Tmbf1_3,Tmbf2_3,Vlsr2_3,Tmbs1_3,Tmbs2_3
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,Tmb_4,Tmbf1_4,Tmbf2_4,Vlsr2_4,Tmbs1_4,Tmbs2_4
if(ar__5(32:41) .EQ. '') Tmbf1_4 = rNULL__
if(ar__5(44:53) .EQ. '') Tmbf2_4 = rNULL__
if(ar__5(56:65) .EQ. '') Vlsr2_4 = rNULL__
if(ar__5(68:77) .EQ. '') Tmbs1_4 = rNULL__
if(ar__5(80:89) .EQ. '') Tmbs2_4 = rNULL__
c ..............Just test output...........
write(6,6)
+ Vlsr_4,Tmb_4,Tmbf1_4,Tmbf2_4,Vlsr2_4,Tmbs1_4,Tmbs2_4
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end