Conversion of standardized ReadMe file for
file /./ftp/cats/VII/115 into FORTRAN code for reading data files line by line.
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.5, on 2013-Jun-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__=-1.e37) ! NULL real number
parameter (iNULL__=-2147483647) ! NULL int number
C=============================================================================
Cat. VII/115 Surface Photometry of ESO-Uppsala Galaxies (Lauberts+ 1989)
*================================================================================
*The Surface Photometry Catalogue of the ESO-Uppsala Galaxies
* Lauberts A., Valentijn E.A.
* <The Surface Photometry Catalogue of the ESO/Uppsala Galaxies, European
* Southern Observatory (1989)>
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'esolv1.dat' ! *The ESO-LV catalog (most important parameters)
integer*4 nr__
parameter (nr__=15467) ! Number of records
character*474 ar__ ! Full-size record
real*8 RA_ ! (deg) Right Ascension B1950
real*8 Dec ! (deg) Declination B1950
integer*4 Seq ! Running number Identification
real*4 TType ! Revised Morphological TType (T) (1)
integer*4 r_TType ! [1-4] Source of TType: 1=visual revised,
* 2=Corwin, 3=old, 4=automated (section 6) (S)
integer*4 RAh ! (h) Right Ascension B1950 (hours) (R_A)
integer*4 RAm ! (min) Right Ascension B1950 (minutes)
integer*4 RAcs ! (cs) Right Ascension B1950 (0.01seconds)
character*1 DE_ ! Declination B1950 (sign) (DEC)
integer*4 DEd ! (deg) Declination B1950 (degrees)
integer*4 DEm ! (arcmin) Declination B1950 (minutes)
integer*4 DEds ! (100mas) Declination B1950 (0.1seconds)
integer*4 Bf ! [-6,2]? Flag on B (B_FLAG) (3)
integer*4 Rf ! [-6,2]? Flag on R (R_FLAG) (3)
integer*4 ESO_LV ! Identification in catalog (ESO) (2)
integer*4 NGC ! ? Number in NGC catalog (NGC)
integer*4 IC ! ? Number in IC catalog (IC)
real*8 GLON ! (deg) Galactic longitude (L_II)
real*8 GLAT ! (deg) Galactic latitude (B_II)
real*4 cz ! (km/s) ? Heliocentric radial velocity (CZ)
real*4 e_cz ! (km/s) ? Mean error on cz (E_CZ)
real*4 BTmag ! (mag) ? Total B magnitude (B_TOT)
real*4 RTmag ! (mag) ? Total R magnitude (R_TOT)
real*4 B26mag ! (mag) ? B mag. in 26mag/arcsec2 isophote (B_B26)
real*4 R26mag ! (mag) ? R mag. in 26mag/arcsec2 isophote (R_B26)
real*4 B25mag ! (mag) ? B mag. in 25mag/arcsec2 isophote (B_B25)
real*4 R25mag ! (mag) ? R mag. in 25mag/arcsec2 isophote (R_B25)
real*4 v_B_R_e ! (mag) ? B-R color in half-light ellipse (BR_E)
real*4 v_B_R_10 ! (mag) ? B-R color in 10" diameter circle (BR_10)
real*4 v_B_R_T ! (mag) ? B-R total color (BR_T)
real*4 mu0_B ! (mag/arcsec2) ? Average blue central surface brightness
* within 10" diameter circular aperture (BS_0)
real*4 D25_B ! (arcsec) ? Isophotal 25mag/arcsec2 Major Diameter (D25)
real*4 D26_B ! (arcsec) ? Isophotal 26mag/arcsec2 Major Diameter (D26)
real*4 D27_B ! (arcsec) ? Isophotal 27mag/arcsec2 Major Diameter (D27)
real*4 a_b_1 ! ? a/b from ellipse fit to B octants (AB_OCT)
real*4 PA_1 ! (deg) ? Position angle of B ellipse fit (PA_OCT)
real*4 a_b_org ! ? a/b from original ESO-Uppsala (AB_ORG)
real*4 a_b_0 ! ? a/b of adopted ellipse (AB_ADO)
real*4 PA_org ! (deg) ? Original position angle (PA_ORG)
real*4 PA_0 ! (deg) ? Adopted Position angle (PA_ADO)
real*4 D70_B ! (arcsec) ? Major diameter at 70% B light (D70)
real*4 D80_B ! (arcsec) ? Major diameter at 80% B light (D80)
real*4 D90_B ! (arcsec) ? Major diameter at 90% B light (D90)
real*4 De_B ! (arcsec) ? Major diameter at half light (D_E)
real*4 Dorg_B ! (arcsec) ? Major diameter, original (DB_EU)
real*4 DN_B ! (arcsec) ? Diameter of circle within which the mean B
* surface brightness is 20.75mag/arcsec2
real*4 Ae_B ! (arcsec) ? Aperture of circle enclosing half total B
* light (A_E)
real*4 Ntot ! (deg-2) ? Total density of galaxies (NG_T)
real*4 mue_B ! (mag/arcsec2) ? B surface brightness at half B total light
* (only if De(B)>5arcsec) (BS_E)
real*4 mue_R ! (mag/arcsec2) ? R surface brightness at half B total light
* (only if De(B)>5arcsec) (RS_E)
real*4 mu0_R ! (mag/arcsec2) ? Average red central surface brightness
* within 10" diameter circular aperture (RS_0)
real*4 Noct_B ! ? deVaucouleurs exponent fit in B (NOCT_B) (5)
real*4 Noct_R ! ? deVaucouleurs exponent fit to R (NOCT_R) (5)
real*4 D_B_R ! ? Difference B-R of regions (outside-inside)
* B effective radius (BR_D)
real*4 Abh ! (mag) ? Burstein-Heiles extinction in B (AEXT)
integer*4 C ! [0,2] Sample parameter (C) (4)
character*52 Comment ! Original classification and comments
*Note (1): TType is a number corresponding to classifications:
* -------------------------------------
* TType Hubble Corwin ESO-Uppsala
* -------------------------------------
* -5 E E E
* -4 cD
* -3 E-S0 L- E-S0
* -2 S0 L S0
* -1 L+
* 0 S0/a S0 S0-a
* 1 Sa S1 Sa
* 2 Sa-b S2 Sa-b
* 3 Sb S3 Sb
* 4 Sb-c S4 Sb-c
* 5 S5 S...
* 6 Sc S6 Sc, Sc-d
* 7 S7 S../Irr
* 8 Sc-Irr S8 Sd
* 9 S9
* 10 Irr I Irr
* 11 unknown
* -------------------------------------
*Note (2): the ESO-LV number is made of 3 digits for the field name,
* 3 digits for the running number in the field,
* and a last digit for companions.
*Note (3): a positive number indicates the number of radii suspected
* to contain saturated central parts; negative values indicate the
* number of radii where photographic data have been replaced with
* photoelectric measurements.
*Note (4): sample parameter is the following:
* 0 = not complete sample;
* 1 = ESO-Uppsala catalog of objects with diameter>1'
* 2 = Vmax subsample of 1 (see section 9)
*Note (5): de Vaucouleurs law assumes the fits in surface brightness
* I = Io exp(-(r/{alpha})^N^)
* The exponent N is contained in the catalog; the FITS version includes
* also the scale length {alpha}
C=============================================================================
C Loading file 'esolv1.dat' ! *The ESO-LV catalog (most important parameters)
C Format for file interpretation
1 format(
+ 4X,I5,3X,F4.1,4X,I1,1X,I2,I2,I4,1X,A1,I2,I2,I3,2X,I2,3X,I2,4X,
+ I7,2X,I5,2X,I5,2X,F7.3,2X,F7.3,2X,F6.0,3X,F4.0,4X,F5.2,4X,
+ F5.2,4X,F5.2,4X,F5.2,5X,F5.2,4X,F5.2,4X,F5.2,4X,F5.2,4X,F5.2,
+ 3X,F5.2,4X,F6.1,3X,F6.1,3X,F6.1,2X,F5.2,4X,F5.1,4X,F5.2,4X,
+ F5.2,4X,F5.1,4X,F5.1,3X,F6.1,3X,F6.1,5X,F6.1,3X,F6.1,3X,F6.1,
+ 4X,F5.1,3X,F6.1,4X,F6.3,4X,F5.2,4X,F5.2,5X,F5.2,4X,F5.3,4X,
+ F5.3,4X,F5.2,3X,F6.3,4X,I1,2X,A52)
C Effective file loading
open(unit=1,file='esolv1.dat', status='old')
write(6,*) '....Loading file: esolv1.dat'
do i__=1,15467
read(1,'(A474)')ar__
read(ar__,1)
+ Seq,TType,r_TType,RAh,RAm,RAcs,DE_,DEd,DEm,DEds,Bf,Rf,ESO_LV,
+ NGC,IC,GLON,GLAT,cz,e_cz,BTmag,RTmag,B26mag,R26mag,B25mag,
+ R25mag,v_B_R_e,v_B_R_10,v_B_R_T,mu0_B,D25_B,D26_B,D27_B,a_b_1,
+ PA_1,a_b_org,a_b_0,PA_org,PA_0,D70_B,D80_B,D90_B,De_B,Dorg_B,
+ DN_B,Ae_B,Ntot,mue_B,mue_R,mu0_R,Noct_B,Noct_R,D_B_R,Abh,C,
+ Comment
if(ar__(42:43) .EQ. '') Bf = iNULL__
if(ar__(47:48) .EQ. '') Rf = iNULL__
if(ar__(62:66) .EQ. '') NGC = iNULL__
if(ar__(69:73) .EQ. '') IC = iNULL__
if(ar__(94:99) .EQ. '') cz = rNULL__
if(ar__(103:106) .EQ. '') e_cz = rNULL__
if(ar__(111:115) .EQ. '') BTmag = rNULL__
if(ar__(120:124) .EQ. '') RTmag = rNULL__
if(ar__(129:133) .EQ. '') B26mag = rNULL__
if(ar__(138:142) .EQ. '') R26mag = rNULL__
if(ar__(148:152) .EQ. '') B25mag = rNULL__
if(ar__(157:161) .EQ. '') R25mag = rNULL__
if(ar__(166:170) .EQ. '') v_B_R_e = rNULL__
if(ar__(175:179) .EQ. '') v_B_R_10 = rNULL__
if(ar__(184:188) .EQ. '') v_B_R_T = rNULL__
if(ar__(192:196) .EQ. '') mu0_B = rNULL__
if(ar__(201:206) .EQ. '') D25_B = rNULL__
if(ar__(210:215) .EQ. '') D26_B = rNULL__
if(ar__(219:224) .EQ. '') D27_B = rNULL__
if(ar__(227:231) .EQ. '') a_b_1 = rNULL__
if(ar__(236:240) .EQ. '') PA_1 = rNULL__
if(ar__(245:249) .EQ. '') a_b_org = rNULL__
if(ar__(254:258) .EQ. '') a_b_0 = rNULL__
if(ar__(263:267) .EQ. '') PA_org = rNULL__
if(ar__(272:276) .EQ. '') PA_0 = rNULL__
if(ar__(280:285) .EQ. '') D70_B = rNULL__
if(ar__(289:294) .EQ. '') D80_B = rNULL__
if(ar__(300:305) .EQ. '') D90_B = rNULL__
if(ar__(309:314) .EQ. '') De_B = rNULL__
if(ar__(318:323) .EQ. '') Dorg_B = rNULL__
if(ar__(328:332) .EQ. '') DN_B = rNULL__
if(ar__(336:341) .EQ. '') Ae_B = rNULL__
if(ar__(346:351) .EQ. '') Ntot = rNULL__
if(ar__(356:360) .EQ. '') mue_B = rNULL__
if(ar__(365:369) .EQ. '') mue_R = rNULL__
if(ar__(375:379) .EQ. '') mu0_R = rNULL__
if(ar__(384:388) .EQ. '') Noct_B = rNULL__
if(ar__(393:397) .EQ. '') Noct_R = rNULL__
if(ar__(402:406) .EQ. '') D_B_R = rNULL__
if(ar__(410:415) .EQ. '') Abh = rNULL__
c Derive coordinates RA_ and Dec from input data
c (RA_ and Dec are set to rNULL_ when unknown)
RA_ = RAh
if(RAh .GE. 0) RA_=RAh*15.
if(RAm .GE. 0) RA_=RA_+RAm/4.
if(RAcs .GE. 0) RA_=RA_+RAcs/24000.
Dec = DEd
if(DEm .GE. 0) Dec=Dec+DEm/60.
if(DEds .GE. 0) Dec=Dec+DEds/36000.
if(DE_.EQ.'-'.AND.Dec.GE.0) Dec=-Dec
c ..............Just test output...........
write(6,1)
+ Seq,TType,r_TType,RAh,RAm,RAcs,DE_,DEd,DEm,DEds,Bf,Rf,ESO_LV,
+ NGC,IC,GLON,GLAT,cz,e_cz,BTmag,RTmag,B26mag,R26mag,B25mag,
+ R25mag,v_B_R_e,v_B_R_10,v_B_R_T,mu0_B,D25_B,D26_B,D27_B,a_b_1,
+ PA_1,a_b_org,a_b_0,PA_org,PA_0,D70_B,D80_B,D90_B,De_B,Dorg_B,
+ DN_B,Ae_B,Ntot,mue_B,mue_R,mu0_R,Noct_B,Noct_R,D_B_R,Abh,C,
+ Comment
write(6,'(6H Pos: 2F8.4)') RA_,Dec
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end