Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/620/A163 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-Sep-12
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/620/A163 Cores in California molecular cloud (Zhang+, 2018)
*================================================================================
*Physical properties and chemical composition of the cores in the California
*molecular cloud.
* Zhang G.-Y., Xu J.-L., Vasyunin A.I., Semenov D.A., Wang J.-J., Dib S.,
* Liu T., Liu S.-Y., Zhang C.-P., Liu X.-L., Wang K., Li D., Wu Z.-Z.,
* Yuan J.-H., Li D.-L., Gao Y.
* <Astron. Astrophys. 620, A163 (2018)>
* =2018A&A...620A.163Z (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'cores.dat' ! Parameters of 300 cores obtained from the
Herschel H_2_ column density map (table 3)
integer*4 nr__
parameter (nr__=300) ! Number of records
character*86 ar__ ! Full-size record
C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg ! (deg) Right Ascension J2000.0
real*8 DEdeg ! (deg) Declination J2000.0
C ---------------------------------- ! (position vector(s) in degrees)
character*15 Core ! Core number (CMCHerschel-NNN)
integer*4 RAh ! (h) Right ascension (J2000.0) (1)
integer*4 RAm ! (min) Right ascension (J2000.0)
real*4 RAs ! (s) Right ascension (J2000.0)
character*1 DE_ ! Declination sign (J2000.0) (1)
integer*4 DEd ! (deg) Declination (J2000.0)
integer*4 DEm ! (arcmin) Declination (J2000.0)
real*4 DEs ! (arcsec) Declination (J2000.0)
real*4 MajAxis ! (arcsec) Major axis of the ellipse (2)
real*4 MinAxis ! (arcsec) Minor axis of the ellipse
real*4 PA ! (deg) Position angle
real*4 Rad ! (pc) Core radius (3)
real*4 Td ! (K) Dust temperature (4)
real*4 n_H2 ! (cm-3) Number density (5)
real*4 M ! (Msun) Core mass
real*4 M_BE ! (Msun) Critical Bonnor-Ebert mass
character*3 Type ! Core type (6)
*Note (1): Right ascension and declination are the center positions of the core
* ellipse shape. The cores are sorted from north to south.
*Note (2): The values of the axes of the ellipse are equal to the FWHMs of the
* equivalent Gaussian.
*Note (3): The core radius is deconvolved to remove the effect of the
* telescope beam.
*Note (4): The average dust temperature in core ellipse shape.
*Note (5): The number density is calculated on the assumption that the core is
* spherical.
*Note (6): Core type as follows:
* USL = gravitationally unbound starless core
* PRE = bound prestellar core
* PRO = protostellar core
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! IRAM 30m observed positions
integer*4 nr__1
parameter (nr__1=30) ! Number of records
character*88 ar__1 ! Full-size record
C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_1 ! (deg) Right Ascension J2000.0
real*8 DEdeg_1 ! (deg) Declination J2000.0
C ---------------------------------- ! (position vector(s) in degrees)
character*6 No ! Observation number (CMC-NN)
character*1 n_No ! [*] Note on No (1)
integer*4 RAh_1 ! (h) Single-pointing observation right ascension
* (J2000.0)
integer*4 RAm_1 ! (min) Single-pointing observation right ascension
* (J2000.0)
real*4 RAs_1 ! (s) Single-pointing observation right ascension
* (J2000.0)
character*1 DE__1 ! Single-pointing observation declination sign
* (J2000.0)
integer*4 DEd_1 ! (deg) Single-pointing observation declination
* (J2000.0)
integer*4 DEm_1 ! (arcmin) Single-pointing observation declination
* (J2000.0)
real*4 DEs_1 ! (arcsec) Single-pointing observation declination
* (J2000.0)
real*4 Td_1 ! (K) Dust average temperature in one beam
* (29", IRAM 30m|86GHz)
real*4 SH2 ! (10+21cm-2) H_2_ average column density in one beam
* ({SIGMA}H_2_)
real*4 R ! (pc) ? Radius
real*4 n_H2_1 ! (10+5cm-3) ? Number density
real*4 M_1 ! (Msun) ? Herschel core mass
real*4 MV ! (Msun) ? Virial mass
real*4 M_BE_1 ! (Msun) ? Critical Bonnor-Ebert mass
character*6 Type_1 ! Core type (2)
character*15 Core_1 ! Core number (CMCHerschel-NNN)
*Note (1): CMC-2 is galaxy 3C 111.
*Note (2): Core type as follows:
* PRE = gravitationally bound prestellar core
* PRO = protostellar core
* REF = reference position that is offset from the cores
* Galaxy = galaxy
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'h13cop.dat' ! Properties of H^13^CO^+^(1-0) (table 4)
integer*4 nr__2
parameter (nr__2=30) ! Number of records
character*62 ar__2 ! Full-size record
character*6 No_1 ! Observation number (CMC-NN)
real*4 Tmb ! (K) ? Corrected main-beam temperature
real*4 e_Tmb ! (K) ? rms uncertainty on Tmb
real*4 I ! (K.km/s) ? Integrated main-beam temperature
real*4 e_I ! (K.km/s) ? rms uncertainty on I
real*4 FWHM ! (km/s) ? Full width at half-maximum of the
* Gaussian fitting profile
real*4 e_FWHM ! (km/s) ? rms uncertainty on FWHM
real*4 Vlsr ! (km/s) ? Local standard of rest velocity
real*4 e_Vlsr ! (km/s) ? rms uncertainty on Vlsr
real*4 depth ! ? Optical depth (1)
real*4 N ! (10+12cm-2) ? Column density
real*4 X ! (10-11) ? Abundance
*Note (1): We assumed a constant abundance ratio of 50 for [C/^13^C] in the CMC.
* The optical depths of H^13^CO^+^ and HCO^+^ and HN^13^C and HNC were obtained
* by comparing the measured brightness temperatures.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hn13c.dat' ! Properties of HN^13^C (1-0) (table 4)
integer*4 nr__3
parameter (nr__3=30) ! Number of records
character*62 ar__3 ! Full-size record
character*6 No_2 ! Observation number (CMC-NN)
real*4 Tmb_1 ! (K) ? Corrected main-beam temperature
real*4 e_Tmb_1 ! (K) ? rms uncertainty on Tmb
real*4 I_1 ! (K.km/s) ? Integrated main-beam temperature
real*4 e_I_1 ! (K.km/s) ? rms uncertainty on I
real*4 FWHM_1 ! (km/s) ? Full width at half-maximum of the
* Gaussian fitting profile
real*4 e_FWHM_1 ! (km/s) ? rms uncertainty on FWHM
real*4 Vlsr_1 ! (km/s) ? Local standard of rest velocity
real*4 e_Vlsr_1 ! (km/s) ? rms uncertainty on Vlsr
real*4 depth_1 ! ? Optical depth (1)
real*4 N_1 ! (10+12cm-2) ? Column density
real*4 X_1 ! (10-11) ? Abundance
*Note (1): We assumed a constant abundance ratio of 50 for [C/^13^C] in the CMC.
* The optical depths of H^13^CO^+^ and HCO^+^ and HN^13^C and HNC were obtained
* by comparing the measured brightness temperatures.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'n2hp.dat' ! Properties of N_2_H^+^ (1-0) (table 4)
integer*4 nr__4
parameter (nr__4=30) ! Number of records
character*69 ar__4 ! Full-size record
character*6 No_3 ! Observation number (CMC-NN)
real*4 Tmb_2 ! (K) ? Corrected main-beam temperature
real*4 e_Tmb_2 ! (K) ? rms uncertainty on Tmb
real*4 I_2 ! (K.km/s) ? Integrated main-beam temperature
real*4 e_I_2 ! (K.km/s) ? rms uncertainty on I
real*4 FWHM_2 ! (km/s) ? Full width at half-maximum of the
* Gaussian fitting profile
real*4 e_FWHM_2 ! (km/s) ? rms uncertainty on FWHM
real*4 Vlsr_2 ! (km/s) ? Local standard of rest velocity
real*4 e_Vlsr_2 ! (km/s) ? rms uncertainty on Vlsr
real*4 depth_2 ! ? Optical depth (1)
real*4 e_depth ! ? rms uncertainty on depth (1)
real*4 N_2 ! (10+13cm-2) ? Column density (1)
real*4 X_2 ! (10-10) ? Abundance
*Note (1): Optical depth and column density for N_2_H^+^ are estimated by its
* component JF_1_F=(101-012).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'c2h.dat' ! Properties of C_2_H (1-0) (table 4)
integer*4 nr__5
parameter (nr__5=30) ! Number of records
character*70 ar__5 ! Full-size record
character*6 No_4 ! Observation number (CMC-NN)
real*4 Tmb_3 ! (K) ? Corrected main-beam temperature
real*4 e_Tmb_3 ! (K) ? rms uncertainty on Tmb
real*4 I_3 ! (K.km/s) ? Integrated main-beam temperature
real*4 e_I_3 ! (K.km/s) ? rms uncertainty on I
real*4 FWHM_3 ! (km/s) ? Full width at half-maximum of the
* Gaussian fitting profile
real*4 e_FWHM_3 ! (km/s) ? rms uncertainty on FWHM
real*4 Vlsr_3 ! (km/s) ? Local standard of rest velocity
real*4 e_Vlsr_3 ! (km/s) ? rms uncertainty on Vlsr
real*4 depth_3 ! ? Optical depth (1)
real*4 e_depth_1 ! ? rms uncertainty on depth
real*4 N_3 ! (10+14cm-2) ? Column density (1)
real*4 X_3 ! (10-9) ? Abundance
*Note (1): Optical depth for C_2_H is estimated by its main component
* JF=(3/2,2-1/2,1) with the HFS method in the CLASS software. Column density
* is also estimated by its main component.
* Optical depth and column density for HCN are estimated by its main component
* JF=(12-01).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hcn.dat' ! Properties of HCN (1-0) (table 4)
integer*4 nr__6
parameter (nr__6=30) ! Number of records
character*70 ar__6 ! Full-size record
character*6 No_5 ! Observation number (CMC-NN)
real*4 Tmb_4 ! (K) ? Corrected main-beam temperature
real*4 e_Tmb_4 ! (K) ? rms uncertainty on Tmb
real*4 I_4 ! (K.km/s) ? Integrated main-beam temperature
real*4 e_I_4 ! (K.km/s) ? rms uncertainty on I
real*4 FWHM_4 ! (km/s) ? Full width at half-maximum of the
* Gaussian fitting profile
real*4 e_FWHM_4 ! (km/s) ? rms uncertainty on FWHM
real*4 Vlsr_4 ! (km/s) ? Local standard of rest velocity
real*4 e_Vlsr_4 ! (km/s) ? rms uncertainty on Vlsr
real*4 depth_4 ! ? Optical depth (1)
real*4 e_depth_2 ! ? rms uncertainty on depth
real*4 N_4 ! (10+14cm-2) ? Column density (1)
real*4 X_4 ! (10-9) ? Abundance
*Note (1): Optical depth for C_2_H is estimated by its main component
* JF=(3/2,2-1/2,1) with the HFS method in the CLASS software. Column density
* is also estimated by its main component.
* Optical depth and column density for HCN are estimated by its main component
* JF=(12-01).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hcop.dat' ! Properties of HCO^+^ (1-0) (table 4)
integer*4 nr__7
parameter (nr__7=30) ! Number of records
character*70 ar__7 ! Full-size record
character*6 No_6 ! Observation number (CMC-NN)
real*4 Tmb_5 ! (K) ? Corrected main-beam temperature
real*4 e_Tmb_5 ! (K) ? rms uncertainty on Tmb
real*4 I_5 ! (K.km/s) ? Integrated main-beam temperature
real*4 e_I_5 ! (K.km/s) ? rms uncertainty on I
real*4 FWHM_5 ! (km/s) ? Full width at half-maximum of the
* Gaussian fitting profile
real*4 e_FWHM_5 ! (km/s) ? rms uncertainty on FWHM
real*4 Vlsr_5 ! (km/s) ? Local standard of rest velocity
real*4 e_Vlsr_5 ! (km/s) ? rms uncertainty on Vlsr
real*4 depth_5 ! ? Optical depth (1)
real*4 e_depth_3 ! ? rms uncertainty on depth
real*4 N_5 ! (10+14cm-2) ? Column density (1)
real*4 X_5 ! (10-9) ? Abundance
*Note (1): Optical depth for C_2_H is estimated by its main component
* JF=(3/2,2-1/2,1) with the HFS method in the CLASS software. Column density
* is also estimated by its main component.
* Optical depth and column density for HCN are estimated by its main component
* JF=(12-01).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hnc.dat' ! Properties of HNC (1-0) (table 4)
integer*4 nr__8
parameter (nr__8=30) ! Number of records
character*70 ar__8 ! Full-size record
character*6 No_7 ! Observation number (CMC-NN)
real*4 Tmb_6 ! (K) ? Corrected main-beam temperature
real*4 e_Tmb_6 ! (K) ? rms uncertainty on Tmb
real*4 I_6 ! (K.km/s) ? Integrated main-beam temperature
real*4 e_I_6 ! (K.km/s) ? rms uncertainty on I
real*4 FWHM_6 ! (km/s) ? Full width at half-maximum of the
* Gaussian fitting profile
real*4 e_FWHM_6 ! (km/s) ? rms uncertainty on FWHM
real*4 Vlsr_6 ! (km/s) ? Local standard of rest velocity
real*4 e_Vlsr_6 ! (km/s) ? rms uncertainty on Vlsr
real*4 depth_6 ! ? Optical depth (1)
real*4 e_depth_4 ! ? rms uncertainty on depth
real*4 N_6 ! (10+14cm-2) ? Column density (1)
real*4 X_6 ! (10-9) ? Abundance
*Note (1): Optical depth for C_2_H is estimated by its main component
* JF=(3/2,2-1/2,1) with the HFS method in the CLASS software. Column density
* is also estimated by its main component.
* Optical depth and column density for HCN are estimated by its main component
* JF=(12-01).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'c18o.dat' ! Properties of C^18^O (1-0) (table 4)
integer*4 nr__9
parameter (nr__9=30) ! Number of records
character*58 ar__9 ! Full-size record
character*6 No_8 ! Observation number (CMC-NN)
real*4 Tmb_7 ! (K) ? Corrected main-beam temperature
real*4 e_Tmb_7 ! (K) ? rms uncertainty on Tmb
real*4 I_7 ! (K.km/s) ? Integrated main-beam temperature
real*4 e_I_7 ! (K.km/s) ? rms uncertainty on I
real*4 FWHM_7 ! (km/s) ? Full width at half-maximum of the
* Gaussian fitting profile
real*4 e_FWHM_7 ! (km/s) ? rms uncertainty on FWHM
real*4 Vlsr_7 ! (km/s) ? Local standard of rest velocity
real*4 e_Vlsr_7 ! (km/s) ? rms uncertainty on Vlsr
real*4 N_7 ! (10+15cm-2) ? Column density
real*4 X_7 ! (10-8) ? Abundance
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for '13co.dat' ! Properties of ^13^CO (1-0) (table 4)
integer*4 nr__10
parameter (nr__10=30) ! Number of records
character*59 ar__10 ! Full-size record
character*6 No_9 ! Observation number (CMC-NN)
real*4 Tmb_8 ! (K) Corrected main-beam temperature
real*4 e_Tmb_8 ! (K) rms uncertainty on Tmb
real*4 I_8 ! (K.km/s) Integrated main-beam temperature
real*4 e_I_8 ! (K.km/s) rms uncertainty on I
real*4 FWHM_8 ! (km/s) Full width at half-maximum of the
* Gaussian fitting profile
real*4 e_FWHM_8 ! (km/s) rms uncertainty on FWHM
real*4 Vlsr_8 ! (km/s) Local standard of rest velocity
real*4 e_Vlsr_8 ! (km/s) rms uncertainty on Vlsr
real*4 N_8 ! (10+16cm-2) Column density
real*4 X_8 ! (10-6) Abundance
C=============================================================================
C Loading file 'cores.dat' ! Parameters of 300 cores obtained from the
* Herschel H_2_ column density map (table 3)
C Format for file interpretation
1 format(
+ A15,1X,I2,1X,I2,1X,F6.3,1X,A1,I2,1X,I2,1X,F5.2,1X,F5.1,1X,
+ F4.1,1X,F5.1,1X,F4.2,1X,F4.1,1X,F4.1,1X,F4.1,1X,F3.1,1X,A3)
C Effective file loading
open(unit=1,status='old',file=
+'cores.dat')
write(6,*) '....Loading file: cores.dat'
do i__=1,300
read(1,'(A86)')ar__
read(ar__,1)
+ Core,RAh,RAm,RAs,DE_,DEd,DEm,DEs,MajAxis,MinAxis,PA,Rad,Td,
+ n_H2,M,M_BE,Type
RAdeg = rNULL__
DEdeg = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh .GT. -180) RAdeg=RAh*15.
if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
if(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
if(DEd .GE. 0) DEdeg=DEd
if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c ..............Just test output...........
write(6,1)
+ Core,RAh,RAm,RAs,DE_,DEd,DEm,DEs,MajAxis,MinAxis,PA,Rad,Td,
+ n_H2,M,M_BE,Type
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table1.dat' ! IRAM 30m observed positions
C Format for file interpretation
2 format(
+ A6,A1,I2,1X,I2,1X,F6.3,1X,A1,I2,1X,I2,1X,F5.2,1X,F4.1,1X,F4.1,
+ 1X,F4.2,1X,F4.1,1X,F4.1,1X,F3.1,1X,F3.1,1X,A6,1X,A15)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,30
read(1,'(A88)')ar__1
read(ar__1,2)
+ No,n_No,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,Td_1,SH2,R,
+ n_H2_1,M_1,MV,M_BE_1,Type_1,Core_1
if(ar__1(44:47) .EQ. '') R = rNULL__
if(ar__1(49:52) .EQ. '') n_H2_1 = rNULL__
if(ar__1(54:57) .EQ. '') M_1 = rNULL__
if(ar__1(59:61) .EQ. '') MV = rNULL__
if(ar__1(63:65) .EQ. '') M_BE_1 = rNULL__
RAdeg_1 = rNULL__
DEdeg_1 = rNULL__
c Derive coordinates RAdeg_1 and DEdeg_1 from input data
c (RAdeg_1 and DEdeg_1 are set to rNULL__ when unknown)
if(RAh_1 .GT. -180) RAdeg_1=RAh_1*15.
if(RAm_1 .GT. -180) RAdeg_1=RAdeg_1+RAm_1/4.
if(RAs_1 .GT. -180) RAdeg_1=RAdeg_1+RAs_1/240.
if(DEd_1 .GE. 0) DEdeg_1=DEd_1
if(DEm_1 .GE. 0) DEdeg_1=DEdeg_1+DEm_1/60.
if(DEs_1 .GE. 0) DEdeg_1=DEdeg_1+DEs_1/3600.
if(DE__1.EQ.'-'.AND.DEdeg_1.GE.0) DEdeg_1=-DEdeg_1
c ..............Just test output...........
write(6,2)
+ No,n_No,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,Td_1,SH2,R,
+ n_H2_1,M_1,MV,M_BE_1,Type_1,Core_1
write(6,'(6H Pos: 2F8.4)') RAdeg_1,DEdeg_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'h13cop.dat' ! Properties of H^13^CO^+^(1-0) (table 4)
C Format for file interpretation
3 format(
+ A6,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,
+ F4.2,1X,F4.2,1X,F4.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'h13cop.dat')
write(6,*) '....Loading file: h13cop.dat'
do i__=1,30
read(1,'(A62)')ar__2
read(ar__2,3)
+ No_1,Tmb,e_Tmb,I,e_I,FWHM,e_FWHM,Vlsr,e_Vlsr,depth,N,X
if(ar__2(8:11) .EQ. '') Tmb = rNULL__
if(ar__2(13:16) .EQ. '') e_Tmb = rNULL__
if(ar__2(18:21) .EQ. '') I = rNULL__
if(ar__2(23:26) .EQ. '') e_I = rNULL__
if(ar__2(28:31) .EQ. '') FWHM = rNULL__
if(ar__2(33:36) .EQ. '') e_FWHM = rNULL__
if(ar__2(38:42) .EQ. '') Vlsr = rNULL__
if(ar__2(44:47) .EQ. '') e_Vlsr = rNULL__
if(ar__2(49:52) .EQ. '') depth = rNULL__
if(ar__2(54:57) .EQ. '') N = rNULL__
if(ar__2(59:62) .EQ. '') X = rNULL__
c ..............Just test output...........
write(6,3)
+ No_1,Tmb,e_Tmb,I,e_I,FWHM,e_FWHM,Vlsr,e_Vlsr,depth,N,X
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'hn13c.dat' ! Properties of HN^13^C (1-0) (table 4)
C Format for file interpretation
4 format(
+ A6,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,
+ F4.2,1X,F4.2,1X,F4.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'hn13c.dat')
write(6,*) '....Loading file: hn13c.dat'
do i__=1,30
read(1,'(A62)')ar__3
read(ar__3,4)
+ No_2,Tmb_1,e_Tmb_1,I_1,e_I_1,FWHM_1,e_FWHM_1,Vlsr_1,e_Vlsr_1,
+ depth_1,N_1,X_1
if(ar__3(8:11) .EQ. '') Tmb_1 = rNULL__
if(ar__3(13:16) .EQ. '') e_Tmb_1 = rNULL__
if(ar__3(18:21) .EQ. '') I_1 = rNULL__
if(ar__3(23:26) .EQ. '') e_I_1 = rNULL__
if(ar__3(28:31) .EQ. '') FWHM_1 = rNULL__
if(ar__3(33:36) .EQ. '') e_FWHM_1 = rNULL__
if(ar__3(38:42) .EQ. '') Vlsr_1 = rNULL__
if(ar__3(44:47) .EQ. '') e_Vlsr_1 = rNULL__
if(ar__3(49:52) .EQ. '') depth_1 = rNULL__
if(ar__3(54:57) .EQ. '') N_1 = rNULL__
if(ar__3(59:62) .EQ. '') X_1 = rNULL__
c ..............Just test output...........
write(6,4)
+ No_2,Tmb_1,e_Tmb_1,I_1,e_I_1,FWHM_1,e_FWHM_1,Vlsr_1,e_Vlsr_1,
+ depth_1,N_1,X_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'n2hp.dat' ! Properties of N_2_H^+^ (1-0) (table 4)
C Format for file interpretation
5 format(
+ A6,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F6.2,1X,
+ F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'n2hp.dat')
write(6,*) '....Loading file: n2hp.dat'
do i__=1,30
read(1,'(A69)')ar__4
read(ar__4,5)
+ No_3,Tmb_2,e_Tmb_2,I_2,e_I_2,FWHM_2,e_FWHM_2,Vlsr_2,e_Vlsr_2,
+ depth_2,e_depth,N_2,X_2
if(ar__4(8:11) .EQ. '') Tmb_2 = rNULL__
if(ar__4(13:16) .EQ. '') e_Tmb_2 = rNULL__
if(ar__4(18:21) .EQ. '') I_2 = rNULL__
if(ar__4(23:26) .EQ. '') e_I_2 = rNULL__
if(ar__4(28:31) .EQ. '') FWHM_2 = rNULL__
if(ar__4(33:37) .EQ. '') e_FWHM_2 = rNULL__
if(ar__4(39:44) .EQ. '') Vlsr_2 = rNULL__
if(ar__4(46:49) .EQ. '') e_Vlsr_2 = rNULL__
if(ar__4(51:54) .EQ. '') depth_2 = rNULL__
if(ar__4(56:59) .EQ. '') e_depth = rNULL__
if(ar__4(61:64) .EQ. '') N_2 = rNULL__
if(ar__4(66:69) .EQ. '') X_2 = rNULL__
c ..............Just test output...........
write(6,5)
+ No_3,Tmb_2,e_Tmb_2,I_2,e_I_2,FWHM_2,e_FWHM_2,Vlsr_2,e_Vlsr_2,
+ depth_2,e_depth,N_2,X_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'c2h.dat' ! Properties of C_2_H (1-0) (table 4)
C Format for file interpretation
6 format(
+ A6,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,
+ F4.2,1X,F5.2,1X,F5.2,1X,F4.2,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'c2h.dat')
write(6,*) '....Loading file: c2h.dat'
do i__=1,30
read(1,'(A70)')ar__5
read(ar__5,6)
+ No_4,Tmb_3,e_Tmb_3,I_3,e_I_3,FWHM_3,e_FWHM_3,Vlsr_3,e_Vlsr_3,
+ depth_3,e_depth_1,N_3,X_3
if(ar__5(8:11) .EQ. '') Tmb_3 = rNULL__
if(ar__5(13:16) .EQ. '') e_Tmb_3 = rNULL__
if(ar__5(18:21) .EQ. '') I_3 = rNULL__
if(ar__5(23:26) .EQ. '') e_I_3 = rNULL__
if(ar__5(28:31) .EQ. '') FWHM_3 = rNULL__
if(ar__5(33:36) .EQ. '') e_FWHM_3 = rNULL__
if(ar__5(38:42) .EQ. '') Vlsr_3 = rNULL__
if(ar__5(44:47) .EQ. '') e_Vlsr_3 = rNULL__
if(ar__5(49:53) .EQ. '') depth_3 = rNULL__
if(ar__5(55:59) .EQ. '') e_depth_1 = rNULL__
if(ar__5(61:64) .EQ. '') N_3 = rNULL__
if(ar__5(66:70) .EQ. '') X_3 = rNULL__
c ..............Just test output...........
write(6,6)
+ No_4,Tmb_3,e_Tmb_3,I_3,e_I_3,FWHM_3,e_FWHM_3,Vlsr_3,e_Vlsr_3,
+ depth_3,e_depth_1,N_3,X_3
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'hcn.dat' ! Properties of HCN (1-0) (table 4)
C Format for file interpretation
7 format(
+ A6,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,
+ F4.2,1X,F5.2,1X,F5.2,1X,F4.2,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'hcn.dat')
write(6,*) '....Loading file: hcn.dat'
do i__=1,30
read(1,'(A70)')ar__6
read(ar__6,7)
+ No_5,Tmb_4,e_Tmb_4,I_4,e_I_4,FWHM_4,e_FWHM_4,Vlsr_4,e_Vlsr_4,
+ depth_4,e_depth_2,N_4,X_4
if(ar__6(8:11) .EQ. '') Tmb_4 = rNULL__
if(ar__6(13:16) .EQ. '') e_Tmb_4 = rNULL__
if(ar__6(18:21) .EQ. '') I_4 = rNULL__
if(ar__6(23:26) .EQ. '') e_I_4 = rNULL__
if(ar__6(28:31) .EQ. '') FWHM_4 = rNULL__
if(ar__6(33:36) .EQ. '') e_FWHM_4 = rNULL__
if(ar__6(38:42) .EQ. '') Vlsr_4 = rNULL__
if(ar__6(44:47) .EQ. '') e_Vlsr_4 = rNULL__
if(ar__6(49:53) .EQ. '') depth_4 = rNULL__
if(ar__6(55:59) .EQ. '') e_depth_2 = rNULL__
if(ar__6(61:64) .EQ. '') N_4 = rNULL__
if(ar__6(66:70) .EQ. '') X_4 = rNULL__
c ..............Just test output...........
write(6,7)
+ No_5,Tmb_4,e_Tmb_4,I_4,e_I_4,FWHM_4,e_FWHM_4,Vlsr_4,e_Vlsr_4,
+ depth_4,e_depth_2,N_4,X_4
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'hcop.dat' ! Properties of HCO^+^ (1-0) (table 4)
C Format for file interpretation
8 format(
+ A6,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,
+ F4.2,1X,F5.2,1X,F5.2,1X,F4.2,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'hcop.dat')
write(6,*) '....Loading file: hcop.dat'
do i__=1,30
read(1,'(A70)')ar__7
read(ar__7,8)
+ No_6,Tmb_5,e_Tmb_5,I_5,e_I_5,FWHM_5,e_FWHM_5,Vlsr_5,e_Vlsr_5,
+ depth_5,e_depth_3,N_5,X_5
if(ar__7(8:11) .EQ. '') Tmb_5 = rNULL__
if(ar__7(13:16) .EQ. '') e_Tmb_5 = rNULL__
if(ar__7(18:21) .EQ. '') I_5 = rNULL__
if(ar__7(23:26) .EQ. '') e_I_5 = rNULL__
if(ar__7(28:31) .EQ. '') FWHM_5 = rNULL__
if(ar__7(33:36) .EQ. '') e_FWHM_5 = rNULL__
if(ar__7(38:42) .EQ. '') Vlsr_5 = rNULL__
if(ar__7(44:47) .EQ. '') e_Vlsr_5 = rNULL__
if(ar__7(49:53) .EQ. '') depth_5 = rNULL__
if(ar__7(55:59) .EQ. '') e_depth_3 = rNULL__
if(ar__7(61:64) .EQ. '') N_5 = rNULL__
if(ar__7(66:70) .EQ. '') X_5 = rNULL__
c ..............Just test output...........
write(6,8)
+ No_6,Tmb_5,e_Tmb_5,I_5,e_I_5,FWHM_5,e_FWHM_5,Vlsr_5,e_Vlsr_5,
+ depth_5,e_depth_3,N_5,X_5
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'hnc.dat' ! Properties of HNC (1-0) (table 4)
C Format for file interpretation
9 format(
+ A6,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,
+ F4.2,1X,F5.2,1X,F5.2,1X,F4.2,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'hnc.dat')
write(6,*) '....Loading file: hnc.dat'
do i__=1,30
read(1,'(A70)')ar__8
read(ar__8,9)
+ No_7,Tmb_6,e_Tmb_6,I_6,e_I_6,FWHM_6,e_FWHM_6,Vlsr_6,e_Vlsr_6,
+ depth_6,e_depth_4,N_6,X_6
if(ar__8(8:11) .EQ. '') Tmb_6 = rNULL__
if(ar__8(13:16) .EQ. '') e_Tmb_6 = rNULL__
if(ar__8(18:21) .EQ. '') I_6 = rNULL__
if(ar__8(23:26) .EQ. '') e_I_6 = rNULL__
if(ar__8(28:31) .EQ. '') FWHM_6 = rNULL__
if(ar__8(33:36) .EQ. '') e_FWHM_6 = rNULL__
if(ar__8(38:42) .EQ. '') Vlsr_6 = rNULL__
if(ar__8(44:47) .EQ. '') e_Vlsr_6 = rNULL__
if(ar__8(49:53) .EQ. '') depth_6 = rNULL__
if(ar__8(55:59) .EQ. '') e_depth_4 = rNULL__
if(ar__8(61:64) .EQ. '') N_6 = rNULL__
if(ar__8(66:70) .EQ. '') X_6 = rNULL__
c ..............Just test output...........
write(6,9)
+ No_7,Tmb_6,e_Tmb_6,I_6,e_I_6,FWHM_6,e_FWHM_6,Vlsr_6,e_Vlsr_6,
+ depth_6,e_depth_4,N_6,X_6
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'c18o.dat' ! Properties of C^18^O (1-0) (table 4)
C Format for file interpretation
10 format(
+ A6,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,
+ F4.2,1X,F4.2,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'c18o.dat')
write(6,*) '....Loading file: c18o.dat'
do i__=1,30
read(1,'(A58)')ar__9
read(ar__9,10)
+ No_8,Tmb_7,e_Tmb_7,I_7,e_I_7,FWHM_7,e_FWHM_7,Vlsr_7,e_Vlsr_7,
+ N_7,X_7
if(ar__9(8:11) .EQ. '') Tmb_7 = rNULL__
if(ar__9(13:16) .EQ. '') e_Tmb_7 = rNULL__
if(ar__9(18:21) .EQ. '') I_7 = rNULL__
if(ar__9(23:26) .EQ. '') e_I_7 = rNULL__
if(ar__9(28:31) .EQ. '') FWHM_7 = rNULL__
if(ar__9(33:36) .EQ. '') e_FWHM_7 = rNULL__
if(ar__9(38:42) .EQ. '') Vlsr_7 = rNULL__
if(ar__9(44:47) .EQ. '') e_Vlsr_7 = rNULL__
if(ar__9(49:52) .EQ. '') N_7 = rNULL__
if(ar__9(54:58) .EQ. '') X_7 = rNULL__
c ..............Just test output...........
write(6,10)
+ No_8,Tmb_7,e_Tmb_7,I_7,e_I_7,FWHM_7,e_FWHM_7,Vlsr_7,e_Vlsr_7,
+ N_7,X_7
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file '13co.dat' ! Properties of ^13^CO (1-0) (table 4)
C Format for file interpretation
11 format(
+ A6,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,
+ F4.2,1X,F4.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'13co.dat')
write(6,*) '....Loading file: 13co.dat'
do i__=1,30
read(1,'(A59)')ar__10
read(ar__10,11)
+ No_9,Tmb_8,e_Tmb_8,I_8,e_I_8,FWHM_8,e_FWHM_8,Vlsr_8,e_Vlsr_8,
+ N_8,X_8
c ..............Just test output...........
write(6,11)
+ No_9,Tmb_8,e_Tmb_8,I_8,e_I_8,FWHM_8,e_FWHM_8,Vlsr_8,e_Vlsr_8,
+ N_8,X_8
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end