Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/620/A163 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-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 (nr__) ! (deg) Right Ascension J2000.0
real*8 DEdeg (nr__) ! (deg) Declination J2000.0
C ---------------------------------- ! (position vector(s) in degrees)
character*15 Core (nr__) ! Core number (CMCHerschel-NNN)
integer*4 RAh (nr__) ! (h) Right ascension (J2000.0) (1)
integer*4 RAm (nr__) ! (min) Right ascension (J2000.0)
real*4 RAs (nr__) ! (s) Right ascension (J2000.0)
character*1 DE_ (nr__) ! Declination sign (J2000.0) (1)
integer*4 DEd (nr__) ! (deg) Declination (J2000.0)
integer*4 DEm (nr__) ! (arcmin) Declination (J2000.0)
real*4 DEs (nr__) ! (arcsec) Declination (J2000.0)
real*4 MajAxis (nr__) ! (arcsec) Major axis of the ellipse (2)
real*4 MinAxis (nr__) ! (arcsec) Minor axis of the ellipse
real*4 PA (nr__) ! (deg) Position angle
real*4 Rad (nr__) ! (pc) Core radius (3)
real*4 Td (nr__) ! (K) Dust temperature (4)
real*4 n_H2 (nr__) ! (cm-3) Number density (5)
real*4 M (nr__) ! (Msun) Core mass
real*4 M_BE (nr__) ! (Msun) Critical Bonnor-Ebert mass
character*3 Type (nr__) ! 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 (nr__1) ! (deg) Right Ascension J2000.0
real*8 DEdeg_1 (nr__1) ! (deg) Declination J2000.0
C ---------------------------------- ! (position vector(s) in degrees)
character*6 No (nr__1) ! Observation number (CMC-NN)
character*1 n_No (nr__1) ! [*] Note on No (1)
integer*4 RAh_1 (nr__1) ! (h) Single-pointing observation right ascension
* (J2000.0)
integer*4 RAm_1 (nr__1) ! (min) Single-pointing observation right ascension
* (J2000.0)
real*4 RAs_1 (nr__1) ! (s) Single-pointing observation right ascension
* (J2000.0)
character*1 DE__1 (nr__1) ! Single-pointing observation declination sign
* (J2000.0)
integer*4 DEd_1 (nr__1) ! (deg) Single-pointing observation declination
* (J2000.0)
integer*4 DEm_1 (nr__1) ! (arcmin) Single-pointing observation declination
* (J2000.0)
real*4 DEs_1 (nr__1) ! (arcsec) Single-pointing observation declination
* (J2000.0)
real*4 Td_1 (nr__1) ! (K) Dust average temperature in one beam
* (29", IRAM 30m|86GHz)
real*4 SH2 (nr__1) ! (10+21cm-2) H_2_ average column density in one beam
* ({SIGMA}H_2_)
real*4 R (nr__1) ! (pc) ? Radius
real*4 n_H2_1 (nr__1) ! (10+5cm-3) ? Number density
real*4 M_1 (nr__1) ! (Msun) ? Herschel core mass
real*4 MV (nr__1) ! (Msun) ? Virial mass
real*4 M_BE_1 (nr__1) ! (Msun) ? Critical Bonnor-Ebert mass
character*6 Type_1 (nr__1) ! Core type (2)
character*15 Core_1 (nr__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 (nr__2) ! Observation number (CMC-NN)
real*4 Tmb (nr__2) ! (K) ? Corrected main-beam temperature
real*4 e_Tmb (nr__2) ! (K) ? rms uncertainty on Tmb
real*4 I (nr__2) ! (K.km/s) ? Integrated main-beam temperature
real*4 e_I (nr__2) ! (K.km/s) ? rms uncertainty on I
real*4 FWHM (nr__2) ! (km/s) ? Full width at half-maximum of the
* Gaussian fitting profile
real*4 e_FWHM (nr__2) ! (km/s) ? rms uncertainty on FWHM
real*4 Vlsr (nr__2) ! (km/s) ? Local standard of rest velocity
real*4 e_Vlsr (nr__2) ! (km/s) ? rms uncertainty on Vlsr
real*4 depth (nr__2) ! ? Optical depth (1)
real*4 N (nr__2) ! (10+12cm-2) ? Column density
real*4 X (nr__2) ! (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 (nr__3) ! Observation number (CMC-NN)
real*4 Tmb_1 (nr__3) ! (K) ? Corrected main-beam temperature
real*4 e_Tmb_1 (nr__3) ! (K) ? rms uncertainty on Tmb
real*4 I_1 (nr__3) ! (K.km/s) ? Integrated main-beam temperature
real*4 e_I_1 (nr__3) ! (K.km/s) ? rms uncertainty on I
real*4 FWHM_1 (nr__3) ! (km/s) ? Full width at half-maximum of the
* Gaussian fitting profile
real*4 e_FWHM_1 (nr__3) ! (km/s) ? rms uncertainty on FWHM
real*4 Vlsr_1 (nr__3) ! (km/s) ? Local standard of rest velocity
real*4 e_Vlsr_1 (nr__3) ! (km/s) ? rms uncertainty on Vlsr
real*4 depth_1 (nr__3) ! ? Optical depth (1)
real*4 N_1 (nr__3) ! (10+12cm-2) ? Column density
real*4 X_1 (nr__3) ! (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 (nr__4) ! Observation number (CMC-NN)
real*4 Tmb_2 (nr__4) ! (K) ? Corrected main-beam temperature
real*4 e_Tmb_2 (nr__4) ! (K) ? rms uncertainty on Tmb
real*4 I_2 (nr__4) ! (K.km/s) ? Integrated main-beam temperature
real*4 e_I_2 (nr__4) ! (K.km/s) ? rms uncertainty on I
real*4 FWHM_2 (nr__4) ! (km/s) ? Full width at half-maximum of the
* Gaussian fitting profile
real*4 e_FWHM_2 (nr__4) ! (km/s) ? rms uncertainty on FWHM
real*4 Vlsr_2 (nr__4) ! (km/s) ? Local standard of rest velocity
real*4 e_Vlsr_2 (nr__4) ! (km/s) ? rms uncertainty on Vlsr
real*4 depth_2 (nr__4) ! ? Optical depth (1)
real*4 e_depth (nr__4) ! ? rms uncertainty on depth (1)
real*4 N_2 (nr__4) ! (10+13cm-2) ? Column density (1)
real*4 X_2 (nr__4) ! (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 (nr__5) ! Observation number (CMC-NN)
real*4 Tmb_3 (nr__5) ! (K) ? Corrected main-beam temperature
real*4 e_Tmb_3 (nr__5) ! (K) ? rms uncertainty on Tmb
real*4 I_3 (nr__5) ! (K.km/s) ? Integrated main-beam temperature
real*4 e_I_3 (nr__5) ! (K.km/s) ? rms uncertainty on I
real*4 FWHM_3 (nr__5) ! (km/s) ? Full width at half-maximum of the
* Gaussian fitting profile
real*4 e_FWHM_3 (nr__5) ! (km/s) ? rms uncertainty on FWHM
real*4 Vlsr_3 (nr__5) ! (km/s) ? Local standard of rest velocity
real*4 e_Vlsr_3 (nr__5) ! (km/s) ? rms uncertainty on Vlsr
real*4 depth_3 (nr__5) ! ? Optical depth (1)
real*4 e_depth_1 (nr__5) ! ? rms uncertainty on depth
real*4 N_3 (nr__5) ! (10+14cm-2) ? Column density (1)
real*4 X_3 (nr__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 '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 (nr__6) ! Observation number (CMC-NN)
real*4 Tmb_4 (nr__6) ! (K) ? Corrected main-beam temperature
real*4 e_Tmb_4 (nr__6) ! (K) ? rms uncertainty on Tmb
real*4 I_4 (nr__6) ! (K.km/s) ? Integrated main-beam temperature
real*4 e_I_4 (nr__6) ! (K.km/s) ? rms uncertainty on I
real*4 FWHM_4 (nr__6) ! (km/s) ? Full width at half-maximum of the
* Gaussian fitting profile
real*4 e_FWHM_4 (nr__6) ! (km/s) ? rms uncertainty on FWHM
real*4 Vlsr_4 (nr__6) ! (km/s) ? Local standard of rest velocity
real*4 e_Vlsr_4 (nr__6) ! (km/s) ? rms uncertainty on Vlsr
real*4 depth_4 (nr__6) ! ? Optical depth (1)
real*4 e_depth_2 (nr__6) ! ? rms uncertainty on depth
real*4 N_4 (nr__6) ! (10+14cm-2) ? Column density (1)
real*4 X_4 (nr__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 '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 (nr__7) ! Observation number (CMC-NN)
real*4 Tmb_5 (nr__7) ! (K) ? Corrected main-beam temperature
real*4 e_Tmb_5 (nr__7) ! (K) ? rms uncertainty on Tmb
real*4 I_5 (nr__7) ! (K.km/s) ? Integrated main-beam temperature
real*4 e_I_5 (nr__7) ! (K.km/s) ? rms uncertainty on I
real*4 FWHM_5 (nr__7) ! (km/s) ? Full width at half-maximum of the
* Gaussian fitting profile
real*4 e_FWHM_5 (nr__7) ! (km/s) ? rms uncertainty on FWHM
real*4 Vlsr_5 (nr__7) ! (km/s) ? Local standard of rest velocity
real*4 e_Vlsr_5 (nr__7) ! (km/s) ? rms uncertainty on Vlsr
real*4 depth_5 (nr__7) ! ? Optical depth (1)
real*4 e_depth_3 (nr__7) ! ? rms uncertainty on depth
real*4 N_5 (nr__7) ! (10+14cm-2) ? Column density (1)
real*4 X_5 (nr__7) ! (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 (nr__8) ! Observation number (CMC-NN)
real*4 Tmb_6 (nr__8) ! (K) ? Corrected main-beam temperature
real*4 e_Tmb_6 (nr__8) ! (K) ? rms uncertainty on Tmb
real*4 I_6 (nr__8) ! (K.km/s) ? Integrated main-beam temperature
real*4 e_I_6 (nr__8) ! (K.km/s) ? rms uncertainty on I
real*4 FWHM_6 (nr__8) ! (km/s) ? Full width at half-maximum of the
* Gaussian fitting profile
real*4 e_FWHM_6 (nr__8) ! (km/s) ? rms uncertainty on FWHM
real*4 Vlsr_6 (nr__8) ! (km/s) ? Local standard of rest velocity
real*4 e_Vlsr_6 (nr__8) ! (km/s) ? rms uncertainty on Vlsr
real*4 depth_6 (nr__8) ! ? Optical depth (1)
real*4 e_depth_4 (nr__8) ! ? rms uncertainty on depth
real*4 N_6 (nr__8) ! (10+14cm-2) ? Column density (1)
real*4 X_6 (nr__8) ! (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 (nr__9) ! Observation number (CMC-NN)
real*4 Tmb_7 (nr__9) ! (K) ? Corrected main-beam temperature
real*4 e_Tmb_7 (nr__9) ! (K) ? rms uncertainty on Tmb
real*4 I_7 (nr__9) ! (K.km/s) ? Integrated main-beam temperature
real*4 e_I_7 (nr__9) ! (K.km/s) ? rms uncertainty on I
real*4 FWHM_7 (nr__9) ! (km/s) ? Full width at half-maximum of the
* Gaussian fitting profile
real*4 e_FWHM_7 (nr__9) ! (km/s) ? rms uncertainty on FWHM
real*4 Vlsr_7 (nr__9) ! (km/s) ? Local standard of rest velocity
real*4 e_Vlsr_7 (nr__9) ! (km/s) ? rms uncertainty on Vlsr
real*4 N_7 (nr__9) ! (10+15cm-2) ? Column density
real*4 X_7 (nr__9) ! (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 (nr__10) ! Observation number (CMC-NN)
real*4 Tmb_8 (nr__10) ! (K) Corrected main-beam temperature
real*4 e_Tmb_8 (nr__10) ! (K) rms uncertainty on Tmb
real*4 I_8 (nr__10) ! (K.km/s) Integrated main-beam temperature
real*4 e_I_8 (nr__10) ! (K.km/s) rms uncertainty on I
real*4 FWHM_8 (nr__10) ! (km/s) Full width at half-maximum of the
* Gaussian fitting profile
real*4 e_FWHM_8 (nr__10) ! (km/s) rms uncertainty on FWHM
real*4 Vlsr_8 (nr__10) ! (km/s) Local standard of rest velocity
real*4 e_Vlsr_8 (nr__10) ! (km/s) rms uncertainty on Vlsr
real*4 N_8 (nr__10) ! (10+16cm-2) Column density
real*4 X_8 (nr__10) ! (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(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),MajAxis(i__),MinAxis(i__),PA(i__),Rad(i__),
+ Td(i__),n_H2(i__),M(i__),M_BE(i__),Type(i__)
RAdeg(i__) = rNULL__
DEdeg(i__) = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c ..............Just test output...........
write(6,1)
+ Core(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),MajAxis(i__),MinAxis(i__),PA(i__),Rad(i__),
+ Td(i__),n_H2(i__),M(i__),M_BE(i__),Type(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
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(i__),n_No(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__),
+ DEd_1(i__),DEm_1(i__),DEs_1(i__),Td_1(i__),SH2(i__),R(i__),
+ n_H2_1(i__),M_1(i__),MV(i__),M_BE_1(i__),Type_1(i__),
+ Core_1(i__)
if(ar__1(44:47) .EQ. '') R(i__) = rNULL__
if(ar__1(49:52) .EQ. '') n_H2_1(i__) = rNULL__
if(ar__1(54:57) .EQ. '') M_1(i__) = rNULL__
if(ar__1(59:61) .EQ. '') MV(i__) = rNULL__
if(ar__1(63:65) .EQ. '') M_BE_1(i__) = rNULL__
RAdeg_1(i__) = rNULL__
DEdeg_1(i__) = 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(i__) .GT. -180) RAdeg_1(i__)=RAh_1(i__)*15.
if(RAm_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAm_1(i__)/4.
if(RAs_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAs_1(i__)/240.
if(DEd_1(i__) .GE. 0) DEdeg_1(i__)=DEd_1(i__)
if(DEm_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEm_1(i__)/60.
if(DEs_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEs_1(i__)/3600.
if(DE__1(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__)
c ..............Just test output...........
write(6,2)
+ No(i__),n_No(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__),
+ DEd_1(i__),DEm_1(i__),DEs_1(i__),Td_1(i__),SH2(i__),R(i__),
+ n_H2_1(i__),M_1(i__),MV(i__),M_BE_1(i__),Type_1(i__),
+ Core_1(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg_1(i__),DEdeg_1(i__)
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(i__),Tmb(i__),e_Tmb(i__),I(i__),e_I(i__),FWHM(i__),
+ e_FWHM(i__),Vlsr(i__),e_Vlsr(i__),depth(i__),N(i__),X(i__)
if(ar__2(8:11) .EQ. '') Tmb(i__) = rNULL__
if(ar__2(13:16) .EQ. '') e_Tmb(i__) = rNULL__
if(ar__2(18:21) .EQ. '') I(i__) = rNULL__
if(ar__2(23:26) .EQ. '') e_I(i__) = rNULL__
if(ar__2(28:31) .EQ. '') FWHM(i__) = rNULL__
if(ar__2(33:36) .EQ. '') e_FWHM(i__) = rNULL__
if(ar__2(38:42) .EQ. '') Vlsr(i__) = rNULL__
if(ar__2(44:47) .EQ. '') e_Vlsr(i__) = rNULL__
if(ar__2(49:52) .EQ. '') depth(i__) = rNULL__
if(ar__2(54:57) .EQ. '') N(i__) = rNULL__
if(ar__2(59:62) .EQ. '') X(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ No_1(i__),Tmb(i__),e_Tmb(i__),I(i__),e_I(i__),FWHM(i__),
+ e_FWHM(i__),Vlsr(i__),e_Vlsr(i__),depth(i__),N(i__),X(i__)
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(i__),Tmb_1(i__),e_Tmb_1(i__),I_1(i__),e_I_1(i__),
+ FWHM_1(i__),e_FWHM_1(i__),Vlsr_1(i__),e_Vlsr_1(i__),
+ depth_1(i__),N_1(i__),X_1(i__)
if(ar__3(8:11) .EQ. '') Tmb_1(i__) = rNULL__
if(ar__3(13:16) .EQ. '') e_Tmb_1(i__) = rNULL__
if(ar__3(18:21) .EQ. '') I_1(i__) = rNULL__
if(ar__3(23:26) .EQ. '') e_I_1(i__) = rNULL__
if(ar__3(28:31) .EQ. '') FWHM_1(i__) = rNULL__
if(ar__3(33:36) .EQ. '') e_FWHM_1(i__) = rNULL__
if(ar__3(38:42) .EQ. '') Vlsr_1(i__) = rNULL__
if(ar__3(44:47) .EQ. '') e_Vlsr_1(i__) = rNULL__
if(ar__3(49:52) .EQ. '') depth_1(i__) = rNULL__
if(ar__3(54:57) .EQ. '') N_1(i__) = rNULL__
if(ar__3(59:62) .EQ. '') X_1(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ No_2(i__),Tmb_1(i__),e_Tmb_1(i__),I_1(i__),e_I_1(i__),
+ FWHM_1(i__),e_FWHM_1(i__),Vlsr_1(i__),e_Vlsr_1(i__),
+ depth_1(i__),N_1(i__),X_1(i__)
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(i__),Tmb_2(i__),e_Tmb_2(i__),I_2(i__),e_I_2(i__),
+ FWHM_2(i__),e_FWHM_2(i__),Vlsr_2(i__),e_Vlsr_2(i__),
+ depth_2(i__),e_depth(i__),N_2(i__),X_2(i__)
if(ar__4(8:11) .EQ. '') Tmb_2(i__) = rNULL__
if(ar__4(13:16) .EQ. '') e_Tmb_2(i__) = rNULL__
if(ar__4(18:21) .EQ. '') I_2(i__) = rNULL__
if(ar__4(23:26) .EQ. '') e_I_2(i__) = rNULL__
if(ar__4(28:31) .EQ. '') FWHM_2(i__) = rNULL__
if(ar__4(33:37) .EQ. '') e_FWHM_2(i__) = rNULL__
if(ar__4(39:44) .EQ. '') Vlsr_2(i__) = rNULL__
if(ar__4(46:49) .EQ. '') e_Vlsr_2(i__) = rNULL__
if(ar__4(51:54) .EQ. '') depth_2(i__) = rNULL__
if(ar__4(56:59) .EQ. '') e_depth(i__) = rNULL__
if(ar__4(61:64) .EQ. '') N_2(i__) = rNULL__
if(ar__4(66:69) .EQ. '') X_2(i__) = rNULL__
c ..............Just test output...........
write(6,5)
+ No_3(i__),Tmb_2(i__),e_Tmb_2(i__),I_2(i__),e_I_2(i__),
+ FWHM_2(i__),e_FWHM_2(i__),Vlsr_2(i__),e_Vlsr_2(i__),
+ depth_2(i__),e_depth(i__),N_2(i__),X_2(i__)
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(i__),Tmb_3(i__),e_Tmb_3(i__),I_3(i__),e_I_3(i__),
+ FWHM_3(i__),e_FWHM_3(i__),Vlsr_3(i__),e_Vlsr_3(i__),
+ depth_3(i__),e_depth_1(i__),N_3(i__),X_3(i__)
if(ar__5(8:11) .EQ. '') Tmb_3(i__) = rNULL__
if(ar__5(13:16) .EQ. '') e_Tmb_3(i__) = rNULL__
if(ar__5(18:21) .EQ. '') I_3(i__) = rNULL__
if(ar__5(23:26) .EQ. '') e_I_3(i__) = rNULL__
if(ar__5(28:31) .EQ. '') FWHM_3(i__) = rNULL__
if(ar__5(33:36) .EQ. '') e_FWHM_3(i__) = rNULL__
if(ar__5(38:42) .EQ. '') Vlsr_3(i__) = rNULL__
if(ar__5(44:47) .EQ. '') e_Vlsr_3(i__) = rNULL__
if(ar__5(49:53) .EQ. '') depth_3(i__) = rNULL__
if(ar__5(55:59) .EQ. '') e_depth_1(i__) = rNULL__
if(ar__5(61:64) .EQ. '') N_3(i__) = rNULL__
if(ar__5(66:70) .EQ. '') X_3(i__) = rNULL__
c ..............Just test output...........
write(6,6)
+ No_4(i__),Tmb_3(i__),e_Tmb_3(i__),I_3(i__),e_I_3(i__),
+ FWHM_3(i__),e_FWHM_3(i__),Vlsr_3(i__),e_Vlsr_3(i__),
+ depth_3(i__),e_depth_1(i__),N_3(i__),X_3(i__)
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(i__),Tmb_4(i__),e_Tmb_4(i__),I_4(i__),e_I_4(i__),
+ FWHM_4(i__),e_FWHM_4(i__),Vlsr_4(i__),e_Vlsr_4(i__),
+ depth_4(i__),e_depth_2(i__),N_4(i__),X_4(i__)
if(ar__6(8:11) .EQ. '') Tmb_4(i__) = rNULL__
if(ar__6(13:16) .EQ. '') e_Tmb_4(i__) = rNULL__
if(ar__6(18:21) .EQ. '') I_4(i__) = rNULL__
if(ar__6(23:26) .EQ. '') e_I_4(i__) = rNULL__
if(ar__6(28:31) .EQ. '') FWHM_4(i__) = rNULL__
if(ar__6(33:36) .EQ. '') e_FWHM_4(i__) = rNULL__
if(ar__6(38:42) .EQ. '') Vlsr_4(i__) = rNULL__
if(ar__6(44:47) .EQ. '') e_Vlsr_4(i__) = rNULL__
if(ar__6(49:53) .EQ. '') depth_4(i__) = rNULL__
if(ar__6(55:59) .EQ. '') e_depth_2(i__) = rNULL__
if(ar__6(61:64) .EQ. '') N_4(i__) = rNULL__
if(ar__6(66:70) .EQ. '') X_4(i__) = rNULL__
c ..............Just test output...........
write(6,7)
+ No_5(i__),Tmb_4(i__),e_Tmb_4(i__),I_4(i__),e_I_4(i__),
+ FWHM_4(i__),e_FWHM_4(i__),Vlsr_4(i__),e_Vlsr_4(i__),
+ depth_4(i__),e_depth_2(i__),N_4(i__),X_4(i__)
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(i__),Tmb_5(i__),e_Tmb_5(i__),I_5(i__),e_I_5(i__),
+ FWHM_5(i__),e_FWHM_5(i__),Vlsr_5(i__),e_Vlsr_5(i__),
+ depth_5(i__),e_depth_3(i__),N_5(i__),X_5(i__)
if(ar__7(8:11) .EQ. '') Tmb_5(i__) = rNULL__
if(ar__7(13:16) .EQ. '') e_Tmb_5(i__) = rNULL__
if(ar__7(18:21) .EQ. '') I_5(i__) = rNULL__
if(ar__7(23:26) .EQ. '') e_I_5(i__) = rNULL__
if(ar__7(28:31) .EQ. '') FWHM_5(i__) = rNULL__
if(ar__7(33:36) .EQ. '') e_FWHM_5(i__) = rNULL__
if(ar__7(38:42) .EQ. '') Vlsr_5(i__) = rNULL__
if(ar__7(44:47) .EQ. '') e_Vlsr_5(i__) = rNULL__
if(ar__7(49:53) .EQ. '') depth_5(i__) = rNULL__
if(ar__7(55:59) .EQ. '') e_depth_3(i__) = rNULL__
if(ar__7(61:64) .EQ. '') N_5(i__) = rNULL__
if(ar__7(66:70) .EQ. '') X_5(i__) = rNULL__
c ..............Just test output...........
write(6,8)
+ No_6(i__),Tmb_5(i__),e_Tmb_5(i__),I_5(i__),e_I_5(i__),
+ FWHM_5(i__),e_FWHM_5(i__),Vlsr_5(i__),e_Vlsr_5(i__),
+ depth_5(i__),e_depth_3(i__),N_5(i__),X_5(i__)
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(i__),Tmb_6(i__),e_Tmb_6(i__),I_6(i__),e_I_6(i__),
+ FWHM_6(i__),e_FWHM_6(i__),Vlsr_6(i__),e_Vlsr_6(i__),
+ depth_6(i__),e_depth_4(i__),N_6(i__),X_6(i__)
if(ar__8(8:11) .EQ. '') Tmb_6(i__) = rNULL__
if(ar__8(13:16) .EQ. '') e_Tmb_6(i__) = rNULL__
if(ar__8(18:21) .EQ. '') I_6(i__) = rNULL__
if(ar__8(23:26) .EQ. '') e_I_6(i__) = rNULL__
if(ar__8(28:31) .EQ. '') FWHM_6(i__) = rNULL__
if(ar__8(33:36) .EQ. '') e_FWHM_6(i__) = rNULL__
if(ar__8(38:42) .EQ. '') Vlsr_6(i__) = rNULL__
if(ar__8(44:47) .EQ. '') e_Vlsr_6(i__) = rNULL__
if(ar__8(49:53) .EQ. '') depth_6(i__) = rNULL__
if(ar__8(55:59) .EQ. '') e_depth_4(i__) = rNULL__
if(ar__8(61:64) .EQ. '') N_6(i__) = rNULL__
if(ar__8(66:70) .EQ. '') X_6(i__) = rNULL__
c ..............Just test output...........
write(6,9)
+ No_7(i__),Tmb_6(i__),e_Tmb_6(i__),I_6(i__),e_I_6(i__),
+ FWHM_6(i__),e_FWHM_6(i__),Vlsr_6(i__),e_Vlsr_6(i__),
+ depth_6(i__),e_depth_4(i__),N_6(i__),X_6(i__)
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(i__),Tmb_7(i__),e_Tmb_7(i__),I_7(i__),e_I_7(i__),
+ FWHM_7(i__),e_FWHM_7(i__),Vlsr_7(i__),e_Vlsr_7(i__),N_7(i__),
+ X_7(i__)
if(ar__9(8:11) .EQ. '') Tmb_7(i__) = rNULL__
if(ar__9(13:16) .EQ. '') e_Tmb_7(i__) = rNULL__
if(ar__9(18:21) .EQ. '') I_7(i__) = rNULL__
if(ar__9(23:26) .EQ. '') e_I_7(i__) = rNULL__
if(ar__9(28:31) .EQ. '') FWHM_7(i__) = rNULL__
if(ar__9(33:36) .EQ. '') e_FWHM_7(i__) = rNULL__
if(ar__9(38:42) .EQ. '') Vlsr_7(i__) = rNULL__
if(ar__9(44:47) .EQ. '') e_Vlsr_7(i__) = rNULL__
if(ar__9(49:52) .EQ. '') N_7(i__) = rNULL__
if(ar__9(54:58) .EQ. '') X_7(i__) = rNULL__
c ..............Just test output...........
write(6,10)
+ No_8(i__),Tmb_7(i__),e_Tmb_7(i__),I_7(i__),e_I_7(i__),
+ FWHM_7(i__),e_FWHM_7(i__),Vlsr_7(i__),e_Vlsr_7(i__),N_7(i__),
+ X_7(i__)
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(i__),Tmb_8(i__),e_Tmb_8(i__),I_8(i__),e_I_8(i__),
+ FWHM_8(i__),e_FWHM_8(i__),Vlsr_8(i__),e_Vlsr_8(i__),N_8(i__),
+ X_8(i__)
c ..............Just test output...........
write(6,11)
+ No_9(i__),Tmb_8(i__),e_Tmb_8(i__),I_8(i__),e_I_8(i__),
+ FWHM_8(i__),e_FWHM_8(i__),Vlsr_8(i__),e_Vlsr_8(i__),N_8(i__),
+ X_8(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end