Conversion of standardized ReadMe file for
file /./ftp/cats/VII/221 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. VII/221 PSCz catalog (Saunders+, 2000)
*================================================================================
*The PSCz catalogue
* Saunders W., Sutherland W.J., Maddox S.J., Keeble O., Oliver S.J.,
* Rowan-Robinson M., McMahon R.G., Efstathiou G.P., Tadros H., White S.D.M.,
* Frenk C.S., Carraminana A., Hawkins M.R.S.
* <Mon. Not. R. Astron. Soc. 317, 55 (2000)>
* =2000MNRAS.317...55S
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'pscz.dat' ! Main catalogue
integer*4 nr__
parameter (nr__=18351) ! Number of records
character*533 ar__ ! Full-size record
C B1950 position composed of: RAdeg DEdeg
character*11 CNAME (nr__) ! *Name e.g. Q12345+4856
integer*4 RAh (nr__) ! (h) Right ascension B1950
integer*4 RAm (nr__) ! (min) Right ascension B1950
integer*4 RAds (nr__) ! (0.1s) Right ascension B1950
character*1 DE_ (nr__) ! Declination sign B1950
integer*4 DEd (nr__) ! (deg) Declination B1950
integer*4 DEm (nr__) ! (arcmin) Declination B1950
integer*4 DEs (nr__) ! (arcsec) Declination B1950
real*8 S12 (nr__) ! (Jy) Flux at 12um
real*8 S25 (nr__) ! (Jy) Flux at 25um
real*8 S60 (nr__) ! (Jy) Flux at 60um
real*8 S100 (nr__) ! (Jy) Flux at 100um
integer*4 q_S12 (nr__) ! *[1,9] Flux origin and quality flag at 12um
integer*4 q_S25 (nr__) ! *[1,9] Flux origin and quality flag at 25um
integer*4 q_S60 (nr__) ! *[1,9] Flux origin and quality flag at 60um
integer*4 q_S100 (nr__) ! *[1,9] Flux origin and quality flag at 100um
integer*4 SES12 (nr__) ! [0,8] Small Extended Source flag at 12um
integer*4 SES25 (nr__) ! [0,8] Small Extended Source flag at 25um
integer*4 SES60 (nr__) ! [0,8] Small Extended Source flag at 60um
integer*4 SES100 (nr__) ! [0,8] Small Extended Source flag at 100um
integer*4 e_S12 (nr__) ! (%) Flux uncertainty at 12um
integer*4 e_S25 (nr__) ! (%) Flux uncertainty at 25um
integer*4 e_S60 (nr__) ! (%) Flux uncertainty at 60um
integer*4 e_S100 (nr__) ! (%) Flux uncertainty at 100um
integer*4 PosMin (nr__) ! (arcsec) Minor axis of 2{sigma} error ellipse
integer*4 PosMax (nr__) ! (arcsec) Major axis of 2{sigma} error ellipse
integer*4 PosPA (nr__) ! (deg) Position angle E of N of major axis
* of error ellipse
integer*4 Cirr1 (nr__) ! *Cirrus flag
integer*4 Cirr2 (nr__) ! *Cirrus flag
integer*4 Cirr3 (nr__) ! (MJy/sr) *Cirrus flag
character*1 CC1 (nr__) ! [A-N X] Correlation Coefficient at 12um
character*1 CC2 (nr__) ! [A-N X] Correlation Coefficient at 25um
character*1 CC3 (nr__) ! [A-M X1] Correlation Coefficient at 60um
character*1 CC4 (nr__) ! [A-N X1] Correlation Coefficient at 100um
character*1 CONF (nr__) ! [0-9 ACDEF] Confusion flag, hex encoded
character*1 DISC (nr__) ! [0-9 A-D] Discrepant flux flag, hex encoded
character*1 HSD (nr__) ! [0-8] High Source Density flag, hex encoded
integer*4 PNEARH (nr__) ! [0-9] Nearby hours-confirmed sources
integer*4 PNEARW (nr__) ! [0-5] Nearby weeks-confirmed sources
integer*4 CES12 (nr__) ! *Nearby seconds confirmed extended sources
* at 12um
integer*4 CES25 (nr__) ! *Nearby seconds confirmed extended sources
* at 25um
integer*4 CES60 (nr__) ! *Nearby seconds confirmed extended sources
* at 60um
integer*4 CES100 (nr__) ! *Nearby seconds confirmed extended sources
* at 100um
integer*4 SNR60 (nr__) ! 10 x signal/noise of weakest detection
* at 60um
integer*4 NHCON (nr__) ! Number of hours-confirmed source detections
integer*4 MHCON (nr__) ! Potential number of hours-confirmed source
* detections
integer*4 FCOR60 (nr__) ! *1000*flux correction factor applied at 60um
character*5 PLATE (nr__) ! Plate number (POSS or SRC)
character*7 NOM (nr__) ! Nominal plate centre in RA,dec
integer*4 NUM (nr__) ! Identifier for source on this plate
integer*4 RAo_h (nr__) ! (h) ? Optical right ascension B1950
integer*4 RAo_m (nr__) ! (min) ? Optical right ascension B1950
real*4 RAo_s (nr__) ! (s) ? Optical right ascension B1950
character*1 DEo__ (nr__) ! ? Optical declination sign B1950
integer*4 DEo_d (nr__) ! (deg) ? Optical declination B1950
integer*4 DEo_m (nr__) ! (arcmin) ? Optical declination B1950
real*4 DEo_s (nr__) ! (arcsec) ? Optical declination B1950
integer*4 ADX (nr__) ! (arcsec) ? North/South Delta Optical/IRAS position
integer*4 ADY (nr__) ! (arcsec) ? East/West Delta Optical/IRAS position
real*4 ALK (nr__) ! ? log_10 likelihood for optical
* identification
integer*4 REL (nr__) ! (%) ? Prob. correct optical identification
character*1 CAND (nr__) ! ? Optical candidate A,B etc
integer*4 IOPT (nr__) ! *? Schmidt plate material used
real*4 Bmag (nr__) ! (mag) ? Estimated optical B magnitude
integer*4 MajAxiO (nr__) ! (arcsec) ? Optical major axis
integer*4 MinAxiO (nr__) ! (arcsec) ? Optical minor axis
integer*4 PAO (nr__) ! (deg) ? Optical major axis position angle E of N
real*8 SuBr (nr__) ! (mag/arcsec+2) ? Average optical surface brightness
integer*4 IDTYPE (nr__) ! Cross-identification type from PSC
character*12 Oname (nr__) ! Optical name
real*4 mag (nr__) ! (mag) Zwicky/RC3 magnitude
real*4 D1X (nr__) ! (arcmin) *RC3/UGC/MCG/ESO D25 diameter from CATX
integer*4 AR (nr__) ! (%) *Axial ratio d25/D25 from CATX
integer*4 Hvel (nr__) ! (km/s) Best available heliocentric velocity (cz)
integer*4 e_Hvel (nr__) ! (km/s) Error on heliocentric velocity
integer*4 r_Hvel (nr__) ! *Reference for heliocentric velocity
integer*4 Hvel2 (nr__) ! (km/s) Heliocentric velocity (cz)
* measured as part of PSCz survey
integer*4 e_Hvel2 (nr__) ! (km/s) Error on PSCz heliocentric velocity
integer*4 r_Hvel2 (nr__) ! *Reference for PSCz heliocentric velocity
integer*4 TType (nr__) ! De Vaucouleurs Type, 99 = unknown
character*8 PGC (nr__) ! Catalogue of Principal Galaxies
* (PGC, Cat. <VII/119>) name
real*8 RAdeg (nr__) ! (deg) Right ascension B1950 (IRAS source)
real*8 DEdeg (nr__) ! (deg) Declination B1950 (IRAS source)
real*8 GLON (nr__) ! (deg) Galactic longitude
real*8 GLAT (nr__) ! (deg) Galactic latitude
integer*4 cz (nr__) ! (km/s) Velocity cz in Local Group frame
real*4 Dist (nr__) ! (Mpc) Distance from Rowan-Robinson 1988
integer*4 IDclus (nr__) ! Cluster identification from
* Rowan-Robinson 1988
real*4 ALI100 (nr__) ! ([MJy/sr]) *log10(100um background in MJy/sr)
real*4 ALI100B (nr__) ! ([MJy/sr]) *log10(100um background in MJy/sr
* for lune bin)
integer*4 NLBIN (nr__) ! Lune bin number
character*2 CLASS (nr__) ! *Optical classification code
character*5 CSEC (nr__) ! PSCz sector, e.g. A/123
integer*4 Xpos (nr__) ! (mm) X position from nominal bottom left corner
* of Schmidt plate
integer*4 Ypos (nr__) ! (mm) Y position from nominal bottom left corner
* of Schmidt plate
character*1 CSTAT (nr__) ! [JORXPCELS?0] Optical material used
character*1 CAREA (nr__) ! h=high|b|, p=other PSCz,
* l=BTP, c=in coverage gaps
character*1 CZSTAT (nr__) ! *[#-1234lpxq] Redshift status
character*1 C12CLAS (nr__) ! Classification in 1.2Jy survey
real*8 F3PSC (nr__) ! (Jy) PSC flux at 60um
real*8 F4PSC (nr__) ! (Jy) PSC flux at 100um
real*4 ABB (nr__) ! (mag) *Estimated average B extinction in lune bin
* derived from 100 micron background
character*1 CPSCZ (nr__) ! *[^&=-+] Source category
integer*4 MHCONB (nr__) ! * 2 or >=3 hours-confirmed coverages
* for this lune bin
character*1 CQDOT (nr__) ! [q] q = in QDOT 1-in-6 sample
integer*4 W20 (nr__) ! (km/s) *HI (21cm) 20% velocity width or flag
integer*4 e_W20 (nr__) ! (km/s) Error on W20
integer*4 W50 (nr__) ! (km/s) *HI (21cm) 50% velocity width or flag
integer*4 e_W50 (nr__) ! (km/s) Error on W50
integer*4 r_W20 (nr__) ! *HI reference
real*8 FP12 (nr__) ! (Jy) Addscan flux at 12um using point
* source template
real*8 FP25 (nr__) ! (Jy) Addscan flux at 25um using point
* source template
real*8 FP60 (nr__) ! (Jy) Addscan flux at 60um using point
* source template
real*8 FP100 (nr__) ! (Jy) Addscan flux at 100um using point
* source template
integer*4 e_FP12 (nr__) ! (mJy) Statistical error in FP12
integer*4 e_FP25 (nr__) ! (mJy) Statistical error in FP25
integer*4 e_FP60 (nr__) ! (mJy) Statistical error in FP60
integer*4 e_FP100 (nr__) ! (mJy) Statistical error in FP100
integer*4 BP12 (nr__) ! Background fit subtracted or not at 12um
integer*4 BP25 (nr__) ! Background fit subtracted or not at 25um
integer*4 BP60 (nr__) ! Background fit subtracted or not at 60um
integer*4 BP100 (nr__) ! Background fit subtracted or not at 100um
real*8 FE12 (nr__) ! (Jy) Coadded or extended addscan flux at 12um
real*8 FE25 (nr__) ! (Jy) Coadded or extended addscan flux at 25um
real*8 FE60 (nr__) ! (Jy) Coadded or extended addscan flux at 60um
real*8 FE100 (nr__) ! (Jy) Coadded or extended addscan flux at 100um
integer*4 e_FE12 (nr__) ! (mJy) Statistical error in FE12
integer*4 e_FE25 (nr__) ! (mJy) Statistical error in FE25
integer*4 e_FE60 (nr__) ! (mJy) Statistical error in FE60
integer*4 e_FE100 (nr__) ! (mJy) Statistical error in FE100
integer*4 BE12 (nr__) ! *[0-2] Background fit at 12un
integer*4 BE25 (nr__) ! *[0-2] Background fit at 25um
integer*4 BE60 (nr__) ! *[0-2] Background fit at 60um
integer*4 BE100 (nr__) ! *[0-2] Background fit at 100um
real*8 WIDTH2 (nr__) ! (arcmin+2) Addscan width squared at 60um
real*8 e_WIDTH2 (nr__) ! (arcmin+2) Error on WIDTH2
*Note on CPSCZ: source category:
* = PSC (or PSC reject file) source, and not a large-diameter optical galaxy.
* + addscanned large optical galaxy added to PSCz.
* - PSC source associated with a "+" addscan object (normally excluded).
* & high latitude PSCz source, but in a bin excluded by psczmask.dat .
* ^ high latitude optical galaxy, in a bin excluded by psczmask.dat .
*Note on CZSTAT: redshift status:
* 1 PSCz redshift, excellent quality
* 2 PSCz redshift, good quality
* 3 PSCz redshift, low s/n but secure
* 4 PSCz redshift, marginal detection, unreliable
* l literature redshift, see r_Hvel for details.
* p private communication, see r_Hvel for details.
* q redshift from QDOT (J/MNRAS/308/897) or QCCOD
* # redshift wanted for PSCz, but still unknown
* x Give up, too optically faint (B or b_J > 19.5 mag)
* - Don't need it: galactic source, duplicate, part of larger galaxy etc.
*Note on CNAME: IRAS Names.
* Digits hhmmm+ddmm as in PSC. The initial letter means:
* Q: from QIGC (Cat. J/MNRAS/253/485)
* P: from PSC (Cat. II/215)
* R: from PSC reject file
* O: from optical catalogue
*Note on q_S12, q_S25, q_S60, q_S100:
* 1 = upper limit, 2 = moderate, 3 = OK, 4 = new HCON merged flux,
* 5 = addscanned flux/1.1, 6 = coadded flux,
* 7 = addscanned flux with point source template, 9 = single HCON.
*Note on CES12, CES25, CES60, CES100:
* 1 = flagged as confirmed extended source,
* 2 = in SSS (Small Scale Structure) catalog (Cat. <VII/73>)
* 3 = in Large Galaxy Catalog,
* 4 = addscan flux,
* 5 = addscan flux, source is double,
* 6 = coadd flux,
* 7 = coadd flux, source is double,
* 8 = addscanned flux, Yahil software used,
* 9 = addscanned point source flux, Yahil software used.
* 2-8 data from QIGC, there should be none in here except class='ds' sources.
*Note on Cirr1, Cirr2, Cirr3:
* The cirrus flag measures the number of sources detected only at
* 100{mu}m in a 1-{deg}^2^ area around the source
*Note on FCOR60: FCOR60 = NINT(FCOR(3) *1000)
* FCOR = correction factors applied to correct for overestimated faint
* fluxes as in Explanatory Supp. For the supplementary single-HCON
* sources, we have applied a factor of 1.121, as used in the PSC for 2
* out of possible 4 HCON detections. For merged single-HCON sources,
* unmerged in PSC due to software errors, we have corrected the
* correction factors accordingly.
*Note on IOPT: Source of Schmidt plate material.
* 1 = POSS-O/APM, 2 = SRC-bJ/APM,
* 3,4,5 = SRC-bJ/Cosmos, 6 = SR/Cosmos,
* 7 = POSS-E/APM, 9 = RC3 diameters .
* The optical information from plate scans is for the candidate with the
* best likelihood, not necessarily that for which the redshift is taken.
* The probability that these are the same depends strongly on latitude;
* at high latitude they are very reliable.
*
* The optical material is heterogeneous and except for IOPT=2 is
* uncalibrated. IOPT=6 magnitudes are estimated red mags, otherwise B or
* b_J. Due to saturation, magnitudes brighter than about 16 mag are very
* unreliable. Where a correct bright galaxy has been misclassified by
* the source-finding software as a star, the magnitude will be
* meaningless.
*Note on D1X, AR:
* CATX is the all-sky optical galaxy compilation based on the CGCG, UGC,
* MGC and the ESO catalogues
*Note on r_Hvel, r_Hvel2: Velocity references
* (negative reference numbers are PSCz or private communication)
* -99 INT/FOS or AAT/FORS
* -98 NOT
* -97 CTIO 60"
* -96 Fairall for PSCz
* -95 Nancay measurements for PSCz
* -94 AAT service data
* -93 INT/IDS500
* -92 Cananea
* -91 Parkes
* -90 WHT/FOSII
* -32 Juraszek
* -31 Visnanathan
* -30 Kraan-Kortweg
* -29 Lu, Pantoja
* -28 Nancay
* -27 Chamaraux, Yamada, Takata, Nakanishi
* -25 WHT Ultraluminous
* -24 QCCOD pole extension
* -23 QCCOD
* -21 Giovanelli + Haynes
* -20 DaCosta
* -10 - 70 as in Huchra 1991, See Cat. <VII/193>
* 83 QDOT, Lawrence et al., Cat. <J/MNRAS/308/897>
* 84 Leech et al., 1988MNRAS.231..977L
* 85 NED, see http://nedwww.ipac.caltech.edu/
* 86 PGC (Cat. <VII/119>), LEDA
* 87 1.2Jy, Cat. <VII/185>
* 88 Allen et al., 1990ApJ...355..130A
* 89 Seeberger
* 90 RC3 HI, Cat. <VII/155>
* 91 RC3 optical, Cat. <VII/155>
* 92 Fairall
* 97 Tully 1987ApJ...321..280T
* 98 Lonsdale and Hacking 1989ApJ...339..712L,
* Lonsdale 1986 (see Eales et al., 1988ApJ...328..530E)
*Note on ALI100, ALI100B, ABB:
* 100um backgrounds from Rowan-Robinson et al., 1991MNRAS.249..729R
* ALI4 is the background in a 0.5 degree beam centered on the source,
* corrected for that point source only.
* ALI4B is the background averaged over lune bin, and is corrected for
* all IGAL sources.
* ABB is based on ALI4B, and has been crudely corrected for
* temperature variation of dust across the galaxy.
*Note on MHCONB:
* MHCONB is the number of HCON coverages, crudely estimated from
* Explanatory Supp Fig I.C.1, for comparing source counts, fluxes etc.
*Note on W20, W50:
* The HI data is extremely heterogeneous in quality, reduction and
* linewidth definition.
* -999 = private,
* -997 = not reduced,
* -996 = to be observed,
* -995 = s/n < 5
* -994 = detected OK but not measured,
* -1 = non-detection,
* -2 = confused,
*Note on r_W20: HI references:
* -5 Parkes (Cat. <VIII/15>)
* -4 Nancay for PSCz
* -3 Lu
* -2 Bottinelli
* -1 Giovanelli and Haynes
* 1 LEDA
* 2 Pantoja 1997, Cat. <J/AJ/113/905>
* 3 Seeberger 1994, Cat. <J/A+A/286/17>
*Note on BE12, BE25, BE60, BE100:
* 0=no 1=yes 2=coadded flux
*Note on CLASS: object classification from optical data etc.
* g : galaxy or galaxies
* gf: faint galaxy, B or b_j > 19.5 mag.
* go: optical galaxy with 1D addscanned flux
* gr: Rice et al galaxy with 2D coadded photometry
* ss: Star, confirmed spectroscopically (by us, 1.2 Jansky or de Grijp)
* bs: Bright star from plate or telescope TV, no spectrum
* is: Star from IRAS colours and plate/TV appearance
* c : Cirrus
* pn: Planetary Nebula
* rn: Reflection Nebula
* h2: HII region
* zz: Other galactic source
* lm: in Large Magellanic Cloud
* sm: in Small Magellanic Cloud
* vb: Visually blank on sky survey plates, no non-stellar id
* kb: K' image blank, no non-stellar id
* rb: Radio blank, no 6 or 20cm id
* de: delete from catalogue, part of extended galaxy
* d1: delete from catalogue, broken up 1HCON source already listed elsewhere
* du: deleted optical galaxy, part of multiple system unresolved by IRAS
* dr: deleted optical galaxy, already in as Rice etal galaxy
* ao: PSC source associated with addscanned optical galaxy
* ar: PSC source associated with Rice etal optical galaxy
*
* The divisions above between Galactic classifications 'pn','rn' etc are
* not rigorous.
*
* The CLASS='go' and 'gr' sources were optically selected (see
* accompanying paper) and their IRAS fluxes determined by addscanning or
* coadding. Their fluxes (F12 to F100) are addscan or coadded fluxes,
* reduced by a factor 1.1 to bring them statistically into line with PSC
* fluxes. Where the resulting flux F(3) is 0.6 Jy or greater, they are
* added to the catalogue and flagged with CPSCZ='+'. There are 1467 such
* sources representing 1402 individual galaxies resolved by IRAS. The
* associated PSC source is retained in the catalogue, but is flagged as
* CPSCZ='-' and will normally be excluded (1290 sources).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'psczvs.dat' ! Short version of main catalogue
integer*4 nr__1
parameter (nr__1=18411) ! Number of records
character*100 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*9 CPNAME (nr__1) ! PSCz name eg A/123/456
character*1 CPSCZ_1 (nr__1) ! Source category
character*1 CZSTAT_1 (nr__1) ! Redshift status
character*2 CLASS_1 (nr__1) ! Optical classification code
integer*4 RAh_1 (nr__1) ! (h) Rigth ascension B1950
integer*4 RAm_1 (nr__1) ! (min) Rigth ascension B1950
integer*4 RAds_1 (nr__1) ! (0.1s) Rigth ascension B1950
character*1 DE__1 (nr__1) ! Declination sign B1950
integer*4 DEd_1 (nr__1) ! (deg) Declination B1950
integer*4 DEm_1 (nr__1) ! (arcmin) Declination B1950
integer*4 DEs_1 (nr__1) ! (arcsec) Declination B1950
real*4 GLON_1 (nr__1) ! (deg) Galactic longitude
real*4 GLAT_1 (nr__1) ! (deg) Galactic latitude
real*8 F60 (nr__1) ! (Jy) Flux at 60um
integer*4 Vhel (nr__1) ! (km/s) Best available heliocentric velocity (cz)
integer*4 e_Vhel (nr__1) ! (km/s) Error on heliocentric velocity
integer*4 r_Vhel (nr__1) ! Reference for heliocentric velocity
character*12 Oname_1 (nr__1) ! Optical name
integer*4 I100B (nr__1) ! (MJy/sr) 100 micron background for that lune bin
real*4 ABB_1 (nr__1) ! (mag) Estimated average B extinction for that
* lune bin
character*1 CHSD (nr__1) ! High Source Density flag, hex encoded
character*1 MHCONB_1 (nr__1) ! 2 or >=3 hours-confirmed coverages for
* this lune bin
real*4 F100_F60 (nr__1) ! Ratio of 100micron/60micron flux,
* (negative means upper limit at 100um)
integer*4 NLBIN_1 (nr__1) ! Lune bin
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'psczcg.dat' ! Sources close to coverage gap
integer*4 nr__2
parameter (nr__2=60) ! Number of records
character*533 ar__2 ! Full-size record
C B1950 position composed of: RAdeg DEdeg
character*11 CNAME_1 (nr__2) ! *Name e.g. Q12345+4856
integer*4 RAh_2 (nr__2) ! (h) Right ascension B1950
integer*4 RAm_2 (nr__2) ! (min) Right ascension B1950
integer*4 RAds_2 (nr__2) ! (0.1s) Right ascension B1950
character*1 DE__2 (nr__2) ! Declination sign B1950
integer*4 DEd_2 (nr__2) ! (deg) Declination B1950
integer*4 DEm_2 (nr__2) ! (arcmin) Declination B1950
integer*4 DEs_2 (nr__2) ! (arcsec) Declination B1950
real*8 S12_1 (nr__2) ! (Jy) Flux at 12um
real*8 S25_1 (nr__2) ! (Jy) Flux at 25um
real*8 S60_1 (nr__2) ! (Jy) Flux at 60um
real*8 S100_1 (nr__2) ! (Jy) Flux at 100um
integer*4 q_S12_1 (nr__2) ! *[1,9] Flux origin and quality flag at 12um
integer*4 q_S25_1 (nr__2) ! *[1,9] Flux origin and quality flag at 25um
integer*4 q_S60_1 (nr__2) ! *[1,9] Flux origin and quality flag at 60um
integer*4 q_S100_1 (nr__2) ! *[1,9] Flux origin and quality flag at 100um
integer*4 SES12_1 (nr__2) ! [0,8] Small Extended Source flag at 12um
integer*4 SES25_1 (nr__2) ! [0,8] Small Extended Source flag at 25um
integer*4 SES60_1 (nr__2) ! [0,8] Small Extended Source flag at 60um
integer*4 SES100_1 (nr__2) ! [0,8] Small Extended Source flag at 100um
integer*4 e_S12_1 (nr__2) ! (%) Flux uncertainty at 12um
integer*4 e_S25_1 (nr__2) ! (%) Flux uncertainty at 25um
integer*4 e_S60_1 (nr__2) ! (%) Flux uncertainty at 60um
integer*4 e_S100_1 (nr__2) ! (%) Flux uncertainty at 100um
integer*4 PosMin_1 (nr__2) ! (arcsec) Minor axis of 2{sigma} error ellipse
integer*4 PosMax_1 (nr__2) ! (arcsec) Major axis of 2{sigma} error ellipse
integer*4 PosPA_1 (nr__2) ! (deg) Position angle E of N of major axis
* of error ellipse
integer*4 Cirr1_1 (nr__2) ! *Cirrus flag
integer*4 Cirr2_1 (nr__2) ! *Cirrus flag
integer*4 Cirr3_1 (nr__2) ! (MJy/sr) *Cirrus flag
character*1 CC1_1 (nr__2) ! [A-N X] Correlation Coefficient at 12um
character*1 CC2_1 (nr__2) ! [A-N X] Correlation Coefficient at 25um
character*1 CC3_1 (nr__2) ! [A-M X1] Correlation Coefficient at 60um
character*1 CC4_1 (nr__2) ! [A-N X1] Correlation Coefficient at 100um
character*1 CONF_1 (nr__2) ! [0-9 ACDEF] Confusion flag, hex encoded
character*1 DISC_1 (nr__2) ! [0-9 A-D] Discrepant flux flag, hex encoded
character*1 HSD_1 (nr__2) ! [0-8] High Source Density flag, hex encoded
integer*4 PNEARH_1 (nr__2) ! [0-9] Nearby hours-confirmed sources
integer*4 PNEARW_1 (nr__2) ! [0-5] Nearby weeks-confirmed sources
integer*4 CES12_1 (nr__2) ! *Nearby seconds confirmed extended sources
* at 12um
integer*4 CES25_1 (nr__2) ! *Nearby seconds confirmed extended sources
* at 25um
integer*4 CES60_1 (nr__2) ! *Nearby seconds confirmed extended sources
* at 60um
integer*4 CES100_1 (nr__2) ! *Nearby seconds confirmed extended sources
* at 100um
integer*4 SNR60_1 (nr__2) ! 10 x signal/noise of weakest detection
* at 60um
integer*4 NHCON_1 (nr__2) ! Number of hours-confirmed source detections
integer*4 MHCON_1 (nr__2) ! Potential number of hours-confirmed source
* detections
integer*4 FCOR60_1 (nr__2) ! *1000*flux correction factor applied at 60um
character*5 PLATE_1 (nr__2) ! Plate number (POSS or SRC)
character*7 NOM_1 (nr__2) ! Nominal plate centre in RA,dec
integer*4 NUM_1 (nr__2) ! Identifier for source on this plate
integer*4 RAo_h_1 (nr__2) ! (h) ? Optical right ascension B1950
integer*4 RAo_m_1 (nr__2) ! (min) ? Optical right ascension B1950
real*4 RAo_s_1 (nr__2) ! (s) ? Optical right ascension B1950
character*1 DEo___1 (nr__2) ! ? Optical declination sign B1950
integer*4 DEo_d_1 (nr__2) ! (deg) ? Optical declination B1950
integer*4 DEo_m_1 (nr__2) ! (arcmin) ? Optical declination B1950
real*4 DEo_s_1 (nr__2) ! (arcsec) ? Optical declination B1950
integer*4 ADX_1 (nr__2) ! (arcsec) ? North/South Delta Optical/IRAS position
integer*4 ADY_1 (nr__2) ! (arcsec) ? East/West Delta Optical/IRAS position
real*4 ALK_1 (nr__2) ! ? log_10 likelihood for optical
* identification
integer*4 REL_1 (nr__2) ! (%) ? Prob. correct optical identification
character*1 CAND_1 (nr__2) ! ? Optical candidate A,B etc
integer*4 IOPT_1 (nr__2) ! *? Schmidt plate material used
real*4 Bmag_1 (nr__2) ! (mag) ? Estimated optical B magnitude
integer*4 MajAxiO_1 (nr__2) ! (arcsec) ? Optical major axis
integer*4 MinAxiO_1 (nr__2) ! (arcsec) ? Optical minor axis
integer*4 PAO_1 (nr__2) ! (deg) ? Optical major axis position angle E of N
real*8 SuBr_1 (nr__2) ! (mag/arcsec+2) ? Average optical surface brightness
integer*4 IDTYPE_1 (nr__2) ! Cross-identification type from PSC
character*12 Oname_2 (nr__2) ! Optical name
real*4 mag_1 (nr__2) ! (mag) Zwicky/RC3 magnitude
real*4 D1X_1 (nr__2) ! (arcmin) *RC3/UGC/MCG/ESO D25 diameter from CATX
integer*4 AR_1 (nr__2) ! (%) *Axial ratio d25/D25 from CATX
integer*4 Hvel_1 (nr__2) ! (km/s) Best available heliocentric velocity (cz)
integer*4 e_Hvel_1 (nr__2) ! (km/s) Error on heliocentric velocity
integer*4 r_Hvel_1 (nr__2) ! *Reference for heliocentric velocity
integer*4 Hvel2_1 (nr__2) ! (km/s) Heliocentric velocity (cz)
* measured as part of PSCz survey
integer*4 e_Hvel2_1 (nr__2) ! (km/s) Error on PSCz heliocentric velocity
integer*4 r_Hvel2_1 (nr__2) ! *Reference for PSCz heliocentric velocity
integer*4 TType_1 (nr__2) ! De Vaucouleurs Type, 99 = unknown
character*8 PGC_1 (nr__2) ! Catalogue of Principal Galaxies
* (PGC, Cat. <VII/119>) name
real*8 RAdeg_2 (nr__2) ! (deg) Right ascension B1950 (IRAS source)
real*8 DEdeg_2 (nr__2) ! (deg) Declination B1950 (IRAS source)
real*8 GLON_2 (nr__2) ! (deg) Galactic longitude
real*8 GLAT_2 (nr__2) ! (deg) Galactic latitude
integer*4 cz_1 (nr__2) ! (km/s) Velocity cz in Local Group frame
real*4 Dist_1 (nr__2) ! (Mpc) Distance from Rowan-Robinson 1988
integer*4 IDclus_1 (nr__2) ! Cluster identification from
* Rowan-Robinson 1988
real*4 ALI100_1 (nr__2) ! ([MJy/sr]) *log10(100um background in MJy/sr)
real*4 ALI100B_1 (nr__2) ! ([MJy/sr]) *log10(100um background in MJy/sr
* for lune bin)
integer*4 NLBIN_2 (nr__2) ! Lune bin number
character*2 CLASS_2 (nr__2) ! *Optical classification code
character*5 CSEC_1 (nr__2) ! PSCz sector, e.g. A/123
integer*4 Xpos_1 (nr__2) ! (mm) X position from nominal bottom left corner
* of Schmidt plate
integer*4 Ypos_1 (nr__2) ! (mm) Y position from nominal bottom left corner
* of Schmidt plate
character*1 CSTAT_1 (nr__2) ! [JORXPCELS?0] Optical material used
character*1 CAREA_1 (nr__2) ! h=high|b|, p=other PSCz,
* l=BTP, c=in coverage gaps
character*1 CZSTAT_2 (nr__2) ! *[#-1234lpxq] Redshift status
character*1 C12CLAS_1 (nr__2) ! Classification in 1.2Jy survey
real*8 F3PSC_1 (nr__2) ! (Jy) PSC flux at 60um
real*8 F4PSC_1 (nr__2) ! (Jy) PSC flux at 100um
real*4 ABB_2 (nr__2) ! (mag) *Estimated average B extinction in lune bin
* derived from 100 micron background
character*1 CPSCZ_2 (nr__2) ! *[^&=-+] Source category
integer*4 MHCONB_2 (nr__2) ! * 2 or >=3 hours-confirmed coverages
* for this lune bin
character*1 CQDOT_1 (nr__2) ! [q] q = in QDOT 1-in-6 sample
integer*4 W20_1 (nr__2) ! (km/s) *HI (21cm) 20% velocity width or flag
integer*4 e_W20_1 (nr__2) ! (km/s) Error on W20
integer*4 W50_1 (nr__2) ! (km/s) *HI (21cm) 50% velocity width or flag
integer*4 e_W50_1 (nr__2) ! (km/s) Error on W50
integer*4 r_W20_1 (nr__2) ! *HI reference
real*8 FP12_1 (nr__2) ! (Jy) Addscan flux at 12um using point
* source template
real*8 FP25_1 (nr__2) ! (Jy) Addscan flux at 25um using point
* source template
real*8 FP60_1 (nr__2) ! (Jy) Addscan flux at 60um using point
* source template
real*8 FP100_1 (nr__2) ! (Jy) Addscan flux at 100um using point
* source template
integer*4 e_FP12_1 (nr__2) ! (mJy) Statistical error in FP12
integer*4 e_FP25_1 (nr__2) ! (mJy) Statistical error in FP25
integer*4 e_FP60_1 (nr__2) ! (mJy) Statistical error in FP60
integer*4 e_FP100_1 (nr__2) ! (mJy) Statistical error in FP100
integer*4 BP12_1 (nr__2) ! Background fit subtracted or not at 12um
integer*4 BP25_1 (nr__2) ! Background fit subtracted or not at 25um
integer*4 BP60_1 (nr__2) ! Background fit subtracted or not at 60um
integer*4 BP100_1 (nr__2) ! Background fit subtracted or not at 100um
real*8 FE12_1 (nr__2) ! (Jy) Coadded or extended addscan flux at 12um
real*8 FE25_1 (nr__2) ! (Jy) Coadded or extended addscan flux at 25um
real*8 FE60_1 (nr__2) ! (Jy) Coadded or extended addscan flux at 60um
real*8 FE100_1 (nr__2) ! (Jy) Coadded or extended addscan flux at 100um
integer*4 e_FE12_1 (nr__2) ! (mJy) Statistical error in FE12
integer*4 e_FE25_1 (nr__2) ! (mJy) Statistical error in FE25
integer*4 e_FE60_1 (nr__2) ! (mJy) Statistical error in FE60
integer*4 e_FE100_1 (nr__2) ! (mJy) Statistical error in FE100
integer*4 BE12_1 (nr__2) ! *[0-2] Background fit at 12un
integer*4 BE25_1 (nr__2) ! *[0-2] Background fit at 25um
integer*4 BE60_1 (nr__2) ! *[0-2] Background fit at 60um
integer*4 BE100_1 (nr__2) ! *[0-2] Background fit at 100um
real*8 WIDTH2_1 (nr__2) ! (arcmin+2) Addscan width squared at 60um
real*8 e_WIDTH2_1 (nr__2) ! (arcmin+2) Error on WIDTH2
*Note on CPSCZ: source category:
* = PSC (or PSC reject file) source, and not a large-diameter optical galaxy.
* + addscanned large optical galaxy added to PSCz.
* - PSC source associated with a "+" addscan object (normally excluded).
* & high latitude PSCz source, but in a bin excluded by psczmask.dat .
* ^ high latitude optical galaxy, in a bin excluded by psczmask.dat .
*Note on CZSTAT: redshift status:
* 1 PSCz redshift, excellent quality
* 2 PSCz redshift, good quality
* 3 PSCz redshift, low s/n but secure
* 4 PSCz redshift, marginal detection, unreliable
* l literature redshift, see r_Hvel for details.
* p private communication, see r_Hvel for details.
* q redshift from QDOT (J/MNRAS/308/897) or QCCOD
* # redshift wanted for PSCz, but still unknown
* x Give up, too optically faint (B or b_J > 19.5 mag)
* - Don't need it: galactic source, duplicate, part of larger galaxy etc.
*Note on CNAME: IRAS Names.
* Digits hhmmm+ddmm as in PSC. The initial letter means:
* Q: from QIGC (Cat. J/MNRAS/253/485)
* P: from PSC (Cat. II/215)
* R: from PSC reject file
* O: from optical catalogue
*Note on q_S12, q_S25, q_S60, q_S100:
* 1 = upper limit, 2 = moderate, 3 = OK, 4 = new HCON merged flux,
* 5 = addscanned flux/1.1, 6 = coadded flux,
* 7 = addscanned flux with point source template, 9 = single HCON.
*Note on CES12, CES25, CES60, CES100:
* 1 = flagged as confirmed extended source,
* 2 = in SSS (Small Scale Structure) catalog (Cat. <VII/73>)
* 3 = in Large Galaxy Catalog,
* 4 = addscan flux,
* 5 = addscan flux, source is double,
* 6 = coadd flux,
* 7 = coadd flux, source is double,
* 8 = addscanned flux, Yahil software used,
* 9 = addscanned point source flux, Yahil software used.
* 2-8 data from QIGC, there should be none in here except class='ds' sources.
*Note on Cirr1, Cirr2, Cirr3:
* The cirrus flag measures the number of sources detected only at
* 100{mu}m in a 1-{deg}^2^ area around the source
*Note on FCOR60: FCOR60 = NINT(FCOR(3) *1000)
* FCOR = correction factors applied to correct for overestimated faint
* fluxes as in Explanatory Supp. For the supplementary single-HCON
* sources, we have applied a factor of 1.121, as used in the PSC for 2
* out of possible 4 HCON detections. For merged single-HCON sources,
* unmerged in PSC due to software errors, we have corrected the
* correction factors accordingly.
*Note on IOPT: Source of Schmidt plate material.
* 1 = POSS-O/APM, 2 = SRC-bJ/APM,
* 3,4,5 = SRC-bJ/Cosmos, 6 = SR/Cosmos,
* 7 = POSS-E/APM, 9 = RC3 diameters .
* The optical information from plate scans is for the candidate with the
* best likelihood, not necessarily that for which the redshift is taken.
* The probability that these are the same depends strongly on latitude;
* at high latitude they are very reliable.
*
* The optical material is heterogeneous and except for IOPT=2 is
* uncalibrated. IOPT=6 magnitudes are estimated red mags, otherwise B or
* b_J. Due to saturation, magnitudes brighter than about 16 mag are very
* unreliable. Where a correct bright galaxy has been misclassified by
* the source-finding software as a star, the magnitude will be
* meaningless.
*Note on D1X, AR:
* CATX is the all-sky optical galaxy compilation based on the CGCG, UGC,
* MGC and the ESO catalogues
*Note on r_Hvel, r_Hvel2: Velocity references
* (negative reference numbers are PSCz or private communication)
* -99 INT/FOS or AAT/FORS
* -98 NOT
* -97 CTIO 60"
* -96 Fairall for PSCz
* -95 Nancay measurements for PSCz
* -94 AAT service data
* -93 INT/IDS500
* -92 Cananea
* -91 Parkes
* -90 WHT/FOSII
* -32 Juraszek
* -31 Visnanathan
* -30 Kraan-Kortweg
* -29 Lu, Pantoja
* -28 Nancay
* -27 Chamaraux, Yamada, Takata, Nakanishi
* -25 WHT Ultraluminous
* -24 QCCOD pole extension
* -23 QCCOD
* -21 Giovanelli + Haynes
* -20 DaCosta
* -10 - 70 as in Huchra 1991, See Cat. <VII/193>
* 83 QDOT, Lawrence et al., Cat. <J/MNRAS/308/897>
* 84 Leech et al., 1988MNRAS.231..977L
* 85 NED, see http://nedwww.ipac.caltech.edu/
* 86 PGC (Cat. <VII/119>), LEDA
* 87 1.2Jy, Cat. <VII/185>
* 88 Allen et al., 1990ApJ...355..130A
* 89 Seeberger
* 90 RC3 HI, Cat. <VII/155>
* 91 RC3 optical, Cat. <VII/155>
* 92 Fairall
* 97 Tully 1987ApJ...321..280T
* 98 Lonsdale and Hacking 1989ApJ...339..712L,
* Lonsdale 1986 (see Eales et al., 1988ApJ...328..530E)
*Note on ALI100, ALI100B, ABB:
* 100um backgrounds from Rowan-Robinson et al., 1991MNRAS.249..729R
* ALI4 is the background in a 0.5 degree beam centered on the source,
* corrected for that point source only.
* ALI4B is the background averaged over lune bin, and is corrected for
* all IGAL sources.
* ABB is based on ALI4B, and has been crudely corrected for
* temperature variation of dust across the galaxy.
*Note on MHCONB:
* MHCONB is the number of HCON coverages, crudely estimated from
* Explanatory Supp Fig I.C.1, for comparing source counts, fluxes etc.
*Note on W20, W50:
* The HI data is extremely heterogeneous in quality, reduction and
* linewidth definition.
* -999 = private,
* -997 = not reduced,
* -996 = to be observed,
* -995 = s/n < 5
* -994 = detected OK but not measured,
* -1 = non-detection,
* -2 = confused,
*Note on r_W20: HI references:
* -5 Parkes (Cat. <VIII/15>)
* -4 Nancay for PSCz
* -3 Lu
* -2 Bottinelli
* -1 Giovanelli and Haynes
* 1 LEDA
* 2 Pantoja 1997, Cat. <J/AJ/113/905>
* 3 Seeberger 1994, Cat. <J/A+A/286/17>
*Note on BE12, BE25, BE60, BE100:
* 0=no 1=yes 2=coadded flux
*Note on CLASS: object classification from optical data etc.
* g : galaxy or galaxies
* gf: faint galaxy, B or b_j > 19.5 mag.
* go: optical galaxy with 1D addscanned flux
* gr: Rice et al galaxy with 2D coadded photometry
* ss: Star, confirmed spectroscopically (by us, 1.2 Jansky or de Grijp)
* bs: Bright star from plate or telescope TV, no spectrum
* is: Star from IRAS colours and plate/TV appearance
* c : Cirrus
* pn: Planetary Nebula
* rn: Reflection Nebula
* h2: HII region
* zz: Other galactic source
* lm: in Large Magellanic Cloud
* sm: in Small Magellanic Cloud
* vb: Visually blank on sky survey plates, no non-stellar id
* kb: K' image blank, no non-stellar id
* rb: Radio blank, no 6 or 20cm id
* de: delete from catalogue, part of extended galaxy
* d1: delete from catalogue, broken up 1HCON source already listed elsewhere
* du: deleted optical galaxy, part of multiple system unresolved by IRAS
* dr: deleted optical galaxy, already in as Rice etal galaxy
* ao: PSC source associated with addscanned optical galaxy
* ar: PSC source associated with Rice etal optical galaxy
*
* The divisions above between Galactic classifications 'pn','rn' etc are
* not rigorous.
*
* The CLASS='go' and 'gr' sources were optically selected (see
* accompanying paper) and their IRAS fluxes determined by addscanning or
* coadding. Their fluxes (F12 to F100) are addscan or coadded fluxes,
* reduced by a factor 1.1 to bring them statistically into line with PSC
* fluxes. Where the resulting flux F(3) is 0.6 Jy or greater, they are
* added to the catalogue and flagged with CPSCZ='+'. There are 1467 such
* sources representing 1402 individual galaxies resolved by IRAS. The
* associated PSC source is retained in the catalogue, but is flagged as
* CPSCZ='-' and will normally be excluded (1290 sources).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'psczcgvs.dat' ! Short version of psczcg.dat
integer*4 nr__3
parameter (nr__3=60) ! Number of records
character*100 ar__3 ! Full-size record
C B1950 position composed of: RAh RAm RAds DE- DEd DEm DEs
real*8 RAdeg_3 (nr__3) ! (deg) Right Ascension B1950
real*8 DEdeg_3 (nr__3) ! (deg) Declination B1950
C ---------------------------------- ! (position vector(s) in degrees)
character*9 CPNAME_1 (nr__3) ! PSCz name eg A/123/456
character*1 CPSCZ_3 (nr__3) ! Source category
character*1 CZSTAT_3 (nr__3) ! Redshift status
character*2 CLASS_3 (nr__3) ! Optical classification code
integer*4 RAh_3 (nr__3) ! (h) Rigth ascension B1950
integer*4 RAm_3 (nr__3) ! (min) Rigth ascension B1950
integer*4 RAds_3 (nr__3) ! (0.1s) Rigth ascension B1950
character*1 DE__3 (nr__3) ! Declination sign B1950
integer*4 DEd_3 (nr__3) ! (deg) Declination B1950
integer*4 DEm_3 (nr__3) ! (arcmin) Declination B1950
integer*4 DEs_3 (nr__3) ! (arcsec) Declination B1950
real*4 GLON_3 (nr__3) ! (deg) Galactic longitude
real*4 GLAT_3 (nr__3) ! (deg) Galactic latitude
real*8 F60_1 (nr__3) ! (Jy) Flux at 60um
integer*4 Vhel_1 (nr__3) ! (km/s) Best available heliocentric velocity (cz)
integer*4 e_Vhel_1 (nr__3) ! (km/s) Error on heliocentric velocity
integer*4 r_Vhel_1 (nr__3) ! Reference for heliocentric velocity
character*12 Oname_3 (nr__3) ! Optical name
integer*4 I100B_1 (nr__3) ! (MJy/sr) 100 micron background for that lune bin
real*4 ABB_3 (nr__3) ! (mag) Estimated average B extinction for that
* lune bin
character*1 CHSD_1 (nr__3) ! High Source Density flag, hex encoded
character*1 MHCONB_3 (nr__3) ! 2 or >=3 hours-confirmed coverages for
* this lune bin
real*4 F100_F60_1 (nr__3) ! Ratio of 100micron/60micron flux,
* (negative means upper limit at 100um)
integer*4 NLBIN_3 (nr__3) ! Lune bin
C=============================================================================
C Loading file 'pscz.dat' ! Main catalogue
C Format for file interpretation
1 format(
+ 1X,A11,1X,I2,I2,I3,1X,A1,I2,I2,I2,F8.2,F8.2,F8.2,F8.2,1X,I1,
+ I1,I1,I1,1X,I1,I1,I1,I1,1X,I3,I3,I3,I3,1X,I3,I3,I3,1X,I1,I1,
+ I3,1X,A1,A1,A1,A1,1X,A1,A1,A1,1X,I1,I1,1X,I1,I1,I1,I1,I5,I2,
+ I2,1X,I4,1X,A5,1X,A7,1X,I3,1X,I2,I2,F5.2,1X,A1,I2,I2,F4.1,1X,
+ I3,I3,F4.1,I3,A1,I1,F6.2,I3,I3,I3,F7.3,I2,A12,F6.2,F5.1,I3,I6,
+ I4,I3,I6,I4,I3,I2,1X,A8,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,I6,
+ F5.1,I4,1X,F5.2,F5.2,I6,1X,A2,1X,A5,I4,I4,1X,A1,A1,A1,1X,A1,
+ F8.2,F8.2,F5.2,1X,A1,I2,1X,A1,I4,I2,I4,I2,I2,F7.3,F9.3,F9.3,
+ F9.3,1X,I4,I4,I4,I4,1X,I1,I1,I1,I1,F9.3,F9.3,F9.3,F9.3,1X,I4,
+ I4,I4,I4,1X,I1,I1,I1,I1,F8.4,F8.4)
C Effective file loading
open(unit=1,status='old',file=
+'pscz.dat')
write(6,*) '....Loading file: pscz.dat'
do i__=1,18351
read(1,'(A533)')ar__
read(ar__,1)
+ CNAME(i__),RAh(i__),RAm(i__),RAds(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),S12(i__),S25(i__),S60(i__),S100(i__),
+ q_S12(i__),q_S25(i__),q_S60(i__),q_S100(i__),SES12(i__),
+ SES25(i__),SES60(i__),SES100(i__),e_S12(i__),e_S25(i__),
+ e_S60(i__),e_S100(i__),PosMin(i__),PosMax(i__),PosPA(i__),
+ Cirr1(i__),Cirr2(i__),Cirr3(i__),CC1(i__),CC2(i__),CC3(i__),
+ CC4(i__),CONF(i__),DISC(i__),HSD(i__),PNEARH(i__),PNEARW(i__),
+ CES12(i__),CES25(i__),CES60(i__),CES100(i__),SNR60(i__),
+ NHCON(i__),MHCON(i__),FCOR60(i__),PLATE(i__),NOM(i__),
+ NUM(i__),RAo_h(i__),RAo_m(i__),RAo_s(i__),DEo__(i__),
+ DEo_d(i__),DEo_m(i__),DEo_s(i__),ADX(i__),ADY(i__),ALK(i__),
+ REL(i__),CAND(i__),IOPT(i__),Bmag(i__),MajAxiO(i__),
+ MinAxiO(i__),PAO(i__),SuBr(i__),IDTYPE(i__),Oname(i__),
+ mag(i__),D1X(i__),AR(i__),Hvel(i__),e_Hvel(i__),r_Hvel(i__),
+ Hvel2(i__),e_Hvel2(i__),r_Hvel2(i__),TType(i__),PGC(i__),
+ RAdeg(i__),DEdeg(i__),GLON(i__),GLAT(i__),cz(i__),Dist(i__),
+ IDclus(i__),ALI100(i__),ALI100B(i__),NLBIN(i__),CLASS(i__),
+ CSEC(i__),Xpos(i__),Ypos(i__),CSTAT(i__),CAREA(i__),
+ CZSTAT(i__),C12CLAS(i__),F3PSC(i__),F4PSC(i__),ABB(i__),
+ CPSCZ(i__),MHCONB(i__),CQDOT(i__),W20(i__),e_W20(i__),
+ W50(i__),e_W50(i__),r_W20(i__),FP12(i__),FP25(i__),FP60(i__),
+ FP100(i__),e_FP12(i__),e_FP25(i__),e_FP60(i__),e_FP100(i__),
+ BP12(i__),BP25(i__),BP60(i__),BP100(i__),FE12(i__),FE25(i__),
+ FE60(i__),FE100(i__),e_FE12(i__),e_FE25(i__),e_FE60(i__),
+ e_FE100(i__),BE12(i__),BE25(i__),BE60(i__),BE100(i__),
+ WIDTH2(i__),e_WIDTH2(i__)
if(ar__(150:151) .EQ. '') RAo_h(i__) = iNULL__
if(ar__(152:153) .EQ. '') RAo_m(i__) = iNULL__
if(ar__(154:158) .EQ. '') RAo_s(i__) = rNULL__
if(ar__(161:162) .EQ. '') DEo_d(i__) = iNULL__
if(ar__(163:164) .EQ. '') DEo_m(i__) = iNULL__
if(ar__(165:168) .EQ. '') DEo_s(i__) = rNULL__
if(ar__(170:172) .EQ. '') ADX(i__) = iNULL__
if(ar__(173:175) .EQ. '') ADY(i__) = iNULL__
if(ar__(176:179) .EQ. '') ALK(i__) = rNULL__
if(ar__(180:182) .EQ. '') REL(i__) = iNULL__
if(ar__(184:184) .EQ. '') IOPT(i__) = iNULL__
if(ar__(185:190) .EQ. '') Bmag(i__) = rNULL__
if(ar__(191:193) .EQ. '') MajAxiO(i__) = iNULL__
if(ar__(194:196) .EQ. '') MinAxiO(i__) = iNULL__
if(ar__(197:199) .EQ. '') PAO(i__) = iNULL__
if(ar__(200:206) .EQ. '') SuBr(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ CNAME(i__),RAh(i__),RAm(i__),RAds(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),S12(i__),S25(i__),S60(i__),S100(i__),
+ q_S12(i__),q_S25(i__),q_S60(i__),q_S100(i__),SES12(i__),
+ SES25(i__),SES60(i__),SES100(i__),e_S12(i__),e_S25(i__),
+ e_S60(i__),e_S100(i__),PosMin(i__),PosMax(i__),PosPA(i__),
+ Cirr1(i__),Cirr2(i__),Cirr3(i__),CC1(i__),CC2(i__),CC3(i__),
+ CC4(i__),CONF(i__),DISC(i__),HSD(i__),PNEARH(i__),PNEARW(i__),
+ CES12(i__),CES25(i__),CES60(i__),CES100(i__),SNR60(i__),
+ NHCON(i__),MHCON(i__),FCOR60(i__),PLATE(i__),NOM(i__),
+ NUM(i__),RAo_h(i__),RAo_m(i__),RAo_s(i__),DEo__(i__),
+ DEo_d(i__),DEo_m(i__),DEo_s(i__),ADX(i__),ADY(i__),ALK(i__),
+ REL(i__),CAND(i__),IOPT(i__),Bmag(i__),MajAxiO(i__),
+ MinAxiO(i__),PAO(i__),SuBr(i__),IDTYPE(i__),Oname(i__),
+ mag(i__),D1X(i__),AR(i__),Hvel(i__),e_Hvel(i__),r_Hvel(i__),
+ Hvel2(i__),e_Hvel2(i__),r_Hvel2(i__),TType(i__),PGC(i__),
+ RAdeg(i__),DEdeg(i__),GLON(i__),GLAT(i__),cz(i__),Dist(i__),
+ IDclus(i__),ALI100(i__),ALI100B(i__),NLBIN(i__),CLASS(i__),
+ CSEC(i__),Xpos(i__),Ypos(i__),CSTAT(i__),CAREA(i__),
+ CZSTAT(i__),C12CLAS(i__),F3PSC(i__),F4PSC(i__),ABB(i__),
+ CPSCZ(i__),MHCONB(i__),CQDOT(i__),W20(i__),e_W20(i__),
+ W50(i__),e_W50(i__),r_W20(i__),FP12(i__),FP25(i__),FP60(i__),
+ FP100(i__),e_FP12(i__),e_FP25(i__),e_FP60(i__),e_FP100(i__),
+ BP12(i__),BP25(i__),BP60(i__),BP100(i__),FE12(i__),FE25(i__),
+ FE60(i__),FE100(i__),e_FE12(i__),e_FE25(i__),e_FE60(i__),
+ e_FE100(i__),BE12(i__),BE25(i__),BE60(i__),BE100(i__),
+ WIDTH2(i__),e_WIDTH2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'psczvs.dat' ! Short version of main catalogue
C Format for file interpretation
2 format(
+ A9,1X,A1,A1,A2,1X,I2,I2,I3,1X,A1,I2,I2,I2,F6.1,F6.1,F8.2,I6,
+ I4,I3,A12,I4,F5.2,1X,A1,1X,A1,F6.2,I6)
C Effective file loading
open(unit=1,status='old',file=
+'psczvs.dat')
write(6,*) '....Loading file: psczvs.dat'
do i__=1,18411
read(1,'(A100)')ar__1
read(ar__1,2)
+ CPNAME(i__),CPSCZ_1(i__),CZSTAT_1(i__),CLASS_1(i__),
+ RAh_1(i__),RAm_1(i__),RAds_1(i__),DE__1(i__),DEd_1(i__),
+ DEm_1(i__),DEs_1(i__),GLON_1(i__),GLAT_1(i__),F60(i__),
+ Vhel(i__),e_Vhel(i__),r_Vhel(i__),Oname_1(i__),I100B(i__),
+ ABB_1(i__),CHSD(i__),MHCONB_1(i__),F100_F60(i__),NLBIN_1(i__)
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_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAds_1(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)
+ CPNAME(i__),CPSCZ_1(i__),CZSTAT_1(i__),CLASS_1(i__),
+ RAh_1(i__),RAm_1(i__),RAds_1(i__),DE__1(i__),DEd_1(i__),
+ DEm_1(i__),DEs_1(i__),GLON_1(i__),GLAT_1(i__),F60(i__),
+ Vhel(i__),e_Vhel(i__),r_Vhel(i__),Oname_1(i__),I100B(i__),
+ ABB_1(i__),CHSD(i__),MHCONB_1(i__),F100_F60(i__),NLBIN_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 'psczcg.dat' ! Sources close to coverage gap
C Format for file interpretation
3 format(
+ 1X,A11,1X,I2,I2,I3,1X,A1,I2,I2,I2,F8.2,F8.2,F8.2,F8.2,1X,I1,
+ I1,I1,I1,1X,I1,I1,I1,I1,1X,I3,I3,I3,I3,1X,I3,I3,I3,1X,I1,I1,
+ I3,1X,A1,A1,A1,A1,1X,A1,A1,A1,1X,I1,I1,1X,I1,I1,I1,I1,I5,I2,
+ I2,1X,I4,1X,A5,1X,A7,1X,I3,1X,I2,I2,F5.2,1X,A1,I2,I2,F4.1,1X,
+ I3,I3,F4.1,I3,A1,I1,F6.2,I3,I3,I3,F7.3,I2,A12,F6.2,F5.1,I3,I6,
+ I4,I3,I6,I4,I3,I2,1X,A8,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,I6,
+ F5.1,I4,1X,F5.2,F5.2,I6,1X,A2,1X,A5,I4,I4,1X,A1,A1,A1,1X,A1,
+ F8.2,F8.2,F5.2,1X,A1,I2,1X,A1,I4,I2,I4,I2,I2,F7.3,F9.3,F9.3,
+ F9.3,1X,I4,I4,I4,I4,1X,I1,I1,I1,I1,F9.3,F9.3,F9.3,F9.3,1X,I4,
+ I4,I4,I4,1X,I1,I1,I1,I1,F8.4,F8.4)
C Effective file loading
open(unit=1,status='old',file=
+'psczcg.dat')
write(6,*) '....Loading file: psczcg.dat'
do i__=1,60
read(1,'(A533)')ar__2
read(ar__2,3)
+ CNAME_1(i__),RAh_2(i__),RAm_2(i__),RAds_2(i__),DE__2(i__),
+ DEd_2(i__),DEm_2(i__),DEs_2(i__),S12_1(i__),S25_1(i__),
+ S60_1(i__),S100_1(i__),q_S12_1(i__),q_S25_1(i__),q_S60_1(i__),
+ q_S100_1(i__),SES12_1(i__),SES25_1(i__),SES60_1(i__),
+ SES100_1(i__),e_S12_1(i__),e_S25_1(i__),e_S60_1(i__),
+ e_S100_1(i__),PosMin_1(i__),PosMax_1(i__),PosPA_1(i__),
+ Cirr1_1(i__),Cirr2_1(i__),Cirr3_1(i__),CC1_1(i__),CC2_1(i__),
+ CC3_1(i__),CC4_1(i__),CONF_1(i__),DISC_1(i__),HSD_1(i__),
+ PNEARH_1(i__),PNEARW_1(i__),CES12_1(i__),CES25_1(i__),
+ CES60_1(i__),CES100_1(i__),SNR60_1(i__),NHCON_1(i__),
+ MHCON_1(i__),FCOR60_1(i__),PLATE_1(i__),NOM_1(i__),NUM_1(i__),
+ RAo_h_1(i__),RAo_m_1(i__),RAo_s_1(i__),DEo___1(i__),
+ DEo_d_1(i__),DEo_m_1(i__),DEo_s_1(i__),ADX_1(i__),ADY_1(i__),
+ ALK_1(i__),REL_1(i__),CAND_1(i__),IOPT_1(i__),Bmag_1(i__),
+ MajAxiO_1(i__),MinAxiO_1(i__),PAO_1(i__),SuBr_1(i__),
+ IDTYPE_1(i__),Oname_2(i__),mag_1(i__),D1X_1(i__),AR_1(i__),
+ Hvel_1(i__),e_Hvel_1(i__),r_Hvel_1(i__),Hvel2_1(i__),
+ e_Hvel2_1(i__),r_Hvel2_1(i__),TType_1(i__),PGC_1(i__),
+ RAdeg_2(i__),DEdeg_2(i__),GLON_2(i__),GLAT_2(i__),cz_1(i__),
+ Dist_1(i__),IDclus_1(i__),ALI100_1(i__),ALI100B_1(i__),
+ NLBIN_2(i__),CLASS_2(i__),CSEC_1(i__),Xpos_1(i__),Ypos_1(i__),
+ CSTAT_1(i__),CAREA_1(i__),CZSTAT_2(i__),C12CLAS_1(i__),
+ F3PSC_1(i__),F4PSC_1(i__),ABB_2(i__),CPSCZ_2(i__),
+ MHCONB_2(i__),CQDOT_1(i__),W20_1(i__),e_W20_1(i__),W50_1(i__),
+ e_W50_1(i__),r_W20_1(i__),FP12_1(i__),FP25_1(i__),FP60_1(i__),
+ FP100_1(i__),e_FP12_1(i__),e_FP25_1(i__),e_FP60_1(i__),
+ e_FP100_1(i__),BP12_1(i__),BP25_1(i__),BP60_1(i__),
+ BP100_1(i__),FE12_1(i__),FE25_1(i__),FE60_1(i__),FE100_1(i__),
+ e_FE12_1(i__),e_FE25_1(i__),e_FE60_1(i__),e_FE100_1(i__),
+ BE12_1(i__),BE25_1(i__),BE60_1(i__),BE100_1(i__),
+ WIDTH2_1(i__),e_WIDTH2_1(i__)
if(ar__2(150:151) .EQ. '') RAo_h_1(i__) = iNULL__
if(ar__2(152:153) .EQ. '') RAo_m_1(i__) = iNULL__
if(ar__2(154:158) .EQ. '') RAo_s_1(i__) = rNULL__
if(ar__2(161:162) .EQ. '') DEo_d_1(i__) = iNULL__
if(ar__2(163:164) .EQ. '') DEo_m_1(i__) = iNULL__
if(ar__2(165:168) .EQ. '') DEo_s_1(i__) = rNULL__
if(ar__2(170:172) .EQ. '') ADX_1(i__) = iNULL__
if(ar__2(173:175) .EQ. '') ADY_1(i__) = iNULL__
if(ar__2(176:179) .EQ. '') ALK_1(i__) = rNULL__
if(ar__2(180:182) .EQ. '') REL_1(i__) = iNULL__
if(ar__2(184:184) .EQ. '') IOPT_1(i__) = iNULL__
if(ar__2(185:190) .EQ. '') Bmag_1(i__) = rNULL__
if(ar__2(191:193) .EQ. '') MajAxiO_1(i__) = iNULL__
if(ar__2(194:196) .EQ. '') MinAxiO_1(i__) = iNULL__
if(ar__2(197:199) .EQ. '') PAO_1(i__) = iNULL__
if(ar__2(200:206) .EQ. '') SuBr_1(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ CNAME_1(i__),RAh_2(i__),RAm_2(i__),RAds_2(i__),DE__2(i__),
+ DEd_2(i__),DEm_2(i__),DEs_2(i__),S12_1(i__),S25_1(i__),
+ S60_1(i__),S100_1(i__),q_S12_1(i__),q_S25_1(i__),q_S60_1(i__),
+ q_S100_1(i__),SES12_1(i__),SES25_1(i__),SES60_1(i__),
+ SES100_1(i__),e_S12_1(i__),e_S25_1(i__),e_S60_1(i__),
+ e_S100_1(i__),PosMin_1(i__),PosMax_1(i__),PosPA_1(i__),
+ Cirr1_1(i__),Cirr2_1(i__),Cirr3_1(i__),CC1_1(i__),CC2_1(i__),
+ CC3_1(i__),CC4_1(i__),CONF_1(i__),DISC_1(i__),HSD_1(i__),
+ PNEARH_1(i__),PNEARW_1(i__),CES12_1(i__),CES25_1(i__),
+ CES60_1(i__),CES100_1(i__),SNR60_1(i__),NHCON_1(i__),
+ MHCON_1(i__),FCOR60_1(i__),PLATE_1(i__),NOM_1(i__),NUM_1(i__),
+ RAo_h_1(i__),RAo_m_1(i__),RAo_s_1(i__),DEo___1(i__),
+ DEo_d_1(i__),DEo_m_1(i__),DEo_s_1(i__),ADX_1(i__),ADY_1(i__),
+ ALK_1(i__),REL_1(i__),CAND_1(i__),IOPT_1(i__),Bmag_1(i__),
+ MajAxiO_1(i__),MinAxiO_1(i__),PAO_1(i__),SuBr_1(i__),
+ IDTYPE_1(i__),Oname_2(i__),mag_1(i__),D1X_1(i__),AR_1(i__),
+ Hvel_1(i__),e_Hvel_1(i__),r_Hvel_1(i__),Hvel2_1(i__),
+ e_Hvel2_1(i__),r_Hvel2_1(i__),TType_1(i__),PGC_1(i__),
+ RAdeg_2(i__),DEdeg_2(i__),GLON_2(i__),GLAT_2(i__),cz_1(i__),
+ Dist_1(i__),IDclus_1(i__),ALI100_1(i__),ALI100B_1(i__),
+ NLBIN_2(i__),CLASS_2(i__),CSEC_1(i__),Xpos_1(i__),Ypos_1(i__),
+ CSTAT_1(i__),CAREA_1(i__),CZSTAT_2(i__),C12CLAS_1(i__),
+ F3PSC_1(i__),F4PSC_1(i__),ABB_2(i__),CPSCZ_2(i__),
+ MHCONB_2(i__),CQDOT_1(i__),W20_1(i__),e_W20_1(i__),W50_1(i__),
+ e_W50_1(i__),r_W20_1(i__),FP12_1(i__),FP25_1(i__),FP60_1(i__),
+ FP100_1(i__),e_FP12_1(i__),e_FP25_1(i__),e_FP60_1(i__),
+ e_FP100_1(i__),BP12_1(i__),BP25_1(i__),BP60_1(i__),
+ BP100_1(i__),FE12_1(i__),FE25_1(i__),FE60_1(i__),FE100_1(i__),
+ e_FE12_1(i__),e_FE25_1(i__),e_FE60_1(i__),e_FE100_1(i__),
+ BE12_1(i__),BE25_1(i__),BE60_1(i__),BE100_1(i__),
+ WIDTH2_1(i__),e_WIDTH2_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'psczcgvs.dat' ! Short version of psczcg.dat
C Format for file interpretation
4 format(
+ A9,1X,A1,A1,A2,1X,I2,I2,I3,1X,A1,I2,I2,I2,F6.1,F6.1,F8.2,I6,
+ I4,I3,A12,I4,F5.2,1X,A1,1X,A1,F6.2,I6)
C Effective file loading
open(unit=1,status='old',file=
+'psczcgvs.dat')
write(6,*) '....Loading file: psczcgvs.dat'
do i__=1,60
read(1,'(A100)')ar__3
read(ar__3,4)
+ CPNAME_1(i__),CPSCZ_3(i__),CZSTAT_3(i__),CLASS_3(i__),
+ RAh_3(i__),RAm_3(i__),RAds_3(i__),DE__3(i__),DEd_3(i__),
+ DEm_3(i__),DEs_3(i__),GLON_3(i__),GLAT_3(i__),F60_1(i__),
+ Vhel_1(i__),e_Vhel_1(i__),r_Vhel_1(i__),Oname_3(i__),
+ I100B_1(i__),ABB_3(i__),CHSD_1(i__),MHCONB_3(i__),
+ F100_F60_1(i__),NLBIN_3(i__)
RAdeg_3(i__) = rNULL__
DEdeg_3(i__) = rNULL__
c Derive coordinates RAdeg_3 and DEdeg_3 from input data
c (RAdeg_3 and DEdeg_3 are set to rNULL__ when unknown)
if(RAh_3(i__) .GT. -180) RAdeg_3(i__)=RAh_3(i__)*15.
if(RAm_3(i__) .GT. -180) RAdeg_3(i__)=RAdeg_3(i__)+RAm_3(i__)/4.
if(RAds_3(i__) .GT. -180) RAdeg_3(i__)=RAdeg_3(i__)+RAds_3(i__)/2400.
if(DEd_3(i__) .GE. 0) DEdeg_3(i__)=DEd_3(i__)
if(DEm_3(i__) .GE. 0) DEdeg_3(i__)=DEdeg_3(i__)+DEm_3(i__)/60.
if(DEs_3(i__) .GE. 0) DEdeg_3(i__)=DEdeg_3(i__)+DEs_3(i__)/3600.
if(DE__3(i__).EQ.'-'.AND.DEdeg_3(i__).GE.0) DEdeg_3(i__)=-DEdeg_3(i__)
c ..............Just test output...........
write(6,4)
+ CPNAME_1(i__),CPSCZ_3(i__),CZSTAT_3(i__),CLASS_3(i__),
+ RAh_3(i__),RAm_3(i__),RAds_3(i__),DE__3(i__),DEd_3(i__),
+ DEm_3(i__),DEs_3(i__),GLON_3(i__),GLAT_3(i__),F60_1(i__),
+ Vhel_1(i__),e_Vhel_1(i__),r_Vhel_1(i__),Oname_3(i__),
+ I100B_1(i__),ABB_3(i__),CHSD_1(i__),MHCONB_3(i__),
+ F100_F60_1(i__),NLBIN_3(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg_3(i__),DEdeg_3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end