Conversion of standardized ReadMe file for
catalog 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.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. VII/61A Catalogue of Markarian Galaxies (Markarian+ 1967-81)
*================================================================================
*Catalogue of Markarian Galaxies
* Markarian, B. E.
* <Astrofizika 3, 55 (1967)>
* <Astrofizika 5, 443 (1969)>
* <Astrofizika 5, 581 (1969)>
* Markarian, B. E., Lipovetskii, V. A.
* <Astrofizika 7, 511 (1971)>
* <Astrofizika 8, 155 (1972)>
* <Astrofizika 9, 487 (1973)>
* <Astrofizika 10, 307 (1974)>
* <Astrofizika 12, 389 (1976)>
* <Astrofizika 12, 657 (1976)>
* Markarian, B. E., Lipovetskii, V. A., Stepanian, Dzh. A.
* <Astrofizika 13, 225 (1977)>
* <Astrofizika 13, 397 (1977)>
* <Astrofizika 15, 201 (1979)>
* <Astrofizika 15, 363 (1979)>
* <Astrofizika 15, 549 (1979)>
* <Astrofizika 17, 619 (1981)>
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'catalog.dat' ! The Catalogue
integer*4 nr__
parameter (nr__=1525) ! Number of records
character*74 ar__ ! Full-size record
C 1950.0 position composed of: RAh RAm DE- DEd DEm
real*8 RAdeg ! (deg) Right Ascension 1950.0
real*8 DEdeg ! (deg) Declination 1950.0
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 Mkn ! [1,1500] Markarian ID
character*1 m_Mkn ! [A] for additional number
character*1 Active ! [QS] 'Q' for QSO, 'S' for Active galaxy
character*14 Name1 ! Cross ID 1: MCG, NGC, IC, or HZ numbers
character*10 Name2 ! Cross ID 2: MCG, NGC, IC, or HZ numbers
integer*4 RAh ! (h) Hours RA, 1950.0
real*4 RAm ! (min) Minutes RA, 1950.0
character*1 DE_ ! Sign Dec, 1950.0
integer*4 DEd ! (deg) Degrees Dec, 1950.0
integer*4 DEm ! (arcmin) Minutes Dec, 1950.0
integer*4 MajAxis ! (arcsec) Major axis
integer*4 MinAxis ! (arcsec) Minor axis
real*8 Mag ! (mag) Estimated mag. from objective prism spectrum
character*6 SpType ! Spectral type, see adc.doc for description
C=============================================================================
C Loading file 'catalog.dat' ! The Catalogue
C Format for file interpretation
1 format(4X,I4,A1,A1,1X,A14,A10,I2,F6.2,A1,I2,I3,I5,I5,F7.2,2X,A6)
C Effective file loading
open(unit=1,status='old',file=
+'catalog.dat')
write(6,*) '....Loading file: catalog.dat'
do i__=1,1525
read(1,'(A74)')ar__
read(ar__,1)
+ Mkn,m_Mkn,Active,Name1,Name2,RAh,RAm,DE_,DEd,DEm,MajAxis,
+ MinAxis,Mag,SpType
RAdeg = rNULL__
DEdeg = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh .GT. -180) RAdeg=RAh*15.
if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
if(DEd .GE. 0) DEdeg=DEd
if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c ..............Just test output...........
write(6,1)
+ Mkn,m_Mkn,Active,Name1,Name2,RAh,RAm,DE_,DEd,DEm,MajAxis,
+ MinAxis,Mag,SpType
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end