Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/175/462 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. J/ApJS/175/462 Spectroscopy of Coma early-type galaxies (Corsini+, 2008)
*================================================================================
*Spatially resolved spectroscopy of Coma cluster early-type galaxies.
*IV. Completing the data set.
* Corsini E.M., Wegner G., Saglia R.P., Thomas J., Bender R., Thomas D.
* <Astrophys. J. Suppl. Ser., 175, 462-484 (2008)>
* =2008ApJS..175..462C
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'galaxies.dat' ! List of the sample galaxies
integer*4 nr__
parameter (nr__=15) ! Number of records
character*42 ar__ ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg (nr__) ! (deg) Right Ascension J2000
real*8 DEdeg (nr__) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*8 Name (nr__) ! Galaxy name (GMP NNNN)
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 (J2000)
integer*4 DEd (nr__) ! (deg) Declination (J2000)
integer*4 DEm (nr__) ! (arcmin) Declination (J2000)
integer*4 DEs (nr__) ! (arcsec) Declination (J2000)
character*8 OName (nr__) ! Other name
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Photometric parameters of the sample galaxies
observed with HST/WFPC2
integer*4 nr__1
parameter (nr__1=251) ! Number of records
character*118 ar__1 ! Full-size record
character*8 Name_1 (nr__1) ! Galaxy name (GMP NNNN)
real*4 a (nr__1) ! (arcsec) Semi-major axis
real*4 e_a (nr__1) ! (arcsec) Uncertainty in a
real*4 muR (nr__1) ! (mag/arcsec2) Azimuthally averaged R band surface
* brightness
real*4 e_muR (nr__1) ! (mag/arcsec2) Uncertainty in muR
real*4 ell (nr__1) ! Ellipticity (1-b/a)
real*4 e_ell (nr__1) ! Uncertainty in ell
real*4 PA (nr__1) ! (deg) Position angle
real*4 e_PA (nr__1) ! (deg) Uncertainty in PA
real*4 xpos (nr__1) ! (arcsec) The X offset from center coordinate
real*4 ypos (nr__1) ! (arcsec) The Y offset from center coordinate
real*4 e_pos (nr__1) ! (arcsec) Positional error (1)
real*4 a3_a (nr__1) ! (%) ? Third cosine Fourier coefficient (2)
real*4 b3_a (nr__1) ! (%) ? Third sine Fourier coefficient (2)
real*4 a4_a (nr__1) ! (%) Fourth cosine Fourier coefficient (2)
real*4 b4_a (nr__1) ! (%) ? Fourth sine Fourier coefficient (2)
real*4 a6_a (nr__1) ! (%) ? Sixth cosine Fourier coefficient (2)
real*4 b6_a (nr__1) ! (%) ? Sixth sine Fourier coefficient (2)
real*4 Err (nr__1) ! Fourier coefficient error (3)
*Note (1): Error on the center coordinates derived from the residual rms of
* the ellipse fit to the isophotes: e_pos=rms_fit_/N^0.5^ with N{<=}128
* the number of fitted points of the isophote.
*Note (2): values times 100.
*Note (3): Error of Fourier coefficients defined as
* Err=(({Sigma}_i=10_^N/2^(a_i_^2^+b_i_^2^))/(N/2-10))^0.5^*100/a.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Stellar kinematics of the sample galaxies
integer*4 nr__2
parameter (nr__2=143) ! Number of records
character*76 ar__2 ! Full-size record
character*8 Name_2 (nr__2) ! Galaxy name (GMP NNNN)
real*4 r (nr__2) ! (arcsec) [-90/90] Radial distance
real*4 RV (nr__2) ! (km/s) Line-of-sight radial velocity (1)
real*4 e_RV (nr__2) ! (km/s) Uncertainty in RV
real*4 sigma (nr__2) ! (km/s) Velocity dispersion
real*4 e_sigma (nr__2) ! (km/s) Uncertainty in sigma
real*4 H3 (nr__2) ! Hermite decomposition H_3_ coefficient
real*4 e_H3 (nr__2) ! Uncertainty in H3
real*4 H4 (nr__2) ! Hermite decomposition H_4_ coefficient
real*4 e_H4 (nr__2) ! Uncertainty in H4
integer*4 PA_1 (nr__2) ! (deg) Position angle
real*4 Off (nr__2) ! (arcsec) Offset from galaxy center
integer*4 Run (nr__2) ! [1/5] Observation run; see Table 4 (2)
*Note (1): Relative to galaxy center.
*Note (2): Table 4 :
* --------------------------------------------------
* Telescope HET 2.4m Hiltner
* --------------------------------------------------
* Spectrograph LRS Modular
* CCD Ford Aerospace "Echelle" SITe
* Pixel number 3072x3072 2048x2048
* Pixel size ({mu}/m2) 15 x 15 24 x 24
* Binning 2 x 2 1 x 1
* RON (e^-^) 5.2 7.9
* Gain (e^-^ADU^-1^) 1.8 2.7
* Scale (arcsec/pixel) 0.47 0.606
* Slit width (arcsec) 1.0 1.9
* --------------------------------------------------
* ------------------------------------------------------------------------------
* Parameter Run 1 Run 2 Run 3 Run 4 Run 5
* ------------------------------------------------------------------------------
* Telescope HET 2.4m 2.4m HET 2.4m
* Grating (gr/mm) Grism G2 600 1200 1200 Echelle E2 316 1200
* Dispersion 1.99 0.91 1.01 0.72 1.01
* ({AA}/pixel)
* Wavelength range 4284-7400 4702-6401 4644-6455 4870-5956 4544-6363
* ({AA})
* Instrumental FWHM 4.7 2.5 2.3 2.1 2.3
* ({AA})
* Instrumental {sigma} 114 61 53 51 53
* (km/s)
* ------------------------------------------------------------------------------
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! Line strength indices of the sample galaxies
integer*4 nr__3
parameter (nr__3=143) ! Number of records
character*86 ar__3 ! Full-size record
character*8 Name_3 (nr__3) ! Galaxy name (GMP NNNN)
real*4 r_1 (nr__3) ! (arcsec) [-90/90] Radial distance
real*4 Hb (nr__3) ! (0.1nm) H{beta} line strength index
real*4 e_Hb (nr__3) ! (0.1nm) Uncertainty in Hb
real*4 v_MgFe_ (nr__3) ! (0.1nm) ? Combinded Magnesium-Iron index
real*4 e__MgFe_ (nr__3) ! (0.1nm) ? Uncertainty in [MgFe]
real*4 v_Fe_ (nr__3) ! (0.1nm) Average Iron index
real*4 e__Fe_ (nr__3) ! (0.1nm) Uncertainty in <Fe>
real*4 Mgb (nr__3) ! (0.1nm) ? Mgb line strength index
real*4 e_Mgb (nr__3) ! (0.1nm) ? Uncertainty in Mgb
real*4 Mg2 (nr__3) ! (mag) ? Mg_2_ line strength index
real*4 e_Mg2 (nr__3) ! (mag) ? Uncertainty in Mg2
integer*4 PA_2 (nr__3) ! (deg) Position angle
real*4 Off_1 (nr__3) ! (arcsec) Offset for galaxy center
integer*4 Run_1 (nr__3) ! [1/5] Observation run; see Table 4 (above)
C=============================================================================
C Loading file 'galaxies.dat' ! List of the sample galaxies
C Format for file interpretation
1 format(A8,4X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,2X,A8)
C Effective file loading
open(unit=1,status='old',file=
+'galaxies.dat')
write(6,*) '....Loading file: galaxies.dat'
do i__=1,15
read(1,'(A42)')ar__
read(ar__,1)
+ Name(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),OName(i__)
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,1)
+ Name(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),OName(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 'table1.dat' ! Photometric parameters of the sample galaxies
* observed with HST/WFPC2
C Format for file interpretation
2 format(
+ A8,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.1,1X,
+ F4.1,1X,F6.3,1X,F6.3,1X,F5.3,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,
+ 1X,F5.2,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,251
read(1,'(A118)')ar__1
read(ar__1,2)
+ Name_1(i__),a(i__),e_a(i__),muR(i__),e_muR(i__),ell(i__),
+ e_ell(i__),PA(i__),e_PA(i__),xpos(i__),ypos(i__),e_pos(i__),
+ a3_a(i__),b3_a(i__),a4_a(i__),b4_a(i__),a6_a(i__),b6_a(i__),
+ Err(i__)
if(ar__1(79:83) .EQ. '') a3_a(i__) = rNULL__
if(ar__1(85:89) .EQ. '') b3_a(i__) = rNULL__
if(ar__1(97:101) .EQ. '') b4_a(i__) = rNULL__
if(ar__1(103:107) .EQ. '') a6_a(i__) = rNULL__
if(ar__1(109:113) .EQ. '') b6_a(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1(i__),a(i__),e_a(i__),muR(i__),e_muR(i__),ell(i__),
+ e_ell(i__),PA(i__),e_PA(i__),xpos(i__),ypos(i__),e_pos(i__),
+ a3_a(i__),b3_a(i__),a4_a(i__),b4_a(i__),a6_a(i__),b6_a(i__),
+ Err(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Stellar kinematics of the sample galaxies
C Format for file interpretation
3 format(
+ A8,1X,F6.2,1X,F6.1,1X,F5.1,1X,F5.1,1X,F4.1,1X,F6.3,1X,F5.3,1X,
+ F6.3,1X,F5.3,1X,I3,1X,F4.1,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,143
read(1,'(A76)')ar__2
read(ar__2,3)
+ Name_2(i__),r(i__),RV(i__),e_RV(i__),sigma(i__),e_sigma(i__),
+ H3(i__),e_H3(i__),H4(i__),e_H4(i__),PA_1(i__),Off(i__),
+ Run(i__)
c ..............Just test output...........
write(6,3)
+ Name_2(i__),r(i__),RV(i__),e_RV(i__),sigma(i__),e_sigma(i__),
+ H3(i__),e_H3(i__),H4(i__),e_H4(i__),PA_1(i__),Off(i__),
+ Run(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! Line strength indices of the sample galaxies
C Format for file interpretation
4 format(
+ A8,1X,F6.2,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,
+ F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,I3,1X,F4.1,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,143
read(1,'(A86)')ar__3
read(ar__3,4)
+ Name_3(i__),r_1(i__),Hb(i__),e_Hb(i__),v_MgFe_(i__),
+ e__MgFe_(i__),v_Fe_(i__),e__Fe_(i__),Mgb(i__),e_Mgb(i__),
+ Mg2(i__),e_Mg2(i__),PA_2(i__),Off_1(i__),Run_1(i__)
if(ar__3(29:33) .EQ. '') v_MgFe_(i__) = rNULL__
if(ar__3(35:39) .EQ. '') e__MgFe_(i__) = rNULL__
if(ar__3(53:57) .EQ. '') Mgb(i__) = rNULL__
if(ar__3(59:63) .EQ. '') e_Mgb(i__) = rNULL__
if(ar__3(65:69) .EQ. '') Mg2(i__) = rNULL__
if(ar__3(71:75) .EQ. '') e_Mg2(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ Name_3(i__),r_1(i__),Hb(i__),e_Hb(i__),v_MgFe_(i__),
+ e__MgFe_(i__),v_Fe_(i__),e__Fe_(i__),Mgb(i__),e_Mgb(i__),
+ Mg2(i__),e_Mg2(i__),PA_2(i__),Off_1(i__),Run_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end