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-25
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/A+AS/134/59 The central Fornax Cluster. I. Galaxy photometry (Hilker+ 1999)
*================================================================================
*The central region of the Fornax cluster. I. A catalog and photometric
*properties of galaxies in selected CCD fields
* Hilker M., Kissler-Patig M., Richtler T., Infante L., Quintana H.
* <Astron. Astrophys. Suppl. Ser. 134, 59 (1999)>
* =1999A&AS..134...59H (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'appendxa.dat' ! Photometric catalog
integer*4 nr__
parameter (nr__=873) ! Number of records
character*96 ar__ ! Full-size record
real*8 RA_ (nr__) ! (deg) Right Ascension J2000
real*8 Dec (nr__) ! (deg) Declination J2000
character*9 CGF (nr__) ! Catalog name of the object (1)
character*1 n_CGF (nr__) ! [*] *: objects that are also listed in
* appendb.dat file
integer*4 RAh (nr__) ! (h) Right ascension (J2000)
integer*4 RAm (nr__) ! (min) Right ascension (J2000)
real*4 RAs (nr__) ! (s) Right ascension (J2000)
character*1 DE_ (nr__) ! Declination sign
integer*4 DEd (nr__) ! (deg) Declination (J2000)
integer*4 DEm (nr__) ! (arcmin) Declination (J2000)
real*4 DEs (nr__) ! (arcsec) Declination (J2000)
real*4 Vtot (nr__) ! (mag) Total V apparent magnitude
character*1 nbo (nr__) ! [n] Neighbouring object (2)
real*4 sbVp (nr__) ! (mag/arcsec2) V peak surface brightness
real*4 V_I (nr__) ! (mag) ? (V-I) color within 3 arcsec aperture
real*4 ell (nr__) ! ? Ellipticity, ell = 1 - b/a
real*4 PA (nr__) ! (deg) [-90/90] Position angle of the major axis (3)
real*4 aphot (nr__) ! (arcsec) Size of the major axis
real*4 d26 (nr__) ! (arcsec) ? Size of major axis at sbV = 26 mag/arcsec2
real*4 aeff (nr__) ! (arcsec) ? Effective semi-major axis
real*4 sbVeff (nr__) ! (mag/arcsec2) V effective surface brightness
*Note (1):
* The catalog names are prefixed by the acronym CGF (Catalog of
* Galaxies in Fornax) followed by a sequence number of the field and the
* sequence number of the galaxy in this field (ordered with decreasing
* magnitude). In the catalog the galaxies are ordered by increasing
* right ascension.
*Note (2):
* The letter "n" indicates that a neighbouring object is present within
* two isophotal radii.
*Note (3):
* 0 degree is in east direction, positive angles towards the south,
* and negative angles towards the north direction
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'appendxb.dat' ! Catalog of profile fit parameters
integer*4 nr__1
parameter (nr__1=192) ! Number of records
character*119 ar__1 ! Full-size record
real*8 RA__1 (nr__1) ! (deg) Right Ascension J2000
real*8 Dec_1 (nr__1) ! (deg) Declination J2000
character*9 CGF_1 (nr__1) ! Catalog name of the object (see Appendix A)
integer*4 RAh_1 (nr__1) ! (h) Right ascension (J2000)
integer*4 RAm_1 (nr__1) ! (min) Right ascension (J2000)
real*4 RAs_1 (nr__1) ! (s) Right ascension (J2000)
character*1 DE__1 (nr__1) ! Declination sign
integer*4 DEd_1 (nr__1) ! (deg) Declination (J2000)
integer*4 DEm_1 (nr__1) ! (arcmin) Declination (J2000)
real*4 DEs_1 (nr__1) ! (arcsec) Declination (J2000)
character*8 MType (nr__1) ! Morphological Type
real*4 Vtot_1 (nr__1) ! (mag) Total V apparent magnitude
character*1 nbo_1 (nr__1) ! Neighboring object (see Appendix A)
real*4 rcap (nr__1) ! (arcsec) Apparent core radius (1)
real*4 rccor (nr__1) ! (arcsec) ? Corrected core radius (using
* Kormendy 1985ApJ...295...73K)
real*4 sbVcor (nr__1) ! (mag/arcsec2) ? Corrected central surface brightness (sb)
real*4 nsers (nr__1) ! ? Exponent of Sersic (1968, Atlas de galaxias
* australes) profile fit
real*4 sbVser (nr__1) ! (mag/arcsec2) ? Central surface brightness of Sersic fit
real*4 r0sers (nr__1) ! (arcsec) ? Scale length of Sersic profile fit
real*4 sbVdV (nr__1) ! (mag/arcsec2) ? Central sb of de Vaucouleurs profile fit
real*4 refdV (nr__1) ! (arcsec) ? Effective radius of de Vaucouleurs profile
real*4 sbVexp (nr__1) ! (mag/arcsec2) ? Central sb of exponential profile fit
real*4 alpha (nr__1) ! (arcsec) ? Scale length of exponential profile fit
character*1 flag (nr__1) ! [bdne] Exponential profile fit flag (2)
*Note (1):
* Radius where the apparent central surface intensity has its half
* value
*Note (2): the meaning of the flag is:
* b: bulge is present
* d: light deficiency in the center
* n: nucleus present
* e: exponential law at all radii
C=============================================================================
C Loading file 'appendxa.dat' ! Photometric catalog
C Format for file interpretation
1 format(
+ A9,A1,2X,I1,1X,I2,1X,F5.2,2X,A1,I2,1X,I2,1X,F4.1,2X,F5.2,1X,
+ A1,2X,F5.2,2X,F5.2,2X,F4.2,2X,F5.1,2X,F4.1,2X,F4.1,2X,F4.1,2X,
+ F5.2)
C Effective file loading
open(unit=1,file='appendxa.dat', status='old')
write(6,*) '....Loading file: appendxa.dat'
do i__=1,873
read(1,'(A96)')ar__
read(ar__,1)
+ CGF(i__),n_CGF(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
+ DEd(i__),DEm(i__),DEs(i__),Vtot(i__),nbo(i__),sbVp(i__),
+ V_I(i__),ell(i__),PA(i__),aphot(i__),d26(i__),aeff(i__),
+ sbVeff(i__)
if(ar__(54:58) .EQ. '') V_I(i__) = rNULL__
if(ar__(61:64) .EQ. '') ell(i__) = rNULL__
if(ar__(80:83) .EQ. '') d26(i__) = rNULL__
if(ar__(86:89) .EQ. '') aeff(i__) = rNULL__
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)
+ CGF(i__),n_CGF(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
+ DEd(i__),DEm(i__),DEs(i__),Vtot(i__),nbo(i__),sbVp(i__),
+ V_I(i__),ell(i__),PA(i__),aphot(i__),d26(i__),aeff(i__),
+ sbVeff(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 'appendxb.dat' ! Catalog of profile fit parameters
C Format for file interpretation
2 format(
+ A9,2X,I1,1X,I2,1X,F5.2,2X,A1,I2,1X,I2,1X,F4.1,2X,A8,2X,F5.2,
+ 1X,A1,2X,F3.1,2X,F3.1,2X,F5.2,2X,F4.2,2X,F5.2,2X,F4.2,2X,F5.2,
+ 2X,F5.2,2X,F5.2,1X,F5.2,2X,A1)
C Effective file loading
open(unit=1,file='appendxb.dat', status='old')
write(6,*) '....Loading file: appendxb.dat'
do i__=1,192
read(1,'(A119)')ar__1
read(ar__1,2)
+ CGF_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__),
+ DEd_1(i__),DEm_1(i__),DEs_1(i__),MType(i__),Vtot_1(i__),
+ nbo_1(i__),rcap(i__),rccor(i__),sbVcor(i__),nsers(i__),
+ sbVser(i__),r0sers(i__),sbVdV(i__),refdV(i__),sbVexp(i__),
+ alpha(i__),flag(i__)
if(ar__1(61:63) .EQ. '') rccor(i__) = rNULL__
if(ar__1(66:70) .EQ. '') sbVcor(i__) = rNULL__
if(ar__1(73:76) .EQ. '') nsers(i__) = rNULL__
if(ar__1(79:83) .EQ. '') sbVser(i__) = rNULL__
if(ar__1(86:89) .EQ. '') r0sers(i__) = rNULL__
if(ar__1(92:96) .EQ. '') sbVdV(i__) = rNULL__
if(ar__1(99:103) .EQ. '') refdV(i__) = rNULL__
if(ar__1(106:110) .EQ. '') sbVexp(i__) = rNULL__
if(ar__1(112:116) .EQ. '') alpha(i__) = rNULL__
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(i__) = RAh_1(i__)
if(RAh_1(i__) .GE. 0) RA__1(i__)=RAh_1(i__)*15.
if(RAm_1(i__) .GE. 0) RA__1(i__)=RA__1(i__)+RAm_1(i__)/4.
if(RAs_1(i__) .GE. 0) RA__1(i__)=RA__1(i__)+RAs_1(i__)/240.
Dec_1(i__) = DEd_1(i__)
if(DEm_1(i__) .GE. 0) Dec_1(i__)=Dec_1(i__)+DEm_1(i__)/60.
if(DEs_1(i__) .GE. 0) Dec_1(i__)=Dec_1(i__)+DEs_1(i__)/3600.
if(DE__1(i__).EQ.'-'.AND.Dec_1(i__).GE.0) Dec_1(i__)=-Dec_1(i__)
c ..............Just test output...........
write(6,2)
+ CGF_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__),
+ DEd_1(i__),DEm_1(i__),DEs_1(i__),MType(i__),Vtot_1(i__),
+ nbo_1(i__),rcap(i__),rccor(i__),sbVcor(i__),nsers(i__),
+ sbVser(i__),r0sers(i__),sbVdV(i__),refdV(i__),sbVexp(i__),
+ alpha(i__),flag(i__)
write(6,'(6H Pos: 2F8.4)') RA__1(i__),Dec_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end