Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/374/399 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.5, on 2013-Jun-19
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__
ter (rNULL__=-1.e37) ! NULL real number
parameter (iNULL__=-2147483647) ! NULL int number
C=============================================================================
Cat. J/MNRAS/374/399 Infrared star clusters with |b|<20{deg} (Froebrich+, 2007)
*================================================================================
*A systematic survey for infrared star clusters with |b|<20{deg} using 2MASS.
* Froebrich D., Scholz A., Raftery C.L.
* <Mon. Not. R. Astron. Soc., 374, 399-408 (2007)>
* =2007MNRAS.374..399F
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea1.dat' ! Properties of the probable clusters in our sample
of newly detected star cluster candidates
integer*4 nr__
parameter (nr__=473) ! Number of records
character*118 ar__ ! Full-size record
real*8 RA_ ! (deg) Right Ascension J2000
real*8 Dec ! (deg) Declination J2000
integer*4 v_FSR2007_ ! Sequential number
real*8 GLON ! (deg) Galactic longitude
real*8 GLAT ! (deg) Galactic latitude
integer*4 RAh ! (h) Right ascension (J2000)
integer*4 RAm ! (min) Right ascension (J2000)
integer*4 RAs ! (s) Right ascension (J2000)
character*1 DE_ ! Declination sign (J2000)
integer*4 DEd ! (deg) Declination (J2000)
integer*4 DEm ! (arcmin) Declination (J2000)
integer*4 DEs ! (arcsec) Declination (J2000)
real*4 rcore ! (deg) Core radius as fitted to the H-band radial
* star density profile
real*4 rtid ! (deg) Tidal radius as fitted to the H-band radial
* star density profile
integer*4 Dc ! (arcmin-2) Central star density in stars/arcmin^2^ as
* fitted to H-band radial star density profile
integer*4 Dcont ! ?=-99 Star density contrast in the H band (1)
integer*4 N ! Number of stars in the cluster as obtained
* using eq. (3) and the star density profile in
* the H band
integer*4 Nc ! Corrected number of stars in the cluster
* determined using eq. (4) at an apparent
* magnitude of 15 in the H band
real*4 logRN ! ?=99.0 Measure logRN as determined in
* Sect. 4.4 (2)
integer*4 Q ! [0,6] Quality flag determined from our analysis
* of the radial density profile fit (3)
integer*4 Qall ! Detailed quality flags (3)
character*31 Ident ! Name of possible known cluster with erroneous
* coordinates in SIMBAD
*Note (1): determined as ratio of the central star density and the square
* root of the background star density. Values of 99 indicate objects
* with a fitted background star density of zero.
*Note (2): RN=(P^Op^/P^Gl^)/(681/86), where P^Op^ is the Open cluster
* probability and P^Gl^ the globular cluster probability. If logRN is
* negative, the object is probably a globular cluster, otherwise it is
* an open cluster. Values of 99.0 indicate cluster candidates with
* P^Gl^=0.
*Note (3): the details are in Sect. 4.2 of the paper.
* Q contains the number of negative tests (i.e. digits set to 1) in Qall.
* Qall is a 6-digit binary number identifying which of the quality
* tests are negative:
* #1 = the cluster is detected automatically (0) or manually (1);
* #2 = the scatter of the three fitted core radii is larger than 30% (1)
* or smaller (0);
* #3 = the scatter of the three fitted tidal radii is larger than 66% (1)
* or smaller (0);
* #4 = the contrast of the central to background star density
* [Ic/sqrt(Ibg)] in all filters is larger than five (0) or not (1);
* #5 = the ratio of tidal to core radius in all three filters is
* between 3 and 45 (0) or not (1);
* #6 = the core radius in all three filters is larger than twice the
* error of the coordinates, i.e. 0.01{deg} (0) or not (1).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tableb1.dat' ! Properties of the possible clusters in our sample
of newly detected star cluster candidates
integer*4 nr__1
parameter (nr__1=548) ! Number of records
character*118 ar__1 ! Full-size record
real*8 RA__1 ! (deg) Right Ascension J2000
real*8 Dec_1 ! (deg) Declination J2000
integer*4 v_FSR2007__1 ! Sequential number
real*8 GLON_1 ! (deg) Galactic longitude
real*8 GLAT_1 ! (deg) Galactic latitude
integer*4 RAh_1 ! (h) Right ascension (J2000)
integer*4 RAm_1 ! (min) Right ascension (J2000)
integer*4 RAs_1 ! (s) Right ascension (J2000)
character*1 DE__1 ! Declination sign (J2000)
integer*4 DEd_1 ! (deg) Declination (J2000)
integer*4 DEm_1 ! (arcmin) Declination (J2000)
integer*4 DEs_1 ! (arcsec) Declination (J2000)
real*4 rcore_1 ! (deg) Core radius as fitted to the H-band radial
* star density profile
real*4 rtid_1 ! (deg) Tidal radius as fitted to the H-band radial
* star density profile
integer*4 Dc_1 ! (arcmin-2) Central star density in stars/arcmin^2^ as
* fitted to H-band radial star density profile
integer*4 Dcont_1 ! ?=-99 Star density contrast in the H band (1)
integer*4 N_1 ! Number of stars in the cluster as obtained
* using eq. (3) and the star density profile in
* the H band
integer*4 Nc_1 ! Corrected number of stars in the cluster
* determined using eq. (4) at an apparent
* magnitude of 15 in the H band
real*4 logRN_1 ! ?=99.0 Measure logRN as determined in
* Sect. 4.4 (2)
integer*4 Q_1 ! [0,6] Quality flag determined from our analysis
* of the radial density profile fit (3)
integer*4 Qall_1 ! Detailed quality flags (3)
character*31 Ident_1 ! Name of possible known cluster with erroneous
* coordinates in SIMBAD
*Note (1): determined as ratio of the central star density and the square
* root of the background star density. Values of 99 indicate objects
* with a fitted background star density of zero.
*Note (2): RN=(P^Op^/P^Gl^)/(681/86), where P^Op^ is the Open cluster
* probability and P^Gl^ the globular cluster probability. If logRN is
* negative, the object is probably a globular cluster, otherwise it is
* an open cluster. Values of 99.0 indicate cluster candidates with
* P^Gl^=0.
*Note (3): the details are in Sect. 4.2 of the paper.
* Q contains the number of negative tests (i.e. digits set to 1) in Qall.
* Qall is a 6-digit binary number identifying which of the quality
* tests are negative:
* #1 = the cluster is detected automatically (0) or manually (1);
* #2 = the scatter of the three fitted core radii is larger than 30% (1)
* or smaller (0);
* #3 = the scatter of the three fitted tidal radii is larger than 66% (1)
* or smaller (0);
* #4 = the contrast of the central to background star density
* [Ic/sqrt(Ibg)] in all filters is larger than five (0) or not (1);
* #5 = the ratio of tidal to core radius in all three filters is
* between 3 and 45 (0) or not (1);
* #6 = the core radius in all three filters is larger than twice the
* error of the coordinates, i.e. 0.01{deg} (0) or not (1).
C=============================================================================
C Loading file 'tablea1.dat' ! Properties of the probable clusters in our sample
* of newly detected star cluster candidates
C Format for file interpretation
1 format(
+ I4,1X,F7.3,1X,F7.3,1X,I2,1X,I2,1X,I2,1X,A1,I2,1X,I2,1X,I2,1X,
+ F5.3,1X,F5.3,1X,I2,1X,I3,1X,I5,1X,I5,1X,F6.3,1X,I1,1X,I6,1X,
+ A31)
C Effective file loading
open(unit=1,file='tablea1.dat', status='old')
write(6,*) '....Loading file: tablea1.dat'
do i__=1,473
read(1,'(A118)')ar__
read(ar__,1)
+ v_FSR2007_,GLON,GLAT,RAh,RAm,RAs,DE_,DEd,DEm,DEs,rcore,rtid,
+ Dc,Dcont,N,Nc,logRN,Q,Qall,Ident
c Derive coordinates RA_ and Dec from input data
c (RA_ and Dec are set to rNULL_ when unknown)
RA_ = RAh
if(RAh .GE. 0) RA_=RAh*15.
if(RAm .GE. 0) RA_=RA_+RAm/4.
if(RAs .GE. 0) RA_=RA_+RAs/240.
Dec = DEd
if(DEm .GE. 0) Dec=Dec+DEm/60.
if(DEs .GE. 0) Dec=Dec+DEs/3600.
if(DE_.EQ.'-'.AND.Dec.GE.0) Dec=-Dec
c ..............Just test output...........
write(6,1)
+ v_FSR2007_,GLON,GLAT,RAh,RAm,RAs,DE_,DEd,DEm,DEs,rcore,rtid,
+ Dc,Dcont,N,Nc,logRN,Q,Qall,Ident
write(6,'(6H Pos: 2F8.4)') RA_,Dec
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tableb1.dat' ! Properties of the possible clusters in our sample
* of newly detected star cluster candidates
C Format for file interpretation
2 format(
+ I4,1X,F7.3,1X,F7.3,1X,I2,1X,I2,1X,I2,1X,A1,I2,1X,I2,1X,I2,1X,
+ F5.3,1X,F5.3,1X,I2,1X,I3,1X,I5,1X,I5,1X,F6.3,1X,I1,1X,I6,1X,
+ A31)
C Effective file loading
open(unit=1,file='tableb1.dat', status='old')
write(6,*) '....Loading file: tableb1.dat'
do i__=1,548
read(1,'(A118)')ar__1
read(ar__1,2)
+ v_FSR2007__1,GLON_1,GLAT_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,
+ DEm_1,DEs_1,rcore_1,rtid_1,Dc_1,Dcont_1,N_1,Nc_1,logRN_1,Q_1,
+ Qall_1,Ident_1
c Derive coordinates RA__1 and Dec_1 from input data
c (RA__1 and Dec_1 are set to rNULL_ when unknown)
RA__1 = RAh_1
if(RAh_1 .GE. 0) RA__1=RAh_1*15.
if(RAm_1 .GE. 0) RA__1=RA__1+RAm_1/4.
if(RAs_1 .GE. 0) RA__1=RA__1+RAs_1/240.
Dec_1 = DEd_1
if(DEm_1 .GE. 0) Dec_1=Dec_1+DEm_1/60.
if(DEs_1 .GE. 0) Dec_1=Dec_1+DEs_1/3600.
if(DE__1.EQ.'-'.AND.Dec_1.GE.0) Dec_1=-Dec_1
c ..............Just test output...........
write(6,2)
+ v_FSR2007__1,GLON_1,GLAT_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,
+ DEm_1,DEs_1,rcore_1,rtid_1,Dc_1,Dcont_1,N_1,Nc_1,logRN_1,Q_1,
+ Qall_1,Ident_1
write(6,'(6H Pos: 2F8.4)') RA__1,Dec_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end