Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/594/A26 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. J/A+A/594/A26 Second Planck Catalogue of Compact Sources (PCCS2) (Planck+, 2016)
*================================================================================
*Planck 2015 results. XXVI. The Second Planck Catalogue of Compact Sources.
* Planck collaboration
* Ade P.A.R., Aghanim N., Argueso F., Arnaud M., Ashdown M., Aumont J.,
* Baccigalupi C., Banday A.J., Barreiro R.B., Bartolo N., Battaner E.,
* Beichman C., Benabed K., Benoit A., Benoit-Levy A., Bernard J.-P.,
* Bersanelli M., Bielewicz P., Bock J.J., Bohringer H., Bonaldi A.,
* Bonavera L., Bond J.R., Borrill J., Bouchet F.R., Boulanger F., Bucher M.,
* Burigana C., Butler R.C., Calabrese E., Cardoso J.-F., Carvalho P.,
* Catalano A., Challinor A., Chamballu A., Chary R.-R., Chiang H.C.,
* Christensen P.R., Clemens M., Clements D.L., Colombi S., Colombo L.P.L.,
* Combet C., Couchot F., Coulais A., Crill B.P., Curto A., Cuttaia F.,
* Danese L., Davies R.D., Davis R.J., De Bernardis P., De Rosa A.,
* De Zotti G., Delabrouille J., Desert F.-X., Dickinson C., Diego J.M.,
* Dole H., Donzelli S., Dore O., Douspis M., Ducout A., Dupac X.,
* Efstathiou G., Elsner F., Ensslin T.A., Eriksen H.K., Falgarone E.,
* Fergusson J., Finelli F., Forni O., Frailis M., Fraisse A.A.,
* Franceschi E., Frejsel A., Galeotta S., Galli S., Ganga K., Giard M.,
* Giraud-Heraud Y., Gjerlow E., Gonzalez-Nuevo J., Gorski K.M., Gratton S.,
* Gregorio A., Gruppuso A., Gudmundsson J.E., Hansen F.K., Hanson D.,
* Harrison D.L., Helou G., Henrot-Versille S., Hernandez-Monteagudo C.,
* Herranz D., Hildebrandt S.R., Hivon E., Hobson M., Holmes W.A.,
* Hornstrup A., Hovest W., Huffenberger K.M., Hurier G., Jaffe A.H.,
* Jaffe T.R., Jones W.C., Juvela M., Keihanen E., Keskitalo R., Kisner T.S.,
* Kneissl R., Knoche J., Kunz M., Kurki-Suonio H., Lagache G.,
* Lahteenmaki A., Lamarre J.-M., Lasenby A., Lattanzi M., Lawrence C.R.,
* Leahy J.P., Leonardi R., Leon-Tavares J., Lesgourgues J., Levrier F.,
* Liguori M., Lilje P.B., Linden-Vornle M., Lopez-Caniego M., Lubin P.M.,
* Macias-Perez J.F., Maggio G., Maino D., Mandolesi N., Mangilli A.,
* Maris M., Marshall D.J., Martin P.G., Martinez-Gonzalez E., Masi S.,
* Matarrese S., Mcgehee P., Meinhold P.R., Melchiorri A., Mendes L.,
* Mennella A., Migliaccio M., Mitra S., Miville-Deschenes M.-A., Moneti A.,
* Montier L., Morgante G., Mortlock D., Moss A., Munshi D., Murphy J.A.,
* Naselsky P., Nati F., Natoli P., Negrello M., Netterfield C.B.,
* Norgaard-Nielsen H.U., Noviello F., Novikov D., Novikov I., Oxborrow C.A.,
* Paci F., Pagano L., Pajot F., Paladini R., Paoletti D., Partridge B.,
* Pasian F., Patanchon G., Pearson T.J., Perdereau O., Perotto L.,
* Perrotta F., Pettorino V., Piacentini F., Piat M., Pierpaoli E.,
* Pietrobon D., Plaszczynski S., Pointecouteau E., Polenta G., Pratt G.W.,
* Prezeau G., Prunet S., Puget J.-L., Rachen J.P., Reach W.T., Rebolo R.,
* Reinecke M., Remazeilles M., Renault C., Renzi A., Ristorcelli I.,
* Rocha G., Rosset C., Rossetti M., Roudier G., Rowan-Robinson M.,
* Rubino-Martin J.A., Rusholme B., Sandri M., Sanghera H.S., Santos D.,
* Savelainen M., Savini G., Scott D., Seiffert M.D., Shellard E.P.S.,
* Spencer L.D., Stolyarov V., Sudiwala R., Sunyaev R., Sutton D.,
* Suur-Uski A.-S., Sygnet J.-F., Tauber J.A., Terenzi L., Toffolatti L.,
* Tomasi M., Tornikoski M., Tristram M., Tucci M., Tuovinen J., Turler M.,
* Umana G., Valenziano L., Valiviita J., Van Tent B., Vielva P., Villa F.,
* Wade L.A., Walter B., Wandelt B.D., Wehus I.K., Yvon D., Zacchei A.,
* Zonca A.
* <Astron. Astrophys., 594, A26 (2016)>
* =2016A&A...594A..26P (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'pccs030.dat' ! Second Catalogue of Compact Sources, 30GHz
integer*4 nr__
parameter (nr__=1560) ! Number of records
character*573 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*23 Name ! Source name, PCCS2 FFF GLLL.ll+BB.bbb
* (NAME) (G1)
real*8 GLON ! (deg) Galactic longitude based on extraction
* algorithm (GLON)
real*8 GLAT ! (deg) Galactic latitude based on extraction
* algorithm (GLAT)
real*8 RAdeg ! (deg) Right ascension (J2000) transformed
* from (GLON,GLAT) (RA)
real*8 DEdeg ! (deg) Declination (J2000) transformed
* from (GLON,GLAT) (DEC)
real*4 DetFlux ! (mJy) Flux density of source as determined by
* detection method (DETFLUX)
real*4 e_DetFlux ! (mJy) Uncertainty (1 sigma) in DetFlux
* (DETFLUX_ERR)
real*4 AperFlux ! (mJy) Flux density of source as determined from
* the aperture photometry (APERFLUX)
real*4 e_AperFlux ! (mJy) Uncertainty (1 sigma) in AperFlux
* (APERFLUX_ERR)
real*4 PSFFlux ! (mJy) Flux density of source as determined from
* PSF fitting (PSFFLUX)
real*4 e_PSFFlux ! (mJy) Uncertainty (1 sigma) in PSFFlux
* (PSFFLUX_ERR)
real*4 GauFlux ! (mJy) Flux density of source as determined from
* 2-D Gaussian fitting (GAUFLUX)
real*4 e_GauFlux ! (mJy) ?=- Uncertainty (1 sigma) in GauFlux
* (GAUFLUX_ERR)
real*8 GauSemi1 ! (arcmin) Gaussian fit along axis 1 (FWHM)
* (GAU_SEMI1) (G2)
real*8 e_GauSemi1 ! (arcmin) ?=- Uncertainty (1 sigma) in GauSemi1
* (GAU_SEMI1_ERR)
real*8 GauSemi2 ! (arcmin) Gaussian fit along axis 2 (FWHM)
* (GAU_SEMI2) (G2)
real*8 e_GauSemi2 ! (arcmin) ?=- Uncertainty (1 sigma) in GauSemi2
* (GAU_SEMI2_ERR)
real*4 GauTheta ! (deg) Gaussian fit orientation angle
* (GAU_THETA) (G2)
real*4 e_GauTheta ! (deg) ?=- Uncertainty (1 sigma) in GauTheta
* (GAU_THETA_ERR)
real*8 GauFWHMEff ! (arcmin) Gaussian fit effective FWHM (GAU_FWHM_EFF)
integer*4 Extended ! ?=-1 Extended source flag (EXTENDED) (G5)
integer*4 ExtVal ! ?=-1 External validation flag (EXT_VAL) (G6)
character*24 ERCSC ! Name of the ERCSC counterpart, if any
* (ERCSC)
character*23 PCCS1 ! Name of the PCCS counterpart, if any (PCCS1)
integer*4 HRel ! ?=-1 Highest reliability catalogue to which
* the source belong
* (HIGHEST_RELIABILITY_CAT) (G7)
*
*---------- only for 100, 143, 217 and 353GHz -----
integer*4 CirrusN ! ?=- Number of sources (S/N>5) detected at
* 857GHz within a 1-degree radius (CIRRUS_N)
real*8 SkyBri ! (MJy/sr) ?=- The mean 857GHz brightness within a
* 2{deg} radius. This may be used as another
* indicator of cirrus contamination
* (SKY_BRIGHTNESS)
*
*------------- Polarization measurements (30-353GHz only) ----------
real*8 P ! (mJy) ?=- Polarization flux density of the sources
* as determined by a matched filter (P) (G3)
real*8 e_P ! (mJy) ?=- Uncertainty (1 sigma) in P (P_ERR)
real*8 AnP ! (deg) ?=- Orientation of polarization with respect
* to NGP (ANGLE_P) (G2) (G3)
real*8 e_AnP ! (deg) ?=- Uncertainty (1 sigma) in AnP
* (ANGLE_P_ERR) (G3)
real*8 ApP ! (mJy) ?=- Orientation of polarization with respect
* to NGP (APER_P) (G2) (G3)
real*8 e_ApP ! (mJy) ?=- Uncertainty (1 sigma) in AperP
* (APER_P_ERR) (G3)
real*8 ApAnP ! (deg) ?=- Orientation of polarization with
* respect to NGP (APER_ANGLE_P) (G2) (G3)
real*8 e_ApAnP ! (deg) ?=- Uncertainty (1 sigma) in ApAneP
* (APER_ANGLE_P_ERR) (G3)
real*8 PUL ! (mJy) ?=- Polarization flux density 99.99% upper
* limit (P_UPPER_LIMIT) (G9)
real*8 ApPUL ! (mJy) ?=- Polarization flux density 99.99% upper
* limit (APER_P_UPPER_LIMIT) (G9)
*
*---------- Marginal polarization measurements (100-353GHz only) --------
integer*4 PStat ! ?=- Polarization detection status
* (P_STAT) (G4)
real*8 PX ! (mJy) ?=- Polarization flux density of the sources
* as determined by a matched filter using
* Bayesian polarization estimator (PX)
real*8 e_PX ! (mJy) ?=- PX uncertainty; lower 95% error bar
* (PX_ERR_LOWER)
real*8 E_PX_1 ! (mJy) ?=- PX uncertainty; upper 95% error bar
* (PX_ERR_UPPER)
real*8 AnPX ! (deg) ?=- Orientation of polarization with respect
* to NGP using Bayesian polarization
* estimator (ANGLE_PX) (G2)
real*8 e_AnPX ! (deg) ?=- AnglePX uncertainty; lower 95% error bar
* (ANGLE_PX_ERR_LOWER)
real*8 E_AnPX_1 ! (deg) ?=- AnglePX uncertainty; upper 95% error bar
* (ANGLE_PX_ERR_UPPER)
*
*------------------Fluxes at 217, 353 and 545 GHz (857GHz only)------------------
real*8 ApFlux217 ! (mJy) ?=- Flux density at 217GHz of source as
* determined from the aperture photometry
* (APERFLUX_227)
real*8 e_ApFlux217 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux217
* (APERFLUX_ERR_227)
real*8 ApFlux353 ! (mJy) ?=- Flux density at 353GHz of source as
* determined from the aperture photometry
* (APERFLUX_353)
real*8 e_ApFlux353 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux353
* (APERFLUX_ERR_353)
real*8 ApFlux545 ! (mJy) ?=- Flux density at 545GHz of source as
* determined from the aperture photometry
* (APERFLUX_545)
real*8 e_ApFlux545 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux545
* (APERFLUX_ERR_545)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'pccs044.dat' ! Second Catalogue of Compact Sources, 44GHz
integer*4 nr__1
parameter (nr__1=934) ! Number of records
character*573 ar__1 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*23 Name_1 ! Source name, PCCS2 FFF GLLL.ll+BB.bbb
* (NAME) (G1)
real*8 GLON_1 ! (deg) Galactic longitude based on extraction
* algorithm (GLON)
real*8 GLAT_1 ! (deg) Galactic latitude based on extraction
* algorithm (GLAT)
real*8 RAdeg_1 ! (deg) Right ascension (J2000) transformed
* from (GLON,GLAT) (RA)
real*8 DEdeg_1 ! (deg) Declination (J2000) transformed
* from (GLON,GLAT) (DEC)
real*4 DetFlux_1 ! (mJy) Flux density of source as determined by
* detection method (DETFLUX)
real*4 e_DetFlux_1 ! (mJy) Uncertainty (1 sigma) in DetFlux
* (DETFLUX_ERR)
real*4 AperFlux_1 ! (mJy) Flux density of source as determined from
* the aperture photometry (APERFLUX)
real*4 e_AperFlux_1 ! (mJy) Uncertainty (1 sigma) in AperFlux
* (APERFLUX_ERR)
real*4 PSFFlux_1 ! (mJy) Flux density of source as determined from
* PSF fitting (PSFFLUX)
real*4 e_PSFFlux_1 ! (mJy) Uncertainty (1 sigma) in PSFFlux
* (PSFFLUX_ERR)
real*4 GauFlux_1 ! (mJy) Flux density of source as determined from
* 2-D Gaussian fitting (GAUFLUX)
real*4 e_GauFlux_1 ! (mJy) ?=- Uncertainty (1 sigma) in GauFlux
* (GAUFLUX_ERR)
real*8 GauSemi1_1 ! (arcmin) Gaussian fit along axis 1 (FWHM)
* (GAU_SEMI1) (G2)
real*8 e_GauSemi1_1 ! (arcmin) ?=- Uncertainty (1 sigma) in GauSemi1
* (GAU_SEMI1_ERR)
real*8 GauSemi2_1 ! (arcmin) Gaussian fit along axis 2 (FWHM)
* (GAU_SEMI2) (G2)
real*8 e_GauSemi2_1 ! (arcmin) ?=- Uncertainty (1 sigma) in GauSemi2
* (GAU_SEMI2_ERR)
real*4 GauTheta_1 ! (deg) Gaussian fit orientation angle
* (GAU_THETA) (G2)
real*4 e_GauTheta_1 ! (deg) ?=- Uncertainty (1 sigma) in GauTheta
* (GAU_THETA_ERR)
real*8 GauFWHMEff_1 ! (arcmin) Gaussian fit effective FWHM (GAU_FWHM_EFF)
integer*4 Extended_1 ! ?=-1 Extended source flag (EXTENDED) (G5)
integer*4 ExtVal_1 ! ?=-1 External validation flag (EXT_VAL) (G6)
character*24 ERCSC_1 ! Name of the ERCSC counterpart, if any
* (ERCSC)
character*23 PCCS1_1 ! Name of the PCCS counterpart, if any (PCCS1)
integer*4 HRel_1 ! ?=-1 Highest reliability catalogue to which
* the source belong
* (HIGHEST_RELIABILITY_CAT) (G7)
*
*---------- only for 100, 143, 217 and 353GHz -----
integer*4 CirrusN_1 ! ?=- Number of sources (S/N>5) detected at
* 857GHz within a 1-degree radius (CIRRUS_N)
real*8 SkyBri_1 ! (MJy/sr) ?=- The mean 857GHz brightness within a
* 2{deg} radius. This may be used as another
* indicator of cirrus contamination
* (SKY_BRIGHTNESS)
*
*------------- Polarization measurements (30-353GHz only) ----------
real*8 P_1 ! (mJy) ?=- Polarization flux density of the sources
* as determined by a matched filter (P) (G3)
real*8 e_P_1 ! (mJy) ?=- Uncertainty (1 sigma) in P (P_ERR)
real*8 AnP_1 ! (deg) ?=- Orientation of polarization with respect
* to NGP (ANGLE_P) (G2) (G3)
real*8 e_AnP_1 ! (deg) ?=- Uncertainty (1 sigma) in AnP
* (ANGLE_P_ERR) (G3)
real*8 ApP_1 ! (mJy) ?=- Orientation of polarization with respect
* to NGP (APER_P) (G2) (G3)
real*8 e_ApP_1 ! (mJy) ?=- Uncertainty (1 sigma) in AperP
* (APER_P_ERR) (G3)
real*8 ApAnP_1 ! (deg) ?=- Orientation of polarization with
* respect to NGP (APER_ANGLE_P) (G2) (G3)
real*8 e_ApAnP_1 ! (deg) ?=- Uncertainty (1 sigma) in ApAneP
* (APER_ANGLE_P_ERR) (G3)
real*8 PUL_1 ! (mJy) ?=- Polarization flux density 99.99% upper
* limit (P_UPPER_LIMIT) (G9)
real*8 ApPUL_1 ! (mJy) ?=- Polarization flux density 99.99% upper
* limit (APER_P_UPPER_LIMIT) (G9)
*
*---------- Marginal polarization measurements (100-353GHz only) --------
integer*4 PStat_1 ! ?=- Polarization detection status
* (P_STAT) (G4)
real*8 PX_1 ! (mJy) ?=- Polarization flux density of the sources
* as determined by a matched filter using
* Bayesian polarization estimator (PX)
real*8 e_PX_2 ! (mJy) ?=- PX uncertainty; lower 95% error bar
* (PX_ERR_LOWER)
real*8 E_PX_3 ! (mJy) ?=- PX uncertainty; upper 95% error bar
* (PX_ERR_UPPER)
real*8 AnPX_1 ! (deg) ?=- Orientation of polarization with respect
* to NGP using Bayesian polarization
* estimator (ANGLE_PX) (G2)
real*8 e_AnPX_2 ! (deg) ?=- AnglePX uncertainty; lower 95% error bar
* (ANGLE_PX_ERR_LOWER)
real*8 E_AnPX_3 ! (deg) ?=- AnglePX uncertainty; upper 95% error bar
* (ANGLE_PX_ERR_UPPER)
*
*------------------Fluxes at 217, 353 and 545 GHz (857GHz only)------------------
real*8 ApFlux217_1 ! (mJy) ?=- Flux density at 217GHz of source as
* determined from the aperture photometry
* (APERFLUX_227)
real*8 e_ApFlux217_1 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux217
* (APERFLUX_ERR_227)
real*8 ApFlux353_1 ! (mJy) ?=- Flux density at 353GHz of source as
* determined from the aperture photometry
* (APERFLUX_353)
real*8 e_ApFlux353_1 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux353
* (APERFLUX_ERR_353)
real*8 ApFlux545_1 ! (mJy) ?=- Flux density at 545GHz of source as
* determined from the aperture photometry
* (APERFLUX_545)
real*8 e_ApFlux545_1 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux545
* (APERFLUX_ERR_545)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'pccs070.dat' ! Second Catalogue of Compact Sources, 70GHz
integer*4 nr__2
parameter (nr__2=1296) ! Number of records
character*573 ar__2 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*23 Name_2 ! Source name, PCCS2 FFF GLLL.ll+BB.bbb
* (NAME) (G1)
real*8 GLON_2 ! (deg) Galactic longitude based on extraction
* algorithm (GLON)
real*8 GLAT_2 ! (deg) Galactic latitude based on extraction
* algorithm (GLAT)
real*8 RAdeg_2 ! (deg) Right ascension (J2000) transformed
* from (GLON,GLAT) (RA)
real*8 DEdeg_2 ! (deg) Declination (J2000) transformed
* from (GLON,GLAT) (DEC)
real*4 DetFlux_2 ! (mJy) Flux density of source as determined by
* detection method (DETFLUX)
real*4 e_DetFlux_2 ! (mJy) Uncertainty (1 sigma) in DetFlux
* (DETFLUX_ERR)
real*4 AperFlux_2 ! (mJy) Flux density of source as determined from
* the aperture photometry (APERFLUX)
real*4 e_AperFlux_2 ! (mJy) Uncertainty (1 sigma) in AperFlux
* (APERFLUX_ERR)
real*4 PSFFlux_2 ! (mJy) Flux density of source as determined from
* PSF fitting (PSFFLUX)
real*4 e_PSFFlux_2 ! (mJy) Uncertainty (1 sigma) in PSFFlux
* (PSFFLUX_ERR)
real*4 GauFlux_2 ! (mJy) Flux density of source as determined from
* 2-D Gaussian fitting (GAUFLUX)
real*4 e_GauFlux_2 ! (mJy) ?=- Uncertainty (1 sigma) in GauFlux
* (GAUFLUX_ERR)
real*8 GauSemi1_2 ! (arcmin) Gaussian fit along axis 1 (FWHM)
* (GAU_SEMI1) (G2)
real*8 e_GauSemi1_2 ! (arcmin) ?=- Uncertainty (1 sigma) in GauSemi1
* (GAU_SEMI1_ERR)
real*8 GauSemi2_2 ! (arcmin) Gaussian fit along axis 2 (FWHM)
* (GAU_SEMI2) (G2)
real*8 e_GauSemi2_2 ! (arcmin) ?=- Uncertainty (1 sigma) in GauSemi2
* (GAU_SEMI2_ERR)
real*4 GauTheta_2 ! (deg) Gaussian fit orientation angle
* (GAU_THETA) (G2)
real*4 e_GauTheta_2 ! (deg) ?=- Uncertainty (1 sigma) in GauTheta
* (GAU_THETA_ERR)
real*8 GauFWHMEff_2 ! (arcmin) Gaussian fit effective FWHM (GAU_FWHM_EFF)
integer*4 Extended_2 ! ?=-1 Extended source flag (EXTENDED) (G5)
integer*4 ExtVal_2 ! ?=-1 External validation flag (EXT_VAL) (G6)
character*24 ERCSC_2 ! Name of the ERCSC counterpart, if any
* (ERCSC)
character*23 PCCS1_2 ! Name of the PCCS counterpart, if any (PCCS1)
integer*4 HRel_2 ! ?=-1 Highest reliability catalogue to which
* the source belong
* (HIGHEST_RELIABILITY_CAT) (G7)
*
*---------- only for 100, 143, 217 and 353GHz -----
integer*4 CirrusN_2 ! ?=- Number of sources (S/N>5) detected at
* 857GHz within a 1-degree radius (CIRRUS_N)
real*8 SkyBri_2 ! (MJy/sr) ?=- The mean 857GHz brightness within a
* 2{deg} radius. This may be used as another
* indicator of cirrus contamination
* (SKY_BRIGHTNESS)
*
*------------- Polarization measurements (30-353GHz only) ----------
real*8 P_2 ! (mJy) ?=- Polarization flux density of the sources
* as determined by a matched filter (P) (G3)
real*8 e_P_2 ! (mJy) ?=- Uncertainty (1 sigma) in P (P_ERR)
real*8 AnP_2 ! (deg) ?=- Orientation of polarization with respect
* to NGP (ANGLE_P) (G2) (G3)
real*8 e_AnP_2 ! (deg) ?=- Uncertainty (1 sigma) in AnP
* (ANGLE_P_ERR) (G3)
real*8 ApP_2 ! (mJy) ?=- Orientation of polarization with respect
* to NGP (APER_P) (G2) (G3)
real*8 e_ApP_2 ! (mJy) ?=- Uncertainty (1 sigma) in AperP
* (APER_P_ERR) (G3)
real*8 ApAnP_2 ! (deg) ?=- Orientation of polarization with
* respect to NGP (APER_ANGLE_P) (G2) (G3)
real*8 e_ApAnP_2 ! (deg) ?=- Uncertainty (1 sigma) in ApAneP
* (APER_ANGLE_P_ERR) (G3)
real*8 PUL_2 ! (mJy) ?=- Polarization flux density 99.99% upper
* limit (P_UPPER_LIMIT) (G9)
real*8 ApPUL_2 ! (mJy) ?=- Polarization flux density 99.99% upper
* limit (APER_P_UPPER_LIMIT) (G9)
*
*---------- Marginal polarization measurements (100-353GHz only) --------
integer*4 PStat_2 ! ?=- Polarization detection status
* (P_STAT) (G4)
real*8 PX_2 ! (mJy) ?=- Polarization flux density of the sources
* as determined by a matched filter using
* Bayesian polarization estimator (PX)
real*8 e_PX_4 ! (mJy) ?=- PX uncertainty; lower 95% error bar
* (PX_ERR_LOWER)
real*8 E_PX_5 ! (mJy) ?=- PX uncertainty; upper 95% error bar
* (PX_ERR_UPPER)
real*8 AnPX_2 ! (deg) ?=- Orientation of polarization with respect
* to NGP using Bayesian polarization
* estimator (ANGLE_PX) (G2)
real*8 e_AnPX_4 ! (deg) ?=- AnglePX uncertainty; lower 95% error bar
* (ANGLE_PX_ERR_LOWER)
real*8 E_AnPX_5 ! (deg) ?=- AnglePX uncertainty; upper 95% error bar
* (ANGLE_PX_ERR_UPPER)
*
*------------------Fluxes at 217, 353 and 545 GHz (857GHz only)------------------
real*8 ApFlux217_2 ! (mJy) ?=- Flux density at 217GHz of source as
* determined from the aperture photometry
* (APERFLUX_227)
real*8 e_ApFlux217_2 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux217
* (APERFLUX_ERR_227)
real*8 ApFlux353_2 ! (mJy) ?=- Flux density at 353GHz of source as
* determined from the aperture photometry
* (APERFLUX_353)
real*8 e_ApFlux353_2 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux353
* (APERFLUX_ERR_353)
real*8 ApFlux545_2 ! (mJy) ?=- Flux density at 545GHz of source as
* determined from the aperture photometry
* (APERFLUX_545)
real*8 e_ApFlux545_2 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux545
* (APERFLUX_ERR_545)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'pccs100.dat' ! Second Catalogue of Compact Sources, 100GHz
integer*4 nr__3
parameter (nr__3=1742) ! Number of records
character*573 ar__3 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*23 Name_3 ! Source name, PCCS2 FFF GLLL.ll+BB.bbb
* (NAME) (G1)
real*8 GLON_3 ! (deg) Galactic longitude based on extraction
* algorithm (GLON)
real*8 GLAT_3 ! (deg) Galactic latitude based on extraction
* algorithm (GLAT)
real*8 RAdeg_3 ! (deg) Right ascension (J2000) transformed
* from (GLON,GLAT) (RA)
real*8 DEdeg_3 ! (deg) Declination (J2000) transformed
* from (GLON,GLAT) (DEC)
real*4 DetFlux_3 ! (mJy) Flux density of source as determined by
* detection method (DETFLUX)
real*4 e_DetFlux_3 ! (mJy) Uncertainty (1 sigma) in DetFlux
* (DETFLUX_ERR)
real*4 AperFlux_3 ! (mJy) Flux density of source as determined from
* the aperture photometry (APERFLUX)
real*4 e_AperFlux_3 ! (mJy) Uncertainty (1 sigma) in AperFlux
* (APERFLUX_ERR)
real*4 PSFFlux_3 ! (mJy) Flux density of source as determined from
* PSF fitting (PSFFLUX)
real*4 e_PSFFlux_3 ! (mJy) Uncertainty (1 sigma) in PSFFlux
* (PSFFLUX_ERR)
real*4 GauFlux_3 ! (mJy) Flux density of source as determined from
* 2-D Gaussian fitting (GAUFLUX)
real*4 e_GauFlux_3 ! (mJy) ?=- Uncertainty (1 sigma) in GauFlux
* (GAUFLUX_ERR)
real*8 GauSemi1_3 ! (arcmin) Gaussian fit along axis 1 (FWHM)
* (GAU_SEMI1) (G2)
real*8 e_GauSemi1_3 ! (arcmin) ?=- Uncertainty (1 sigma) in GauSemi1
* (GAU_SEMI1_ERR)
real*8 GauSemi2_3 ! (arcmin) Gaussian fit along axis 2 (FWHM)
* (GAU_SEMI2) (G2)
real*8 e_GauSemi2_3 ! (arcmin) ?=- Uncertainty (1 sigma) in GauSemi2
* (GAU_SEMI2_ERR)
real*4 GauTheta_3 ! (deg) Gaussian fit orientation angle
* (GAU_THETA) (G2)
real*4 e_GauTheta_3 ! (deg) ?=- Uncertainty (1 sigma) in GauTheta
* (GAU_THETA_ERR)
real*8 GauFWHMEff_3 ! (arcmin) Gaussian fit effective FWHM (GAU_FWHM_EFF)
integer*4 Extended_3 ! ?=-1 Extended source flag (EXTENDED) (G5)
integer*4 ExtVal_3 ! ?=-1 External validation flag (EXT_VAL) (G6)
character*24 ERCSC_3 ! Name of the ERCSC counterpart, if any
* (ERCSC)
character*23 PCCS1_3 ! Name of the PCCS counterpart, if any (PCCS1)
integer*4 HRel_3 ! ?=-1 Highest reliability catalogue to which
* the source belong
* (HIGHEST_RELIABILITY_CAT) (G7)
*
*---------- only for 100, 143, 217 and 353GHz -----
integer*4 CirrusN_3 ! ?=- Number of sources (S/N>5) detected at
* 857GHz within a 1-degree radius (CIRRUS_N)
real*8 SkyBri_3 ! (MJy/sr) ?=- The mean 857GHz brightness within a
* 2{deg} radius. This may be used as another
* indicator of cirrus contamination
* (SKY_BRIGHTNESS)
*
*------------- Polarization measurements (30-353GHz only) ----------
real*8 P_3 ! (mJy) ?=- Polarization flux density of the sources
* as determined by a matched filter (P) (G3)
real*8 e_P_3 ! (mJy) ?=- Uncertainty (1 sigma) in P (P_ERR)
real*8 AnP_3 ! (deg) ?=- Orientation of polarization with respect
* to NGP (ANGLE_P) (G2) (G3)
real*8 e_AnP_3 ! (deg) ?=- Uncertainty (1 sigma) in AnP
* (ANGLE_P_ERR) (G3)
real*8 ApP_3 ! (mJy) ?=- Orientation of polarization with respect
* to NGP (APER_P) (G2) (G3)
real*8 e_ApP_3 ! (mJy) ?=- Uncertainty (1 sigma) in AperP
* (APER_P_ERR) (G3)
real*8 ApAnP_3 ! (deg) ?=- Orientation of polarization with
* respect to NGP (APER_ANGLE_P) (G2) (G3)
real*8 e_ApAnP_3 ! (deg) ?=- Uncertainty (1 sigma) in ApAneP
* (APER_ANGLE_P_ERR) (G3)
real*8 PUL_3 ! (mJy) ?=- Polarization flux density 99.99% upper
* limit (P_UPPER_LIMIT) (G9)
real*8 ApPUL_3 ! (mJy) ?=- Polarization flux density 99.99% upper
* limit (APER_P_UPPER_LIMIT) (G9)
*
*---------- Marginal polarization measurements (100-353GHz only) --------
integer*4 PStat_3 ! ?=- Polarization detection status
* (P_STAT) (G4)
real*8 PX_3 ! (mJy) ?=- Polarization flux density of the sources
* as determined by a matched filter using
* Bayesian polarization estimator (PX)
real*8 e_PX_6 ! (mJy) ?=- PX uncertainty; lower 95% error bar
* (PX_ERR_LOWER)
real*8 E_PX_7 ! (mJy) ?=- PX uncertainty; upper 95% error bar
* (PX_ERR_UPPER)
real*8 AnPX_3 ! (deg) ?=- Orientation of polarization with respect
* to NGP using Bayesian polarization
* estimator (ANGLE_PX) (G2)
real*8 e_AnPX_6 ! (deg) ?=- AnglePX uncertainty; lower 95% error bar
* (ANGLE_PX_ERR_LOWER)
real*8 E_AnPX_7 ! (deg) ?=- AnglePX uncertainty; upper 95% error bar
* (ANGLE_PX_ERR_UPPER)
*
*------------------Fluxes at 217, 353 and 545 GHz (857GHz only)------------------
real*8 ApFlux217_3 ! (mJy) ?=- Flux density at 217GHz of source as
* determined from the aperture photometry
* (APERFLUX_227)
real*8 e_ApFlux217_3 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux217
* (APERFLUX_ERR_227)
real*8 ApFlux353_3 ! (mJy) ?=- Flux density at 353GHz of source as
* determined from the aperture photometry
* (APERFLUX_353)
real*8 e_ApFlux353_3 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux353
* (APERFLUX_ERR_353)
real*8 ApFlux545_3 ! (mJy) ?=- Flux density at 545GHz of source as
* determined from the aperture photometry
* (APERFLUX_545)
real*8 e_ApFlux545_3 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux545
* (APERFLUX_ERR_545)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'pccs143.dat' ! Second Catalogue of Compact Sources, 143GHz
integer*4 nr__4
parameter (nr__4=2160) ! Number of records
character*573 ar__4 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*23 Name_4 ! Source name, PCCS2 FFF GLLL.ll+BB.bbb
* (NAME) (G1)
real*8 GLON_4 ! (deg) Galactic longitude based on extraction
* algorithm (GLON)
real*8 GLAT_4 ! (deg) Galactic latitude based on extraction
* algorithm (GLAT)
real*8 RAdeg_4 ! (deg) Right ascension (J2000) transformed
* from (GLON,GLAT) (RA)
real*8 DEdeg_4 ! (deg) Declination (J2000) transformed
* from (GLON,GLAT) (DEC)
real*4 DetFlux_4 ! (mJy) Flux density of source as determined by
* detection method (DETFLUX)
real*4 e_DetFlux_4 ! (mJy) Uncertainty (1 sigma) in DetFlux
* (DETFLUX_ERR)
real*4 AperFlux_4 ! (mJy) Flux density of source as determined from
* the aperture photometry (APERFLUX)
real*4 e_AperFlux_4 ! (mJy) Uncertainty (1 sigma) in AperFlux
* (APERFLUX_ERR)
real*4 PSFFlux_4 ! (mJy) Flux density of source as determined from
* PSF fitting (PSFFLUX)
real*4 e_PSFFlux_4 ! (mJy) Uncertainty (1 sigma) in PSFFlux
* (PSFFLUX_ERR)
real*4 GauFlux_4 ! (mJy) Flux density of source as determined from
* 2-D Gaussian fitting (GAUFLUX)
real*4 e_GauFlux_4 ! (mJy) ?=- Uncertainty (1 sigma) in GauFlux
* (GAUFLUX_ERR)
real*8 GauSemi1_4 ! (arcmin) Gaussian fit along axis 1 (FWHM)
* (GAU_SEMI1) (G2)
real*8 e_GauSemi1_4 ! (arcmin) ?=- Uncertainty (1 sigma) in GauSemi1
* (GAU_SEMI1_ERR)
real*8 GauSemi2_4 ! (arcmin) Gaussian fit along axis 2 (FWHM)
* (GAU_SEMI2) (G2)
real*8 e_GauSemi2_4 ! (arcmin) ?=- Uncertainty (1 sigma) in GauSemi2
* (GAU_SEMI2_ERR)
real*4 GauTheta_4 ! (deg) Gaussian fit orientation angle
* (GAU_THETA) (G2)
real*4 e_GauTheta_4 ! (deg) ?=- Uncertainty (1 sigma) in GauTheta
* (GAU_THETA_ERR)
real*8 GauFWHMEff_4 ! (arcmin) Gaussian fit effective FWHM (GAU_FWHM_EFF)
integer*4 Extended_4 ! ?=-1 Extended source flag (EXTENDED) (G5)
integer*4 ExtVal_4 ! ?=-1 External validation flag (EXT_VAL) (G6)
character*24 ERCSC_4 ! Name of the ERCSC counterpart, if any
* (ERCSC)
character*23 PCCS1_4 ! Name of the PCCS counterpart, if any (PCCS1)
integer*4 HRel_4 ! ?=-1 Highest reliability catalogue to which
* the source belong
* (HIGHEST_RELIABILITY_CAT) (G7)
*
*---------- only for 100, 143, 217 and 353GHz -----
integer*4 CirrusN_4 ! ?=- Number of sources (S/N>5) detected at
* 857GHz within a 1-degree radius (CIRRUS_N)
real*8 SkyBri_4 ! (MJy/sr) ?=- The mean 857GHz brightness within a
* 2{deg} radius. This may be used as another
* indicator of cirrus contamination
* (SKY_BRIGHTNESS)
*
*------------- Polarization measurements (30-353GHz only) ----------
real*8 P_4 ! (mJy) ?=- Polarization flux density of the sources
* as determined by a matched filter (P) (G3)
real*8 e_P_4 ! (mJy) ?=- Uncertainty (1 sigma) in P (P_ERR)
real*8 AnP_4 ! (deg) ?=- Orientation of polarization with respect
* to NGP (ANGLE_P) (G2) (G3)
real*8 e_AnP_4 ! (deg) ?=- Uncertainty (1 sigma) in AnP
* (ANGLE_P_ERR) (G3)
real*8 ApP_4 ! (mJy) ?=- Orientation of polarization with respect
* to NGP (APER_P) (G2) (G3)
real*8 e_ApP_4 ! (mJy) ?=- Uncertainty (1 sigma) in AperP
* (APER_P_ERR) (G3)
real*8 ApAnP_4 ! (deg) ?=- Orientation of polarization with
* respect to NGP (APER_ANGLE_P) (G2) (G3)
real*8 e_ApAnP_4 ! (deg) ?=- Uncertainty (1 sigma) in ApAneP
* (APER_ANGLE_P_ERR) (G3)
real*8 PUL_4 ! (mJy) ?=- Polarization flux density 99.99% upper
* limit (P_UPPER_LIMIT) (G9)
real*8 ApPUL_4 ! (mJy) ?=- Polarization flux density 99.99% upper
* limit (APER_P_UPPER_LIMIT) (G9)
*
*---------- Marginal polarization measurements (100-353GHz only) --------
integer*4 PStat_4 ! ?=- Polarization detection status
* (P_STAT) (G4)
real*8 PX_4 ! (mJy) ?=- Polarization flux density of the sources
* as determined by a matched filter using
* Bayesian polarization estimator (PX)
real*8 e_PX_8 ! (mJy) ?=- PX uncertainty; lower 95% error bar
* (PX_ERR_LOWER)
real*8 E_PX_9 ! (mJy) ?=- PX uncertainty; upper 95% error bar
* (PX_ERR_UPPER)
real*8 AnPX_4 ! (deg) ?=- Orientation of polarization with respect
* to NGP using Bayesian polarization
* estimator (ANGLE_PX) (G2)
real*8 e_AnPX_8 ! (deg) ?=- AnglePX uncertainty; lower 95% error bar
* (ANGLE_PX_ERR_LOWER)
real*8 E_AnPX_9 ! (deg) ?=- AnglePX uncertainty; upper 95% error bar
* (ANGLE_PX_ERR_UPPER)
*
*------------------Fluxes at 217, 353 and 545 GHz (857GHz only)------------------
real*8 ApFlux217_4 ! (mJy) ?=- Flux density at 217GHz of source as
* determined from the aperture photometry
* (APERFLUX_227)
real*8 e_ApFlux217_4 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux217
* (APERFLUX_ERR_227)
real*8 ApFlux353_4 ! (mJy) ?=- Flux density at 353GHz of source as
* determined from the aperture photometry
* (APERFLUX_353)
real*8 e_ApFlux353_4 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux353
* (APERFLUX_ERR_353)
real*8 ApFlux545_4 ! (mJy) ?=- Flux density at 545GHz of source as
* determined from the aperture photometry
* (APERFLUX_545)
real*8 e_ApFlux545_4 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux545
* (APERFLUX_ERR_545)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'pccs217.dat' ! Second Catalogue of Compact Sources, 217GHz
integer*4 nr__5
parameter (nr__5=2135) ! Number of records
character*573 ar__5 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*23 Name_5 ! Source name, PCCS2 FFF GLLL.ll+BB.bbb
* (NAME) (G1)
real*8 GLON_5 ! (deg) Galactic longitude based on extraction
* algorithm (GLON)
real*8 GLAT_5 ! (deg) Galactic latitude based on extraction
* algorithm (GLAT)
real*8 RAdeg_5 ! (deg) Right ascension (J2000) transformed
* from (GLON,GLAT) (RA)
real*8 DEdeg_5 ! (deg) Declination (J2000) transformed
* from (GLON,GLAT) (DEC)
real*4 DetFlux_5 ! (mJy) Flux density of source as determined by
* detection method (DETFLUX)
real*4 e_DetFlux_5 ! (mJy) Uncertainty (1 sigma) in DetFlux
* (DETFLUX_ERR)
real*4 AperFlux_5 ! (mJy) Flux density of source as determined from
* the aperture photometry (APERFLUX)
real*4 e_AperFlux_5 ! (mJy) Uncertainty (1 sigma) in AperFlux
* (APERFLUX_ERR)
real*4 PSFFlux_5 ! (mJy) Flux density of source as determined from
* PSF fitting (PSFFLUX)
real*4 e_PSFFlux_5 ! (mJy) Uncertainty (1 sigma) in PSFFlux
* (PSFFLUX_ERR)
real*4 GauFlux_5 ! (mJy) Flux density of source as determined from
* 2-D Gaussian fitting (GAUFLUX)
real*4 e_GauFlux_5 ! (mJy) ?=- Uncertainty (1 sigma) in GauFlux
* (GAUFLUX_ERR)
real*8 GauSemi1_5 ! (arcmin) Gaussian fit along axis 1 (FWHM)
* (GAU_SEMI1) (G2)
real*8 e_GauSemi1_5 ! (arcmin) ?=- Uncertainty (1 sigma) in GauSemi1
* (GAU_SEMI1_ERR)
real*8 GauSemi2_5 ! (arcmin) Gaussian fit along axis 2 (FWHM)
* (GAU_SEMI2) (G2)
real*8 e_GauSemi2_5 ! (arcmin) ?=- Uncertainty (1 sigma) in GauSemi2
* (GAU_SEMI2_ERR)
real*4 GauTheta_5 ! (deg) Gaussian fit orientation angle
* (GAU_THETA) (G2)
real*4 e_GauTheta_5 ! (deg) ?=- Uncertainty (1 sigma) in GauTheta
* (GAU_THETA_ERR)
real*8 GauFWHMEff_5 ! (arcmin) Gaussian fit effective FWHM (GAU_FWHM_EFF)
integer*4 Extended_5 ! ?=-1 Extended source flag (EXTENDED) (G5)
integer*4 ExtVal_5 ! ?=-1 External validation flag (EXT_VAL) (G6)
character*24 ERCSC_5 ! Name of the ERCSC counterpart, if any
* (ERCSC)
character*23 PCCS1_5 ! Name of the PCCS counterpart, if any (PCCS1)
integer*4 HRel_5 ! ?=-1 Highest reliability catalogue to which
* the source belong
* (HIGHEST_RELIABILITY_CAT) (G7)
*
*---------- only for 100, 143, 217 and 353GHz -----
integer*4 CirrusN_5 ! ?=- Number of sources (S/N>5) detected at
* 857GHz within a 1-degree radius (CIRRUS_N)
real*8 SkyBri_5 ! (MJy/sr) ?=- The mean 857GHz brightness within a
* 2{deg} radius. This may be used as another
* indicator of cirrus contamination
* (SKY_BRIGHTNESS)
*
*------------- Polarization measurements (30-353GHz only) ----------
real*8 P_5 ! (mJy) ?=- Polarization flux density of the sources
* as determined by a matched filter (P) (G3)
real*8 e_P_5 ! (mJy) ?=- Uncertainty (1 sigma) in P (P_ERR)
real*8 AnP_5 ! (deg) ?=- Orientation of polarization with respect
* to NGP (ANGLE_P) (G2) (G3)
real*8 e_AnP_5 ! (deg) ?=- Uncertainty (1 sigma) in AnP
* (ANGLE_P_ERR) (G3)
real*8 ApP_5 ! (mJy) ?=- Orientation of polarization with respect
* to NGP (APER_P) (G2) (G3)
real*8 e_ApP_5 ! (mJy) ?=- Uncertainty (1 sigma) in AperP
* (APER_P_ERR) (G3)
real*8 ApAnP_5 ! (deg) ?=- Orientation of polarization with
* respect to NGP (APER_ANGLE_P) (G2) (G3)
real*8 e_ApAnP_5 ! (deg) ?=- Uncertainty (1 sigma) in ApAneP
* (APER_ANGLE_P_ERR) (G3)
real*8 PUL_5 ! (mJy) ?=- Polarization flux density 99.99% upper
* limit (P_UPPER_LIMIT) (G9)
real*8 ApPUL_5 ! (mJy) ?=- Polarization flux density 99.99% upper
* limit (APER_P_UPPER_LIMIT) (G9)
*
*---------- Marginal polarization measurements (100-353GHz only) --------
integer*4 PStat_5 ! ?=- Polarization detection status
* (P_STAT) (G4)
real*8 PX_5 ! (mJy) ?=- Polarization flux density of the sources
* as determined by a matched filter using
* Bayesian polarization estimator (PX)
real*8 e_PX_10 ! (mJy) ?=- PX uncertainty; lower 95% error bar
* (PX_ERR_LOWER)
real*8 E_PX_11 ! (mJy) ?=- PX uncertainty; upper 95% error bar
* (PX_ERR_UPPER)
real*8 AnPX_5 ! (deg) ?=- Orientation of polarization with respect
* to NGP using Bayesian polarization
* estimator (ANGLE_PX) (G2)
real*8 e_AnPX_10 ! (deg) ?=- AnglePX uncertainty; lower 95% error bar
* (ANGLE_PX_ERR_LOWER)
real*8 E_AnPX_11 ! (deg) ?=- AnglePX uncertainty; upper 95% error bar
* (ANGLE_PX_ERR_UPPER)
*
*------------------Fluxes at 217, 353 and 545 GHz (857GHz only)------------------
real*8 ApFlux217_5 ! (mJy) ?=- Flux density at 217GHz of source as
* determined from the aperture photometry
* (APERFLUX_227)
real*8 e_ApFlux217_5 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux217
* (APERFLUX_ERR_227)
real*8 ApFlux353_5 ! (mJy) ?=- Flux density at 353GHz of source as
* determined from the aperture photometry
* (APERFLUX_353)
real*8 e_ApFlux353_5 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux353
* (APERFLUX_ERR_353)
real*8 ApFlux545_5 ! (mJy) ?=- Flux density at 545GHz of source as
* determined from the aperture photometry
* (APERFLUX_545)
real*8 e_ApFlux545_5 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux545
* (APERFLUX_ERR_545)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'pccs353.dat' ! Second Catalogue of Compact Sources, 353GHz
integer*4 nr__6
parameter (nr__6=1344) ! Number of records
character*573 ar__6 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*23 Name_6 ! Source name, PCCS2 FFF GLLL.ll+BB.bbb
* (NAME) (G1)
real*8 GLON_6 ! (deg) Galactic longitude based on extraction
* algorithm (GLON)
real*8 GLAT_6 ! (deg) Galactic latitude based on extraction
* algorithm (GLAT)
real*8 RAdeg_6 ! (deg) Right ascension (J2000) transformed
* from (GLON,GLAT) (RA)
real*8 DEdeg_6 ! (deg) Declination (J2000) transformed
* from (GLON,GLAT) (DEC)
real*4 DetFlux_6 ! (mJy) Flux density of source as determined by
* detection method (DETFLUX)
real*4 e_DetFlux_6 ! (mJy) Uncertainty (1 sigma) in DetFlux
* (DETFLUX_ERR)
real*4 AperFlux_6 ! (mJy) Flux density of source as determined from
* the aperture photometry (APERFLUX)
real*4 e_AperFlux_6 ! (mJy) Uncertainty (1 sigma) in AperFlux
* (APERFLUX_ERR)
real*4 PSFFlux_6 ! (mJy) Flux density of source as determined from
* PSF fitting (PSFFLUX)
real*4 e_PSFFlux_6 ! (mJy) Uncertainty (1 sigma) in PSFFlux
* (PSFFLUX_ERR)
real*4 GauFlux_6 ! (mJy) Flux density of source as determined from
* 2-D Gaussian fitting (GAUFLUX)
real*4 e_GauFlux_6 ! (mJy) ?=- Uncertainty (1 sigma) in GauFlux
* (GAUFLUX_ERR)
real*8 GauSemi1_6 ! (arcmin) Gaussian fit along axis 1 (FWHM)
* (GAU_SEMI1) (G2)
real*8 e_GauSemi1_6 ! (arcmin) ?=- Uncertainty (1 sigma) in GauSemi1
* (GAU_SEMI1_ERR)
real*8 GauSemi2_6 ! (arcmin) Gaussian fit along axis 2 (FWHM)
* (GAU_SEMI2) (G2)
real*8 e_GauSemi2_6 ! (arcmin) ?=- Uncertainty (1 sigma) in GauSemi2
* (GAU_SEMI2_ERR)
real*4 GauTheta_6 ! (deg) Gaussian fit orientation angle
* (GAU_THETA) (G2)
real*4 e_GauTheta_6 ! (deg) ?=- Uncertainty (1 sigma) in GauTheta
* (GAU_THETA_ERR)
real*8 GauFWHMEff_6 ! (arcmin) Gaussian fit effective FWHM (GAU_FWHM_EFF)
integer*4 Extended_6 ! ?=-1 Extended source flag (EXTENDED) (G5)
integer*4 ExtVal_6 ! ?=-1 External validation flag (EXT_VAL) (G6)
character*24 ERCSC_6 ! Name of the ERCSC counterpart, if any
* (ERCSC)
character*23 PCCS1_6 ! Name of the PCCS counterpart, if any (PCCS1)
integer*4 HRel_6 ! ?=-1 Highest reliability catalogue to which
* the source belong
* (HIGHEST_RELIABILITY_CAT) (G7)
*
*---------- only for 100, 143, 217 and 353GHz -----
integer*4 CirrusN_6 ! ?=- Number of sources (S/N>5) detected at
* 857GHz within a 1-degree radius (CIRRUS_N)
real*8 SkyBri_6 ! (MJy/sr) ?=- The mean 857GHz brightness within a
* 2{deg} radius. This may be used as another
* indicator of cirrus contamination
* (SKY_BRIGHTNESS)
*
*------------- Polarization measurements (30-353GHz only) ----------
real*8 P_6 ! (mJy) ?=- Polarization flux density of the sources
* as determined by a matched filter (P) (G3)
real*8 e_P_6 ! (mJy) ?=- Uncertainty (1 sigma) in P (P_ERR)
real*8 AnP_6 ! (deg) ?=- Orientation of polarization with respect
* to NGP (ANGLE_P) (G2) (G3)
real*8 e_AnP_6 ! (deg) ?=- Uncertainty (1 sigma) in AnP
* (ANGLE_P_ERR) (G3)
real*8 ApP_6 ! (mJy) ?=- Orientation of polarization with respect
* to NGP (APER_P) (G2) (G3)
real*8 e_ApP_6 ! (mJy) ?=- Uncertainty (1 sigma) in AperP
* (APER_P_ERR) (G3)
real*8 ApAnP_6 ! (deg) ?=- Orientation of polarization with
* respect to NGP (APER_ANGLE_P) (G2) (G3)
real*8 e_ApAnP_6 ! (deg) ?=- Uncertainty (1 sigma) in ApAneP
* (APER_ANGLE_P_ERR) (G3)
real*8 PUL_6 ! (mJy) ?=- Polarization flux density 99.99% upper
* limit (P_UPPER_LIMIT) (G9)
real*8 ApPUL_6 ! (mJy) ?=- Polarization flux density 99.99% upper
* limit (APER_P_UPPER_LIMIT) (G9)
*
*---------- Marginal polarization measurements (100-353GHz only) --------
integer*4 PStat_6 ! ?=- Polarization detection status
* (P_STAT) (G4)
real*8 PX_6 ! (mJy) ?=- Polarization flux density of the sources
* as determined by a matched filter using
* Bayesian polarization estimator (PX)
real*8 e_PX_12 ! (mJy) ?=- PX uncertainty; lower 95% error bar
* (PX_ERR_LOWER)
real*8 E_PX_13 ! (mJy) ?=- PX uncertainty; upper 95% error bar
* (PX_ERR_UPPER)
real*8 AnPX_6 ! (deg) ?=- Orientation of polarization with respect
* to NGP using Bayesian polarization
* estimator (ANGLE_PX) (G2)
real*8 e_AnPX_12 ! (deg) ?=- AnglePX uncertainty; lower 95% error bar
* (ANGLE_PX_ERR_LOWER)
real*8 E_AnPX_13 ! (deg) ?=- AnglePX uncertainty; upper 95% error bar
* (ANGLE_PX_ERR_UPPER)
*
*------------------Fluxes at 217, 353 and 545 GHz (857GHz only)------------------
real*8 ApFlux217_6 ! (mJy) ?=- Flux density at 217GHz of source as
* determined from the aperture photometry
* (APERFLUX_227)
real*8 e_ApFlux217_6 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux217
* (APERFLUX_ERR_227)
real*8 ApFlux353_6 ! (mJy) ?=- Flux density at 353GHz of source as
* determined from the aperture photometry
* (APERFLUX_353)
real*8 e_ApFlux353_6 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux353
* (APERFLUX_ERR_353)
real*8 ApFlux545_6 ! (mJy) ?=- Flux density at 545GHz of source as
* determined from the aperture photometry
* (APERFLUX_545)
real*8 e_ApFlux545_6 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux545
* (APERFLUX_ERR_545)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'pccs545.dat' ! Second Catalogue of Compact Sources, 545GHz
integer*4 nr__7
parameter (nr__7=1694) ! Number of records
character*573 ar__7 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*23 Name_7 ! Source name, PCCS2 FFF GLLL.ll+BB.bbb
* (NAME) (G1)
real*8 GLON_7 ! (deg) Galactic longitude based on extraction
* algorithm (GLON)
real*8 GLAT_7 ! (deg) Galactic latitude based on extraction
* algorithm (GLAT)
real*8 RAdeg_7 ! (deg) Right ascension (J2000) transformed
* from (GLON,GLAT) (RA)
real*8 DEdeg_7 ! (deg) Declination (J2000) transformed
* from (GLON,GLAT) (DEC)
real*4 DetFlux_7 ! (mJy) Flux density of source as determined by
* detection method (DETFLUX)
real*4 e_DetFlux_7 ! (mJy) Uncertainty (1 sigma) in DetFlux
* (DETFLUX_ERR)
real*4 AperFlux_7 ! (mJy) Flux density of source as determined from
* the aperture photometry (APERFLUX)
real*4 e_AperFlux_7 ! (mJy) Uncertainty (1 sigma) in AperFlux
* (APERFLUX_ERR)
real*4 PSFFlux_7 ! (mJy) Flux density of source as determined from
* PSF fitting (PSFFLUX)
real*4 e_PSFFlux_7 ! (mJy) Uncertainty (1 sigma) in PSFFlux
* (PSFFLUX_ERR)
real*4 GauFlux_7 ! (mJy) Flux density of source as determined from
* 2-D Gaussian fitting (GAUFLUX)
real*4 e_GauFlux_7 ! (mJy) ?=- Uncertainty (1 sigma) in GauFlux
* (GAUFLUX_ERR)
real*8 GauSemi1_7 ! (arcmin) Gaussian fit along axis 1 (FWHM)
* (GAU_SEMI1) (G2)
real*8 e_GauSemi1_7 ! (arcmin) ?=- Uncertainty (1 sigma) in GauSemi1
* (GAU_SEMI1_ERR)
real*8 GauSemi2_7 ! (arcmin) Gaussian fit along axis 2 (FWHM)
* (GAU_SEMI2) (G2)
real*8 e_GauSemi2_7 ! (arcmin) ?=- Uncertainty (1 sigma) in GauSemi2
* (GAU_SEMI2_ERR)
real*4 GauTheta_7 ! (deg) Gaussian fit orientation angle
* (GAU_THETA) (G2)
real*4 e_GauTheta_7 ! (deg) ?=- Uncertainty (1 sigma) in GauTheta
* (GAU_THETA_ERR)
real*8 GauFWHMEff_7 ! (arcmin) Gaussian fit effective FWHM (GAU_FWHM_EFF)
integer*4 Extended_7 ! ?=-1 Extended source flag (EXTENDED) (G5)
integer*4 ExtVal_7 ! ?=-1 External validation flag (EXT_VAL) (G6)
character*24 ERCSC_7 ! Name of the ERCSC counterpart, if any
* (ERCSC)
character*23 PCCS1_7 ! Name of the PCCS counterpart, if any (PCCS1)
integer*4 HRel_7 ! ?=-1 Highest reliability catalogue to which
* the source belong
* (HIGHEST_RELIABILITY_CAT) (G7)
*
*---------- only for 100, 143, 217 and 353GHz -----
integer*4 CirrusN_7 ! ?=- Number of sources (S/N>5) detected at
* 857GHz within a 1-degree radius (CIRRUS_N)
real*8 SkyBri_7 ! (MJy/sr) ?=- The mean 857GHz brightness within a
* 2{deg} radius. This may be used as another
* indicator of cirrus contamination
* (SKY_BRIGHTNESS)
*
*------------- Polarization measurements (30-353GHz only) ----------
real*8 P_7 ! (mJy) ?=- Polarization flux density of the sources
* as determined by a matched filter (P) (G3)
real*8 e_P_7 ! (mJy) ?=- Uncertainty (1 sigma) in P (P_ERR)
real*8 AnP_7 ! (deg) ?=- Orientation of polarization with respect
* to NGP (ANGLE_P) (G2) (G3)
real*8 e_AnP_7 ! (deg) ?=- Uncertainty (1 sigma) in AnP
* (ANGLE_P_ERR) (G3)
real*8 ApP_7 ! (mJy) ?=- Orientation of polarization with respect
* to NGP (APER_P) (G2) (G3)
real*8 e_ApP_7 ! (mJy) ?=- Uncertainty (1 sigma) in AperP
* (APER_P_ERR) (G3)
real*8 ApAnP_7 ! (deg) ?=- Orientation of polarization with
* respect to NGP (APER_ANGLE_P) (G2) (G3)
real*8 e_ApAnP_7 ! (deg) ?=- Uncertainty (1 sigma) in ApAneP
* (APER_ANGLE_P_ERR) (G3)
real*8 PUL_7 ! (mJy) ?=- Polarization flux density 99.99% upper
* limit (P_UPPER_LIMIT) (G9)
real*8 ApPUL_7 ! (mJy) ?=- Polarization flux density 99.99% upper
* limit (APER_P_UPPER_LIMIT) (G9)
*
*---------- Marginal polarization measurements (100-353GHz only) --------
integer*4 PStat_7 ! ?=- Polarization detection status
* (P_STAT) (G4)
real*8 PX_7 ! (mJy) ?=- Polarization flux density of the sources
* as determined by a matched filter using
* Bayesian polarization estimator (PX)
real*8 e_PX_14 ! (mJy) ?=- PX uncertainty; lower 95% error bar
* (PX_ERR_LOWER)
real*8 E_PX_15 ! (mJy) ?=- PX uncertainty; upper 95% error bar
* (PX_ERR_UPPER)
real*8 AnPX_7 ! (deg) ?=- Orientation of polarization with respect
* to NGP using Bayesian polarization
* estimator (ANGLE_PX) (G2)
real*8 e_AnPX_14 ! (deg) ?=- AnglePX uncertainty; lower 95% error bar
* (ANGLE_PX_ERR_LOWER)
real*8 E_AnPX_15 ! (deg) ?=- AnglePX uncertainty; upper 95% error bar
* (ANGLE_PX_ERR_UPPER)
*
*------------------Fluxes at 217, 353 and 545 GHz (857GHz only)------------------
real*8 ApFlux217_7 ! (mJy) ?=- Flux density at 217GHz of source as
* determined from the aperture photometry
* (APERFLUX_227)
real*8 e_ApFlux217_7 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux217
* (APERFLUX_ERR_227)
real*8 ApFlux353_7 ! (mJy) ?=- Flux density at 353GHz of source as
* determined from the aperture photometry
* (APERFLUX_353)
real*8 e_ApFlux353_7 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux353
* (APERFLUX_ERR_353)
real*8 ApFlux545_7 ! (mJy) ?=- Flux density at 545GHz of source as
* determined from the aperture photometry
* (APERFLUX_545)
real*8 e_ApFlux545_7 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux545
* (APERFLUX_ERR_545)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'pccs857.dat' ! Second Catalogue of Compact Sources, 857GHz
integer*4 nr__8
parameter (nr__8=4891) ! Number of records
character*573 ar__8 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*23 Name_8 ! Source name, PCCS2 FFF GLLL.ll+BB.bbb
* (NAME) (G1)
real*8 GLON_8 ! (deg) Galactic longitude based on extraction
* algorithm (GLON)
real*8 GLAT_8 ! (deg) Galactic latitude based on extraction
* algorithm (GLAT)
real*8 RAdeg_8 ! (deg) Right ascension (J2000) transformed
* from (GLON,GLAT) (RA)
real*8 DEdeg_8 ! (deg) Declination (J2000) transformed
* from (GLON,GLAT) (DEC)
real*4 DetFlux_8 ! (mJy) Flux density of source as determined by
* detection method (DETFLUX)
real*4 e_DetFlux_8 ! (mJy) Uncertainty (1 sigma) in DetFlux
* (DETFLUX_ERR)
real*4 AperFlux_8 ! (mJy) Flux density of source as determined from
* the aperture photometry (APERFLUX)
real*4 e_AperFlux_8 ! (mJy) Uncertainty (1 sigma) in AperFlux
* (APERFLUX_ERR)
real*4 PSFFlux_8 ! (mJy) Flux density of source as determined from
* PSF fitting (PSFFLUX)
real*4 e_PSFFlux_8 ! (mJy) Uncertainty (1 sigma) in PSFFlux
* (PSFFLUX_ERR)
real*4 GauFlux_8 ! (mJy) Flux density of source as determined from
* 2-D Gaussian fitting (GAUFLUX)
real*4 e_GauFlux_8 ! (mJy) ?=- Uncertainty (1 sigma) in GauFlux
* (GAUFLUX_ERR)
real*8 GauSemi1_8 ! (arcmin) Gaussian fit along axis 1 (FWHM)
* (GAU_SEMI1) (G2)
real*8 e_GauSemi1_8 ! (arcmin) ?=- Uncertainty (1 sigma) in GauSemi1
* (GAU_SEMI1_ERR)
real*8 GauSemi2_8 ! (arcmin) Gaussian fit along axis 2 (FWHM)
* (GAU_SEMI2) (G2)
real*8 e_GauSemi2_8 ! (arcmin) ?=- Uncertainty (1 sigma) in GauSemi2
* (GAU_SEMI2_ERR)
real*4 GauTheta_8 ! (deg) Gaussian fit orientation angle
* (GAU_THETA) (G2)
real*4 e_GauTheta_8 ! (deg) ?=- Uncertainty (1 sigma) in GauTheta
* (GAU_THETA_ERR)
real*8 GauFWHMEff_8 ! (arcmin) Gaussian fit effective FWHM (GAU_FWHM_EFF)
integer*4 Extended_8 ! ?=-1 Extended source flag (EXTENDED) (G5)
integer*4 ExtVal_8 ! ?=-1 External validation flag (EXT_VAL) (G6)
character*24 ERCSC_8 ! Name of the ERCSC counterpart, if any
* (ERCSC)
character*23 PCCS1_8 ! Name of the PCCS counterpart, if any (PCCS1)
integer*4 HRel_8 ! ?=-1 Highest reliability catalogue to which
* the source belong
* (HIGHEST_RELIABILITY_CAT) (G7)
*
*---------- only for 100, 143, 217 and 353GHz -----
integer*4 CirrusN_8 ! ?=- Number of sources (S/N>5) detected at
* 857GHz within a 1-degree radius (CIRRUS_N)
real*8 SkyBri_8 ! (MJy/sr) ?=- The mean 857GHz brightness within a
* 2{deg} radius. This may be used as another
* indicator of cirrus contamination
* (SKY_BRIGHTNESS)
*
*------------- Polarization measurements (30-353GHz only) ----------
real*8 P_8 ! (mJy) ?=- Polarization flux density of the sources
* as determined by a matched filter (P) (G3)
real*8 e_P_8 ! (mJy) ?=- Uncertainty (1 sigma) in P (P_ERR)
real*8 AnP_8 ! (deg) ?=- Orientation of polarization with respect
* to NGP (ANGLE_P) (G2) (G3)
real*8 e_AnP_8 ! (deg) ?=- Uncertainty (1 sigma) in AnP
* (ANGLE_P_ERR) (G3)
real*8 ApP_8 ! (mJy) ?=- Orientation of polarization with respect
* to NGP (APER_P) (G2) (G3)
real*8 e_ApP_8 ! (mJy) ?=- Uncertainty (1 sigma) in AperP
* (APER_P_ERR) (G3)
real*8 ApAnP_8 ! (deg) ?=- Orientation of polarization with
* respect to NGP (APER_ANGLE_P) (G2) (G3)
real*8 e_ApAnP_8 ! (deg) ?=- Uncertainty (1 sigma) in ApAneP
* (APER_ANGLE_P_ERR) (G3)
real*8 PUL_8 ! (mJy) ?=- Polarization flux density 99.99% upper
* limit (P_UPPER_LIMIT) (G9)
real*8 ApPUL_8 ! (mJy) ?=- Polarization flux density 99.99% upper
* limit (APER_P_UPPER_LIMIT) (G9)
*
*---------- Marginal polarization measurements (100-353GHz only) --------
integer*4 PStat_8 ! ?=- Polarization detection status
* (P_STAT) (G4)
real*8 PX_8 ! (mJy) ?=- Polarization flux density of the sources
* as determined by a matched filter using
* Bayesian polarization estimator (PX)
real*8 e_PX_16 ! (mJy) ?=- PX uncertainty; lower 95% error bar
* (PX_ERR_LOWER)
real*8 E_PX_17 ! (mJy) ?=- PX uncertainty; upper 95% error bar
* (PX_ERR_UPPER)
real*8 AnPX_8 ! (deg) ?=- Orientation of polarization with respect
* to NGP using Bayesian polarization
* estimator (ANGLE_PX) (G2)
real*8 e_AnPX_16 ! (deg) ?=- AnglePX uncertainty; lower 95% error bar
* (ANGLE_PX_ERR_LOWER)
real*8 E_AnPX_17 ! (deg) ?=- AnglePX uncertainty; upper 95% error bar
* (ANGLE_PX_ERR_UPPER)
*
*------------------Fluxes at 217, 353 and 545 GHz (857GHz only)------------------
real*8 ApFlux217_8 ! (mJy) ?=- Flux density at 217GHz of source as
* determined from the aperture photometry
* (APERFLUX_227)
real*8 e_ApFlux217_8 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux217
* (APERFLUX_ERR_227)
real*8 ApFlux353_8 ! (mJy) ?=- Flux density at 353GHz of source as
* determined from the aperture photometry
* (APERFLUX_353)
real*8 e_ApFlux353_8 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux353
* (APERFLUX_ERR_353)
real*8 ApFlux545_8 ! (mJy) ?=- Flux density at 545GHz of source as
* determined from the aperture photometry
* (APERFLUX_545)
real*8 e_ApFlux545_8 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux545
* (APERFLUX_ERR_545)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'pccs100e.dat' ! Second Catalogue of Compact Sources, 100GHz,
unknown reliability
integer*4 nr__9
parameter (nr__9=2487) ! Number of records
character*586 ar__9 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*24 Name_9 ! Source name, PCCS2E FFF GLLL.ll+BB.bbb
* (NAME) (G1)
real*8 GLON_9 ! (deg) Galactic longitude based on extraction
* algorithm (GLON)
real*8 GLAT_9 ! (deg) Galactic latitude based on extraction
* algorithm (GLAT)
real*8 RAdeg_9 ! (deg) Right ascension (J2000) transformed
* from (GLON,GLAT) (RA)
real*8 DEdeg_9 ! (deg) Declination (J2000) transformed
* from (GLON,GLAT) (DEC)
real*4 DetFlux_9 ! (mJy) Flux density of source as determined by
* detection method (DETFLUX)
real*4 e_DetFlux_9 ! (mJy) Uncertainty (1 sigma) in DetFlux
* (DETFLUX_ERR)
real*4 AperFlux_9 ! (mJy) Flux density of source as determined from
* the aperture photometry (APERFLUX)
real*4 e_AperFlux_9 ! (mJy) Uncertainty (1 sigma) in AperFlux
* (APERFLUX_ERR)
real*4 PSFFlux_9 ! (mJy) Flux density of source as determined from
* PSF fitting (PSFFLUX)
real*4 e_PSFFlux_9 ! (mJy) Uncertainty (1 sigma) in PSFFlux
* (PSFFLUX_ERR)
real*4 GauFlux_9 ! (mJy) Flux density of source as determined from
* 2-D Gaussian fitting (GAUFLUX)
real*4 e_GauFlux_9 ! (mJy) Uncertainty (1 sigma) in GauFlux
* (GAUFLUX_ERR)
real*8 GauSemi1_9 ! (arcmin) Gaussian fit along axis 1 (FWHM)
* (GAU_SEMI1) (G2)
real*8 e_GauSemi1_9 ! (arcmin) ?=- Uncertainty (1 sigma) in GauSemi1
* (GAU_SEMI1_ERR)
real*8 GauSemi2_9 ! (arcmin) Gaussian fit along axis 2 (FWHM)
* (GAU_SEMI2) (G2)
real*8 e_GauSemi2_9 ! (arcmin) ?=- Uncertainty (1 sigma) in GauSemi2
* (GAU_SEMI2_ERR)
real*4 GauTheta_9 ! (deg) Gaussian fit orientation angle
* (GAU_THETA) (G2)
real*4 e_GauTheta_9 ! (deg) Uncertainty (1 sigma) in GauTheta
* (GAU_THETA_ERR)
real*8 GauFWHMEff_9 ! (arcmin) Gaussian fit effective FWHM (GAU_FWHM_EFF)
*
*------------- Polarization measurements (30-353GHz only) ----------
real*8 P_9 ! (mJy) ?=- Polarization flux density of the sources
* as determined by a matched filter (P) (G3)
real*8 e_P_9 ! (mJy) ?=- Uncertainty (1 sigma) in P (P_ERR)
real*8 AnP_9 ! (deg) ?=- Orientation of polarization with respect
* to NGP (ANGLE_P) (G2) (G3)
real*8 e_AnP_9 ! (deg) ?=- Uncertainty (1 sigma) in AnP
* (ANGLE_P_ERR) (G3)
real*8 ApP_9 ! (mJy) ?=- Orientation of polarization with respect
* to NGP (APER_P) (G2) (G3)
real*8 e_ApP_9 ! (mJy) ?=- Uncertainty (1 sigma) in AperP
* (APER_P_ERR) (G3)
real*8 ApAnP_9 ! (deg) ?=- Orientation of polarization with
* respect to NGP (APER_ANGLE_P) (G2) (G3)
real*8 e_ApAnP_9 ! (deg) ?=- Uncertainty (1 sigma) in ApAneP
* (APER_ANGLE_P_ERR) (G3)
real*8 PUL_9 ! (mJy) ?=- Polarization flux density 99.99% upper
* limit (P_UPPER_LIMIT) (G9)
real*8 ApPUL_9 ! (mJy) ?=- Polarization flux density 99.99% upper
* limit (APER_P_UPPER_LIMIT) (G9)
*
*---------- Marginal polarization measurements (100-353GHz only) --------
integer*4 PStat_9 ! ?=- Polarization detection status
* (P_STAT) (G4)
real*8 PX_9 ! (mJy) ?=- Polarization flux density of the sources
* as determined by a matched filter using
* Bayesian polarization estimator (PX)
real*8 e_PX_18 ! (mJy) ?=- PX uncertainty; lower 95% error bar
* (PX_ERR_LOWER)
real*8 E_PX_19 ! (mJy) ?=- PX uncertainty; upper 95% error bar
* (PX_ERR_UPPER)
real*8 AnPX_9 ! (deg) ?=- Orientation of polarization with respect
* to NGP using Bayesian polarization
* estimator (ANGLE_PX) (G2)
real*8 e_AnPX_18 ! (deg) ?=- AnglePX uncertainty; lower 95% error bar
* (ANGLE_PX_ERR_LOWER)
real*8 E_AnPX_19 ! (deg) ?=- AnglePX uncertainty; upper 95% error bar
* (ANGLE_PX_ERR_UPPER)
*--------------------------------------------------------------------------------
integer*4 Extended_9 ! ?=-1 Extended source flag (EXTENDED) (G5)
integer*4 ExtVal_9 ! ?=-1 External validation flag (EXT_VAL) (G6)
character*24 ERCSC_9 ! Name of the ERCSC counterpart, if any
* (ERCSC)
character*23 PCCS1_9 ! Name of the PCCS counterpart, if any (PCCS1)
integer*4 HRel_9 ! ?=- Highest reliability catalogue to which
* the source belong
* (HIGHEST_RELIABILITY_CAT) (G7)
integer*4 CirrusN_9 ! ?=- Number of sources (S/N>5) detected at
* 857GHz within a 1-degree radius (CIRRUS_N)
real*8 SkyBri_9 ! (MJy/sr) ?=- The mean 857GHz brightness within a
* 2{deg} radius. This may be used as another
* indicator of cirrus contamination
* (SKY_BRIGHTNESS)
integer*4 Zone ! [1/3]? Which zone ?
* (pccs100e.dat and pccs857e.dat only)
* (WHICH_ZONE) (G8)
*------------------Fluxes at 217, 353 and 545 GHz (857GHz only)------------------
real*4 ApFlux217_9 ! (mJy) ?=- Flux density at 217GHz of source as
* determined from the aperture photometry
* (APERFLUX_227)
real*4 e_ApFlux217_9 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux217
* (APERFLUX_ERR_227)
real*4 ApFlux353_9 ! (mJy) ?=- Flux density at 353GHz of source as
* determined from the aperture photometry
* (APERFLUX_353)
real*4 e_ApFlux353_9 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux353
* (APERFLUX_ERR_353)
real*4 ApFlux545_9 ! (mJy) ?=- Flux density at 545GHz of source as
* determined from the aperture photometry
* (APERFLUX_545)
real*4 e_ApFlux545_9 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux545
* (APERFLUX_ERR_545)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'pccs143e.dat' ! Second Catalogue of Compact Sources, 143GHz,
unknown reliability
integer*4 nr__10
parameter (nr__10=4139) ! Number of records
character*586 ar__10 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*24 Name_10 ! Source name, PCCS2E FFF GLLL.ll+BB.bbb
* (NAME) (G1)
real*8 GLON_10 ! (deg) Galactic longitude based on extraction
* algorithm (GLON)
real*8 GLAT_10 ! (deg) Galactic latitude based on extraction
* algorithm (GLAT)
real*8 RAdeg_10 ! (deg) Right ascension (J2000) transformed
* from (GLON,GLAT) (RA)
real*8 DEdeg_10 ! (deg) Declination (J2000) transformed
* from (GLON,GLAT) (DEC)
real*4 DetFlux_10 ! (mJy) Flux density of source as determined by
* detection method (DETFLUX)
real*4 e_DetFlux_10 ! (mJy) Uncertainty (1 sigma) in DetFlux
* (DETFLUX_ERR)
real*4 AperFlux_10 ! (mJy) Flux density of source as determined from
* the aperture photometry (APERFLUX)
real*4 e_AperFlux_10 ! (mJy) Uncertainty (1 sigma) in AperFlux
* (APERFLUX_ERR)
real*4 PSFFlux_10 ! (mJy) Flux density of source as determined from
* PSF fitting (PSFFLUX)
real*4 e_PSFFlux_10 ! (mJy) Uncertainty (1 sigma) in PSFFlux
* (PSFFLUX_ERR)
real*4 GauFlux_10 ! (mJy) Flux density of source as determined from
* 2-D Gaussian fitting (GAUFLUX)
real*4 e_GauFlux_10 ! (mJy) Uncertainty (1 sigma) in GauFlux
* (GAUFLUX_ERR)
real*8 GauSemi1_10 ! (arcmin) Gaussian fit along axis 1 (FWHM)
* (GAU_SEMI1) (G2)
real*8 e_GauSemi1_10 ! (arcmin) ?=- Uncertainty (1 sigma) in GauSemi1
* (GAU_SEMI1_ERR)
real*8 GauSemi2_10 ! (arcmin) Gaussian fit along axis 2 (FWHM)
* (GAU_SEMI2) (G2)
real*8 e_GauSemi2_10 ! (arcmin) ?=- Uncertainty (1 sigma) in GauSemi2
* (GAU_SEMI2_ERR)
real*4 GauTheta_10 ! (deg) Gaussian fit orientation angle
* (GAU_THETA) (G2)
real*4 e_GauTheta_10 ! (deg) Uncertainty (1 sigma) in GauTheta
* (GAU_THETA_ERR)
real*8 GauFWHMEff_10 ! (arcmin) Gaussian fit effective FWHM (GAU_FWHM_EFF)
*
*------------- Polarization measurements (30-353GHz only) ----------
real*8 P_10 ! (mJy) ?=- Polarization flux density of the sources
* as determined by a matched filter (P) (G3)
real*8 e_P_10 ! (mJy) ?=- Uncertainty (1 sigma) in P (P_ERR)
real*8 AnP_10 ! (deg) ?=- Orientation of polarization with respect
* to NGP (ANGLE_P) (G2) (G3)
real*8 e_AnP_10 ! (deg) ?=- Uncertainty (1 sigma) in AnP
* (ANGLE_P_ERR) (G3)
real*8 ApP_10 ! (mJy) ?=- Orientation of polarization with respect
* to NGP (APER_P) (G2) (G3)
real*8 e_ApP_10 ! (mJy) ?=- Uncertainty (1 sigma) in AperP
* (APER_P_ERR) (G3)
real*8 ApAnP_10 ! (deg) ?=- Orientation of polarization with
* respect to NGP (APER_ANGLE_P) (G2) (G3)
real*8 e_ApAnP_10 ! (deg) ?=- Uncertainty (1 sigma) in ApAneP
* (APER_ANGLE_P_ERR) (G3)
real*8 PUL_10 ! (mJy) ?=- Polarization flux density 99.99% upper
* limit (P_UPPER_LIMIT) (G9)
real*8 ApPUL_10 ! (mJy) ?=- Polarization flux density 99.99% upper
* limit (APER_P_UPPER_LIMIT) (G9)
*
*---------- Marginal polarization measurements (100-353GHz only) --------
integer*4 PStat_10 ! ?=- Polarization detection status
* (P_STAT) (G4)
real*8 PX_10 ! (mJy) ?=- Polarization flux density of the sources
* as determined by a matched filter using
* Bayesian polarization estimator (PX)
real*8 e_PX_20 ! (mJy) ?=- PX uncertainty; lower 95% error bar
* (PX_ERR_LOWER)
real*8 E_PX_21 ! (mJy) ?=- PX uncertainty; upper 95% error bar
* (PX_ERR_UPPER)
real*8 AnPX_10 ! (deg) ?=- Orientation of polarization with respect
* to NGP using Bayesian polarization
* estimator (ANGLE_PX) (G2)
real*8 e_AnPX_20 ! (deg) ?=- AnglePX uncertainty; lower 95% error bar
* (ANGLE_PX_ERR_LOWER)
real*8 E_AnPX_21 ! (deg) ?=- AnglePX uncertainty; upper 95% error bar
* (ANGLE_PX_ERR_UPPER)
*--------------------------------------------------------------------------------
integer*4 Extended_10 ! ?=-1 Extended source flag (EXTENDED) (G5)
integer*4 ExtVal_10 ! ?=-1 External validation flag (EXT_VAL) (G6)
character*24 ERCSC_10 ! Name of the ERCSC counterpart, if any
* (ERCSC)
character*23 PCCS1_10 ! Name of the PCCS counterpart, if any (PCCS1)
integer*4 HRel_10 ! ?=- Highest reliability catalogue to which
* the source belong
* (HIGHEST_RELIABILITY_CAT) (G7)
integer*4 CirrusN_10 ! ?=- Number of sources (S/N>5) detected at
* 857GHz within a 1-degree radius (CIRRUS_N)
real*8 SkyBri_10 ! (MJy/sr) ?=- The mean 857GHz brightness within a
* 2{deg} radius. This may be used as another
* indicator of cirrus contamination
* (SKY_BRIGHTNESS)
integer*4 Zone_1 ! [1/3]? Which zone ?
* (pccs100e.dat and pccs857e.dat only)
* (WHICH_ZONE) (G8)
*------------------Fluxes at 217, 353 and 545 GHz (857GHz only)------------------
real*4 ApFlux217_10 ! (mJy) ?=- Flux density at 217GHz of source as
* determined from the aperture photometry
* (APERFLUX_227)
real*4 e_ApFlux217_10 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux217
* (APERFLUX_ERR_227)
real*4 ApFlux353_10 ! (mJy) ?=- Flux density at 353GHz of source as
* determined from the aperture photometry
* (APERFLUX_353)
real*4 e_ApFlux353_10 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux353
* (APERFLUX_ERR_353)
real*4 ApFlux545_10 ! (mJy) ?=- Flux density at 545GHz of source as
* determined from the aperture photometry
* (APERFLUX_545)
real*4 e_ApFlux545_10 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux545
* (APERFLUX_ERR_545)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'pccs217e.dat' ! Second Catalogue of Compact Sources, 217GHz,
unknown reliability
integer*4 nr__11
parameter (nr__11=16842) ! Number of records
character*586 ar__11 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*24 Name_11 ! Source name, PCCS2E FFF GLLL.ll+BB.bbb
* (NAME) (G1)
real*8 GLON_11 ! (deg) Galactic longitude based on extraction
* algorithm (GLON)
real*8 GLAT_11 ! (deg) Galactic latitude based on extraction
* algorithm (GLAT)
real*8 RAdeg_11 ! (deg) Right ascension (J2000) transformed
* from (GLON,GLAT) (RA)
real*8 DEdeg_11 ! (deg) Declination (J2000) transformed
* from (GLON,GLAT) (DEC)
real*4 DetFlux_11 ! (mJy) Flux density of source as determined by
* detection method (DETFLUX)
real*4 e_DetFlux_11 ! (mJy) Uncertainty (1 sigma) in DetFlux
* (DETFLUX_ERR)
real*4 AperFlux_11 ! (mJy) Flux density of source as determined from
* the aperture photometry (APERFLUX)
real*4 e_AperFlux_11 ! (mJy) Uncertainty (1 sigma) in AperFlux
* (APERFLUX_ERR)
real*4 PSFFlux_11 ! (mJy) Flux density of source as determined from
* PSF fitting (PSFFLUX)
real*4 e_PSFFlux_11 ! (mJy) Uncertainty (1 sigma) in PSFFlux
* (PSFFLUX_ERR)
real*4 GauFlux_11 ! (mJy) Flux density of source as determined from
* 2-D Gaussian fitting (GAUFLUX)
real*4 e_GauFlux_11 ! (mJy) Uncertainty (1 sigma) in GauFlux
* (GAUFLUX_ERR)
real*8 GauSemi1_11 ! (arcmin) Gaussian fit along axis 1 (FWHM)
* (GAU_SEMI1) (G2)
real*8 e_GauSemi1_11 ! (arcmin) ?=- Uncertainty (1 sigma) in GauSemi1
* (GAU_SEMI1_ERR)
real*8 GauSemi2_11 ! (arcmin) Gaussian fit along axis 2 (FWHM)
* (GAU_SEMI2) (G2)
real*8 e_GauSemi2_11 ! (arcmin) ?=- Uncertainty (1 sigma) in GauSemi2
* (GAU_SEMI2_ERR)
real*4 GauTheta_11 ! (deg) Gaussian fit orientation angle
* (GAU_THETA) (G2)
real*4 e_GauTheta_11 ! (deg) Uncertainty (1 sigma) in GauTheta
* (GAU_THETA_ERR)
real*8 GauFWHMEff_11 ! (arcmin) Gaussian fit effective FWHM (GAU_FWHM_EFF)
*
*------------- Polarization measurements (30-353GHz only) ----------
real*8 P_11 ! (mJy) ?=- Polarization flux density of the sources
* as determined by a matched filter (P) (G3)
real*8 e_P_11 ! (mJy) ?=- Uncertainty (1 sigma) in P (P_ERR)
real*8 AnP_11 ! (deg) ?=- Orientation of polarization with respect
* to NGP (ANGLE_P) (G2) (G3)
real*8 e_AnP_11 ! (deg) ?=- Uncertainty (1 sigma) in AnP
* (ANGLE_P_ERR) (G3)
real*8 ApP_11 ! (mJy) ?=- Orientation of polarization with respect
* to NGP (APER_P) (G2) (G3)
real*8 e_ApP_11 ! (mJy) ?=- Uncertainty (1 sigma) in AperP
* (APER_P_ERR) (G3)
real*8 ApAnP_11 ! (deg) ?=- Orientation of polarization with
* respect to NGP (APER_ANGLE_P) (G2) (G3)
real*8 e_ApAnP_11 ! (deg) ?=- Uncertainty (1 sigma) in ApAneP
* (APER_ANGLE_P_ERR) (G3)
real*8 PUL_11 ! (mJy) ?=- Polarization flux density 99.99% upper
* limit (P_UPPER_LIMIT) (G9)
real*8 ApPUL_11 ! (mJy) ?=- Polarization flux density 99.99% upper
* limit (APER_P_UPPER_LIMIT) (G9)
*
*---------- Marginal polarization measurements (100-353GHz only) --------
integer*4 PStat_11 ! ?=- Polarization detection status
* (P_STAT) (G4)
real*8 PX_11 ! (mJy) ?=- Polarization flux density of the sources
* as determined by a matched filter using
* Bayesian polarization estimator (PX)
real*8 e_PX_22 ! (mJy) ?=- PX uncertainty; lower 95% error bar
* (PX_ERR_LOWER)
real*8 E_PX_23 ! (mJy) ?=- PX uncertainty; upper 95% error bar
* (PX_ERR_UPPER)
real*8 AnPX_11 ! (deg) ?=- Orientation of polarization with respect
* to NGP using Bayesian polarization
* estimator (ANGLE_PX) (G2)
real*8 e_AnPX_22 ! (deg) ?=- AnglePX uncertainty; lower 95% error bar
* (ANGLE_PX_ERR_LOWER)
real*8 E_AnPX_23 ! (deg) ?=- AnglePX uncertainty; upper 95% error bar
* (ANGLE_PX_ERR_UPPER)
*--------------------------------------------------------------------------------
integer*4 Extended_11 ! ?=-1 Extended source flag (EXTENDED) (G5)
integer*4 ExtVal_11 ! ?=-1 External validation flag (EXT_VAL) (G6)
character*24 ERCSC_11 ! Name of the ERCSC counterpart, if any
* (ERCSC)
character*23 PCCS1_11 ! Name of the PCCS counterpart, if any (PCCS1)
integer*4 HRel_11 ! ?=- Highest reliability catalogue to which
* the source belong
* (HIGHEST_RELIABILITY_CAT) (G7)
integer*4 CirrusN_11 ! ?=- Number of sources (S/N>5) detected at
* 857GHz within a 1-degree radius (CIRRUS_N)
real*8 SkyBri_11 ! (MJy/sr) ?=- The mean 857GHz brightness within a
* 2{deg} radius. This may be used as another
* indicator of cirrus contamination
* (SKY_BRIGHTNESS)
integer*4 Zone_2 ! [1/3]? Which zone ?
* (pccs100e.dat and pccs857e.dat only)
* (WHICH_ZONE) (G8)
*------------------Fluxes at 217, 353 and 545 GHz (857GHz only)------------------
real*4 ApFlux217_11 ! (mJy) ?=- Flux density at 217GHz of source as
* determined from the aperture photometry
* (APERFLUX_227)
real*4 e_ApFlux217_11 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux217
* (APERFLUX_ERR_227)
real*4 ApFlux353_11 ! (mJy) ?=- Flux density at 353GHz of source as
* determined from the aperture photometry
* (APERFLUX_353)
real*4 e_ApFlux353_11 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux353
* (APERFLUX_ERR_353)
real*4 ApFlux545_11 ! (mJy) ?=- Flux density at 545GHz of source as
* determined from the aperture photometry
* (APERFLUX_545)
real*4 e_ApFlux545_11 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux545
* (APERFLUX_ERR_545)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'pccs353e.dat' ! Second Catalogue of Compact Sources, 353GHz,
unknown reliability
integer*4 nr__12
parameter (nr__12=22665) ! Number of records
character*586 ar__12 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*24 Name_12 ! Source name, PCCS2E FFF GLLL.ll+BB.bbb
* (NAME) (G1)
real*8 GLON_12 ! (deg) Galactic longitude based on extraction
* algorithm (GLON)
real*8 GLAT_12 ! (deg) Galactic latitude based on extraction
* algorithm (GLAT)
real*8 RAdeg_12 ! (deg) Right ascension (J2000) transformed
* from (GLON,GLAT) (RA)
real*8 DEdeg_12 ! (deg) Declination (J2000) transformed
* from (GLON,GLAT) (DEC)
real*4 DetFlux_12 ! (mJy) Flux density of source as determined by
* detection method (DETFLUX)
real*4 e_DetFlux_12 ! (mJy) Uncertainty (1 sigma) in DetFlux
* (DETFLUX_ERR)
real*4 AperFlux_12 ! (mJy) Flux density of source as determined from
* the aperture photometry (APERFLUX)
real*4 e_AperFlux_12 ! (mJy) Uncertainty (1 sigma) in AperFlux
* (APERFLUX_ERR)
real*4 PSFFlux_12 ! (mJy) Flux density of source as determined from
* PSF fitting (PSFFLUX)
real*4 e_PSFFlux_12 ! (mJy) Uncertainty (1 sigma) in PSFFlux
* (PSFFLUX_ERR)
real*4 GauFlux_12 ! (mJy) Flux density of source as determined from
* 2-D Gaussian fitting (GAUFLUX)
real*4 e_GauFlux_12 ! (mJy) Uncertainty (1 sigma) in GauFlux
* (GAUFLUX_ERR)
real*8 GauSemi1_12 ! (arcmin) Gaussian fit along axis 1 (FWHM)
* (GAU_SEMI1) (G2)
real*8 e_GauSemi1_12 ! (arcmin) ?=- Uncertainty (1 sigma) in GauSemi1
* (GAU_SEMI1_ERR)
real*8 GauSemi2_12 ! (arcmin) Gaussian fit along axis 2 (FWHM)
* (GAU_SEMI2) (G2)
real*8 e_GauSemi2_12 ! (arcmin) ?=- Uncertainty (1 sigma) in GauSemi2
* (GAU_SEMI2_ERR)
real*4 GauTheta_12 ! (deg) Gaussian fit orientation angle
* (GAU_THETA) (G2)
real*4 e_GauTheta_12 ! (deg) Uncertainty (1 sigma) in GauTheta
* (GAU_THETA_ERR)
real*8 GauFWHMEff_12 ! (arcmin) Gaussian fit effective FWHM (GAU_FWHM_EFF)
*
*------------- Polarization measurements (30-353GHz only) ----------
real*8 P_12 ! (mJy) ?=- Polarization flux density of the sources
* as determined by a matched filter (P) (G3)
real*8 e_P_12 ! (mJy) ?=- Uncertainty (1 sigma) in P (P_ERR)
real*8 AnP_12 ! (deg) ?=- Orientation of polarization with respect
* to NGP (ANGLE_P) (G2) (G3)
real*8 e_AnP_12 ! (deg) ?=- Uncertainty (1 sigma) in AnP
* (ANGLE_P_ERR) (G3)
real*8 ApP_12 ! (mJy) ?=- Orientation of polarization with respect
* to NGP (APER_P) (G2) (G3)
real*8 e_ApP_12 ! (mJy) ?=- Uncertainty (1 sigma) in AperP
* (APER_P_ERR) (G3)
real*8 ApAnP_12 ! (deg) ?=- Orientation of polarization with
* respect to NGP (APER_ANGLE_P) (G2) (G3)
real*8 e_ApAnP_12 ! (deg) ?=- Uncertainty (1 sigma) in ApAneP
* (APER_ANGLE_P_ERR) (G3)
real*8 PUL_12 ! (mJy) ?=- Polarization flux density 99.99% upper
* limit (P_UPPER_LIMIT) (G9)
real*8 ApPUL_12 ! (mJy) ?=- Polarization flux density 99.99% upper
* limit (APER_P_UPPER_LIMIT) (G9)
*
*---------- Marginal polarization measurements (100-353GHz only) --------
integer*4 PStat_12 ! ?=- Polarization detection status
* (P_STAT) (G4)
real*8 PX_12 ! (mJy) ?=- Polarization flux density of the sources
* as determined by a matched filter using
* Bayesian polarization estimator (PX)
real*8 e_PX_24 ! (mJy) ?=- PX uncertainty; lower 95% error bar
* (PX_ERR_LOWER)
real*8 E_PX_25 ! (mJy) ?=- PX uncertainty; upper 95% error bar
* (PX_ERR_UPPER)
real*8 AnPX_12 ! (deg) ?=- Orientation of polarization with respect
* to NGP using Bayesian polarization
* estimator (ANGLE_PX) (G2)
real*8 e_AnPX_24 ! (deg) ?=- AnglePX uncertainty; lower 95% error bar
* (ANGLE_PX_ERR_LOWER)
real*8 E_AnPX_25 ! (deg) ?=- AnglePX uncertainty; upper 95% error bar
* (ANGLE_PX_ERR_UPPER)
*--------------------------------------------------------------------------------
integer*4 Extended_12 ! ?=-1 Extended source flag (EXTENDED) (G5)
integer*4 ExtVal_12 ! ?=-1 External validation flag (EXT_VAL) (G6)
character*24 ERCSC_12 ! Name of the ERCSC counterpart, if any
* (ERCSC)
character*23 PCCS1_12 ! Name of the PCCS counterpart, if any (PCCS1)
integer*4 HRel_12 ! ?=- Highest reliability catalogue to which
* the source belong
* (HIGHEST_RELIABILITY_CAT) (G7)
integer*4 CirrusN_12 ! ?=- Number of sources (S/N>5) detected at
* 857GHz within a 1-degree radius (CIRRUS_N)
real*8 SkyBri_12 ! (MJy/sr) ?=- The mean 857GHz brightness within a
* 2{deg} radius. This may be used as another
* indicator of cirrus contamination
* (SKY_BRIGHTNESS)
integer*4 Zone_3 ! [1/3]? Which zone ?
* (pccs100e.dat and pccs857e.dat only)
* (WHICH_ZONE) (G8)
*------------------Fluxes at 217, 353 and 545 GHz (857GHz only)------------------
real*4 ApFlux217_12 ! (mJy) ?=- Flux density at 217GHz of source as
* determined from the aperture photometry
* (APERFLUX_227)
real*4 e_ApFlux217_12 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux217
* (APERFLUX_ERR_227)
real*4 ApFlux353_12 ! (mJy) ?=- Flux density at 353GHz of source as
* determined from the aperture photometry
* (APERFLUX_353)
real*4 e_ApFlux353_12 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux353
* (APERFLUX_ERR_353)
real*4 ApFlux545_12 ! (mJy) ?=- Flux density at 545GHz of source as
* determined from the aperture photometry
* (APERFLUX_545)
real*4 e_ApFlux545_12 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux545
* (APERFLUX_ERR_545)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'pccs545e.dat' ! Second Catalogue of Compact Sources, 545GHz,
unknown reliability
integer*4 nr__13
parameter (nr__13=31068) ! Number of records
character*586 ar__13 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*24 Name_13 ! Source name, PCCS2E FFF GLLL.ll+BB.bbb
* (NAME) (G1)
real*8 GLON_13 ! (deg) Galactic longitude based on extraction
* algorithm (GLON)
real*8 GLAT_13 ! (deg) Galactic latitude based on extraction
* algorithm (GLAT)
real*8 RAdeg_13 ! (deg) Right ascension (J2000) transformed
* from (GLON,GLAT) (RA)
real*8 DEdeg_13 ! (deg) Declination (J2000) transformed
* from (GLON,GLAT) (DEC)
real*4 DetFlux_13 ! (mJy) Flux density of source as determined by
* detection method (DETFLUX)
real*4 e_DetFlux_13 ! (mJy) Uncertainty (1 sigma) in DetFlux
* (DETFLUX_ERR)
real*4 AperFlux_13 ! (mJy) Flux density of source as determined from
* the aperture photometry (APERFLUX)
real*4 e_AperFlux_13 ! (mJy) Uncertainty (1 sigma) in AperFlux
* (APERFLUX_ERR)
real*4 PSFFlux_13 ! (mJy) Flux density of source as determined from
* PSF fitting (PSFFLUX)
real*4 e_PSFFlux_13 ! (mJy) Uncertainty (1 sigma) in PSFFlux
* (PSFFLUX_ERR)
real*4 GauFlux_13 ! (mJy) Flux density of source as determined from
* 2-D Gaussian fitting (GAUFLUX)
real*4 e_GauFlux_13 ! (mJy) Uncertainty (1 sigma) in GauFlux
* (GAUFLUX_ERR)
real*8 GauSemi1_13 ! (arcmin) Gaussian fit along axis 1 (FWHM)
* (GAU_SEMI1) (G2)
real*8 e_GauSemi1_13 ! (arcmin) ?=- Uncertainty (1 sigma) in GauSemi1
* (GAU_SEMI1_ERR)
real*8 GauSemi2_13 ! (arcmin) Gaussian fit along axis 2 (FWHM)
* (GAU_SEMI2) (G2)
real*8 e_GauSemi2_13 ! (arcmin) ?=- Uncertainty (1 sigma) in GauSemi2
* (GAU_SEMI2_ERR)
real*4 GauTheta_13 ! (deg) Gaussian fit orientation angle
* (GAU_THETA) (G2)
real*4 e_GauTheta_13 ! (deg) Uncertainty (1 sigma) in GauTheta
* (GAU_THETA_ERR)
real*8 GauFWHMEff_13 ! (arcmin) Gaussian fit effective FWHM (GAU_FWHM_EFF)
*
*------------- Polarization measurements (30-353GHz only) ----------
real*8 P_13 ! (mJy) ?=- Polarization flux density of the sources
* as determined by a matched filter (P) (G3)
real*8 e_P_13 ! (mJy) ?=- Uncertainty (1 sigma) in P (P_ERR)
real*8 AnP_13 ! (deg) ?=- Orientation of polarization with respect
* to NGP (ANGLE_P) (G2) (G3)
real*8 e_AnP_13 ! (deg) ?=- Uncertainty (1 sigma) in AnP
* (ANGLE_P_ERR) (G3)
real*8 ApP_13 ! (mJy) ?=- Orientation of polarization with respect
* to NGP (APER_P) (G2) (G3)
real*8 e_ApP_13 ! (mJy) ?=- Uncertainty (1 sigma) in AperP
* (APER_P_ERR) (G3)
real*8 ApAnP_13 ! (deg) ?=- Orientation of polarization with
* respect to NGP (APER_ANGLE_P) (G2) (G3)
real*8 e_ApAnP_13 ! (deg) ?=- Uncertainty (1 sigma) in ApAneP
* (APER_ANGLE_P_ERR) (G3)
real*8 PUL_13 ! (mJy) ?=- Polarization flux density 99.99% upper
* limit (P_UPPER_LIMIT) (G9)
real*8 ApPUL_13 ! (mJy) ?=- Polarization flux density 99.99% upper
* limit (APER_P_UPPER_LIMIT) (G9)
*
*---------- Marginal polarization measurements (100-353GHz only) --------
integer*4 PStat_13 ! ?=- Polarization detection status
* (P_STAT) (G4)
real*8 PX_13 ! (mJy) ?=- Polarization flux density of the sources
* as determined by a matched filter using
* Bayesian polarization estimator (PX)
real*8 e_PX_26 ! (mJy) ?=- PX uncertainty; lower 95% error bar
* (PX_ERR_LOWER)
real*8 E_PX_27 ! (mJy) ?=- PX uncertainty; upper 95% error bar
* (PX_ERR_UPPER)
real*8 AnPX_13 ! (deg) ?=- Orientation of polarization with respect
* to NGP using Bayesian polarization
* estimator (ANGLE_PX) (G2)
real*8 e_AnPX_26 ! (deg) ?=- AnglePX uncertainty; lower 95% error bar
* (ANGLE_PX_ERR_LOWER)
real*8 E_AnPX_27 ! (deg) ?=- AnglePX uncertainty; upper 95% error bar
* (ANGLE_PX_ERR_UPPER)
*--------------------------------------------------------------------------------
integer*4 Extended_13 ! ?=-1 Extended source flag (EXTENDED) (G5)
integer*4 ExtVal_13 ! ?=-1 External validation flag (EXT_VAL) (G6)
character*24 ERCSC_13 ! Name of the ERCSC counterpart, if any
* (ERCSC)
character*23 PCCS1_13 ! Name of the PCCS counterpart, if any (PCCS1)
integer*4 HRel_13 ! ?=- Highest reliability catalogue to which
* the source belong
* (HIGHEST_RELIABILITY_CAT) (G7)
integer*4 CirrusN_13 ! ?=- Number of sources (S/N>5) detected at
* 857GHz within a 1-degree radius (CIRRUS_N)
real*8 SkyBri_13 ! (MJy/sr) ?=- The mean 857GHz brightness within a
* 2{deg} radius. This may be used as another
* indicator of cirrus contamination
* (SKY_BRIGHTNESS)
integer*4 Zone_4 ! [1/3]? Which zone ?
* (pccs100e.dat and pccs857e.dat only)
* (WHICH_ZONE) (G8)
*------------------Fluxes at 217, 353 and 545 GHz (857GHz only)------------------
real*4 ApFlux217_13 ! (mJy) ?=- Flux density at 217GHz of source as
* determined from the aperture photometry
* (APERFLUX_227)
real*4 e_ApFlux217_13 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux217
* (APERFLUX_ERR_227)
real*4 ApFlux353_13 ! (mJy) ?=- Flux density at 353GHz of source as
* determined from the aperture photometry
* (APERFLUX_353)
real*4 e_ApFlux353_13 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux353
* (APERFLUX_ERR_353)
real*4 ApFlux545_13 ! (mJy) ?=- Flux density at 545GHz of source as
* determined from the aperture photometry
* (APERFLUX_545)
real*4 e_ApFlux545_13 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux545
* (APERFLUX_ERR_545)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'pccs857e.dat' ! Second Catalogue of Compact Sources, 857GHz,
unknown reliability (correct file, 08-Dec-2023)
integer*4 nr__14
parameter (nr__14=43290) ! Number of records
character*586 ar__14 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*24 Name_14 ! Source name, PCCS2E FFF GLLL.ll+BB.bbb
* (NAME) (G1)
real*8 GLON_14 ! (deg) Galactic longitude based on extraction
* algorithm (GLON)
real*8 GLAT_14 ! (deg) Galactic latitude based on extraction
* algorithm (GLAT)
real*8 RAdeg_14 ! (deg) Right ascension (J2000) transformed
* from (GLON,GLAT) (RA)
real*8 DEdeg_14 ! (deg) Declination (J2000) transformed
* from (GLON,GLAT) (DEC)
real*4 DetFlux_14 ! (mJy) Flux density of source as determined by
* detection method (DETFLUX)
real*4 e_DetFlux_14 ! (mJy) Uncertainty (1 sigma) in DetFlux
* (DETFLUX_ERR)
real*4 AperFlux_14 ! (mJy) Flux density of source as determined from
* the aperture photometry (APERFLUX)
real*4 e_AperFlux_14 ! (mJy) Uncertainty (1 sigma) in AperFlux
* (APERFLUX_ERR)
real*4 PSFFlux_14 ! (mJy) Flux density of source as determined from
* PSF fitting (PSFFLUX)
real*4 e_PSFFlux_14 ! (mJy) Uncertainty (1 sigma) in PSFFlux
* (PSFFLUX_ERR)
real*4 GauFlux_14 ! (mJy) Flux density of source as determined from
* 2-D Gaussian fitting (GAUFLUX)
real*4 e_GauFlux_14 ! (mJy) Uncertainty (1 sigma) in GauFlux
* (GAUFLUX_ERR)
real*8 GauSemi1_14 ! (arcmin) Gaussian fit along axis 1 (FWHM)
* (GAU_SEMI1) (G2)
real*8 e_GauSemi1_14 ! (arcmin) ?=- Uncertainty (1 sigma) in GauSemi1
* (GAU_SEMI1_ERR)
real*8 GauSemi2_14 ! (arcmin) Gaussian fit along axis 2 (FWHM)
* (GAU_SEMI2) (G2)
real*8 e_GauSemi2_14 ! (arcmin) ?=- Uncertainty (1 sigma) in GauSemi2
* (GAU_SEMI2_ERR)
real*4 GauTheta_14 ! (deg) Gaussian fit orientation angle
* (GAU_THETA) (G2)
real*4 e_GauTheta_14 ! (deg) Uncertainty (1 sigma) in GauTheta
* (GAU_THETA_ERR)
real*8 GauFWHMEff_14 ! (arcmin) Gaussian fit effective FWHM (GAU_FWHM_EFF)
*
*------------- Polarization measurements (30-353GHz only) ----------
real*8 P_14 ! (mJy) ?=- Polarization flux density of the sources
* as determined by a matched filter (P) (G3)
real*8 e_P_14 ! (mJy) ?=- Uncertainty (1 sigma) in P (P_ERR)
real*8 AnP_14 ! (deg) ?=- Orientation of polarization with respect
* to NGP (ANGLE_P) (G2) (G3)
real*8 e_AnP_14 ! (deg) ?=- Uncertainty (1 sigma) in AnP
* (ANGLE_P_ERR) (G3)
real*8 ApP_14 ! (mJy) ?=- Orientation of polarization with respect
* to NGP (APER_P) (G2) (G3)
real*8 e_ApP_14 ! (mJy) ?=- Uncertainty (1 sigma) in AperP
* (APER_P_ERR) (G3)
real*8 ApAnP_14 ! (deg) ?=- Orientation of polarization with
* respect to NGP (APER_ANGLE_P) (G2) (G3)
real*8 e_ApAnP_14 ! (deg) ?=- Uncertainty (1 sigma) in ApAneP
* (APER_ANGLE_P_ERR) (G3)
real*8 PUL_14 ! (mJy) ?=- Polarization flux density 99.99% upper
* limit (P_UPPER_LIMIT) (G9)
real*8 ApPUL_14 ! (mJy) ?=- Polarization flux density 99.99% upper
* limit (APER_P_UPPER_LIMIT) (G9)
*
*---------- Marginal polarization measurements (100-353GHz only) --------
integer*4 PStat_14 ! ?=- Polarization detection status
* (P_STAT) (G4)
real*8 PX_14 ! (mJy) ?=- Polarization flux density of the sources
* as determined by a matched filter using
* Bayesian polarization estimator (PX)
real*8 e_PX_28 ! (mJy) ?=- PX uncertainty; lower 95% error bar
* (PX_ERR_LOWER)
real*8 E_PX_29 ! (mJy) ?=- PX uncertainty; upper 95% error bar
* (PX_ERR_UPPER)
real*8 AnPX_14 ! (deg) ?=- Orientation of polarization with respect
* to NGP using Bayesian polarization
* estimator (ANGLE_PX) (G2)
real*8 e_AnPX_28 ! (deg) ?=- AnglePX uncertainty; lower 95% error bar
* (ANGLE_PX_ERR_LOWER)
real*8 E_AnPX_29 ! (deg) ?=- AnglePX uncertainty; upper 95% error bar
* (ANGLE_PX_ERR_UPPER)
*--------------------------------------------------------------------------------
integer*4 Extended_14 ! ?=-1 Extended source flag (EXTENDED) (G5)
integer*4 ExtVal_14 ! ?=-1 External validation flag (EXT_VAL) (G6)
character*24 ERCSC_14 ! Name of the ERCSC counterpart, if any
* (ERCSC)
character*23 PCCS1_14 ! Name of the PCCS counterpart, if any (PCCS1)
integer*4 HRel_14 ! ?=- Highest reliability catalogue to which
* the source belong
* (HIGHEST_RELIABILITY_CAT) (G7)
integer*4 CirrusN_14 ! ?=- Number of sources (S/N>5) detected at
* 857GHz within a 1-degree radius (CIRRUS_N)
real*8 SkyBri_14 ! (MJy/sr) ?=- The mean 857GHz brightness within a
* 2{deg} radius. This may be used as another
* indicator of cirrus contamination
* (SKY_BRIGHTNESS)
integer*4 Zone_5 ! [1/3]? Which zone ?
* (pccs100e.dat and pccs857e.dat only)
* (WHICH_ZONE) (G8)
*------------------Fluxes at 217, 353 and 545 GHz (857GHz only)------------------
real*4 ApFlux217_14 ! (mJy) ?=- Flux density at 217GHz of source as
* determined from the aperture photometry
* (APERFLUX_227)
real*4 e_ApFlux217_14 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux217
* (APERFLUX_ERR_227)
real*4 ApFlux353_14 ! (mJy) ?=- Flux density at 353GHz of source as
* determined from the aperture photometry
* (APERFLUX_353)
real*4 e_ApFlux353_14 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux353
* (APERFLUX_ERR_353)
real*4 ApFlux545_14 ! (mJy) ?=- Flux density at 545GHz of source as
* determined from the aperture photometry
* (APERFLUX_545)
real*4 e_ApFlux545_14 ! (mJy) ?=- Uncertainty (1 sigma) in ApFlux545
* (APERFLUX_ERR_545)
C=============================================================================
C Loading file 'pccs030.dat' ! Second Catalogue of Compact Sources, 30GHz
C Format for file interpretation
1 format(
+ A23,1X,F11.7,1X,F11.7,1X,F11.7,1X,F11.7,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,F10.6,
+ 1X,F11.6,1X,F11.6,1X,F11.6,1X,F5.1,1X,F4.1,1X,F10.6,1X,I1,1X,
+ I1,1X,A24,1X,A23,1X,I2,1X,I2,1X,F10.6,1X,F12.6,1X,F11.6,1X,
+ F10.6,1X,F10.6,1X,F12.6,1X,F12.6,1X,F12.6,1X,F10.6,1X,F11.6,
+ 1X,F10.6,1X,I1,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F9.6,1X,
+ F9.6,F12.6,1X,F10.6,1X,F12.6,1X,F11.6,1X,F12.6,1X,F11.6)
C Effective file loading
open(unit=1,status='old',file=
+'pccs030.dat')
write(6,*) '....Loading file: pccs030.dat'
do i__=1,1560
read(1,'(A573)')ar__
read(ar__,1)
+ Name,GLON,GLAT,RAdeg,DEdeg,DetFlux,e_DetFlux,AperFlux,
+ e_AperFlux,PSFFlux,e_PSFFlux,GauFlux,e_GauFlux,GauSemi1,
+ e_GauSemi1,GauSemi2,e_GauSemi2,GauTheta,e_GauTheta,GauFWHMEff,
+ Extended,ExtVal,ERCSC,PCCS1,HRel,CirrusN,SkyBri,P,e_P,AnP,
+ e_AnP,ApP,e_ApP,ApAnP,e_ApAnP,PUL,ApPUL,PStat,PX,e_PX,E_PX_1,
+ AnPX,e_AnPX,E_AnPX_1,ApFlux217,e_ApFlux217,ApFlux353,
+ e_ApFlux353,ApFlux545,e_ApFlux545
if (idig(ar__(164:175)).EQ.0) e_GauFlux = rNULL__
if (idig(ar__(188:198)).EQ.0) e_GauSemi1 = rNULL__
if (idig(ar__(212:222)).EQ.0) e_GauSemi2 = rNULL__
if (idig(ar__(230:233)).EQ.0) e_GauTheta = rNULL__
if (CirrusN .EQ. 45) CirrusN = iNULL__
if (idig(ar__(305:314)).EQ.0) SkyBri = rNULL__
if (idig(ar__(316:327)).EQ.0) P = rNULL__
if (idig(ar__(329:339)).EQ.0) e_P = rNULL__
if (idig(ar__(341:350)).EQ.0) AnP = rNULL__
if (idig(ar__(352:361)).EQ.0) e_AnP = rNULL__
if (idig(ar__(363:374)).EQ.0) ApP = rNULL__
if (idig(ar__(376:387)).EQ.0) e_ApP = rNULL__
if (idig(ar__(389:400)).EQ.0) ApAnP = rNULL__
if (idig(ar__(402:411)).EQ.0) e_ApAnP = rNULL__
if (idig(ar__(413:423)).EQ.0) PUL = rNULL__
if (idig(ar__(425:434)).EQ.0) ApPUL = rNULL__
if (PStat .EQ. 45) PStat = iNULL__
if (idig(ar__(438:447)).EQ.0) PX = rNULL__
if (idig(ar__(449:458)).EQ.0) e_PX = rNULL__
if (idig(ar__(460:469)).EQ.0) E_PX_1 = rNULL__
if (idig(ar__(471:480)).EQ.0) AnPX = rNULL__
if (idig(ar__(482:490)).EQ.0) e_AnPX = rNULL__
if (idig(ar__(492:500)).EQ.0) E_AnPX_1 = rNULL__
if (idig(ar__(501:512)).EQ.0) ApFlux217 = rNULL__
if (idig(ar__(514:523)).EQ.0) e_ApFlux217 = rNULL__
if (idig(ar__(525:536)).EQ.0) ApFlux353 = rNULL__
if (idig(ar__(538:548)).EQ.0) e_ApFlux353 = rNULL__
if (idig(ar__(550:561)).EQ.0) ApFlux545 = rNULL__
if (idig(ar__(563:573)).EQ.0) e_ApFlux545 = rNULL__
c ..............Just test output...........
write(6,1)
+ Name,GLON,GLAT,RAdeg,DEdeg,DetFlux,e_DetFlux,AperFlux,
+ e_AperFlux,PSFFlux,e_PSFFlux,GauFlux,e_GauFlux,GauSemi1,
+ e_GauSemi1,GauSemi2,e_GauSemi2,GauTheta,e_GauTheta,GauFWHMEff,
+ Extended,ExtVal,ERCSC,PCCS1,HRel,CirrusN,SkyBri,P,e_P,AnP,
+ e_AnP,ApP,e_ApP,ApAnP,e_ApAnP,PUL,ApPUL,PStat,PX,e_PX,E_PX_1,
+ AnPX,e_AnPX,E_AnPX_1,ApFlux217,e_ApFlux217,ApFlux353,
+ e_ApFlux353,ApFlux545,e_ApFlux545
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'pccs044.dat' ! Second Catalogue of Compact Sources, 44GHz
C Format for file interpretation
2 format(
+ A23,1X,F11.7,1X,F11.7,1X,F11.7,1X,F11.7,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,F10.6,
+ 1X,F11.6,1X,F11.6,1X,F11.6,1X,F5.1,1X,F4.1,1X,F10.6,1X,I1,1X,
+ I1,1X,A24,1X,A23,1X,I2,1X,I2,1X,F10.6,1X,F12.6,1X,F11.6,1X,
+ F10.6,1X,F10.6,1X,F12.6,1X,F12.6,1X,F12.6,1X,F10.6,1X,F11.6,
+ 1X,F10.6,1X,I1,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F9.6,1X,
+ F9.6,F12.6,1X,F10.6,1X,F12.6,1X,F11.6,1X,F12.6,1X,F11.6)
C Effective file loading
open(unit=1,status='old',file=
+'pccs044.dat')
write(6,*) '....Loading file: pccs044.dat'
do i__=1,934
read(1,'(A573)')ar__1
read(ar__1,2)
+ Name_1,GLON_1,GLAT_1,RAdeg_1,DEdeg_1,DetFlux_1,e_DetFlux_1,
+ AperFlux_1,e_AperFlux_1,PSFFlux_1,e_PSFFlux_1,GauFlux_1,
+ e_GauFlux_1,GauSemi1_1,e_GauSemi1_1,GauSemi2_1,e_GauSemi2_1,
+ GauTheta_1,e_GauTheta_1,GauFWHMEff_1,Extended_1,ExtVal_1,
+ ERCSC_1,PCCS1_1,HRel_1,CirrusN_1,SkyBri_1,P_1,e_P_1,AnP_1,
+ e_AnP_1,ApP_1,e_ApP_1,ApAnP_1,e_ApAnP_1,PUL_1,ApPUL_1,PStat_1,
+ PX_1,e_PX_2,E_PX_3,AnPX_1,e_AnPX_2,E_AnPX_3,ApFlux217_1,
+ e_ApFlux217_1,ApFlux353_1,e_ApFlux353_1,ApFlux545_1,
+ e_ApFlux545_1
if (idig(ar__1(164:175)).EQ.0) e_GauFlux_1 = rNULL__
if (idig(ar__1(188:198)).EQ.0) e_GauSemi1_1 = rNULL__
if (idig(ar__1(212:222)).EQ.0) e_GauSemi2_1 = rNULL__
if (idig(ar__1(230:233)).EQ.0) e_GauTheta_1 = rNULL__
if (CirrusN_1 .EQ. 45) CirrusN_1 = iNULL__
if (idig(ar__1(305:314)).EQ.0) SkyBri_1 = rNULL__
if (idig(ar__1(316:327)).EQ.0) P_1 = rNULL__
if (idig(ar__1(329:339)).EQ.0) e_P_1 = rNULL__
if (idig(ar__1(341:350)).EQ.0) AnP_1 = rNULL__
if (idig(ar__1(352:361)).EQ.0) e_AnP_1 = rNULL__
if (idig(ar__1(363:374)).EQ.0) ApP_1 = rNULL__
if (idig(ar__1(376:387)).EQ.0) e_ApP_1 = rNULL__
if (idig(ar__1(389:400)).EQ.0) ApAnP_1 = rNULL__
if (idig(ar__1(402:411)).EQ.0) e_ApAnP_1 = rNULL__
if (idig(ar__1(413:423)).EQ.0) PUL_1 = rNULL__
if (idig(ar__1(425:434)).EQ.0) ApPUL_1 = rNULL__
if (PStat_1 .EQ. 45) PStat_1 = iNULL__
if (idig(ar__1(438:447)).EQ.0) PX_1 = rNULL__
if (idig(ar__1(449:458)).EQ.0) e_PX_2 = rNULL__
if (idig(ar__1(460:469)).EQ.0) E_PX_3 = rNULL__
if (idig(ar__1(471:480)).EQ.0) AnPX_1 = rNULL__
if (idig(ar__1(482:490)).EQ.0) e_AnPX_2 = rNULL__
if (idig(ar__1(492:500)).EQ.0) E_AnPX_3 = rNULL__
if (idig(ar__1(501:512)).EQ.0) ApFlux217_1 = rNULL__
if (idig(ar__1(514:523)).EQ.0) e_ApFlux217_1 = rNULL__
if (idig(ar__1(525:536)).EQ.0) ApFlux353_1 = rNULL__
if (idig(ar__1(538:548)).EQ.0) e_ApFlux353_1 = rNULL__
if (idig(ar__1(550:561)).EQ.0) ApFlux545_1 = rNULL__
if (idig(ar__1(563:573)).EQ.0) e_ApFlux545_1 = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1,GLON_1,GLAT_1,RAdeg_1,DEdeg_1,DetFlux_1,e_DetFlux_1,
+ AperFlux_1,e_AperFlux_1,PSFFlux_1,e_PSFFlux_1,GauFlux_1,
+ e_GauFlux_1,GauSemi1_1,e_GauSemi1_1,GauSemi2_1,e_GauSemi2_1,
+ GauTheta_1,e_GauTheta_1,GauFWHMEff_1,Extended_1,ExtVal_1,
+ ERCSC_1,PCCS1_1,HRel_1,CirrusN_1,SkyBri_1,P_1,e_P_1,AnP_1,
+ e_AnP_1,ApP_1,e_ApP_1,ApAnP_1,e_ApAnP_1,PUL_1,ApPUL_1,PStat_1,
+ PX_1,e_PX_2,E_PX_3,AnPX_1,e_AnPX_2,E_AnPX_3,ApFlux217_1,
+ e_ApFlux217_1,ApFlux353_1,e_ApFlux353_1,ApFlux545_1,
+ e_ApFlux545_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'pccs070.dat' ! Second Catalogue of Compact Sources, 70GHz
C Format for file interpretation
3 format(
+ A23,1X,F11.7,1X,F11.7,1X,F11.7,1X,F11.7,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,F10.6,
+ 1X,F11.6,1X,F11.6,1X,F11.6,1X,F5.1,1X,F4.1,1X,F10.6,1X,I1,1X,
+ I1,1X,A24,1X,A23,1X,I2,1X,I2,1X,F10.6,1X,F12.6,1X,F11.6,1X,
+ F10.6,1X,F10.6,1X,F12.6,1X,F12.6,1X,F12.6,1X,F10.6,1X,F11.6,
+ 1X,F10.6,1X,I1,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F9.6,1X,
+ F9.6,F12.6,1X,F10.6,1X,F12.6,1X,F11.6,1X,F12.6,1X,F11.6)
C Effective file loading
open(unit=1,status='old',file=
+'pccs070.dat')
write(6,*) '....Loading file: pccs070.dat'
do i__=1,1296
read(1,'(A573)')ar__2
read(ar__2,3)
+ Name_2,GLON_2,GLAT_2,RAdeg_2,DEdeg_2,DetFlux_2,e_DetFlux_2,
+ AperFlux_2,e_AperFlux_2,PSFFlux_2,e_PSFFlux_2,GauFlux_2,
+ e_GauFlux_2,GauSemi1_2,e_GauSemi1_2,GauSemi2_2,e_GauSemi2_2,
+ GauTheta_2,e_GauTheta_2,GauFWHMEff_2,Extended_2,ExtVal_2,
+ ERCSC_2,PCCS1_2,HRel_2,CirrusN_2,SkyBri_2,P_2,e_P_2,AnP_2,
+ e_AnP_2,ApP_2,e_ApP_2,ApAnP_2,e_ApAnP_2,PUL_2,ApPUL_2,PStat_2,
+ PX_2,e_PX_4,E_PX_5,AnPX_2,e_AnPX_4,E_AnPX_5,ApFlux217_2,
+ e_ApFlux217_2,ApFlux353_2,e_ApFlux353_2,ApFlux545_2,
+ e_ApFlux545_2
if (idig(ar__2(164:175)).EQ.0) e_GauFlux_2 = rNULL__
if (idig(ar__2(188:198)).EQ.0) e_GauSemi1_2 = rNULL__
if (idig(ar__2(212:222)).EQ.0) e_GauSemi2_2 = rNULL__
if (idig(ar__2(230:233)).EQ.0) e_GauTheta_2 = rNULL__
if (CirrusN_2 .EQ. 45) CirrusN_2 = iNULL__
if (idig(ar__2(305:314)).EQ.0) SkyBri_2 = rNULL__
if (idig(ar__2(316:327)).EQ.0) P_2 = rNULL__
if (idig(ar__2(329:339)).EQ.0) e_P_2 = rNULL__
if (idig(ar__2(341:350)).EQ.0) AnP_2 = rNULL__
if (idig(ar__2(352:361)).EQ.0) e_AnP_2 = rNULL__
if (idig(ar__2(363:374)).EQ.0) ApP_2 = rNULL__
if (idig(ar__2(376:387)).EQ.0) e_ApP_2 = rNULL__
if (idig(ar__2(389:400)).EQ.0) ApAnP_2 = rNULL__
if (idig(ar__2(402:411)).EQ.0) e_ApAnP_2 = rNULL__
if (idig(ar__2(413:423)).EQ.0) PUL_2 = rNULL__
if (idig(ar__2(425:434)).EQ.0) ApPUL_2 = rNULL__
if (PStat_2 .EQ. 45) PStat_2 = iNULL__
if (idig(ar__2(438:447)).EQ.0) PX_2 = rNULL__
if (idig(ar__2(449:458)).EQ.0) e_PX_4 = rNULL__
if (idig(ar__2(460:469)).EQ.0) E_PX_5 = rNULL__
if (idig(ar__2(471:480)).EQ.0) AnPX_2 = rNULL__
if (idig(ar__2(482:490)).EQ.0) e_AnPX_4 = rNULL__
if (idig(ar__2(492:500)).EQ.0) E_AnPX_5 = rNULL__
if (idig(ar__2(501:512)).EQ.0) ApFlux217_2 = rNULL__
if (idig(ar__2(514:523)).EQ.0) e_ApFlux217_2 = rNULL__
if (idig(ar__2(525:536)).EQ.0) ApFlux353_2 = rNULL__
if (idig(ar__2(538:548)).EQ.0) e_ApFlux353_2 = rNULL__
if (idig(ar__2(550:561)).EQ.0) ApFlux545_2 = rNULL__
if (idig(ar__2(563:573)).EQ.0) e_ApFlux545_2 = rNULL__
c ..............Just test output...........
write(6,3)
+ Name_2,GLON_2,GLAT_2,RAdeg_2,DEdeg_2,DetFlux_2,e_DetFlux_2,
+ AperFlux_2,e_AperFlux_2,PSFFlux_2,e_PSFFlux_2,GauFlux_2,
+ e_GauFlux_2,GauSemi1_2,e_GauSemi1_2,GauSemi2_2,e_GauSemi2_2,
+ GauTheta_2,e_GauTheta_2,GauFWHMEff_2,Extended_2,ExtVal_2,
+ ERCSC_2,PCCS1_2,HRel_2,CirrusN_2,SkyBri_2,P_2,e_P_2,AnP_2,
+ e_AnP_2,ApP_2,e_ApP_2,ApAnP_2,e_ApAnP_2,PUL_2,ApPUL_2,PStat_2,
+ PX_2,e_PX_4,E_PX_5,AnPX_2,e_AnPX_4,E_AnPX_5,ApFlux217_2,
+ e_ApFlux217_2,ApFlux353_2,e_ApFlux353_2,ApFlux545_2,
+ e_ApFlux545_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'pccs100.dat' ! Second Catalogue of Compact Sources, 100GHz
C Format for file interpretation
4 format(
+ A23,1X,F11.7,1X,F11.7,1X,F11.7,1X,F11.7,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,F10.6,
+ 1X,F11.6,1X,F11.6,1X,F11.6,1X,F5.1,1X,F4.1,1X,F10.6,1X,I1,1X,
+ I1,1X,A24,1X,A23,1X,I2,1X,I2,1X,F10.6,1X,F12.6,1X,F11.6,1X,
+ F10.6,1X,F10.6,1X,F12.6,1X,F12.6,1X,F12.6,1X,F10.6,1X,F11.6,
+ 1X,F10.6,1X,I1,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F9.6,1X,
+ F9.6,F12.6,1X,F10.6,1X,F12.6,1X,F11.6,1X,F12.6,1X,F11.6)
C Effective file loading
open(unit=1,status='old',file=
+'pccs100.dat')
write(6,*) '....Loading file: pccs100.dat'
do i__=1,1742
read(1,'(A573)')ar__3
read(ar__3,4)
+ Name_3,GLON_3,GLAT_3,RAdeg_3,DEdeg_3,DetFlux_3,e_DetFlux_3,
+ AperFlux_3,e_AperFlux_3,PSFFlux_3,e_PSFFlux_3,GauFlux_3,
+ e_GauFlux_3,GauSemi1_3,e_GauSemi1_3,GauSemi2_3,e_GauSemi2_3,
+ GauTheta_3,e_GauTheta_3,GauFWHMEff_3,Extended_3,ExtVal_3,
+ ERCSC_3,PCCS1_3,HRel_3,CirrusN_3,SkyBri_3,P_3,e_P_3,AnP_3,
+ e_AnP_3,ApP_3,e_ApP_3,ApAnP_3,e_ApAnP_3,PUL_3,ApPUL_3,PStat_3,
+ PX_3,e_PX_6,E_PX_7,AnPX_3,e_AnPX_6,E_AnPX_7,ApFlux217_3,
+ e_ApFlux217_3,ApFlux353_3,e_ApFlux353_3,ApFlux545_3,
+ e_ApFlux545_3
if (idig(ar__3(164:175)).EQ.0) e_GauFlux_3 = rNULL__
if (idig(ar__3(188:198)).EQ.0) e_GauSemi1_3 = rNULL__
if (idig(ar__3(212:222)).EQ.0) e_GauSemi2_3 = rNULL__
if (idig(ar__3(230:233)).EQ.0) e_GauTheta_3 = rNULL__
if (CirrusN_3 .EQ. 45) CirrusN_3 = iNULL__
if (idig(ar__3(305:314)).EQ.0) SkyBri_3 = rNULL__
if (idig(ar__3(316:327)).EQ.0) P_3 = rNULL__
if (idig(ar__3(329:339)).EQ.0) e_P_3 = rNULL__
if (idig(ar__3(341:350)).EQ.0) AnP_3 = rNULL__
if (idig(ar__3(352:361)).EQ.0) e_AnP_3 = rNULL__
if (idig(ar__3(363:374)).EQ.0) ApP_3 = rNULL__
if (idig(ar__3(376:387)).EQ.0) e_ApP_3 = rNULL__
if (idig(ar__3(389:400)).EQ.0) ApAnP_3 = rNULL__
if (idig(ar__3(402:411)).EQ.0) e_ApAnP_3 = rNULL__
if (idig(ar__3(413:423)).EQ.0) PUL_3 = rNULL__
if (idig(ar__3(425:434)).EQ.0) ApPUL_3 = rNULL__
if (PStat_3 .EQ. 45) PStat_3 = iNULL__
if (idig(ar__3(438:447)).EQ.0) PX_3 = rNULL__
if (idig(ar__3(449:458)).EQ.0) e_PX_6 = rNULL__
if (idig(ar__3(460:469)).EQ.0) E_PX_7 = rNULL__
if (idig(ar__3(471:480)).EQ.0) AnPX_3 = rNULL__
if (idig(ar__3(482:490)).EQ.0) e_AnPX_6 = rNULL__
if (idig(ar__3(492:500)).EQ.0) E_AnPX_7 = rNULL__
if (idig(ar__3(501:512)).EQ.0) ApFlux217_3 = rNULL__
if (idig(ar__3(514:523)).EQ.0) e_ApFlux217_3 = rNULL__
if (idig(ar__3(525:536)).EQ.0) ApFlux353_3 = rNULL__
if (idig(ar__3(538:548)).EQ.0) e_ApFlux353_3 = rNULL__
if (idig(ar__3(550:561)).EQ.0) ApFlux545_3 = rNULL__
if (idig(ar__3(563:573)).EQ.0) e_ApFlux545_3 = rNULL__
c ..............Just test output...........
write(6,4)
+ Name_3,GLON_3,GLAT_3,RAdeg_3,DEdeg_3,DetFlux_3,e_DetFlux_3,
+ AperFlux_3,e_AperFlux_3,PSFFlux_3,e_PSFFlux_3,GauFlux_3,
+ e_GauFlux_3,GauSemi1_3,e_GauSemi1_3,GauSemi2_3,e_GauSemi2_3,
+ GauTheta_3,e_GauTheta_3,GauFWHMEff_3,Extended_3,ExtVal_3,
+ ERCSC_3,PCCS1_3,HRel_3,CirrusN_3,SkyBri_3,P_3,e_P_3,AnP_3,
+ e_AnP_3,ApP_3,e_ApP_3,ApAnP_3,e_ApAnP_3,PUL_3,ApPUL_3,PStat_3,
+ PX_3,e_PX_6,E_PX_7,AnPX_3,e_AnPX_6,E_AnPX_7,ApFlux217_3,
+ e_ApFlux217_3,ApFlux353_3,e_ApFlux353_3,ApFlux545_3,
+ e_ApFlux545_3
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'pccs143.dat' ! Second Catalogue of Compact Sources, 143GHz
C Format for file interpretation
5 format(
+ A23,1X,F11.7,1X,F11.7,1X,F11.7,1X,F11.7,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,F10.6,
+ 1X,F11.6,1X,F11.6,1X,F11.6,1X,F5.1,1X,F4.1,1X,F10.6,1X,I1,1X,
+ I1,1X,A24,1X,A23,1X,I2,1X,I2,1X,F10.6,1X,F12.6,1X,F11.6,1X,
+ F10.6,1X,F10.6,1X,F12.6,1X,F12.6,1X,F12.6,1X,F10.6,1X,F11.6,
+ 1X,F10.6,1X,I1,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F9.6,1X,
+ F9.6,F12.6,1X,F10.6,1X,F12.6,1X,F11.6,1X,F12.6,1X,F11.6)
C Effective file loading
open(unit=1,status='old',file=
+'pccs143.dat')
write(6,*) '....Loading file: pccs143.dat'
do i__=1,2160
read(1,'(A573)')ar__4
read(ar__4,5)
+ Name_4,GLON_4,GLAT_4,RAdeg_4,DEdeg_4,DetFlux_4,e_DetFlux_4,
+ AperFlux_4,e_AperFlux_4,PSFFlux_4,e_PSFFlux_4,GauFlux_4,
+ e_GauFlux_4,GauSemi1_4,e_GauSemi1_4,GauSemi2_4,e_GauSemi2_4,
+ GauTheta_4,e_GauTheta_4,GauFWHMEff_4,Extended_4,ExtVal_4,
+ ERCSC_4,PCCS1_4,HRel_4,CirrusN_4,SkyBri_4,P_4,e_P_4,AnP_4,
+ e_AnP_4,ApP_4,e_ApP_4,ApAnP_4,e_ApAnP_4,PUL_4,ApPUL_4,PStat_4,
+ PX_4,e_PX_8,E_PX_9,AnPX_4,e_AnPX_8,E_AnPX_9,ApFlux217_4,
+ e_ApFlux217_4,ApFlux353_4,e_ApFlux353_4,ApFlux545_4,
+ e_ApFlux545_4
if (idig(ar__4(164:175)).EQ.0) e_GauFlux_4 = rNULL__
if (idig(ar__4(188:198)).EQ.0) e_GauSemi1_4 = rNULL__
if (idig(ar__4(212:222)).EQ.0) e_GauSemi2_4 = rNULL__
if (idig(ar__4(230:233)).EQ.0) e_GauTheta_4 = rNULL__
if (CirrusN_4 .EQ. 45) CirrusN_4 = iNULL__
if (idig(ar__4(305:314)).EQ.0) SkyBri_4 = rNULL__
if (idig(ar__4(316:327)).EQ.0) P_4 = rNULL__
if (idig(ar__4(329:339)).EQ.0) e_P_4 = rNULL__
if (idig(ar__4(341:350)).EQ.0) AnP_4 = rNULL__
if (idig(ar__4(352:361)).EQ.0) e_AnP_4 = rNULL__
if (idig(ar__4(363:374)).EQ.0) ApP_4 = rNULL__
if (idig(ar__4(376:387)).EQ.0) e_ApP_4 = rNULL__
if (idig(ar__4(389:400)).EQ.0) ApAnP_4 = rNULL__
if (idig(ar__4(402:411)).EQ.0) e_ApAnP_4 = rNULL__
if (idig(ar__4(413:423)).EQ.0) PUL_4 = rNULL__
if (idig(ar__4(425:434)).EQ.0) ApPUL_4 = rNULL__
if (PStat_4 .EQ. 45) PStat_4 = iNULL__
if (idig(ar__4(438:447)).EQ.0) PX_4 = rNULL__
if (idig(ar__4(449:458)).EQ.0) e_PX_8 = rNULL__
if (idig(ar__4(460:469)).EQ.0) E_PX_9 = rNULL__
if (idig(ar__4(471:480)).EQ.0) AnPX_4 = rNULL__
if (idig(ar__4(482:490)).EQ.0) e_AnPX_8 = rNULL__
if (idig(ar__4(492:500)).EQ.0) E_AnPX_9 = rNULL__
if (idig(ar__4(501:512)).EQ.0) ApFlux217_4 = rNULL__
if (idig(ar__4(514:523)).EQ.0) e_ApFlux217_4 = rNULL__
if (idig(ar__4(525:536)).EQ.0) ApFlux353_4 = rNULL__
if (idig(ar__4(538:548)).EQ.0) e_ApFlux353_4 = rNULL__
if (idig(ar__4(550:561)).EQ.0) ApFlux545_4 = rNULL__
if (idig(ar__4(563:573)).EQ.0) e_ApFlux545_4 = rNULL__
c ..............Just test output...........
write(6,5)
+ Name_4,GLON_4,GLAT_4,RAdeg_4,DEdeg_4,DetFlux_4,e_DetFlux_4,
+ AperFlux_4,e_AperFlux_4,PSFFlux_4,e_PSFFlux_4,GauFlux_4,
+ e_GauFlux_4,GauSemi1_4,e_GauSemi1_4,GauSemi2_4,e_GauSemi2_4,
+ GauTheta_4,e_GauTheta_4,GauFWHMEff_4,Extended_4,ExtVal_4,
+ ERCSC_4,PCCS1_4,HRel_4,CirrusN_4,SkyBri_4,P_4,e_P_4,AnP_4,
+ e_AnP_4,ApP_4,e_ApP_4,ApAnP_4,e_ApAnP_4,PUL_4,ApPUL_4,PStat_4,
+ PX_4,e_PX_8,E_PX_9,AnPX_4,e_AnPX_8,E_AnPX_9,ApFlux217_4,
+ e_ApFlux217_4,ApFlux353_4,e_ApFlux353_4,ApFlux545_4,
+ e_ApFlux545_4
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'pccs217.dat' ! Second Catalogue of Compact Sources, 217GHz
C Format for file interpretation
6 format(
+ A23,1X,F11.7,1X,F11.7,1X,F11.7,1X,F11.7,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,F10.6,
+ 1X,F11.6,1X,F11.6,1X,F11.6,1X,F5.1,1X,F4.1,1X,F10.6,1X,I1,1X,
+ I1,1X,A24,1X,A23,1X,I2,1X,I2,1X,F10.6,1X,F12.6,1X,F11.6,1X,
+ F10.6,1X,F10.6,1X,F12.6,1X,F12.6,1X,F12.6,1X,F10.6,1X,F11.6,
+ 1X,F10.6,1X,I1,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F9.6,1X,
+ F9.6,F12.6,1X,F10.6,1X,F12.6,1X,F11.6,1X,F12.6,1X,F11.6)
C Effective file loading
open(unit=1,status='old',file=
+'pccs217.dat')
write(6,*) '....Loading file: pccs217.dat'
do i__=1,2135
read(1,'(A573)')ar__5
read(ar__5,6)
+ Name_5,GLON_5,GLAT_5,RAdeg_5,DEdeg_5,DetFlux_5,e_DetFlux_5,
+ AperFlux_5,e_AperFlux_5,PSFFlux_5,e_PSFFlux_5,GauFlux_5,
+ e_GauFlux_5,GauSemi1_5,e_GauSemi1_5,GauSemi2_5,e_GauSemi2_5,
+ GauTheta_5,e_GauTheta_5,GauFWHMEff_5,Extended_5,ExtVal_5,
+ ERCSC_5,PCCS1_5,HRel_5,CirrusN_5,SkyBri_5,P_5,e_P_5,AnP_5,
+ e_AnP_5,ApP_5,e_ApP_5,ApAnP_5,e_ApAnP_5,PUL_5,ApPUL_5,PStat_5,
+ PX_5,e_PX_10,E_PX_11,AnPX_5,e_AnPX_10,E_AnPX_11,ApFlux217_5,
+ e_ApFlux217_5,ApFlux353_5,e_ApFlux353_5,ApFlux545_5,
+ e_ApFlux545_5
if (idig(ar__5(164:175)).EQ.0) e_GauFlux_5 = rNULL__
if (idig(ar__5(188:198)).EQ.0) e_GauSemi1_5 = rNULL__
if (idig(ar__5(212:222)).EQ.0) e_GauSemi2_5 = rNULL__
if (idig(ar__5(230:233)).EQ.0) e_GauTheta_5 = rNULL__
if (CirrusN_5 .EQ. 45) CirrusN_5 = iNULL__
if (idig(ar__5(305:314)).EQ.0) SkyBri_5 = rNULL__
if (idig(ar__5(316:327)).EQ.0) P_5 = rNULL__
if (idig(ar__5(329:339)).EQ.0) e_P_5 = rNULL__
if (idig(ar__5(341:350)).EQ.0) AnP_5 = rNULL__
if (idig(ar__5(352:361)).EQ.0) e_AnP_5 = rNULL__
if (idig(ar__5(363:374)).EQ.0) ApP_5 = rNULL__
if (idig(ar__5(376:387)).EQ.0) e_ApP_5 = rNULL__
if (idig(ar__5(389:400)).EQ.0) ApAnP_5 = rNULL__
if (idig(ar__5(402:411)).EQ.0) e_ApAnP_5 = rNULL__
if (idig(ar__5(413:423)).EQ.0) PUL_5 = rNULL__
if (idig(ar__5(425:434)).EQ.0) ApPUL_5 = rNULL__
if (PStat_5 .EQ. 45) PStat_5 = iNULL__
if (idig(ar__5(438:447)).EQ.0) PX_5 = rNULL__
if (idig(ar__5(449:458)).EQ.0) e_PX_10 = rNULL__
if (idig(ar__5(460:469)).EQ.0) E_PX_11 = rNULL__
if (idig(ar__5(471:480)).EQ.0) AnPX_5 = rNULL__
if (idig(ar__5(482:490)).EQ.0) e_AnPX_10 = rNULL__
if (idig(ar__5(492:500)).EQ.0) E_AnPX_11 = rNULL__
if (idig(ar__5(501:512)).EQ.0) ApFlux217_5 = rNULL__
if (idig(ar__5(514:523)).EQ.0) e_ApFlux217_5 = rNULL__
if (idig(ar__5(525:536)).EQ.0) ApFlux353_5 = rNULL__
if (idig(ar__5(538:548)).EQ.0) e_ApFlux353_5 = rNULL__
if (idig(ar__5(550:561)).EQ.0) ApFlux545_5 = rNULL__
if (idig(ar__5(563:573)).EQ.0) e_ApFlux545_5 = rNULL__
c ..............Just test output...........
write(6,6)
+ Name_5,GLON_5,GLAT_5,RAdeg_5,DEdeg_5,DetFlux_5,e_DetFlux_5,
+ AperFlux_5,e_AperFlux_5,PSFFlux_5,e_PSFFlux_5,GauFlux_5,
+ e_GauFlux_5,GauSemi1_5,e_GauSemi1_5,GauSemi2_5,e_GauSemi2_5,
+ GauTheta_5,e_GauTheta_5,GauFWHMEff_5,Extended_5,ExtVal_5,
+ ERCSC_5,PCCS1_5,HRel_5,CirrusN_5,SkyBri_5,P_5,e_P_5,AnP_5,
+ e_AnP_5,ApP_5,e_ApP_5,ApAnP_5,e_ApAnP_5,PUL_5,ApPUL_5,PStat_5,
+ PX_5,e_PX_10,E_PX_11,AnPX_5,e_AnPX_10,E_AnPX_11,ApFlux217_5,
+ e_ApFlux217_5,ApFlux353_5,e_ApFlux353_5,ApFlux545_5,
+ e_ApFlux545_5
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'pccs353.dat' ! Second Catalogue of Compact Sources, 353GHz
C Format for file interpretation
7 format(
+ A23,1X,F11.7,1X,F11.7,1X,F11.7,1X,F11.7,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,F10.6,
+ 1X,F11.6,1X,F11.6,1X,F11.6,1X,F5.1,1X,F4.1,1X,F10.6,1X,I1,1X,
+ I1,1X,A24,1X,A23,1X,I2,1X,I2,1X,F10.6,1X,F12.6,1X,F11.6,1X,
+ F10.6,1X,F10.6,1X,F12.6,1X,F12.6,1X,F12.6,1X,F10.6,1X,F11.6,
+ 1X,F10.6,1X,I1,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F9.6,1X,
+ F9.6,F12.6,1X,F10.6,1X,F12.6,1X,F11.6,1X,F12.6,1X,F11.6)
C Effective file loading
open(unit=1,status='old',file=
+'pccs353.dat')
write(6,*) '....Loading file: pccs353.dat'
do i__=1,1344
read(1,'(A573)')ar__6
read(ar__6,7)
+ Name_6,GLON_6,GLAT_6,RAdeg_6,DEdeg_6,DetFlux_6,e_DetFlux_6,
+ AperFlux_6,e_AperFlux_6,PSFFlux_6,e_PSFFlux_6,GauFlux_6,
+ e_GauFlux_6,GauSemi1_6,e_GauSemi1_6,GauSemi2_6,e_GauSemi2_6,
+ GauTheta_6,e_GauTheta_6,GauFWHMEff_6,Extended_6,ExtVal_6,
+ ERCSC_6,PCCS1_6,HRel_6,CirrusN_6,SkyBri_6,P_6,e_P_6,AnP_6,
+ e_AnP_6,ApP_6,e_ApP_6,ApAnP_6,e_ApAnP_6,PUL_6,ApPUL_6,PStat_6,
+ PX_6,e_PX_12,E_PX_13,AnPX_6,e_AnPX_12,E_AnPX_13,ApFlux217_6,
+ e_ApFlux217_6,ApFlux353_6,e_ApFlux353_6,ApFlux545_6,
+ e_ApFlux545_6
if (idig(ar__6(164:175)).EQ.0) e_GauFlux_6 = rNULL__
if (idig(ar__6(188:198)).EQ.0) e_GauSemi1_6 = rNULL__
if (idig(ar__6(212:222)).EQ.0) e_GauSemi2_6 = rNULL__
if (idig(ar__6(230:233)).EQ.0) e_GauTheta_6 = rNULL__
if (CirrusN_6 .EQ. 45) CirrusN_6 = iNULL__
if (idig(ar__6(305:314)).EQ.0) SkyBri_6 = rNULL__
if (idig(ar__6(316:327)).EQ.0) P_6 = rNULL__
if (idig(ar__6(329:339)).EQ.0) e_P_6 = rNULL__
if (idig(ar__6(341:350)).EQ.0) AnP_6 = rNULL__
if (idig(ar__6(352:361)).EQ.0) e_AnP_6 = rNULL__
if (idig(ar__6(363:374)).EQ.0) ApP_6 = rNULL__
if (idig(ar__6(376:387)).EQ.0) e_ApP_6 = rNULL__
if (idig(ar__6(389:400)).EQ.0) ApAnP_6 = rNULL__
if (idig(ar__6(402:411)).EQ.0) e_ApAnP_6 = rNULL__
if (idig(ar__6(413:423)).EQ.0) PUL_6 = rNULL__
if (idig(ar__6(425:434)).EQ.0) ApPUL_6 = rNULL__
if (PStat_6 .EQ. 45) PStat_6 = iNULL__
if (idig(ar__6(438:447)).EQ.0) PX_6 = rNULL__
if (idig(ar__6(449:458)).EQ.0) e_PX_12 = rNULL__
if (idig(ar__6(460:469)).EQ.0) E_PX_13 = rNULL__
if (idig(ar__6(471:480)).EQ.0) AnPX_6 = rNULL__
if (idig(ar__6(482:490)).EQ.0) e_AnPX_12 = rNULL__
if (idig(ar__6(492:500)).EQ.0) E_AnPX_13 = rNULL__
if (idig(ar__6(501:512)).EQ.0) ApFlux217_6 = rNULL__
if (idig(ar__6(514:523)).EQ.0) e_ApFlux217_6 = rNULL__
if (idig(ar__6(525:536)).EQ.0) ApFlux353_6 = rNULL__
if (idig(ar__6(538:548)).EQ.0) e_ApFlux353_6 = rNULL__
if (idig(ar__6(550:561)).EQ.0) ApFlux545_6 = rNULL__
if (idig(ar__6(563:573)).EQ.0) e_ApFlux545_6 = rNULL__
c ..............Just test output...........
write(6,7)
+ Name_6,GLON_6,GLAT_6,RAdeg_6,DEdeg_6,DetFlux_6,e_DetFlux_6,
+ AperFlux_6,e_AperFlux_6,PSFFlux_6,e_PSFFlux_6,GauFlux_6,
+ e_GauFlux_6,GauSemi1_6,e_GauSemi1_6,GauSemi2_6,e_GauSemi2_6,
+ GauTheta_6,e_GauTheta_6,GauFWHMEff_6,Extended_6,ExtVal_6,
+ ERCSC_6,PCCS1_6,HRel_6,CirrusN_6,SkyBri_6,P_6,e_P_6,AnP_6,
+ e_AnP_6,ApP_6,e_ApP_6,ApAnP_6,e_ApAnP_6,PUL_6,ApPUL_6,PStat_6,
+ PX_6,e_PX_12,E_PX_13,AnPX_6,e_AnPX_12,E_AnPX_13,ApFlux217_6,
+ e_ApFlux217_6,ApFlux353_6,e_ApFlux353_6,ApFlux545_6,
+ e_ApFlux545_6
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'pccs545.dat' ! Second Catalogue of Compact Sources, 545GHz
C Format for file interpretation
8 format(
+ A23,1X,F11.7,1X,F11.7,1X,F11.7,1X,F11.7,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,F10.6,
+ 1X,F11.6,1X,F11.6,1X,F11.6,1X,F5.1,1X,F4.1,1X,F10.6,1X,I1,1X,
+ I1,1X,A24,1X,A23,1X,I2,1X,I2,1X,F10.6,1X,F12.6,1X,F11.6,1X,
+ F10.6,1X,F10.6,1X,F12.6,1X,F12.6,1X,F12.6,1X,F10.6,1X,F11.6,
+ 1X,F10.6,1X,I1,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F9.6,1X,
+ F9.6,F12.6,1X,F10.6,1X,F12.6,1X,F11.6,1X,F12.6,1X,F11.6)
C Effective file loading
open(unit=1,status='old',file=
+'pccs545.dat')
write(6,*) '....Loading file: pccs545.dat'
do i__=1,1694
read(1,'(A573)')ar__7
read(ar__7,8)
+ Name_7,GLON_7,GLAT_7,RAdeg_7,DEdeg_7,DetFlux_7,e_DetFlux_7,
+ AperFlux_7,e_AperFlux_7,PSFFlux_7,e_PSFFlux_7,GauFlux_7,
+ e_GauFlux_7,GauSemi1_7,e_GauSemi1_7,GauSemi2_7,e_GauSemi2_7,
+ GauTheta_7,e_GauTheta_7,GauFWHMEff_7,Extended_7,ExtVal_7,
+ ERCSC_7,PCCS1_7,HRel_7,CirrusN_7,SkyBri_7,P_7,e_P_7,AnP_7,
+ e_AnP_7,ApP_7,e_ApP_7,ApAnP_7,e_ApAnP_7,PUL_7,ApPUL_7,PStat_7,
+ PX_7,e_PX_14,E_PX_15,AnPX_7,e_AnPX_14,E_AnPX_15,ApFlux217_7,
+ e_ApFlux217_7,ApFlux353_7,e_ApFlux353_7,ApFlux545_7,
+ e_ApFlux545_7
if (idig(ar__7(164:175)).EQ.0) e_GauFlux_7 = rNULL__
if (idig(ar__7(188:198)).EQ.0) e_GauSemi1_7 = rNULL__
if (idig(ar__7(212:222)).EQ.0) e_GauSemi2_7 = rNULL__
if (idig(ar__7(230:233)).EQ.0) e_GauTheta_7 = rNULL__
if (CirrusN_7 .EQ. 45) CirrusN_7 = iNULL__
if (idig(ar__7(305:314)).EQ.0) SkyBri_7 = rNULL__
if (idig(ar__7(316:327)).EQ.0) P_7 = rNULL__
if (idig(ar__7(329:339)).EQ.0) e_P_7 = rNULL__
if (idig(ar__7(341:350)).EQ.0) AnP_7 = rNULL__
if (idig(ar__7(352:361)).EQ.0) e_AnP_7 = rNULL__
if (idig(ar__7(363:374)).EQ.0) ApP_7 = rNULL__
if (idig(ar__7(376:387)).EQ.0) e_ApP_7 = rNULL__
if (idig(ar__7(389:400)).EQ.0) ApAnP_7 = rNULL__
if (idig(ar__7(402:411)).EQ.0) e_ApAnP_7 = rNULL__
if (idig(ar__7(413:423)).EQ.0) PUL_7 = rNULL__
if (idig(ar__7(425:434)).EQ.0) ApPUL_7 = rNULL__
if (PStat_7 .EQ. 45) PStat_7 = iNULL__
if (idig(ar__7(438:447)).EQ.0) PX_7 = rNULL__
if (idig(ar__7(449:458)).EQ.0) e_PX_14 = rNULL__
if (idig(ar__7(460:469)).EQ.0) E_PX_15 = rNULL__
if (idig(ar__7(471:480)).EQ.0) AnPX_7 = rNULL__
if (idig(ar__7(482:490)).EQ.0) e_AnPX_14 = rNULL__
if (idig(ar__7(492:500)).EQ.0) E_AnPX_15 = rNULL__
if (idig(ar__7(501:512)).EQ.0) ApFlux217_7 = rNULL__
if (idig(ar__7(514:523)).EQ.0) e_ApFlux217_7 = rNULL__
if (idig(ar__7(525:536)).EQ.0) ApFlux353_7 = rNULL__
if (idig(ar__7(538:548)).EQ.0) e_ApFlux353_7 = rNULL__
if (idig(ar__7(550:561)).EQ.0) ApFlux545_7 = rNULL__
if (idig(ar__7(563:573)).EQ.0) e_ApFlux545_7 = rNULL__
c ..............Just test output...........
write(6,8)
+ Name_7,GLON_7,GLAT_7,RAdeg_7,DEdeg_7,DetFlux_7,e_DetFlux_7,
+ AperFlux_7,e_AperFlux_7,PSFFlux_7,e_PSFFlux_7,GauFlux_7,
+ e_GauFlux_7,GauSemi1_7,e_GauSemi1_7,GauSemi2_7,e_GauSemi2_7,
+ GauTheta_7,e_GauTheta_7,GauFWHMEff_7,Extended_7,ExtVal_7,
+ ERCSC_7,PCCS1_7,HRel_7,CirrusN_7,SkyBri_7,P_7,e_P_7,AnP_7,
+ e_AnP_7,ApP_7,e_ApP_7,ApAnP_7,e_ApAnP_7,PUL_7,ApPUL_7,PStat_7,
+ PX_7,e_PX_14,E_PX_15,AnPX_7,e_AnPX_14,E_AnPX_15,ApFlux217_7,
+ e_ApFlux217_7,ApFlux353_7,e_ApFlux353_7,ApFlux545_7,
+ e_ApFlux545_7
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'pccs857.dat' ! Second Catalogue of Compact Sources, 857GHz
C Format for file interpretation
9 format(
+ A23,1X,F11.7,1X,F11.7,1X,F11.7,1X,F11.7,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,F10.6,
+ 1X,F11.6,1X,F11.6,1X,F11.6,1X,F5.1,1X,F4.1,1X,F10.6,1X,I1,1X,
+ I1,1X,A24,1X,A23,1X,I2,1X,I2,1X,F10.6,1X,F12.6,1X,F11.6,1X,
+ F10.6,1X,F10.6,1X,F12.6,1X,F12.6,1X,F12.6,1X,F10.6,1X,F11.6,
+ 1X,F10.6,1X,I1,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F9.6,1X,
+ F9.6,F12.6,1X,F10.6,1X,F12.6,1X,F11.6,1X,F12.6,1X,F11.6)
C Effective file loading
open(unit=1,status='old',file=
+'pccs857.dat')
write(6,*) '....Loading file: pccs857.dat'
do i__=1,4891
read(1,'(A573)')ar__8
read(ar__8,9)
+ Name_8,GLON_8,GLAT_8,RAdeg_8,DEdeg_8,DetFlux_8,e_DetFlux_8,
+ AperFlux_8,e_AperFlux_8,PSFFlux_8,e_PSFFlux_8,GauFlux_8,
+ e_GauFlux_8,GauSemi1_8,e_GauSemi1_8,GauSemi2_8,e_GauSemi2_8,
+ GauTheta_8,e_GauTheta_8,GauFWHMEff_8,Extended_8,ExtVal_8,
+ ERCSC_8,PCCS1_8,HRel_8,CirrusN_8,SkyBri_8,P_8,e_P_8,AnP_8,
+ e_AnP_8,ApP_8,e_ApP_8,ApAnP_8,e_ApAnP_8,PUL_8,ApPUL_8,PStat_8,
+ PX_8,e_PX_16,E_PX_17,AnPX_8,e_AnPX_16,E_AnPX_17,ApFlux217_8,
+ e_ApFlux217_8,ApFlux353_8,e_ApFlux353_8,ApFlux545_8,
+ e_ApFlux545_8
if (idig(ar__8(164:175)).EQ.0) e_GauFlux_8 = rNULL__
if (idig(ar__8(188:198)).EQ.0) e_GauSemi1_8 = rNULL__
if (idig(ar__8(212:222)).EQ.0) e_GauSemi2_8 = rNULL__
if (idig(ar__8(230:233)).EQ.0) e_GauTheta_8 = rNULL__
if (CirrusN_8 .EQ. 45) CirrusN_8 = iNULL__
if (idig(ar__8(305:314)).EQ.0) SkyBri_8 = rNULL__
if (idig(ar__8(316:327)).EQ.0) P_8 = rNULL__
if (idig(ar__8(329:339)).EQ.0) e_P_8 = rNULL__
if (idig(ar__8(341:350)).EQ.0) AnP_8 = rNULL__
if (idig(ar__8(352:361)).EQ.0) e_AnP_8 = rNULL__
if (idig(ar__8(363:374)).EQ.0) ApP_8 = rNULL__
if (idig(ar__8(376:387)).EQ.0) e_ApP_8 = rNULL__
if (idig(ar__8(389:400)).EQ.0) ApAnP_8 = rNULL__
if (idig(ar__8(402:411)).EQ.0) e_ApAnP_8 = rNULL__
if (idig(ar__8(413:423)).EQ.0) PUL_8 = rNULL__
if (idig(ar__8(425:434)).EQ.0) ApPUL_8 = rNULL__
if (PStat_8 .EQ. 45) PStat_8 = iNULL__
if (idig(ar__8(438:447)).EQ.0) PX_8 = rNULL__
if (idig(ar__8(449:458)).EQ.0) e_PX_16 = rNULL__
if (idig(ar__8(460:469)).EQ.0) E_PX_17 = rNULL__
if (idig(ar__8(471:480)).EQ.0) AnPX_8 = rNULL__
if (idig(ar__8(482:490)).EQ.0) e_AnPX_16 = rNULL__
if (idig(ar__8(492:500)).EQ.0) E_AnPX_17 = rNULL__
if (idig(ar__8(501:512)).EQ.0) ApFlux217_8 = rNULL__
if (idig(ar__8(514:523)).EQ.0) e_ApFlux217_8 = rNULL__
if (idig(ar__8(525:536)).EQ.0) ApFlux353_8 = rNULL__
if (idig(ar__8(538:548)).EQ.0) e_ApFlux353_8 = rNULL__
if (idig(ar__8(550:561)).EQ.0) ApFlux545_8 = rNULL__
if (idig(ar__8(563:573)).EQ.0) e_ApFlux545_8 = rNULL__
c ..............Just test output...........
write(6,9)
+ Name_8,GLON_8,GLAT_8,RAdeg_8,DEdeg_8,DetFlux_8,e_DetFlux_8,
+ AperFlux_8,e_AperFlux_8,PSFFlux_8,e_PSFFlux_8,GauFlux_8,
+ e_GauFlux_8,GauSemi1_8,e_GauSemi1_8,GauSemi2_8,e_GauSemi2_8,
+ GauTheta_8,e_GauTheta_8,GauFWHMEff_8,Extended_8,ExtVal_8,
+ ERCSC_8,PCCS1_8,HRel_8,CirrusN_8,SkyBri_8,P_8,e_P_8,AnP_8,
+ e_AnP_8,ApP_8,e_ApP_8,ApAnP_8,e_ApAnP_8,PUL_8,ApPUL_8,PStat_8,
+ PX_8,e_PX_16,E_PX_17,AnPX_8,e_AnPX_16,E_AnPX_17,ApFlux217_8,
+ e_ApFlux217_8,ApFlux353_8,e_ApFlux353_8,ApFlux545_8,
+ e_ApFlux545_8
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'pccs100e.dat' ! Second Catalogue of Compact Sources, 100GHz,
* unknown reliability
C Format for file interpretation
10 format(
+ A24,1X,F11.7,1X,F11.7,1X,F11.7,1X,F11.7,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,F10.6,
+ 1X,F12.6,1X,F11.6,F12.6,1X,F5.1,1X,F4.1,1X,F10.6,1X,F12.6,1X,
+ F11.6,1X,F10.6,1X,F10.6,1X,F12.6,1X,F12.6,1X,F12.6,1X,F10.6,
+ 1X,F11.6,1X,F10.6,1X,I1,1X,F11.6,1X,F11.6,1X,F11.6,1X,F10.6,
+ 1X,F9.6,1X,F9.6,1X,I1,1X,I1,1X,A24,1X,A23,1X,I2,1X,I3,1X,
+ F10.6,1X,I1,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
+ E12.6)
C Effective file loading
open(unit=1,status='old',file=
+'pccs100e.dat')
write(6,*) '....Loading file: pccs100e.dat'
do i__=1,2487
read(1,'(A586)')ar__9
read(ar__9,10)
+ Name_9,GLON_9,GLAT_9,RAdeg_9,DEdeg_9,DetFlux_9,e_DetFlux_9,
+ AperFlux_9,e_AperFlux_9,PSFFlux_9,e_PSFFlux_9,GauFlux_9,
+ e_GauFlux_9,GauSemi1_9,e_GauSemi1_9,GauSemi2_9,e_GauSemi2_9,
+ GauTheta_9,e_GauTheta_9,GauFWHMEff_9,P_9,e_P_9,AnP_9,e_AnP_9,
+ ApP_9,e_ApP_9,ApAnP_9,e_ApAnP_9,PUL_9,ApPUL_9,PStat_9,PX_9,
+ e_PX_18,E_PX_19,AnPX_9,e_AnPX_18,E_AnPX_19,Extended_9,
+ ExtVal_9,ERCSC_9,PCCS1_9,HRel_9,CirrusN_9,SkyBri_9,Zone,
+ ApFlux217_9,e_ApFlux217_9,ApFlux353_9,e_ApFlux353_9,
+ ApFlux545_9,e_ApFlux545_9
if (idig(ar__9(189:200)).EQ.0) e_GauSemi1_9 = rNULL__
if (idig(ar__9(213:224)).EQ.0) e_GauSemi2_9 = rNULL__
if (idig(ar__9(248:259)).EQ.0) P_9 = rNULL__
if (idig(ar__9(261:271)).EQ.0) e_P_9 = rNULL__
if (idig(ar__9(273:282)).EQ.0) AnP_9 = rNULL__
if (idig(ar__9(284:293)).EQ.0) e_AnP_9 = rNULL__
if (idig(ar__9(295:306)).EQ.0) ApP_9 = rNULL__
if (idig(ar__9(308:319)).EQ.0) e_ApP_9 = rNULL__
if (idig(ar__9(321:332)).EQ.0) ApAnP_9 = rNULL__
if (idig(ar__9(334:343)).EQ.0) e_ApAnP_9 = rNULL__
if (idig(ar__9(345:355)).EQ.0) PUL_9 = rNULL__
if (idig(ar__9(357:366)).EQ.0) ApPUL_9 = rNULL__
if (PStat_9 .EQ. 45) PStat_9 = iNULL__
if (idig(ar__9(370:380)).EQ.0) PX_9 = rNULL__
if (idig(ar__9(382:392)).EQ.0) e_PX_18 = rNULL__
if (idig(ar__9(394:404)).EQ.0) E_PX_19 = rNULL__
if (idig(ar__9(406:415)).EQ.0) AnPX_9 = rNULL__
if (idig(ar__9(417:425)).EQ.0) e_AnPX_18 = rNULL__
if (idig(ar__9(427:435)).EQ.0) E_AnPX_19 = rNULL__
if (HRel_9 .EQ. 45) HRel_9 = iNULL__
if (CirrusN_9 .EQ. 45) CirrusN_9 = iNULL__
if (idig(ar__9(497:506)).EQ.0) SkyBri_9 = rNULL__
if(ar__9(508:508) .EQ. '') Zone = iNULL__
if (idig(ar__9(510:521)).EQ.0) ApFlux217_9 = rNULL__
if (idig(ar__9(523:534)).EQ.0) e_ApFlux217_9 = rNULL__
if (idig(ar__9(536:547)).EQ.0) ApFlux353_9 = rNULL__
if (idig(ar__9(549:560)).EQ.0) e_ApFlux353_9 = rNULL__
if (idig(ar__9(562:573)).EQ.0) ApFlux545_9 = rNULL__
if (idig(ar__9(575:586)).EQ.0) e_ApFlux545_9 = rNULL__
c ..............Just test output...........
write(6,10)
+ Name_9,GLON_9,GLAT_9,RAdeg_9,DEdeg_9,DetFlux_9,e_DetFlux_9,
+ AperFlux_9,e_AperFlux_9,PSFFlux_9,e_PSFFlux_9,GauFlux_9,
+ e_GauFlux_9,GauSemi1_9,e_GauSemi1_9,GauSemi2_9,e_GauSemi2_9,
+ GauTheta_9,e_GauTheta_9,GauFWHMEff_9,P_9,e_P_9,AnP_9,e_AnP_9,
+ ApP_9,e_ApP_9,ApAnP_9,e_ApAnP_9,PUL_9,ApPUL_9,PStat_9,PX_9,
+ e_PX_18,E_PX_19,AnPX_9,e_AnPX_18,E_AnPX_19,Extended_9,
+ ExtVal_9,ERCSC_9,PCCS1_9,HRel_9,CirrusN_9,SkyBri_9,Zone,
+ ApFlux217_9,e_ApFlux217_9,ApFlux353_9,e_ApFlux353_9,
+ ApFlux545_9,e_ApFlux545_9
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'pccs143e.dat' ! Second Catalogue of Compact Sources, 143GHz,
* unknown reliability
C Format for file interpretation
11 format(
+ A24,1X,F11.7,1X,F11.7,1X,F11.7,1X,F11.7,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,F10.6,
+ 1X,F12.6,1X,F11.6,F12.6,1X,F5.1,1X,F4.1,1X,F10.6,1X,F12.6,1X,
+ F11.6,1X,F10.6,1X,F10.6,1X,F12.6,1X,F12.6,1X,F12.6,1X,F10.6,
+ 1X,F11.6,1X,F10.6,1X,I1,1X,F11.6,1X,F11.6,1X,F11.6,1X,F10.6,
+ 1X,F9.6,1X,F9.6,1X,I1,1X,I1,1X,A24,1X,A23,1X,I2,1X,I3,1X,
+ F10.6,1X,I1,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
+ E12.6)
C Effective file loading
open(unit=1,status='old',file=
+'pccs143e.dat')
write(6,*) '....Loading file: pccs143e.dat'
do i__=1,4139
read(1,'(A586)')ar__10
read(ar__10,11)
+ Name_10,GLON_10,GLAT_10,RAdeg_10,DEdeg_10,DetFlux_10,
+ e_DetFlux_10,AperFlux_10,e_AperFlux_10,PSFFlux_10,
+ e_PSFFlux_10,GauFlux_10,e_GauFlux_10,GauSemi1_10,
+ e_GauSemi1_10,GauSemi2_10,e_GauSemi2_10,GauTheta_10,
+ e_GauTheta_10,GauFWHMEff_10,P_10,e_P_10,AnP_10,e_AnP_10,
+ ApP_10,e_ApP_10,ApAnP_10,e_ApAnP_10,PUL_10,ApPUL_10,PStat_10,
+ PX_10,e_PX_20,E_PX_21,AnPX_10,e_AnPX_20,E_AnPX_21,Extended_10,
+ ExtVal_10,ERCSC_10,PCCS1_10,HRel_10,CirrusN_10,SkyBri_10,
+ Zone_1,ApFlux217_10,e_ApFlux217_10,ApFlux353_10,
+ e_ApFlux353_10,ApFlux545_10,e_ApFlux545_10
if (idig(ar__10(189:200)).EQ.0) e_GauSemi1_10 = rNULL__
if (idig(ar__10(213:224)).EQ.0) e_GauSemi2_10 = rNULL__
if (idig(ar__10(248:259)).EQ.0) P_10 = rNULL__
if (idig(ar__10(261:271)).EQ.0) e_P_10 = rNULL__
if (idig(ar__10(273:282)).EQ.0) AnP_10 = rNULL__
if (idig(ar__10(284:293)).EQ.0) e_AnP_10 = rNULL__
if (idig(ar__10(295:306)).EQ.0) ApP_10 = rNULL__
if (idig(ar__10(308:319)).EQ.0) e_ApP_10 = rNULL__
if (idig(ar__10(321:332)).EQ.0) ApAnP_10 = rNULL__
if (idig(ar__10(334:343)).EQ.0) e_ApAnP_10 = rNULL__
if (idig(ar__10(345:355)).EQ.0) PUL_10 = rNULL__
if (idig(ar__10(357:366)).EQ.0) ApPUL_10 = rNULL__
if (PStat_10 .EQ. 45) PStat_10 = iNULL__
if (idig(ar__10(370:380)).EQ.0) PX_10 = rNULL__
if (idig(ar__10(382:392)).EQ.0) e_PX_20 = rNULL__
if (idig(ar__10(394:404)).EQ.0) E_PX_21 = rNULL__
if (idig(ar__10(406:415)).EQ.0) AnPX_10 = rNULL__
if (idig(ar__10(417:425)).EQ.0) e_AnPX_20 = rNULL__
if (idig(ar__10(427:435)).EQ.0) E_AnPX_21 = rNULL__
if (HRel_10 .EQ. 45) HRel_10 = iNULL__
if (CirrusN_10 .EQ. 45) CirrusN_10 = iNULL__
if (idig(ar__10(497:506)).EQ.0) SkyBri_10 = rNULL__
if(ar__10(508:508) .EQ. '') Zone_1 = iNULL__
if (idig(ar__10(510:521)).EQ.0) ApFlux217_10 = rNULL__
if (idig(ar__10(523:534)).EQ.0) e_ApFlux217_10 = rNULL__
if (idig(ar__10(536:547)).EQ.0) ApFlux353_10 = rNULL__
if (idig(ar__10(549:560)).EQ.0) e_ApFlux353_10 = rNULL__
if (idig(ar__10(562:573)).EQ.0) ApFlux545_10 = rNULL__
if (idig(ar__10(575:586)).EQ.0) e_ApFlux545_10 = rNULL__
c ..............Just test output...........
write(6,11)
+ Name_10,GLON_10,GLAT_10,RAdeg_10,DEdeg_10,DetFlux_10,
+ e_DetFlux_10,AperFlux_10,e_AperFlux_10,PSFFlux_10,
+ e_PSFFlux_10,GauFlux_10,e_GauFlux_10,GauSemi1_10,
+ e_GauSemi1_10,GauSemi2_10,e_GauSemi2_10,GauTheta_10,
+ e_GauTheta_10,GauFWHMEff_10,P_10,e_P_10,AnP_10,e_AnP_10,
+ ApP_10,e_ApP_10,ApAnP_10,e_ApAnP_10,PUL_10,ApPUL_10,PStat_10,
+ PX_10,e_PX_20,E_PX_21,AnPX_10,e_AnPX_20,E_AnPX_21,Extended_10,
+ ExtVal_10,ERCSC_10,PCCS1_10,HRel_10,CirrusN_10,SkyBri_10,
+ Zone_1,ApFlux217_10,e_ApFlux217_10,ApFlux353_10,
+ e_ApFlux353_10,ApFlux545_10,e_ApFlux545_10
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'pccs217e.dat' ! Second Catalogue of Compact Sources, 217GHz,
* unknown reliability
C Format for file interpretation
12 format(
+ A24,1X,F11.7,1X,F11.7,1X,F11.7,1X,F11.7,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,F10.6,
+ 1X,F12.6,1X,F11.6,F12.6,1X,F5.1,1X,F4.1,1X,F10.6,1X,F12.6,1X,
+ F11.6,1X,F10.6,1X,F10.6,1X,F12.6,1X,F12.6,1X,F12.6,1X,F10.6,
+ 1X,F11.6,1X,F10.6,1X,I1,1X,F11.6,1X,F11.6,1X,F11.6,1X,F10.6,
+ 1X,F9.6,1X,F9.6,1X,I1,1X,I1,1X,A24,1X,A23,1X,I2,1X,I3,1X,
+ F10.6,1X,I1,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
+ E12.6)
C Effective file loading
open(unit=1,status='old',file=
+'pccs217e.dat')
write(6,*) '....Loading file: pccs217e.dat'
do i__=1,16842
read(1,'(A586)')ar__11
read(ar__11,12)
+ Name_11,GLON_11,GLAT_11,RAdeg_11,DEdeg_11,DetFlux_11,
+ e_DetFlux_11,AperFlux_11,e_AperFlux_11,PSFFlux_11,
+ e_PSFFlux_11,GauFlux_11,e_GauFlux_11,GauSemi1_11,
+ e_GauSemi1_11,GauSemi2_11,e_GauSemi2_11,GauTheta_11,
+ e_GauTheta_11,GauFWHMEff_11,P_11,e_P_11,AnP_11,e_AnP_11,
+ ApP_11,e_ApP_11,ApAnP_11,e_ApAnP_11,PUL_11,ApPUL_11,PStat_11,
+ PX_11,e_PX_22,E_PX_23,AnPX_11,e_AnPX_22,E_AnPX_23,Extended_11,
+ ExtVal_11,ERCSC_11,PCCS1_11,HRel_11,CirrusN_11,SkyBri_11,
+ Zone_2,ApFlux217_11,e_ApFlux217_11,ApFlux353_11,
+ e_ApFlux353_11,ApFlux545_11,e_ApFlux545_11
if (idig(ar__11(189:200)).EQ.0) e_GauSemi1_11 = rNULL__
if (idig(ar__11(213:224)).EQ.0) e_GauSemi2_11 = rNULL__
if (idig(ar__11(248:259)).EQ.0) P_11 = rNULL__
if (idig(ar__11(261:271)).EQ.0) e_P_11 = rNULL__
if (idig(ar__11(273:282)).EQ.0) AnP_11 = rNULL__
if (idig(ar__11(284:293)).EQ.0) e_AnP_11 = rNULL__
if (idig(ar__11(295:306)).EQ.0) ApP_11 = rNULL__
if (idig(ar__11(308:319)).EQ.0) e_ApP_11 = rNULL__
if (idig(ar__11(321:332)).EQ.0) ApAnP_11 = rNULL__
if (idig(ar__11(334:343)).EQ.0) e_ApAnP_11 = rNULL__
if (idig(ar__11(345:355)).EQ.0) PUL_11 = rNULL__
if (idig(ar__11(357:366)).EQ.0) ApPUL_11 = rNULL__
if (PStat_11 .EQ. 45) PStat_11 = iNULL__
if (idig(ar__11(370:380)).EQ.0) PX_11 = rNULL__
if (idig(ar__11(382:392)).EQ.0) e_PX_22 = rNULL__
if (idig(ar__11(394:404)).EQ.0) E_PX_23 = rNULL__
if (idig(ar__11(406:415)).EQ.0) AnPX_11 = rNULL__
if (idig(ar__11(417:425)).EQ.0) e_AnPX_22 = rNULL__
if (idig(ar__11(427:435)).EQ.0) E_AnPX_23 = rNULL__
if (HRel_11 .EQ. 45) HRel_11 = iNULL__
if (CirrusN_11 .EQ. 45) CirrusN_11 = iNULL__
if (idig(ar__11(497:506)).EQ.0) SkyBri_11 = rNULL__
if(ar__11(508:508) .EQ. '') Zone_2 = iNULL__
if (idig(ar__11(510:521)).EQ.0) ApFlux217_11 = rNULL__
if (idig(ar__11(523:534)).EQ.0) e_ApFlux217_11 = rNULL__
if (idig(ar__11(536:547)).EQ.0) ApFlux353_11 = rNULL__
if (idig(ar__11(549:560)).EQ.0) e_ApFlux353_11 = rNULL__
if (idig(ar__11(562:573)).EQ.0) ApFlux545_11 = rNULL__
if (idig(ar__11(575:586)).EQ.0) e_ApFlux545_11 = rNULL__
c ..............Just test output...........
write(6,12)
+ Name_11,GLON_11,GLAT_11,RAdeg_11,DEdeg_11,DetFlux_11,
+ e_DetFlux_11,AperFlux_11,e_AperFlux_11,PSFFlux_11,
+ e_PSFFlux_11,GauFlux_11,e_GauFlux_11,GauSemi1_11,
+ e_GauSemi1_11,GauSemi2_11,e_GauSemi2_11,GauTheta_11,
+ e_GauTheta_11,GauFWHMEff_11,P_11,e_P_11,AnP_11,e_AnP_11,
+ ApP_11,e_ApP_11,ApAnP_11,e_ApAnP_11,PUL_11,ApPUL_11,PStat_11,
+ PX_11,e_PX_22,E_PX_23,AnPX_11,e_AnPX_22,E_AnPX_23,Extended_11,
+ ExtVal_11,ERCSC_11,PCCS1_11,HRel_11,CirrusN_11,SkyBri_11,
+ Zone_2,ApFlux217_11,e_ApFlux217_11,ApFlux353_11,
+ e_ApFlux353_11,ApFlux545_11,e_ApFlux545_11
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'pccs353e.dat' ! Second Catalogue of Compact Sources, 353GHz,
* unknown reliability
C Format for file interpretation
13 format(
+ A24,1X,F11.7,1X,F11.7,1X,F11.7,1X,F11.7,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,F10.6,
+ 1X,F12.6,1X,F11.6,F12.6,1X,F5.1,1X,F4.1,1X,F10.6,1X,F12.6,1X,
+ F11.6,1X,F10.6,1X,F10.6,1X,F12.6,1X,F12.6,1X,F12.6,1X,F10.6,
+ 1X,F11.6,1X,F10.6,1X,I1,1X,F11.6,1X,F11.6,1X,F11.6,1X,F10.6,
+ 1X,F9.6,1X,F9.6,1X,I1,1X,I1,1X,A24,1X,A23,1X,I2,1X,I3,1X,
+ F10.6,1X,I1,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
+ E12.6)
C Effective file loading
open(unit=1,status='old',file=
+'pccs353e.dat')
write(6,*) '....Loading file: pccs353e.dat'
do i__=1,22665
read(1,'(A586)')ar__12
read(ar__12,13)
+ Name_12,GLON_12,GLAT_12,RAdeg_12,DEdeg_12,DetFlux_12,
+ e_DetFlux_12,AperFlux_12,e_AperFlux_12,PSFFlux_12,
+ e_PSFFlux_12,GauFlux_12,e_GauFlux_12,GauSemi1_12,
+ e_GauSemi1_12,GauSemi2_12,e_GauSemi2_12,GauTheta_12,
+ e_GauTheta_12,GauFWHMEff_12,P_12,e_P_12,AnP_12,e_AnP_12,
+ ApP_12,e_ApP_12,ApAnP_12,e_ApAnP_12,PUL_12,ApPUL_12,PStat_12,
+ PX_12,e_PX_24,E_PX_25,AnPX_12,e_AnPX_24,E_AnPX_25,Extended_12,
+ ExtVal_12,ERCSC_12,PCCS1_12,HRel_12,CirrusN_12,SkyBri_12,
+ Zone_3,ApFlux217_12,e_ApFlux217_12,ApFlux353_12,
+ e_ApFlux353_12,ApFlux545_12,e_ApFlux545_12
if (idig(ar__12(189:200)).EQ.0) e_GauSemi1_12 = rNULL__
if (idig(ar__12(213:224)).EQ.0) e_GauSemi2_12 = rNULL__
if (idig(ar__12(248:259)).EQ.0) P_12 = rNULL__
if (idig(ar__12(261:271)).EQ.0) e_P_12 = rNULL__
if (idig(ar__12(273:282)).EQ.0) AnP_12 = rNULL__
if (idig(ar__12(284:293)).EQ.0) e_AnP_12 = rNULL__
if (idig(ar__12(295:306)).EQ.0) ApP_12 = rNULL__
if (idig(ar__12(308:319)).EQ.0) e_ApP_12 = rNULL__
if (idig(ar__12(321:332)).EQ.0) ApAnP_12 = rNULL__
if (idig(ar__12(334:343)).EQ.0) e_ApAnP_12 = rNULL__
if (idig(ar__12(345:355)).EQ.0) PUL_12 = rNULL__
if (idig(ar__12(357:366)).EQ.0) ApPUL_12 = rNULL__
if (PStat_12 .EQ. 45) PStat_12 = iNULL__
if (idig(ar__12(370:380)).EQ.0) PX_12 = rNULL__
if (idig(ar__12(382:392)).EQ.0) e_PX_24 = rNULL__
if (idig(ar__12(394:404)).EQ.0) E_PX_25 = rNULL__
if (idig(ar__12(406:415)).EQ.0) AnPX_12 = rNULL__
if (idig(ar__12(417:425)).EQ.0) e_AnPX_24 = rNULL__
if (idig(ar__12(427:435)).EQ.0) E_AnPX_25 = rNULL__
if (HRel_12 .EQ. 45) HRel_12 = iNULL__
if (CirrusN_12 .EQ. 45) CirrusN_12 = iNULL__
if (idig(ar__12(497:506)).EQ.0) SkyBri_12 = rNULL__
if(ar__12(508:508) .EQ. '') Zone_3 = iNULL__
if (idig(ar__12(510:521)).EQ.0) ApFlux217_12 = rNULL__
if (idig(ar__12(523:534)).EQ.0) e_ApFlux217_12 = rNULL__
if (idig(ar__12(536:547)).EQ.0) ApFlux353_12 = rNULL__
if (idig(ar__12(549:560)).EQ.0) e_ApFlux353_12 = rNULL__
if (idig(ar__12(562:573)).EQ.0) ApFlux545_12 = rNULL__
if (idig(ar__12(575:586)).EQ.0) e_ApFlux545_12 = rNULL__
c ..............Just test output...........
write(6,13)
+ Name_12,GLON_12,GLAT_12,RAdeg_12,DEdeg_12,DetFlux_12,
+ e_DetFlux_12,AperFlux_12,e_AperFlux_12,PSFFlux_12,
+ e_PSFFlux_12,GauFlux_12,e_GauFlux_12,GauSemi1_12,
+ e_GauSemi1_12,GauSemi2_12,e_GauSemi2_12,GauTheta_12,
+ e_GauTheta_12,GauFWHMEff_12,P_12,e_P_12,AnP_12,e_AnP_12,
+ ApP_12,e_ApP_12,ApAnP_12,e_ApAnP_12,PUL_12,ApPUL_12,PStat_12,
+ PX_12,e_PX_24,E_PX_25,AnPX_12,e_AnPX_24,E_AnPX_25,Extended_12,
+ ExtVal_12,ERCSC_12,PCCS1_12,HRel_12,CirrusN_12,SkyBri_12,
+ Zone_3,ApFlux217_12,e_ApFlux217_12,ApFlux353_12,
+ e_ApFlux353_12,ApFlux545_12,e_ApFlux545_12
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'pccs545e.dat' ! Second Catalogue of Compact Sources, 545GHz,
* unknown reliability
C Format for file interpretation
14 format(
+ A24,1X,F11.7,1X,F11.7,1X,F11.7,1X,F11.7,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,F10.6,
+ 1X,F12.6,1X,F11.6,F12.6,1X,F5.1,1X,F4.1,1X,F10.6,1X,F12.6,1X,
+ F11.6,1X,F10.6,1X,F10.6,1X,F12.6,1X,F12.6,1X,F12.6,1X,F10.6,
+ 1X,F11.6,1X,F10.6,1X,I1,1X,F11.6,1X,F11.6,1X,F11.6,1X,F10.6,
+ 1X,F9.6,1X,F9.6,1X,I1,1X,I1,1X,A24,1X,A23,1X,I2,1X,I3,1X,
+ F10.6,1X,I1,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
+ E12.6)
C Effective file loading
open(unit=1,status='old',file=
+'pccs545e.dat')
write(6,*) '....Loading file: pccs545e.dat'
do i__=1,31068
read(1,'(A586)')ar__13
read(ar__13,14)
+ Name_13,GLON_13,GLAT_13,RAdeg_13,DEdeg_13,DetFlux_13,
+ e_DetFlux_13,AperFlux_13,e_AperFlux_13,PSFFlux_13,
+ e_PSFFlux_13,GauFlux_13,e_GauFlux_13,GauSemi1_13,
+ e_GauSemi1_13,GauSemi2_13,e_GauSemi2_13,GauTheta_13,
+ e_GauTheta_13,GauFWHMEff_13,P_13,e_P_13,AnP_13,e_AnP_13,
+ ApP_13,e_ApP_13,ApAnP_13,e_ApAnP_13,PUL_13,ApPUL_13,PStat_13,
+ PX_13,e_PX_26,E_PX_27,AnPX_13,e_AnPX_26,E_AnPX_27,Extended_13,
+ ExtVal_13,ERCSC_13,PCCS1_13,HRel_13,CirrusN_13,SkyBri_13,
+ Zone_4,ApFlux217_13,e_ApFlux217_13,ApFlux353_13,
+ e_ApFlux353_13,ApFlux545_13,e_ApFlux545_13
if (idig(ar__13(189:200)).EQ.0) e_GauSemi1_13 = rNULL__
if (idig(ar__13(213:224)).EQ.0) e_GauSemi2_13 = rNULL__
if (idig(ar__13(248:259)).EQ.0) P_13 = rNULL__
if (idig(ar__13(261:271)).EQ.0) e_P_13 = rNULL__
if (idig(ar__13(273:282)).EQ.0) AnP_13 = rNULL__
if (idig(ar__13(284:293)).EQ.0) e_AnP_13 = rNULL__
if (idig(ar__13(295:306)).EQ.0) ApP_13 = rNULL__
if (idig(ar__13(308:319)).EQ.0) e_ApP_13 = rNULL__
if (idig(ar__13(321:332)).EQ.0) ApAnP_13 = rNULL__
if (idig(ar__13(334:343)).EQ.0) e_ApAnP_13 = rNULL__
if (idig(ar__13(345:355)).EQ.0) PUL_13 = rNULL__
if (idig(ar__13(357:366)).EQ.0) ApPUL_13 = rNULL__
if (PStat_13 .EQ. 45) PStat_13 = iNULL__
if (idig(ar__13(370:380)).EQ.0) PX_13 = rNULL__
if (idig(ar__13(382:392)).EQ.0) e_PX_26 = rNULL__
if (idig(ar__13(394:404)).EQ.0) E_PX_27 = rNULL__
if (idig(ar__13(406:415)).EQ.0) AnPX_13 = rNULL__
if (idig(ar__13(417:425)).EQ.0) e_AnPX_26 = rNULL__
if (idig(ar__13(427:435)).EQ.0) E_AnPX_27 = rNULL__
if (HRel_13 .EQ. 45) HRel_13 = iNULL__
if (CirrusN_13 .EQ. 45) CirrusN_13 = iNULL__
if (idig(ar__13(497:506)).EQ.0) SkyBri_13 = rNULL__
if(ar__13(508:508) .EQ. '') Zone_4 = iNULL__
if (idig(ar__13(510:521)).EQ.0) ApFlux217_13 = rNULL__
if (idig(ar__13(523:534)).EQ.0) e_ApFlux217_13 = rNULL__
if (idig(ar__13(536:547)).EQ.0) ApFlux353_13 = rNULL__
if (idig(ar__13(549:560)).EQ.0) e_ApFlux353_13 = rNULL__
if (idig(ar__13(562:573)).EQ.0) ApFlux545_13 = rNULL__
if (idig(ar__13(575:586)).EQ.0) e_ApFlux545_13 = rNULL__
c ..............Just test output...........
write(6,14)
+ Name_13,GLON_13,GLAT_13,RAdeg_13,DEdeg_13,DetFlux_13,
+ e_DetFlux_13,AperFlux_13,e_AperFlux_13,PSFFlux_13,
+ e_PSFFlux_13,GauFlux_13,e_GauFlux_13,GauSemi1_13,
+ e_GauSemi1_13,GauSemi2_13,e_GauSemi2_13,GauTheta_13,
+ e_GauTheta_13,GauFWHMEff_13,P_13,e_P_13,AnP_13,e_AnP_13,
+ ApP_13,e_ApP_13,ApAnP_13,e_ApAnP_13,PUL_13,ApPUL_13,PStat_13,
+ PX_13,e_PX_26,E_PX_27,AnPX_13,e_AnPX_26,E_AnPX_27,Extended_13,
+ ExtVal_13,ERCSC_13,PCCS1_13,HRel_13,CirrusN_13,SkyBri_13,
+ Zone_4,ApFlux217_13,e_ApFlux217_13,ApFlux353_13,
+ e_ApFlux353_13,ApFlux545_13,e_ApFlux545_13
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'pccs857e.dat' ! Second Catalogue of Compact Sources, 857GHz,
* unknown reliability (correct file, 08-Dec-2023)
C Format for file interpretation
15 format(
+ A24,1X,F11.7,1X,F11.7,1X,F11.7,1X,F11.7,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,F10.6,
+ 1X,F12.6,1X,F11.6,F12.6,1X,F5.1,1X,F4.1,1X,F10.6,1X,F12.6,1X,
+ F11.6,1X,F10.6,1X,F10.6,1X,F12.6,1X,F12.6,1X,F12.6,1X,F10.6,
+ 1X,F11.6,1X,F10.6,1X,I1,1X,F11.6,1X,F11.6,1X,F11.6,1X,F10.6,
+ 1X,F9.6,1X,F9.6,1X,I1,1X,I1,1X,A24,1X,A23,1X,I2,1X,I3,1X,
+ F10.6,1X,I1,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
+ E12.6)
C Effective file loading
open(unit=1,status='old',file=
+'pccs857e.dat')
write(6,*) '....Loading file: pccs857e.dat'
do i__=1,43290
read(1,'(A586)')ar__14
read(ar__14,15)
+ Name_14,GLON_14,GLAT_14,RAdeg_14,DEdeg_14,DetFlux_14,
+ e_DetFlux_14,AperFlux_14,e_AperFlux_14,PSFFlux_14,
+ e_PSFFlux_14,GauFlux_14,e_GauFlux_14,GauSemi1_14,
+ e_GauSemi1_14,GauSemi2_14,e_GauSemi2_14,GauTheta_14,
+ e_GauTheta_14,GauFWHMEff_14,P_14,e_P_14,AnP_14,e_AnP_14,
+ ApP_14,e_ApP_14,ApAnP_14,e_ApAnP_14,PUL_14,ApPUL_14,PStat_14,
+ PX_14,e_PX_28,E_PX_29,AnPX_14,e_AnPX_28,E_AnPX_29,Extended_14,
+ ExtVal_14,ERCSC_14,PCCS1_14,HRel_14,CirrusN_14,SkyBri_14,
+ Zone_5,ApFlux217_14,e_ApFlux217_14,ApFlux353_14,
+ e_ApFlux353_14,ApFlux545_14,e_ApFlux545_14
if (idig(ar__14(189:200)).EQ.0) e_GauSemi1_14 = rNULL__
if (idig(ar__14(213:224)).EQ.0) e_GauSemi2_14 = rNULL__
if (idig(ar__14(248:259)).EQ.0) P_14 = rNULL__
if (idig(ar__14(261:271)).EQ.0) e_P_14 = rNULL__
if (idig(ar__14(273:282)).EQ.0) AnP_14 = rNULL__
if (idig(ar__14(284:293)).EQ.0) e_AnP_14 = rNULL__
if (idig(ar__14(295:306)).EQ.0) ApP_14 = rNULL__
if (idig(ar__14(308:319)).EQ.0) e_ApP_14 = rNULL__
if (idig(ar__14(321:332)).EQ.0) ApAnP_14 = rNULL__
if (idig(ar__14(334:343)).EQ.0) e_ApAnP_14 = rNULL__
if (idig(ar__14(345:355)).EQ.0) PUL_14 = rNULL__
if (idig(ar__14(357:366)).EQ.0) ApPUL_14 = rNULL__
if (PStat_14 .EQ. 45) PStat_14 = iNULL__
if (idig(ar__14(370:380)).EQ.0) PX_14 = rNULL__
if (idig(ar__14(382:392)).EQ.0) e_PX_28 = rNULL__
if (idig(ar__14(394:404)).EQ.0) E_PX_29 = rNULL__
if (idig(ar__14(406:415)).EQ.0) AnPX_14 = rNULL__
if (idig(ar__14(417:425)).EQ.0) e_AnPX_28 = rNULL__
if (idig(ar__14(427:435)).EQ.0) E_AnPX_29 = rNULL__
if (HRel_14 .EQ. 45) HRel_14 = iNULL__
if (CirrusN_14 .EQ. 45) CirrusN_14 = iNULL__
if (idig(ar__14(497:506)).EQ.0) SkyBri_14 = rNULL__
if(ar__14(508:508) .EQ. '') Zone_5 = iNULL__
if (idig(ar__14(510:521)).EQ.0) ApFlux217_14 = rNULL__
if (idig(ar__14(523:534)).EQ.0) e_ApFlux217_14 = rNULL__
if (idig(ar__14(536:547)).EQ.0) ApFlux353_14 = rNULL__
if (idig(ar__14(549:560)).EQ.0) e_ApFlux353_14 = rNULL__
if (idig(ar__14(562:573)).EQ.0) ApFlux545_14 = rNULL__
if (idig(ar__14(575:586)).EQ.0) e_ApFlux545_14 = rNULL__
c ..............Just test output...........
write(6,15)
+ Name_14,GLON_14,GLAT_14,RAdeg_14,DEdeg_14,DetFlux_14,
+ e_DetFlux_14,AperFlux_14,e_AperFlux_14,PSFFlux_14,
+ e_PSFFlux_14,GauFlux_14,e_GauFlux_14,GauSemi1_14,
+ e_GauSemi1_14,GauSemi2_14,e_GauSemi2_14,GauTheta_14,
+ e_GauTheta_14,GauFWHMEff_14,P_14,e_P_14,AnP_14,e_AnP_14,
+ ApP_14,e_ApP_14,ApAnP_14,e_ApAnP_14,PUL_14,ApPUL_14,PStat_14,
+ PX_14,e_PX_28,E_PX_29,AnPX_14,e_AnPX_28,E_AnPX_29,Extended_14,
+ ExtVal_14,ERCSC_14,PCCS1_14,HRel_14,CirrusN_14,SkyBri_14,
+ Zone_5,ApFlux217_14,e_ApFlux217_14,ApFlux353_14,
+ e_ApFlux353_14,ApFlux545_14,e_ApFlux545_14
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