Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/232/3 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/ApJS/232/3 Spectral line surveys of 30 regions (Widicus Weaver+, 2017)
*================================================================================
*Deep, broadband spectral line surveys of molecule-rich interstellar clouds.
* Widicus Weaver S.L., Laas J.C., Zou L., Kroll J.A., Rad M.L., Hays B.M.,
* Sanders J.L., Lis D.C., Cross T.N., Wehres N., McGuire B.A., Sumner M.C.
* <Astrophys. J. Suppl. Ser., 232, 3-3 (2017)>
* =2017ApJS..232....3W (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! The 30 star-forming region positions, and
velocities used in the observations, and
source physical parameters
integer*4 nr__
parameter (nr__=30) ! Number of records
character*102 ar__ ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg (nr__) ! (deg) Right Ascension J2000
real*8 DEdeg (nr__) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*16 Name (nr__) ! Source identifier (as in table 3)
character*1 f_Name (nr__) ! [a] a: this source has not been well-characterized
character*17 Type (nr__) ! Source type
integer*4 RAh (nr__) ! (h) Hour of right ascension (J2000)
integer*4 RAm (nr__) ! (min) Minute of right ascension (J2000)
real*4 RAs (nr__) ! (s) Second of right ascension (J2000)
character*1 DE_ (nr__) ! Sign of declination (J2000)
integer*4 DEd (nr__) ! (deg) Degree of declination (J2000)
integer*4 DEm (nr__) ! (arcmin) Arcminute of declination (J2000)
real*4 DEs (nr__) ! (arcsec) Arcsecond of declination (J2000)
real*4 Vlsr (nr__) ! (km/s) [-57/114] LSR velocity
integer*4 Dist (nr__) ! (pc) [230/13600]? Distance
character*1 l_Lum (nr__) ! [~><] Limit or uncertainty flag on Lum
character*6 Lum (nr__) ! (Lsun) Luminosity
character*1 l_Mass (nr__) ! [~><] Limit or uncertainty flag on Mass
real*8 Mass (nr__) ! (Msun) [0.1/14000]? Mass (gas and dust)
character*10 Ref (nr__) ! Reference(s) (see refs.dat file)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3a.dat' ! Detected molecules in 30 sources, their
LTE parameters, and
associated 1{sigma} uncertainties
integer*4 nr__1
parameter (nr__1=991) ! Number of records
character*104 ar__1 ! Full-size record
character*16 Name_1 (nr__1) ! Source identifier
character*9 Mol (nr__1) ! Molecule identifier
real*4 logNT (nr__1) ! ([cm-2]) [11/18.3]? log of the column density
real*4 e_logNT (nr__1) ! ([cm-2]) [0/510]? The 1{sigma} uncertainty in logNT
real*4 NT (nr__1) ! (cm-2) ? Column density (1)
real*4 e_NT (nr__1) ! (cm-2) ? The 1{sigma} uncertainty in NT (1)
real*8 Trot (nr__1) ! (K) [2.7/1141]? Rotational temperature (1)
real*8 e_Trot (nr__1) ! (K) [0.04/183360]? Trot 1{sigma} uncertainty (1)
real*4 FWHM (nr__1) ! (km/s) [0.5/40] Full-Width at Half-Maximum of line
real*4 e_FWHM (nr__1) ! (km/s) [0.01/461]? The 1{sigma} uncertainty in FWHM (1)
real*4 delV (nr__1) ! (km/s) [-11.4/15] Velocity shift (2)
real*4 e_delV (nr__1) ! (km/s) [0.01/188]? The 1{sigma} uncertainty in delV (1)
*Note (1): Only molecules with more than three unblended lines at different
* upper state energies can be fit to a reliable NT and Trot.
* In table3a.dat, other small molecules are fit with arbitrary NT and
* Trot, and their associated uncertainties may be non-physically large.
* In table3b.dat, other detected small molecules are reported with
* non-LTE parameters as blanks.
*Note (2): With respect to the Local Standard of Rest velocity.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3b.dat' ! Detected molecules in 30 sources, their
LTE parameters, and
associated 1{sigma} uncertainties with
non-LTE parameters as blanks
integer*4 nr__2
parameter (nr__2=991) ! Number of records
character*104 ar__2 ! Full-size record
character*16 Name_2 (nr__2) ! Source identifier
character*9 Mol_1 (nr__2) ! Molecule identifier
real*4 logNT_1 (nr__2) ! ([cm-2]) [11/18.3]? log of the column density
real*4 e_logNT_1 (nr__2) ! ([cm-2]) [0/510]? The 1{sigma} uncertainty in logNT
real*4 NT_1 (nr__2) ! (cm-2) ? Column density (1)
real*4 e_NT_1 (nr__2) ! (cm-2) ? The 1{sigma} uncertainty in NT (1)
real*8 Trot_1 (nr__2) ! (K) [2.7/1141]? Rotational temperature (1)
real*8 e_Trot_1 (nr__2) ! (K) [0.04/183360]? Trot 1{sigma} uncertainty (1)
real*4 FWHM_1 (nr__2) ! (km/s) [0.5/40] Full-Width at Half-Maximum of line
real*4 e_FWHM_1 (nr__2) ! (km/s) [0.01/461]? The 1{sigma} uncertainty in FWHM (1)
real*4 delV_1 (nr__2) ! (km/s) [-11.4/15] Velocity shift (2)
real*4 e_delV_1 (nr__2) ! (km/s) [0.01/188]? The 1{sigma} uncertainty in delV (1)
*Note (1): Only molecules with more than three unblended lines at different
* upper state energies can be fit to a reliable NT and Trot.
* In table3a.dat, other small molecules are fit with arbitrary NT and
* Trot, and their associated uncertainties may be non-physically large.
* In table3b.dat, other detected small molecules are reported with
* non-LTE parameters as blanks.
*Note (2): With respect to the Local Standard of Rest velocity.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'figure1.dat' ! Baseline-subtracted, deconvolved spectra for
the 30 star-forming regions
integer*4 nr__3
parameter (nr__3=1714794) ! Number of records
character*42 ar__3 ! Full-size record
character*16 Name_3 (nr__3) ! Source identifier (as in table 3)
real*4 Freq (nr__3) ! (MHz) [212808/268574] Frequency
real*4 Tb (nr__3) ! (K) [-9.6/100.7] Brightness temperature
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'refs.dat' ! References
integer*4 nr__4
parameter (nr__4=43) ! Number of records
character*70 ar__4 ! Full-size record
integer*4 Ref_1 (nr__4) ! Reference code
character*19 BibCode (nr__4) ! Bibcode of the reference
character*26 Auth (nr__4) ! First author's name(s)
character*20 Comm (nr__4) ! Comment
C=============================================================================
C Loading file 'table1.dat' ! The 30 star-forming region positions, and
* velocities used in the observations, and
* source physical parameters
C Format for file interpretation
1 format(
+ A16,1X,A1,1X,A17,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,
+ 1X,F6.2,1X,I5,1X,A1,A6,1X,A1,F8.2,1X,A10)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,30
read(1,'(A102)')ar__
read(ar__,1)
+ Name(i__),f_Name(i__),Type(i__),RAh(i__),RAm(i__),RAs(i__),
+ DE_(i__),DEd(i__),DEm(i__),DEs(i__),Vlsr(i__),Dist(i__),
+ l_Lum(i__),Lum(i__),l_Mass(i__),Mass(i__),Ref(i__)
if(ar__(69:73) .EQ. '') Dist(i__) = iNULL__
if(ar__(84:91) .EQ. '') Mass(i__) = rNULL__
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)
+ Name(i__),f_Name(i__),Type(i__),RAh(i__),RAm(i__),RAs(i__),
+ DE_(i__),DEd(i__),DEm(i__),DEs(i__),Vlsr(i__),Dist(i__),
+ l_Lum(i__),Lum(i__),l_Mass(i__),Mass(i__),Ref(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 'table3a.dat' ! Detected molecules in 30 sources, their
* LTE parameters, and
* associated 1{sigma} uncertainties
C Format for file interpretation
2 format(
+ A16,1X,A9,1X,F5.2,1X,F6.2,1X,E9.3,1X,E9.3,1X,F7.2,1X,F9.2,1X,
+ F5.2,1X,F6.2,1X,F6.2,1X,F6.2)
C Effective file loading
open(unit=1,status='old',file=
+'table3a.dat')
write(6,*) '....Loading file: table3a.dat'
do i__=1,991
read(1,'(A104)')ar__1
read(ar__1,2)
+ Name_1(i__),Mol(i__),logNT(i__),e_logNT(i__),NT(i__),
+ e_NT(i__),Trot(i__),e_Trot(i__),FWHM(i__),e_FWHM(i__),
+ delV(i__),e_delV(i__)
if(ar__1(28:32) .EQ. '') logNT(i__) = rNULL__
if(ar__1(34:39) .EQ. '') e_logNT(i__) = rNULL__
if(ar__1(41:49) .EQ. '') NT(i__) = rNULL__
if(ar__1(51:59) .EQ. '') e_NT(i__) = rNULL__
if(ar__1(61:67) .EQ. '') Trot(i__) = rNULL__
if(ar__1(69:77) .EQ. '') e_Trot(i__) = rNULL__
if(ar__1(85:90) .EQ. '') e_FWHM(i__) = rNULL__
if(ar__1(99:104) .EQ. '') e_delV(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1(i__),Mol(i__),logNT(i__),e_logNT(i__),NT(i__),
+ e_NT(i__),Trot(i__),e_Trot(i__),FWHM(i__),e_FWHM(i__),
+ delV(i__),e_delV(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3b.dat' ! Detected molecules in 30 sources, their
* LTE parameters, and
* associated 1{sigma} uncertainties with
* non-LTE parameters as blanks
C Format for file interpretation
3 format(
+ A16,1X,A9,1X,F5.2,1X,F6.2,1X,E9.3,1X,E9.3,1X,F7.2,1X,F9.2,1X,
+ F5.2,1X,F6.2,1X,F6.2,1X,F6.2)
C Effective file loading
open(unit=1,status='old',file=
+'table3b.dat')
write(6,*) '....Loading file: table3b.dat'
do i__=1,991
read(1,'(A104)')ar__2
read(ar__2,3)
+ Name_2(i__),Mol_1(i__),logNT_1(i__),e_logNT_1(i__),NT_1(i__),
+ e_NT_1(i__),Trot_1(i__),e_Trot_1(i__),FWHM_1(i__),
+ e_FWHM_1(i__),delV_1(i__),e_delV_1(i__)
if(ar__2(28:32) .EQ. '') logNT_1(i__) = rNULL__
if(ar__2(34:39) .EQ. '') e_logNT_1(i__) = rNULL__
if(ar__2(41:49) .EQ. '') NT_1(i__) = rNULL__
if(ar__2(51:59) .EQ. '') e_NT_1(i__) = rNULL__
if(ar__2(61:67) .EQ. '') Trot_1(i__) = rNULL__
if(ar__2(69:77) .EQ. '') e_Trot_1(i__) = rNULL__
if(ar__2(85:90) .EQ. '') e_FWHM_1(i__) = rNULL__
if(ar__2(99:104) .EQ. '') e_delV_1(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ Name_2(i__),Mol_1(i__),logNT_1(i__),e_logNT_1(i__),NT_1(i__),
+ e_NT_1(i__),Trot_1(i__),e_Trot_1(i__),FWHM_1(i__),
+ e_FWHM_1(i__),delV_1(i__),e_delV_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'figure1.dat' ! Baseline-subtracted, deconvolved spectra for
* the 30 star-forming regions
C Format for file interpretation
4 format(A16,2X,E12.6,1X,E11.4)
C Effective file loading
open(unit=1,status='old',file=
+'figure1.dat')
write(6,*) '....Loading file: figure1.dat'
do i__=1,1714794
read(1,'(A42)')ar__3
read(ar__3,4)Name_3(i__),Freq(i__),Tb(i__)
c ..............Just test output...........
write(6,4)Name_3(i__),Freq(i__),Tb(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'refs.dat' ! References
C Format for file interpretation
5 format(I2,1X,A19,1X,A26,1X,A20)
C Effective file loading
open(unit=1,status='old',file=
+'refs.dat')
write(6,*) '....Loading file: refs.dat'
do i__=1,43
read(1,'(A70)')ar__4
read(ar__4,5)Ref_1(i__),BibCode(i__),Auth(i__),Comm(i__)
c ..............Just test output...........
write(6,5)Ref_1(i__),BibCode(i__),Auth(i__),Comm(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end