Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/519/A40 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-16
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/A+A/519/A40 Spectroscopic study of 65 nearby galaxies (Annibali+, 2010)
*================================================================================
*Nearby early-type galaxies with ionized gas.
*IV. Origin and powering mechanism of the ionized gas.
* Annibali F., Bressan A., Rampazzo R., Zeilinger W.W., Vega O., Panuzzo P.
* <Astron. Astrophys. 519, A40 (2010)>
* =2010A&A...519A..40A
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Sample overview
integer*4 nr__
parameter (nr__=65) ! Number of records
character*89 ar__ ! Full-size record
character*7 Galaxy (nr__) ! Galaxy name
character*13 MType1 (nr__) ! RSA (Cat. VII/51) morphological type
character*13 MType2 (nr__) ! RC3 (Cat. VII/155) morphological type
integer*4 sigmac (nr__) ! (km/s) ?=- Central velocity dispersion
* (from HYPERLEDA, Cat. VII/237)
real*4 re (nr__) ! (arcsec) Effective radius from RC3 (Cat. VII/155)
integer*4 HV (nr__) ! (km/s) Galaxy systemic velocity
real*4 rho (nr__) ! (Mpc-3) ?=- Richness parameter rho_xyz_,
* in gal/Mpc^3^ (1)
real*4 Age (nr__) ! (Gyr) ? Age (from paper II, Cat. J/A+A/445/79)
real*4 e_Age (nr__) ! (Gyr) ? rms uncertainty on Age
real*4 Z (nr__) ! ? Metallicity (from paper II,
* Cat. J/A+A/445/79)
real*4 e_Z (nr__) ! ? rms uncertainty on Z
real*4 v_a_Fe_ (nr__) ! ([-]) ? {alpha}/Fe abundance ratio (from paper II,
* Cat. J/A+A/445/79)
real*4 e__a_Fe_ (nr__) ! ([-]) ? rms uncertainty on [a/Fe]
character*1 Data (nr__) ! [*] * when data in tables 2 and 3
*Note (1): It represents the density of galaxies brighter than -16 apparent
* B-mag in the vicinity of the entry.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Emission line fluxes
integer*4 nr__1
parameter (nr__1=212) ! Number of records
character*154 ar__1 ! Full-size record
character*7 Galaxy_1 (nr__1) ! Galaxy identifier
integer*4 Reg (nr__1) ! [1/4] Annulus code (G1)
real*4 E_B_V (nr__1) ! Reddening
real*4 e_E_B_V (nr__1) ! Error on reddening
real*4 FOII (nr__1) ! (10-19W/m2/arcsec2) ?=99.99 OII 3727 flux (G2)
real*4 e_FOII (nr__1) ! (10-19W/m2/arcsec2) ?=9.99 Error on FOII 3727
real*4 FNeIII (nr__1) ! (10-19W/m2/arcsec2) ?=99.99 NeIII 3869 flux (G2)
real*4 e_FNeIII (nr__1) ! (10-19W/m2/arcsec2) ?=9.99 Error on FNeIII 3869
real*4 FHb (nr__1) ! (10-19W/m2/arcsec2) ?=99.99 Hbeta flux (G2)
real*4 e_FHb (nr__1) ! (10-19W/m2/arcsec2) ?=9.99 Error on FHbeta
real*4 FOIII59 (nr__1) ! (10-19W/m2/arcsec2) ?=99.99 OIII 4959 flux (G2)
real*4 e_FOIII59 (nr__1) ! (10-19W/m2/arcsec2) ?=9.99 Error on FOIII 4959
real*4 FOIII07 (nr__1) ! (10-19W/m2/arcsec2) ?=99.99 OIII 5007 flux (G2)
real*4 e_FOIII07 (nr__1) ! (10-19W/m2/arcsec2) ?=9.99 Error on FOIII 5007
real*4 FOI (nr__1) ! (10-19W/m2/arcsec2) ?=99.99 OI 6300 flux (G2)
real*4 e_FOI (nr__1) ! (10-19W/m2/arcsec2) ?=9.99 Error on FOI 6300
real*4 FNII48 (nr__1) ! (10-19W/m2/arcsec2) ?=99.99 NII 6548 flux (G2)
real*4 e_FNII48 (nr__1) ! (10-19W/m2/arcsec2) ?=9.99 Error on FNII 6548
real*4 FHa (nr__1) ! (10-19W/m2/arcsec2) ?=99.99 Halpha flux (G2)
real*4 e_FHa (nr__1) ! (10-19W/m2/arcsec2) ?=9.99 Error on FHalpha
real*4 FNII84 (nr__1) ! (10-19W/m2/arcsec2) ?=99.99 NII 6584 flux (G2)
real*4 e_FNII84 (nr__1) ! (10-19W/m2/arcsec2) ?=9.99 Error on FNII 6584
real*4 FSII17 (nr__1) ! (10-19W/m2/arcsec2) ?=99.99 SII 6717 flux (G2)
real*4 e_FSII17 (nr__1) ! (10-19W/m2/arcsec2) ?=9.99 Error on FSII 6717
real*4 FSII31 (nr__1) ! (10-19W/m2/arcsec2) ?=99.99 SII 6731 flux (G2)
real*4 e_FSII31 (nr__1) ! (10-19W/m2/arcsec2) ?=9.99 Error on FSII 6731
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Emission line equivalent widths
integer*4 nr__2
parameter (nr__2=212) ! Number of records
character*141 ar__2 ! Full-size record
character*7 Galaxy_2 (nr__2) ! Galaxy identifier
integer*4 Reg_1 (nr__2) ! [1/4] Annulus code (G1)
real*4 WOII (nr__2) ! (0.1nm) ?=99.99 OII 3727 equivalent width (G2)
real*4 e_WOII (nr__2) ! (0.1nm) ?=9.99 Error on WOII 3727
real*4 WNeIII (nr__2) ! (0.1nm) ?=99.99 NeIII 3869 equivalent width (G2)
real*4 e_WNeIII (nr__2) ! (0.1nm) ?=9.99 Error on WNeIII 3869
real*4 WHb (nr__2) ! (0.1nm) ?=99.99 Hbeta equivalent width (G2)
real*4 e_WHb (nr__2) ! (0.1nm) ?=9.99 Error on WHbeta
real*4 WOIII59 (nr__2) ! (0.1nm) ?=99.99 OIII 4959 equivalent width (G2)
real*4 e_WOIII59 (nr__2) ! (0.1nm) ?=9.99 Error on WOIII 4959
real*4 WOIII07 (nr__2) ! (0.1nm) ?=99.99 OIII 5007 equivalent width (G2)
real*4 e_WOIII07 (nr__2) ! (0.1nm) ?=9.99 Error on WOIII 5007
real*4 WOI (nr__2) ! (0.1nm) ?=99.99 OI 6300 equivalent width (G2)
real*4 e_WOI (nr__2) ! (0.1nm) ?=9.99 Error on WOI 6300
real*4 WNII48 (nr__2) ! (0.1nm) ?=99.99 NII 6548 equivalent width (G2)
real*4 e_WNII48 (nr__2) ! (0.1nm) ?=9.99 Error on WNII 6548
real*4 WHa (nr__2) ! (0.1nm) ?=99.99 Halpha equivalent width (G2)
real*4 e_WHa (nr__2) ! (0.1nm) ?=9.99 Error on WHalpha
real*4 WNII84 (nr__2) ! (0.1nm) ?=99.99 NII 6584 equivalent width (G2)
real*4 e_WNII84 (nr__2) ! (0.1nm) ?=9.99 Error on WNII 6584
real*4 WSII17 (nr__2) ! (0.1nm) ?=99.99 SII 6717 equivalent width (G2)
real*4 e_WSII17 (nr__2) ! (0.1nm) ?=9.99 Error on WSII 6717
real*4 WSII31 (nr__2) ! (0.1nm) ?=99.99 SII 6731 equivalent width (G2)
real*4 e_WSII31 (nr__2) ! (0.1nm) ?=9.99 Error on WSII 6731
C=============================================================================
C Loading file 'table1.dat' ! Sample overview
C Format for file interpretation
1 format(
+ A7,1X,A13,1X,A13,2X,I3,1X,F4.1,1X,I4,1X,F4.2,1X,F4.1,1X,F4.1,
+ 1X,F5.3,1X,F5.3,1X,F4.2,1X,F4.2,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,65
read(1,'(A89)')ar__
read(ar__,1)
+ Galaxy(i__),MType1(i__),MType2(i__),sigmac(i__),re(i__),
+ HV(i__),rho(i__),Age(i__),e_Age(i__),Z(i__),e_Z(i__),
+ v_a_Fe_(i__),e__a_Fe_(i__),Data(i__)
if (sigmac(i__) .EQ. 45) sigmac(i__) = iNULL__
if (idig(ar__(52:55)).EQ.0) rho(i__) = rNULL__
if(ar__(57:60) .EQ. '') Age(i__) = rNULL__
if(ar__(62:65) .EQ. '') e_Age(i__) = rNULL__
if(ar__(67:71) .EQ. '') Z(i__) = rNULL__
if(ar__(73:77) .EQ. '') e_Z(i__) = rNULL__
if(ar__(79:82) .EQ. '') v_a_Fe_(i__) = rNULL__
if(ar__(84:87) .EQ. '') e__a_Fe_(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Galaxy(i__),MType1(i__),MType2(i__),sigmac(i__),re(i__),
+ HV(i__),rho(i__),Age(i__),e_Age(i__),Z(i__),e_Z(i__),
+ v_a_Fe_(i__),e__a_Fe_(i__),Data(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Emission line fluxes
C Format for file interpretation
2 format(
+ A7,1X,I1,1X,F6.2,1X,F5.2,1X,F6.2,1X,F4.2,1X,F6.2,1X,F4.2,1X,
+ F6.2,1X,F4.2,1X,F6.2,1X,F4.2,1X,F6.2,1X,F4.2,1X,F6.2,1X,F4.2,
+ 1X,F6.2,1X,F4.2,1X,F6.2,1X,F4.2,1X,F6.2,1X,F4.2,1X,F6.2,1X,
+ F4.2,1X,F6.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,212
read(1,'(A154)')ar__1
read(ar__1,2)
+ Galaxy_1(i__),Reg(i__),E_B_V(i__),e_E_B_V(i__),FOII(i__),
+ e_FOII(i__),FNeIII(i__),e_FNeIII(i__),FHb(i__),e_FHb(i__),
+ FOIII59(i__),e_FOIII59(i__),FOIII07(i__),e_FOIII07(i__),
+ FOI(i__),e_FOI(i__),FNII48(i__),e_FNII48(i__),FHa(i__),
+ e_FHa(i__),FNII84(i__),e_FNII84(i__),FSII17(i__),
+ e_FSII17(i__),FSII31(i__),e_FSII31(i__)
c ..............Just test output...........
write(6,2)
+ Galaxy_1(i__),Reg(i__),E_B_V(i__),e_E_B_V(i__),FOII(i__),
+ e_FOII(i__),FNeIII(i__),e_FNeIII(i__),FHb(i__),e_FHb(i__),
+ FOIII59(i__),e_FOIII59(i__),FOIII07(i__),e_FOIII07(i__),
+ FOI(i__),e_FOI(i__),FNII48(i__),e_FNII48(i__),FHa(i__),
+ e_FHa(i__),FNII84(i__),e_FNII84(i__),FSII17(i__),
+ e_FSII17(i__),FSII31(i__),e_FSII31(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Emission line equivalent widths
C Format for file interpretation
3 format(
+ A7,1X,I1,1X,F6.2,1X,F4.2,1X,F6.2,1X,F4.2,1X,F6.2,1X,F4.2,1X,
+ F6.2,1X,F4.2,1X,F6.2,1X,F4.2,1X,F6.2,1X,F4.2,1X,F6.2,1X,F4.2,
+ 1X,F6.2,1X,F4.2,1X,F6.2,1X,F4.2,1X,F6.2,1X,F4.2,1X,F6.2,1X,
+ F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,212
read(1,'(A141)')ar__2
read(ar__2,3)
+ Galaxy_2(i__),Reg_1(i__),WOII(i__),e_WOII(i__),WNeIII(i__),
+ e_WNeIII(i__),WHb(i__),e_WHb(i__),WOIII59(i__),e_WOIII59(i__),
+ WOIII07(i__),e_WOIII07(i__),WOI(i__),e_WOI(i__),WNII48(i__),
+ e_WNII48(i__),WHa(i__),e_WHa(i__),WNII84(i__),e_WNII84(i__),
+ WSII17(i__),e_WSII17(i__),WSII31(i__),e_WSII31(i__)
c ..............Just test output...........
write(6,3)
+ Galaxy_2(i__),Reg_1(i__),WOII(i__),e_WOII(i__),WNeIII(i__),
+ e_WNeIII(i__),WHb(i__),e_WHb(i__),WOIII59(i__),e_WOIII59(i__),
+ WOIII07(i__),e_WOIII07(i__),WOI(i__),e_WOI(i__),WNII48(i__),
+ e_WNII48(i__),WHa(i__),e_WHa(i__),WNII84(i__),e_WNII84(i__),
+ WSII17(i__),e_WSII17(i__),WSII31(i__),e_WSII31(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end
C Locate position of first digit in string; or return 0
integer function idig(c)
character*(*) c
character*1 c1
integer lc,i
lc=len(c)
idig=0
do i=1,lc
if(c(i:i).ne.' ') go to 1
end do
1 if(i.gt.lc) return
c1=c(i:i)
if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1
if(i.gt.lc) return
c1=c(i:i)
if(c1.ge.'0'.and.c1.le.'9') idig=i
return
end