Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/422/25 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 2026-Sep-10
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/MNRAS/422/25 SDSS DR7 groups, clusters and filaments (Smith+, 2012)
*================================================================================
*Multiscale probability mapping: groups, clusters and an algorithmic search for
*filaments in SDSS.
* Smith A.G., Hopkins A.M., Hunstead R.W., Pimbblet K.A.
* <Mon. Not. R. Astron. Soc., 422, 25-43 (2012)>
* =2012MNRAS.422...25S
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! *Catalogue of MSPM groups and clusters in SDSS DR7
integer*4 nr__
parameter (nr__=10443) ! Number of records
character*70 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 ID (nr__) ! [1/10443] Identification number of MSPM group
* or cluster (MultiScale Probability Mapping)
real*8 RAdeg (nr__) ! (deg) Right Ascension (J2000) (deg)
real*8 DEdeg (nr__) ! (deg) Declination (J2000) (deg)
real*8 z (nr__) ! [0.025/0.24] Redshift
real*4 P (nr__) ! [0.5/1] Overdensity probability at peak P-S (1)
integer*4 N (nr__) ! [4/157] Count of galaxies with r<17.77
real*4 R (nr__) ! (Mpc) Radius in Mpc/h enclosing region with P-S>0.5 (1)
real*4 LDCa (nr__) ! ? Local density contrast (LDC) with radii
* [0.4,2]Mpc/h
character*1 n_LDCa (nr__) ! [i] 'i' for infinity
real*4 LDCb (nr__) ! ? Local density contrast (LDC) with radii
* [1,2]Mpc/h
character*1 n_LDCb (nr__) ! [i] 'i' for infinity
integer*4 sigma (nr__) ! (km/s) Velocity dispersion {sigma}
real*4 Dens (nr__) ! ?=-1 Galaxy density within 0.4Mpc/h in units
* of the background density
*Note (1):
* * P is the overdensity probability, and P>0.5 designates a location
* denser than average
* * S is the shape statistic, and S<0.5 designates a location denser than
* surrounding locations
* * P-S>0.5 designates therefore locations denser than average and denser than
* surrounding locations, guaranteeing P > 0.5 and S < 0.5
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! *Catalogue of MSPM filaments in SDSS DR7
integer*4 nr__1
parameter (nr__1=53) ! Number of records
character*42 ar__1 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 ID_1 (nr__1) ! [1/10443] Identification number
* (from central group or cluster)
real*8 RAdeg_1 (nr__1) ! (deg) Right Ascension (J2000) (deg)
real*8 DEdeg_1 (nr__1) ! (deg) Declination (J2000) (deg)
real*8 z_1 (nr__1) ! Redshift
real*4 Pe (nr__1) ! Highest elongation probability within 10Mpc/h
integer*4 N_1 (nr__1) ! [3/15] Count of groups and clusters within a
* 10Mpc/h radius on the sky and 10Mpc/h in the
* line of sight (cylindrical aperture)
character*2 Mph (nr__1) ! [I II V] Morphology of the filament (1)
*Note (1): the morphology class is based on the appearance in two
* orthogonal planes: the sky, and a plane perpendicular to the
* sky and aligned with the the filament:
* I = Type I (straight). Both are straight
* II = Type II (warped). At least one is curved, with neither being
* uniform or irregular
* V = Type IV (uniform). Both are uniform
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'void.dat' ! *Catalogue of MSPM voids in SDSS DR7
integer*4 nr__2
parameter (nr__2=47) ! Number of records
character*30 ar__2 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 ID_2 (nr__2) ! Identification number (from a lattice of
* sampling positions formed by spectroscopic
* plate centres and a uniform redshift interval)
real*8 RAdeg_2 (nr__2) ! (deg) Centre Right Ascension (J2000)
real*8 DEdeg_2 (nr__2) ! (deg) Centre Declination (J2000)
real*8 z_2 (nr__2) ! Centre Redshift
C=============================================================================
C Loading file 'table1.dat' ! *Catalogue of MSPM groups and clusters in SDSS DR7
C Format for file interpretation
1 format(
+ I5,1X,F8.4,1X,F8.4,1X,F7.5,1X,F5.3,1X,I3,1X,F5.3,1X,F5.1,A1,
+ F4.1,A1,I4,1X,F6.1)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,10443
read(1,'(A70)')ar__
read(ar__,1)
+ ID(i__),RAdeg(i__),DEdeg(i__),z(i__),P(i__),N(i__),R(i__),
+ LDCa(i__),n_LDCa(i__),LDCb(i__),n_LDCb(i__),sigma(i__),
+ Dens(i__)
if(ar__(49:53) .EQ. '') LDCa(i__) = rNULL__
if(ar__(55:58) .EQ. '') LDCb(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ ID(i__),RAdeg(i__),DEdeg(i__),z(i__),P(i__),N(i__),R(i__),
+ LDCa(i__),n_LDCa(i__),LDCb(i__),n_LDCb(i__),sigma(i__),
+ Dens(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! *Catalogue of MSPM filaments in SDSS DR7
C Format for file interpretation
2 format(I5,1X,F8.4,1X,F7.4,1X,F7.5,1X,F5.3,1X,I2,1X,A2)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,53
read(1,'(A42)')ar__1
read(ar__1,2)
+ ID_1(i__),RAdeg_1(i__),DEdeg_1(i__),z_1(i__),Pe(i__),N_1(i__),
+ Mph(i__)
c ..............Just test output...........
write(6,2)
+ ID_1(i__),RAdeg_1(i__),DEdeg_1(i__),z_1(i__),Pe(i__),N_1(i__),
+ Mph(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'void.dat' ! *Catalogue of MSPM voids in SDSS DR7
C Format for file interpretation
3 format(I5,1X,F8.4,1X,F7.4,1X,F7.5)
C Effective file loading
open(unit=1,status='old',file=
+'void.dat')
write(6,*) '....Loading file: void.dat'
do i__=1,47
read(1,'(A30)')ar__2
read(ar__2,3)ID_2(i__),RAdeg_2(i__),DEdeg_2(i__),z_2(i__)
c ..............Just test output...........
write(6,3)ID_2(i__),RAdeg_2(i__),DEdeg_2(i__),z_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end