Conversion of standardized ReadMe file for
file /./ftp/cats/VII/42A 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-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__
parameter (rNULL__=--2147483648.) ! NULL real number
parameter (iNULL__=(-2147483647-1)) ! NULL int number
integer idig ! testing NULL number
C=============================================================================
Cat. VII/42A GMP catalogue: galaxies in the Coma cluster (Godwin+ 1983)
*================================================================================
*The Coma Cluster - I. A catalogue of magnitudes, colours, ellipticities
*and position angles for 6724 galaxies in the field of the Coma cluster
* Godwin J.G., Metcalfe N., Peach J.V.
* <Mon. Not. R. Astron. Soc. 202, 113 (1983)>
* =1983MNRAS.202..113G
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'catalog.dat' ! The catalogue of galaxies, published as
Appendix Microfiche MN202/1 and 2
integer*4 nr__
parameter (nr__=6724) ! Number of records
character*104 ar__ ! Full-size record
integer*4 GMP (nr__) ! [1/6724]+ Identification number
character*1 n_GMP (nr__) ! [?] "?" indicates a probable galaxy (1)
integer*4 Xpos (nr__) ! (arcsec) [-4900/4900]+= increasing westwards (2)
integer*4 Ypos (nr__) ! (arcsec) [-4900/4900] increasing northwards (2)
character*1 n_b26_5 (nr__) ! [*] indicates an overlapped image in b-band
real*4 b26_5 (nr__) ! (mag) b-band integrated magnitude (3)
real*4 r_b26_5 (nr__) ! (arcsec) Isophotal radius in b; rms is 0.5arcsec (4)
character*1 n_r24_75 (nr__) ! [*] indicates an overlapped image in r-band
real*4 r24_75 (nr__) ! (mag) ?=99.99 r-band integrated magnitude (5)
real*4 r_r24_75 (nr__) ! (arcsec) ?=999.9 Isophotal radius in r;
* rms is 0.7arcsec. (4)
real*4 v_b_r (nr__) ! (mag) ?=9.99 (b-r) color index (6)
real*4 r_b_r (nr__) ! (arcsec) ?=99.9 Isophotal radius in (b-r) (4)
real*4 eps (nr__) ! ?=9.9 Ellipticity (7)
integer*4 PosAngle (nr__) ! (deg) ?=999 Position angle of the major axis (8)
real*4 alpha_b (nr__) ! (mag/arcsec) ?=9.99 Slope parameter in b-band (9)
real*4 gamma_b (nr__) ! (mag) ?=9.9 Profile indicator (10)
integer*4 N_points (nr__) ! ?=99 Number of points used to form
* alpha(b) and gamma(b)
integer*4 GP (nr__) ! [1/923]? Number in Godwin and Peach (1977)
integer*4 v_RB67_ (nr__) ! [1/271]? Number in Rood and Baum (1967,1968)
character*19 Com (nr__) ! Comments (11)
*Note (1): Metcalfe and Peach subjectively discriminated between
* galaxies and stars by examination of 14270 images brighter than
* b26.5 = 21.0. The final list of objects includes 6205 galaxies,
* 146 double galaxies, and 2 triple galaxies together with 350 probable
* galaxies and 21 probable double galaxies for a total of 6724 objects.
* An extensive discussion of the completeness and purity of the search
* is given in Section 5 of the paper
*Note (2): X and Y positions in arc seconds relative to the approximate
* center of the cluster, increasing westwards and northwards
* respectively. The origin of these coordinates is at
* RA=12h57.3 Dec=+28d14.4'(1950.0).
* The galaxies are ordered by increasing Xpos.
*Note (3): The b-band magnitude integrated within the 26.5b mag/arsec2
* isophote. This limiting isophote, about 3 percent over the background
* sky brightness, was chosen for optimum signal-to-noise ratio.
* The rms deviation in b26.5, based on the repeatability of measurements
* using different scans of the same plates, is 0.13 mag for images
* brighter than b26.5 = 20, rising to 0.2 mag in the range
* 20 <= b26.5 < 21. The external errors are likely to be of order
* 0.15 mag for images brighter than b26.5 = 20.
*Note (4): The isophotal radius is defined by
* r = sqrt(A/pi)
* where A is the area enclosed within the isophote.
*Note (5): The r-band magnitude integrated within the 24.75r mag/arcsec2
* isophote. This limiting isophote, about 3 percent over the background
* sky brightness, was chosen for optimum signal-to-noise ratio.
* The limiting isophotes for the r-band and b-band magnitudes are
* separated by roughly the mean galaxy color, so that they refer to
* roughly equal metric areas in a given galaxy.
* The rms deviation in r24.75 is 0.2 mag for images brighter than
* r24.75 = 18. The external errors are likely to be of order 0.15 mag
* for images brighter than b26.5 = 20.
* A value of 99.99 indicates the r-band magnitude was omitted because of
* the limited range of surface brightness measurement.
*Note (6): The external error in the (b-r) color is likely to be of order
* 0.15 mag for images brighter than b26.5 = 20.
* A value of 9.99 indicates that the color was not determined.
*Note (7): The ellipticity is defined by
* eps = (1-b/a) , where b/a is the axial ratio of the ellipse best
* approximating the 24r mag/arcsec2 isophote, or, in the event of an
* overlap, the lowest undisturbed isophote. The typical accuracy is
* +/-0.1 for images with isophotal diameter r(r24) >= 5 arcsec.
* A value of 9.9 indicates the ellipticity was not determined.
*Note (8): The position angle, theta, of the major axis of the ellipse
* is measured counter-clockwise on the sky from north. This angle
* is measured to within a few degrees for most images.
* A value of 999 indicates that the position angle was not determined.
*Note (9): The slope parameter is defined as the slope of the logarithmic
* brightness residuals of each image against a polynomial representing
* the standard stellar profile. Its rms dispersion, based on
* measurements from overlapping scan areas, is about 0.05mag/arcsec.
* A value of 9.99 indicates the slope parameter was not determined.
*Note (10): this parameter specifies the amount by which the central
* surface brightness of the b-band image fell short of that which would
* have been expected for a stellar image with the same number of points
* measured.
* A value of 9.9 indicates that this quantity was not determined.
*Note (11): in cases where the object is marked as "double" or "triple"
* the photometric data refer to the combined image of its components.
C=============================================================================
C Loading file 'catalog.dat' ! The catalogue of galaxies, published as
* Appendix Microfiche MN202/1 and 2
C Format for file interpretation
1 format(
+ I4,A1,I6,I6,1X,A1,F5.2,1X,F5.1,1X,A1,F5.2,1X,F5.1,2X,F4.2,1X,
+ F4.1,1X,F3.1,1X,I3,1X,F5.2,1X,F4.1,1X,I2,1X,I3,1X,I3,1X,A19)
C Effective file loading
open(unit=1,status='old',file=
+'catalog.dat')
write(6,*) '....Loading file: catalog.dat'
do i__=1,6724
read(1,'(A104)')ar__
read(ar__,1)
+ GMP(i__),n_GMP(i__),Xpos(i__),Ypos(i__),n_b26_5(i__),
+ b26_5(i__),r_b26_5(i__),n_r24_75(i__),r24_75(i__),
+ r_r24_75(i__),v_b_r(i__),r_b_r(i__),eps(i__),PosAngle(i__),
+ alpha_b(i__),gamma_b(i__),N_points(i__),GP(i__),v_RB67_(i__),
+ Com(i__)
if(ar__(78:80) .EQ. '') GP(i__) = iNULL__
if(ar__(82:84) .EQ. '') v_RB67_(i__) = iNULL__
c ..............Just test output...........
write(6,1)
+ GMP(i__),n_GMP(i__),Xpos(i__),Ypos(i__),n_b26_5(i__),
+ b26_5(i__),r_b26_5(i__),n_r24_75(i__),r24_75(i__),
+ r_r24_75(i__),v_b_r(i__),r_b_r(i__),eps(i__),PosAngle(i__),
+ alpha_b(i__),gamma_b(i__),N_points(i__),GP(i__),v_RB67_(i__),
+ Com(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end