Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/128/2642 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-20
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/AJ/128/2642 Aquarius superclusters. II. (Caretta+, 2004)
*================================================================================
*The Aquarius superclusters. II. Spectroscopic and photometric data.
* Caretta C.A., Maia M.A.G., Willmer C.N.A.
* <Astron. J., 128, 2642-2651 (2004)>
* =2004AJ....128.2642C
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Observing runs
integer*4 nr__
parameter (nr__=25) ! Number of records
character*70 ar__ ! Full-size record
integer*4 Run (nr__) ! Run number of the observation
character*4 Setup (nr__) ! Instrumental setup (1)
integer*4 ObsDate_Y (nr__) ! (yr) Date of observation (Year)
character*3 ObsDate_M (nr__) ! ("month") Date of observation (Month)
integer*4 Grating1 (nr__) ! (mm-1) Diffraction grating for 1k pixel CCDs, in
* units of lines/mm
integer*4 Grating2 (nr__) ! (mm-1) ? Diffraction grating for 2k detectors, in
* units of lines/mm
integer*4 Slit1 (nr__) ! (um) Slit width
integer*4 Slit2 (nr__) ! (um) ? Second value of slit width
character*3 Detector (nr__) ! Detector Identification (2)
real*4 Disp (nr__) ! (0.1nm/pix) Spectral Dispersion
integer*4 Res (nr__) ! (0.1nm) Spectral Resolution
integer*4 Lambda1 (nr__) ! (0.1nm) Wavelength range (lower value)
integer*4 Lambda2 (nr__) ! (0.1nm) Wavelength range (upper value)
integer*4 Ntarg (nr__) ! Number of pointed targets (3)
integer*4 Nobj (nr__) ! Number of observed objects (3)
*Note (1): Instrumental setup defined as follows:
* C215 = CASLEO (Complejo Astronomico El Leoncito) 2.15m with REOSC;
* O160 = OPD (Observatorio do Pico dos Dias) 1.60m with B&C (Boller & Chivens
* spectrograph);
* E152 = ESO 1.52m with B&C (Boller & Chivens spectrograph).
*Note (2): Characteristics of the detectors and instrumental setup used:
* -----------------------------------------------------------------------------
* Tel. CCD Origin Dimensions Pixel Gain Noise Scale Efficiency(%)
* (pixels) (um) (e-/ADU) (e- rms) (arcsec/um) (4000-8000)
* -----------------------------------------------------------------------------
* C215 Tek 1024x1024 24 2.0 7.4 3.96/100 95-70
* O160 48 EEV 770x1152 22.5 3.3 6.6 1.04/100 15-45
* O160 101 SITe 1024x1024 24 5.0 5.5 1.04/100 40-65
* O160 106 SITe 1024x1024 24 5.0 4.1 1.25/100 55-75
* E152 39 Loral 2048x2048 15 1.2 5.0 0.82/100 75-95
* E152 38 Loral 2688x512 15 1.6 7.1 0.82/100 75-90
* -----------------------------------------------------------------------------
*Note (3): Number of pointed "targets" and final number of observed "objects"
* with successful redshift measurement (because of multiobjects on the slit)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Aquarius galaxy catalog
integer*4 nr__1
parameter (nr__1=920) ! Number of records
character*73 ar__1 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg (nr__1) ! (deg) Right Ascension J2000
real*8 DEdeg (nr__1) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 v_CMW2004_ (nr__1) ! Galaxy sequential number
integer*4 RAh (nr__1) ! (h) Hour of right ascension (J2000)
integer*4 RAm (nr__1) ! (min) Minute of right ascension (J2000)
real*4 RAs (nr__1) ! (s) Second of right ascension (J2000)
character*1 DE_ (nr__1) ! Sign of the declination (J2000)
integer*4 DEd (nr__1) ! (deg) Degree of declination (J2000)
integer*4 DEm (nr__1) ! (arcmin) Arcminute of declination (J2000)
real*4 DEs (nr__1) ! (arcsec) Arcsecond of declination (J2000)
real*4 bJmag (nr__1) ! (mag) ? The b_J_ band magnitude
real*4 Rcmag (nr__1) ! (mag) ? The R-band magnitude (Cousins)
character*7 Cluster (nr__1) ! Cluster identification
character*1 u_Cluster (nr__1) ! [?] Uncertainty indicator
integer*4 cz (nr__1) ! (km/s) Heliocentric velocity
integer*4 e_cz (nr__1) ! (km/s) Error in cz
character*3 Code (nr__1) ! Code identifications (1)
real*4 R (nr__1) ! Signal-to-noise ratio of the correlation
* peak (R parameter)
integer*4 Run_1 (nr__1) ! Observation run number (in table1.dat)
*Note (1): Code identifications, defined as follows:
* x = Cross-correlation;
* e = Emission lines;
* c = Combined cross-correlation and emission-line redshift;
* A = Quality flag A ("Excellent" - absorption and/or emission lines
* are well defined, and the R parameter is much greater than 3);
* B = Quality flag B ("Good" - absorption and/or emission lines
* confirmed by visual inspection or similar radial velocities
* for most of the templates, but small R (between 2 and 4));
* C = Quality flag C ("Marginal" - value of R greater than 2, but with
* features either in absorption and/or emission too weak to be firmly
* confirmed).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! Cluster data
integer*4 nr__2
parameter (nr__2=74) ! Number of records
character*105 ar__2 ! Full-size record
C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_1 (nr__2) ! (deg) Right Ascension J2000.0
real*8 DEdeg_1 (nr__2) ! (deg) Declination J2000.0
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 AqrCC (nr__2) ! AqrCC number
character*10 Cluster_1 (nr__2) ! Other identification
integer*4 RAh_1 (nr__2) ! (h) Hour of Right Ascension (J2000.0)
integer*4 RAm_1 (nr__2) ! (min) Minute of Right Ascension (J2000.0)
real*4 RAs_1 (nr__2) ! (s) Second of Right Ascension (J2000.0)
character*1 DE__1 (nr__2) ! Sign of the Declination (J2000.0)
integer*4 DEd_1 (nr__2) ! (deg) Degree of Declination (J2000.0)
integer*4 DEm_1 (nr__2) ! (arcmin) Arcminute of Declination (J2000.0)
integer*4 DEs_1 (nr__2) ! (arcsec) Arcsecond of Declination (J2000.0)
character*2 Center (nr__2) ! Type of projected center (1)
integer*4 Ngal (nr__2) ! ? Number of galaxies identified as members
* from the Aquarius Galaxy Catalog
integer*4 Ngalz (nr__2) ! ? Number of these galaxies with redshifts from
* the literature
character*15 r_Ngalz (nr__2) ! References for literature redshifts
integer*4 Nmem (nr__2) ! Total number of member galaxies whose
* velocities were used to calculate the
* systemic velocity
integer*4 Vsys (nr__2) ! (km/s) Systemic velocity
integer*4 E_Vsys (nr__2) ! (km/s) Upper error limit on V
integer*4 e_Vsys_1 (nr__2) ! (km/s) Lower error limit on V
integer*4 sigma (nr__2) ! (km/s) Velocity dispersion
integer*4 E_sigma (nr__2) ! (km/s) Upper error limit on sigma
integer*4 e_sigma_1 (nr__2) ! (km/s) Lower error limit on sigma
*Note (1): The type of projected center is defined as follows:
* B# = position of the brightest cluster galaxy (where # is the number of
* giant elliptical galaxies considered as BCGs of the cluster),
* B? = the BCG is not clearly defined;
* G = geometric mean of the positions of galaxies confirmed as cluster
* members.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'refs.dat' ! References
integer*4 nr__3
parameter (nr__3=15) ! Number of records
character*73 ar__3 ! Full-size record
integer*4 Ref (nr__3) ! Reference number
character*19 Bibcode (nr__3) ! Bibcode
character*21 Aut (nr__3) ! Author's name
character*23 Com (nr__3) ! Comments
C=============================================================================
C Loading file 'table1.dat' ! Observing runs
C Format for file interpretation
1 format(
+ I2,2X,A4,2X,I4,1X,A3,2X,I3,1X,I3,2X,I3,1X,I3,2X,A3,2X,F3.1,2X,
+ I2,2X,I4,1X,I4,2X,I2,2X,I3)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,25
read(1,'(A70)')ar__
read(ar__,1)
+ Run(i__),Setup(i__),ObsDate_Y(i__),ObsDate_M(i__),
+ Grating1(i__),Grating2(i__),Slit1(i__),Slit2(i__),
+ Detector(i__),Disp(i__),Res(i__),Lambda1(i__),Lambda2(i__),
+ Ntarg(i__),Nobj(i__)
if(ar__(25:27) .EQ. '') Grating2(i__) = iNULL__
if(ar__(34:36) .EQ. '') Slit2(i__) = iNULL__
c ..............Just test output...........
write(6,1)
+ Run(i__),Setup(i__),ObsDate_Y(i__),ObsDate_M(i__),
+ Grating1(i__),Grating2(i__),Slit1(i__),Slit2(i__),
+ Detector(i__),Disp(i__),Res(i__),Lambda1(i__),Lambda2(i__),
+ Ntarg(i__),Nobj(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Aquarius galaxy catalog
C Format for file interpretation
2 format(
+ I3,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F5.2,1X,F5.2,
+ 1X,A7,A1,1X,I6,1X,I3,1X,A3,1X,F5.2,1X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,920
read(1,'(A73)')ar__1
read(ar__1,2)
+ v_CMW2004_(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),bJmag(i__),Rcmag(i__),Cluster(i__),
+ u_Cluster(i__),cz(i__),e_cz(i__),Code(i__),R(i__),Run_1(i__)
if(ar__1(29:33) .EQ. '') bJmag(i__) = rNULL__
if(ar__1(35:39) .EQ. '') Rcmag(i__) = rNULL__
RAdeg(i__) = rNULL__
DEdeg(i__) = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c ..............Just test output...........
write(6,2)
+ v_CMW2004_(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),bJmag(i__),Rcmag(i__),Cluster(i__),
+ u_Cluster(i__),cz(i__),e_cz(i__),Code(i__),R(i__),Run_1(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! Cluster data
C Format for file interpretation
3 format(
+ I3,2X,A10,2X,I2,1X,I2,1X,F4.1,2X,A1,I2,1X,I2,1X,I2,2X,A2,2X,
+ I2,2X,I2,2X,A15,2X,I2,2X,I5,1X,I4,1X,I4,2X,I4,1X,I4,1X,I4)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,74
read(1,'(A105)')ar__2
read(ar__2,3)
+ AqrCC(i__),Cluster_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),
+ DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),Center(i__),
+ Ngal(i__),Ngalz(i__),r_Ngalz(i__),Nmem(i__),Vsys(i__),
+ E_Vsys(i__),e_Vsys_1(i__),sigma(i__),E_sigma(i__),
+ e_sigma_1(i__)
if(ar__2(45:46) .EQ. '') Ngal(i__) = iNULL__
if(ar__2(49:50) .EQ. '') Ngalz(i__) = iNULL__
RAdeg_1(i__) = rNULL__
DEdeg_1(i__) = rNULL__
c Derive coordinates RAdeg_1 and DEdeg_1 from input data
c (RAdeg_1 and DEdeg_1 are set to rNULL__ when unknown)
if(RAh_1(i__) .GT. -180) RAdeg_1(i__)=RAh_1(i__)*15.
if(RAm_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAm_1(i__)/4.
if(RAs_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAs_1(i__)/240.
if(DEd_1(i__) .GE. 0) DEdeg_1(i__)=DEd_1(i__)
if(DEm_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEm_1(i__)/60.
if(DEs_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEs_1(i__)/3600.
if(DE__1(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__)
c ..............Just test output...........
write(6,3)
+ AqrCC(i__),Cluster_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),
+ DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),Center(i__),
+ Ngal(i__),Ngalz(i__),r_Ngalz(i__),Nmem(i__),Vsys(i__),
+ E_Vsys(i__),e_Vsys_1(i__),sigma(i__),E_sigma(i__),
+ e_sigma_1(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg_1(i__),DEdeg_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'refs.dat' ! References
C Format for file interpretation
4 format(I2,2X,A19,2X,A21,2X,A23)
C Effective file loading
open(unit=1,status='old',file=
+'refs.dat')
write(6,*) '....Loading file: refs.dat'
do i__=1,15
read(1,'(A73)')ar__3
read(ar__3,4)Ref(i__),Bibcode(i__),Aut(i__),Com(i__)
c ..............Just test output...........
write(6,4)Ref(i__),Bibcode(i__),Aut(i__),Com(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end