Conversion of standardized ReadMe file for
catalog 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.5, on 2013-May-21
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__=-1.e37) ! NULL real number
parameter (iNULL__=-2147483647) ! NULL int number
C=============================================================================
Cat. J/PASJ/53/1081 CO molecular clouds in Lupus (Tachihara+, 2001)
*================================================================================
*^12^CO molecular cloud survey and global star formation in Lupus
* Tachihara K., Toyoda S., Onishi T., Mizuno A., Fukui Y., Neuhaeuser R.
* <Publ. Astron. Soc. Jap., 53, 1081 (2001)>
* =2001PASJ...53.1081T
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Observed properties of the clouds
integer*4 nr__
parameter (nr__=113) ! Number of records
character*78 ar__ ! Full-size record
real*8 RA_ (nr__) ! (deg) Right Ascension 1950
real*8 Dec (nr__) ! (deg) Declination 1950
integer*4 v_TTO2001_ (nr__) ! Sequential number
character*12 Name (nr__) ! Cloud name (lll.ll+bb.bb)
character*1 n_Name (nr__) ! [a-d] Previously identified clouds (1)
real*8 GLON (nr__) ! (deg) Galactic longitude (2)
real*4 GLAT (nr__) ! (deg) Galactic latitude (2)
integer*4 RAh (nr__) ! (h) Right ascension (1950) (2)
integer*4 RAm (nr__) ! (min) Right ascension (1950) (2)
real*4 RAs (nr__) ! (s) Right ascension (1950) (2)
character*1 DE_ (nr__) ! Declination sign (1950) (2)
integer*4 DEd (nr__) ! (deg) Declination (1950) (2)
integer*4 DEm (nr__) ! (arcmin) Declination (1950) (2)
real*4 DEs (nr__) ! (arcsec) Declination (1950) (2)
real*4 TR_ (nr__) ! (K) ^12^CO Absolute peak antenna temperature (3)
real*4 FWHM (nr__) ! (km/s) ^12^CO FWHM line-width (3)
real*4 Vlsr (nr__) ! (km/s) ^12^CO centre velocity with respect to LSR (3)
real*4 W_CO (nr__) ! (K.km/s) ^12^CO integrated intensities at the peak
* positions of the CO clouds (3)
*Note (1): These clouds correspond to previously identified clouds as follows:
* a = Lupus 1
* b = Lupus 2
* c = Lupus 9
* d = This cloud consists of the known clouds of Lupus 3, 4, 5, and 6.
*Note (2): Position of the ^12^CO peak integrated intensity of the cloud
*Note (3): If the ^12^CO emission at the peak intensity position shows
* apparent multi-velocity components, the parameters are given
* for each component. The peak intensity position of the cloud
* is defined by a Gaussian fitting.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Physical properties of the clouds
integer*4 nr__1
parameter (nr__1=105) ! Number of records
character*38 ar__1 ! Full-size record
integer*4 v_TTO2001__1(nr__1) ! Sequential name
character*12 Name_1 (nr__1) ! Cloud name (lll.ll+bb.bb)
real*4 N_H2 (nr__1) ! (10+20cm-2) H_2_ column density
real*4 Mass (nr__1) ! (solMass) Mass
real*4 Rad (nr__1) ! (pc) Radius
real*4 n_H2_1 (nr__1) ! (10+2cm-3) Average number density
C=============================================================================
C Loading file 'table1.dat' ! Observed properties of the clouds
C Format for file interpretation
1 format(
+ I3,1X,A12,A1,1X,F7.3,1X,F6.3,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,
+ I2,1X,F4.1,1X,F5.2,1X,F4.2,1X,F5.2,1X,F5.2)
C Effective file loading
open(unit=1,file='table1.dat', status='old')
write(6,*) '....Loading file: table1.dat'
do i__=1,113
read(1,'(A78)')ar__
read(ar__,1)
+ v_TTO2001_(i__),Name(i__),n_Name(i__),GLON(i__),GLAT(i__),
+ RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__),
+ DEs(i__),TR_(i__),FWHM(i__),Vlsr(i__),W_CO(i__)
c Derive coordinates RA_ and Dec from input data
c (RA_ and Dec are set to rNULL_ when unknown)
RA_(i__) = RAh(i__)
if(RAh(i__) .GE. 0) RA_(i__)=RAh(i__)*15.
if(RAm(i__) .GE. 0) RA_(i__)=RA_(i__)+RAm(i__)/4.
if(RAs(i__) .GE. 0) RA_(i__)=RA_(i__)+RAs(i__)/240.
Dec(i__) = DEd(i__)
if(DEm(i__) .GE. 0) Dec(i__)=Dec(i__)+DEm(i__)/60.
if(DEs(i__) .GE. 0) Dec(i__)=Dec(i__)+DEs(i__)/3600.
if(DE_(i__).EQ.'-'.AND.Dec(i__).GE.0) Dec(i__)=-Dec(i__)
c ..............Just test output...........
write(6,1)
+ v_TTO2001_(i__),Name(i__),n_Name(i__),GLON(i__),GLAT(i__),
+ RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__),
+ DEs(i__),TR_(i__),FWHM(i__),Vlsr(i__),W_CO(i__)
write(6,'(6H Pos: 2F8.4)') RA_(i__),Dec(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Physical properties of the clouds
C Format for file interpretation
2 format(I3,1X,A12,1X,F4.1,1X,F6.1,1X,F4.2,1X,F4.1)
C Effective file loading
open(unit=1,file='table2.dat', status='old')
write(6,*) '....Loading file: table2.dat'
do i__=1,105
read(1,'(A38)')ar__1
read(ar__1,2)
+ v_TTO2001__1(i__),Name_1(i__),N_H2(i__),Mass(i__),Rad(i__),
+ n_H2_1(i__)
c ..............Just test output...........
write(6,2)
+ v_TTO2001__1(i__),Name_1(i__),N_H2(i__),Mass(i__),Rad(i__),
+ n_H2_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end