Conversion of standardized ReadMe file for
file /./ftp/cats/VII/169 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. VII/169 Optical Identifications of IRAS Point Sources (Wang+ 1986,87,91)
*================================================================================
*The identifications of IRAS point sources - I. A 304 deg^2 field
*centred on the South Galactic Pole
* Wolstencroft R.D., Savage A., Clowes R.G., MacGillivray H.T.,
* Leggett S.K., Kalafi M.
* <Mon. Not. R. Astron. Soc. 223, 279 (1986)>
* =1986MNRAS.223..279W
*An infrared-optical study of IRAS point sources in the Virgo region
* Leggett S.K., Clowes R.G., Kalafi M., MacGillivray H.T., Puxley P.J.,
* Savage A., Wolstencroft R.D.
* <Mon. Not. R. Astron. Soc. 227, 563 (1987)>
* =1987MNRAS.227..563L
*Optical Identifications of IRAS Point Sources
* Wang G., Clowes R.G., Leggett S.K., MacGillivray H.T., Savage A.
* <MNRAS, 248, 112 (1991)>
* =1991MNRAS.248..112W
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'sgp.dat' ! optical identifications in the
South Galactic Pole region
integer*4 nr__
parameter (nr__=312) ! Number of records
character*140 ar__ ! Full-size record
C B1950 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg (nr__) ! (deg) Right Ascension B1950
real*8 DEdeg (nr__) ! (deg) Declination B1950
C ---------------------------------- ! (position vector(s) in degrees)
character*10 IRAS (nr__) ! [ 0-9-]! IRAS name by IAU format
integer*4 RAh (nr__) ! (h) IRAS right ascension (hours) (B1950)
integer*4 RAm (nr__) ! (min) IRAS right ascension (minutes)
real*4 RAs (nr__) ! (s) IRAS right ascension (seconds)
character*1 DE_ (nr__) ! [-]! declination sign
integer*4 DEd (nr__) ! (deg) IRAS declination (degrees) (B1950)
integer*4 DEm (nr__) ! (arcmin) IRAS declination (minutes)
integer*4 DEs (nr__) ! (arcsec) IRAS declination (seconds)
integer*4 SMAJ (nr__) ! (arcsec) IRAS position uncertainty semi-major axis (1)
integer*4 SMIN (nr__) ! (arcsec) semi-minor axis of above (1)
integer*4 THETA (nr__) ! (deg) [0,180[ orientation of SMAJ, east of
* equatorial north (1)
real*8 F12 (nr__) ! (Jy) [0,] IRAS flux density at 12micron (1)
character*1 l_F12 (nr__) ! [L] upper limit flag of F12
real*8 F25 (nr__) ! (Jy) [0,] IRAS flux density at 25micron (1)
character*1 l_F25 (nr__) ! [L] upper limit flag of F25
real*8 F60 (nr__) ! (Jy) [0,] IRAS flux density at 60micron (1)
character*1 l_F60 (nr__) ! [L] upper limit flag of F60
real*8 F100 (nr__) ! (Jy) [0,] IRAS flux density at 100micron (1)
character*1 l_F100 (nr__) ! [L] upper limit flag of F100
character*11 ANAME (nr__) ! association name in IRAS (1) (2)
integer*4 OFF_RA (nr__) ! (arcsec) ? offset (optical - IRAS) for east of RA
integer*4 OFF_DE (nr__) ! (arcsec) ? offset (optical - IRAS) for north of DE
character*10 TYPE (nr__) ! ! morphological type for galaxy, or
* spectral class for star (3)
real*4 BJmag (nr__) ! (mag) ? for galaxy, B_J magnitude integrated
* above the 25 mag of COSMOS isophoto,
* for star, V magnitude (3)
character*1 f1_BJcos (nr__) ! [(] flag for COBJ (not explained in paper)
real*4 BJcos (nr__) ! (mag) ? B_J magnitude measured by COSMOS
character*1 f2_BJcos (nr__) ! [)] flag for COBJ (not explained in paper)
character*13 NOTES (nr__) ! additional notes (4)
*Note (1): See IRAS Point Source Catalog.
*Note (2): See IRAS Explanatory Supplement V.H.9.
*Note (3): Spectral class and magnitude of the star is taken from SAO or
* Bright Star Catalogues.
*Note (4): For abbreviations, see paper p293. 'R' denotes extra remarks on
* individual sources at the end of the Table in the paper.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'virgo.dat' ! optical identifications in the
Virgo cluster region
integer*4 nr__1
parameter (nr__1=206) ! Number of records
character*130 ar__1 ! Full-size record
C B1950 position composed of: RAh RAm RAds DE- DEd DEm DEs
real*8 RAdeg_1 (nr__1) ! (deg) Right Ascension B1950
real*8 DEdeg_1 (nr__1) ! (deg) Declination B1950
C ---------------------------------- ! (position vector(s) in degrees)
character*10 IRAS_1 (nr__1) ! [ 0-9+]! IRAS name by IAU format
character*1 b1_RADE (nr__1) ! [(] open bracket for IRAS position (1)
integer*4 RAh_1 (nr__1) ! (h) += COSMOS right ascension (h) (B1950) (1)
integer*4 RAm_1 (nr__1) ! (min) COSMOS right ascension (minutes) (1)
integer*4 RAds (nr__1) ! (0.1s) [0,600[ COSMOS right ascension (0.1s) (1)
character*1 DE__1 (nr__1) ! [ ] declination sign (1)
integer*4 DEd_1 (nr__1) ! (deg) COSMOS declination (degrees) (B1950) (1)
integer*4 DEm_1 (nr__1) ! (arcmin) COSMOS declination (minutes) (1)
integer*4 DEs_1 (nr__1) ! (arcsec) COSMOS declination (seconds) (1)
character*1 b2_RADE (nr__1) ! [)] close bracket for IRAS position (1)
integer*4 OFF_RA_1 (nr__1) ! (arcsec) ? offset (optical - IRAS) for RA
integer*4 OFF_DE_1 (nr__1) ! (arcsec) ? offset (optical - IRAS) for DE
character*13 NAME (nr__1) ! name of optical source
character*10 TYPE_1 (nr__1) ! ! morphological classification: galaxies
* by eye following de Vaucouleurs (1959);
* stars from (2)
integer*4 TCODE (nr__1) ! [1,14] type code, see p569 of the paper
real*4 BJmag_1 (nr__1) ! (mag) ? for galaxy, B_J magnitude either from
* catalogues or estimated off the plate
* above the 25 mag of COSMOS isophoto;
* for star, V magnitude (2)
real*4 BJcos_1 (nr__1) ! (mag) ? B_J magnitude measured by COSMOS
real*4 F12_1 (nr__1) ! (Jy) [0,] IRAS flux density at 12micron (3)
character*1 l_F12_1 (nr__1) ! [L:] quality flag of F12 (4)
real*4 F25_1 (nr__1) ! (Jy) [0,] IRAS flux density at 25micron (3)
character*1 l_F25_1 (nr__1) ! [L:] quality flag of F25 (4)
real*4 F60_1 (nr__1) ! (Jy) [0,] IRAS flux density at 60micron (3)
character*1 l_F60_1 (nr__1) ! [L:] quality flag of F60 (4)
real*4 F100_1 (nr__1) ! (Jy) [0,] IRAS flux density at 100micron (3)
character*1 l_F100_1 (nr__1) ! [L:] quality flag of F100 (4)
real*4 logIROP (nr__1) ! ? log. of ratio of far-IR luminosity to
* optical luminosity (5)
real*4 logFIR (nr__1) ! ? log. of far-IR luminosity in units of
* solar lum. for those galaxies with
* measured redshifts and better than
* upper limit detections at both
* 60 and 100microns
character*15 NOTES_1 (nr__1) ! notes of flags given in IRAS Point Source
* Catalogue, and for further remarks (6)
*Note (1): For empty fields the IRAS position is given in brackets.
*Note (2): Spectral class and magnitude of the star is taken from SAO or
* The Bright Star Catalogues. For the few uncatalogued stars,
* B magnitude is estimated from the diffraction spike length on the
* plate.
*Note (3): Not color corrected.
*Note (4): 'L' denotes values that are upper limits only,
* ':' denotes moderate quality fluxes.
*Note (5): Ratios are given for those galaxies with better than upper limit
* detections at both 60 and 100microns, and that have COSMOS
* B-magnitudes.
*Note (6): For abbreviations, see paper p569. 'R' denotes extra remarks on
* individual sources at the end of the Table in the paper, p569.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'fohycom.dat' ! optical identifications in the
Fornax, Hydra I, Coma cluster regions
integer*4 nr__2
parameter (nr__2=241) ! Number of records
character*118 ar__2 ! Full-size record
C B1950 position composed of: RAh RAm RAds DE- DEd DEm DEs
real*8 RAdeg_2 (nr__2) ! (deg) Right Ascension B1950
real*8 DEdeg_2 (nr__2) ! (deg) Declination B1950
C ---------------------------------- ! (position vector(s) in degrees)
character*10 IRAS_2 (nr__2) ! [ 0-9x+-]! IRAS name by IAU format (in
* Coma part, "x1247+257" is written once)
integer*4 RAh_2 (nr__2) ! (h) += COSMOS right ascension (h) (B1950) (1)
integer*4 RAm_2 (nr__2) ! (min) COSMOS right ascension (minutes) (1)
integer*4 RAds_1 (nr__2) ! (0.1s) [0,600[ COSMOS right ascension (0.1s) (1)
character*1 DE__2 (nr__2) ! [-+]! declination sign (1)
integer*4 DEd_2 (nr__2) ! (deg) COSMOS declination (degrees) (B1950) (1)
integer*4 DEm_2 (nr__2) ! (arcmin) COSMOS declination (minutes) (1)
integer*4 DEs_2 (nr__2) ! (arcsec) COSMOS declination (seconds) (1)
integer*4 OFF_RA_2 (nr__2) ! (arcsec) ? offset (optical - IRAS) for RA
integer*4 OFF_DE_2 (nr__2) ! (arcsec) ? offset (optical - IRAS) for DE
real*4 F12_2 (nr__2) ! (Jy) [0,] IRAS flux density at 12micron
character*1 l_F12_2 (nr__2) ! [L:] quality flag of F12 (2)
real*4 F25_2 (nr__2) ! (Jy) [0,] IRAS flux density at 25micron
character*1 l_F25_2 (nr__2) ! [L:] quality flag of F25 (2)
real*4 F60_2 (nr__2) ! (Jy) [0,] IRAS flux density at 60micron
character*1 l_F60_2 (nr__2) ! [L:] quality flag of F60 (2)
real*4 F100_2 (nr__2) ! (Jy) [0,] IRAS flux density at 100micron
character*1 l_F100_2 (nr__2) ! [L:] quality flag of F100 (2)
character*8 TYPE_2 (nr__2) ! ! for galaxy, morphological
* classification followed de Vaucouleurs
* (1959);
* for star, spectral classification
integer*4 TCODE_1 (nr__2) ! [1-14] type code, see p113R of the paper
real*4 BJmag_2 (nr__2) ! (mag) ? for galaxy, COSMOS B_J magnitude;
* for star, B magnitude (see paper p113R)
character*1 q_BJmag (nr__2) ! [e] 'e' denotes that B magnitude for
* star is estimated from the catalogued
* V and the B-V (see paper p113R)
real*4 Bmag_cat (nr__2) ! (mag) ? catalogue magnitude:
* for galaxy, the catalogued B, or B_J
* estimated from plates by comparing
* the images with known mag.;
* for star, catalogued (Bright Star or SAO)
* V magnitude
character*11 NAME_1 (nr__2) ! name in optical catalogue, such as NGC,
* ESO, SAO, BS, etc.
character*17 NOTES_2 (nr__2) ! notes of flags given in IRAS Point Source
* Catalogue, and for further remarks (3)
*Note (1): COSMOS positions are measured from UKST plates (see paper
* p113L). For some objects on the edges of plates, catalogue positions
* are given. For empty fields the IRAS Point Source Cat. position is
* given.
*Note (2): 'L' denotes values that are upper limits only,
* ':' denotes moderate quality fluxes.
*Note (3): For abbreviations, see paper p113R. 'R' denotes extra remarks on
* individual sources at the end of the Table in the paper, p117.
C=============================================================================
C Loading file 'sgp.dat' ! optical identifications in the
* South Galactic Pole region
C Format for file interpretation
1 format(
+ A10,2X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,I3,1X,I2,1X,
+ I2,1X,F9.4,A1,F9.4,A1,F9.4,A1,F9.4,A1,1X,A11,I3,1X,I3,2X,A10,
+ F5.2,A1,F5.2,A1,1X,A13)
C Effective file loading
open(unit=1,status='old',file=
+'sgp.dat')
write(6,*) '....Loading file: sgp.dat'
do i__=1,312
read(1,'(A140)')ar__
read(ar__,1)
+ IRAS(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),SMAJ(i__),SMIN(i__),THETA(i__),F12(i__),
+ l_F12(i__),F25(i__),l_F25(i__),F60(i__),l_F60(i__),F100(i__),
+ l_F100(i__),ANAME(i__),OFF_RA(i__),OFF_DE(i__),TYPE(i__),
+ BJmag(i__),f1_BJcos(i__),BJcos(i__),f2_BJcos(i__),NOTES(i__)
if(ar__(96:98) .EQ. '') OFF_RA(i__) = iNULL__
if(ar__(100:102) .EQ. '') OFF_DE(i__) = iNULL__
if(ar__(115:119) .EQ. '') BJmag(i__) = rNULL__
if(ar__(121:125) .EQ. '') BJcos(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,1)
+ IRAS(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),SMAJ(i__),SMIN(i__),THETA(i__),F12(i__),
+ l_F12(i__),F25(i__),l_F25(i__),F60(i__),l_F60(i__),F100(i__),
+ l_F100(i__),ANAME(i__),OFF_RA(i__),OFF_DE(i__),TYPE(i__),
+ BJmag(i__),f1_BJcos(i__),BJcos(i__),f2_BJcos(i__),NOTES(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 'virgo.dat' ! optical identifications in the
* Virgo cluster region
C Format for file interpretation
2 format(
+ A10,1X,A1,I2,I2,I3,A1,I2,I2,I2,A1,I3,1X,I3,1X,A13,1X,A10,1X,
+ I2,2X,F5.2,1X,F5.2,1X,F5.2,A1,1X,F5.2,A1,1X,F5.2,A1,1X,F5.2,
+ A1,1X,F5.2,1X,F4.1,1X,A15)
C Effective file loading
open(unit=1,status='old',file=
+'virgo.dat')
write(6,*) '....Loading file: virgo.dat'
do i__=1,206
read(1,'(A130)')ar__1
read(ar__1,2)
+ IRAS_1(i__),b1_RADE(i__),RAh_1(i__),RAm_1(i__),RAds(i__),
+ DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),b2_RADE(i__),
+ OFF_RA_1(i__),OFF_DE_1(i__),NAME(i__),TYPE_1(i__),TCODE(i__),
+ BJmag_1(i__),BJcos_1(i__),F12_1(i__),l_F12_1(i__),F25_1(i__),
+ l_F25_1(i__),F60_1(i__),l_F60_1(i__),F100_1(i__),
+ l_F100_1(i__),logIROP(i__),logFIR(i__),NOTES_1(i__)
if(ar__1(28:30) .EQ. '') OFF_RA_1(i__) = iNULL__
if(ar__1(32:34) .EQ. '') OFF_DE_1(i__) = iNULL__
if(ar__1(65:69) .EQ. '') BJmag_1(i__) = rNULL__
if(ar__1(71:75) .EQ. '') BJcos_1(i__) = rNULL__
if(ar__1(105:109) .EQ. '') logIROP(i__) = rNULL__
if(ar__1(111:114) .EQ. '') logFIR(i__) = rNULL__
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(RAds(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAds(i__)/2400.
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,2)
+ IRAS_1(i__),b1_RADE(i__),RAh_1(i__),RAm_1(i__),RAds(i__),
+ DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),b2_RADE(i__),
+ OFF_RA_1(i__),OFF_DE_1(i__),NAME(i__),TYPE_1(i__),TCODE(i__),
+ BJmag_1(i__),BJcos_1(i__),F12_1(i__),l_F12_1(i__),F25_1(i__),
+ l_F25_1(i__),F60_1(i__),l_F60_1(i__),F100_1(i__),
+ l_F100_1(i__),logIROP(i__),logFIR(i__),NOTES_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 'fohycom.dat' ! optical identifications in the
* Fornax, Hydra I, Coma cluster regions
C Format for file interpretation
3 format(
+ A10,1X,I2,I2,I3,1X,A1,I2,I2,I2,1X,I3,1X,I3,2X,F5.2,A1,1X,F5.2,
+ A1,1X,F5.2,A1,F6.2,A1,1X,A8,1X,I2,1X,F5.2,A1,1X,F5.2,1X,A11,
+ 1X,A17)
C Effective file loading
open(unit=1,status='old',file=
+'fohycom.dat')
write(6,*) '....Loading file: fohycom.dat'
do i__=1,241
read(1,'(A118)')ar__2
read(ar__2,3)
+ IRAS_2(i__),RAh_2(i__),RAm_2(i__),RAds_1(i__),DE__2(i__),
+ DEd_2(i__),DEm_2(i__),DEs_2(i__),OFF_RA_2(i__),OFF_DE_2(i__),
+ F12_2(i__),l_F12_2(i__),F25_2(i__),l_F25_2(i__),F60_2(i__),
+ l_F60_2(i__),F100_2(i__),l_F100_2(i__),TYPE_2(i__),
+ TCODE_1(i__),BJmag_2(i__),q_BJmag(i__),Bmag_cat(i__),
+ NAME_1(i__),NOTES_2(i__)
if(ar__2(28:30) .EQ. '') OFF_RA_2(i__) = iNULL__
if(ar__2(32:34) .EQ. '') OFF_DE_2(i__) = iNULL__
if(ar__2(77:81) .EQ. '') BJmag_2(i__) = rNULL__
if(ar__2(84:88) .EQ. '') Bmag_cat(i__) = rNULL__
RAdeg_2(i__) = rNULL__
DEdeg_2(i__) = rNULL__
c Derive coordinates RAdeg_2 and DEdeg_2 from input data
c (RAdeg_2 and DEdeg_2 are set to rNULL__ when unknown)
if(RAh_2(i__) .GT. -180) RAdeg_2(i__)=RAh_2(i__)*15.
if(RAm_2(i__) .GT. -180) RAdeg_2(i__)=RAdeg_2(i__)+RAm_2(i__)/4.
if(RAds_1(i__) .GT. -180) RAdeg_2(i__)=RAdeg_2(i__)+RAds_1(i__)/2400.
if(DEd_2(i__) .GE. 0) DEdeg_2(i__)=DEd_2(i__)
if(DEm_2(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEm_2(i__)/60.
if(DEs_2(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEs_2(i__)/3600.
if(DE__2(i__).EQ.'-'.AND.DEdeg_2(i__).GE.0) DEdeg_2(i__)=-DEdeg_2(i__)
c ..............Just test output...........
write(6,3)
+ IRAS_2(i__),RAh_2(i__),RAm_2(i__),RAds_1(i__),DE__2(i__),
+ DEd_2(i__),DEm_2(i__),DEs_2(i__),OFF_RA_2(i__),OFF_DE_2(i__),
+ F12_2(i__),l_F12_2(i__),F25_2(i__),l_F25_2(i__),F60_2(i__),
+ l_F60_2(i__),F100_2(i__),l_F100_2(i__),TYPE_2(i__),
+ TCODE_1(i__),BJmag_2(i__),q_BJmag(i__),Bmag_cat(i__),
+ NAME_1(i__),NOTES_2(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg_2(i__),DEdeg_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end