Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/415/2336 into FORTRAN code for loading all data files into arrays.
Note that special values are assigned to unknown or unspecified
numbers (also called NULL numbers);
when necessary, the coordinate components making up the right ascension
and declination are converted into floating-point numbers
representing these angles in degrees.
program load_ReadMe
C=============================================================================
C F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2024-Apr-16
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/MNRAS/415/2336 Herschel-ATLAS Science Demonstration Catalog (Rigby+, 2011)
*================================================================================
*Herschel-ATLAS: first data release of the Science Demonstration Phase source
*catalogues.
* Rigby E.E., Maddox S.J., Dunne L., Negrello M., Smith D.J.B.,
* Gonzalez-Nuevo J., Herranz D., Lopez-Caniego M., Auld R., Buttiglione S.,
* Baes M., Cava A., Cooray A., Clements D.L., Dariush A., De Zotti G.,
* Dye S., Eales S., Frayer D., Fritz J., Hopwood R., Ibar E., Ivison R.J.,
* Jarvis M., Panuzzo P., Pascale E., Pohlen M., Rodighiero G., Serjeant S.,
* Temi P., Thompson M.A.
* <Mon. Not. R. Astron. Soc., 415, 2336-2348 (2011)>
* =2011MNRAS.415.2336R
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hatlas.dat' ! Spire catalogue
integer*4 nr__
parameter (nr__=6876) ! Number of records
character*1164 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*16 HATLAS (nr__) ! HATLAS IAU identifier from the 250um
* position (JHHMMSS.s+DDMMSS)
integer*4 SDP (nr__) ! [1/23791] SDP catalogue identifier
real*8 RAdeg (nr__) ! (deg) Right ascension (J2000), determined from
* the 250um data
real*8 DEdeg (nr__) ! (deg) Declination (J2000), determined from
* the 250um data
real*4 F250P (nr__) ! (Jy) Herschel/SPIRE 250um point source flux
real*4 F350P (nr__) ! (Jy) Herschel/SPIRE 350um point source flux
real*4 F500P (nr__) ! (Jy) Herschel/SPIRE 500um point source flux
real*4 e_F250P (nr__) ! (Jy) F250P 1-{sigma} error
real*4 e_F350P (nr__) ! (Jy) F350P 1-{sigma} error
real*4 e_F500P (nr__) ! (Jy) F500P 1-{sigma} error
real*4 F250B (nr__) ! (Jy) Herschel/SPIRE best 250um flux (5)
real*4 F350B (nr__) ! (Jy) Herschel/SPIRE best 350um flux (5)
real*4 F500B (nr__) ! (Jy) Herschel/SPIRE best 500um flux (5)
real*4 e_F250B (nr__) ! (Jy) F250B 1-{sigma} flux error
real*4 e_F350B (nr__) ! (Jy) F350B 1-{sigma} flux error
real*4 e_F500B (nr__) ! (Jy) F500B 1-{sigma} flux error
real*8 a_F250B (nr__) ! (arcsec) ?=-99 Aperture radius used for F250B (5)
real*8 a_F350B (nr__) ! (arcsec) ?=-99 Aperture radius used for F350B (5)
real*4 a_F500B (nr__) ! (arcsec) ?=-99 Aperture radius used for F500B (5)
real*4 F100P (nr__) ! (Jy) ?=-99 Herschel/PACS 100um point source flux,
* aperture corrected in 10" radius aperture (1)
real*4 F160P (nr__) ! (Jy) ?=-99 Herschel/PACS 160um point source flux,
* aperture corrected in 15" radius aperture (1)
real*4 e_F100P (nr__) ! (Jy) ?=-1 F100P 1-{sigma} flux (1)
real*4 e_F160P (nr__) ! (Jy) ?=-1 F100P 1-{sigma} flux (1)
real*4 F100B (nr__) ! (Jy) ?=-99 Herschel/PACS best 100um flux (5)
real*4 F160B (nr__) ! (Jy) ?=-99 Herschel/PACS best 160um flux (5)
real*4 e_F100B (nr__) ! (Jy) ?=-1 Best 100um 1 {sigma} flux error (1)
real*4 e_F160B (nr__) ! (Jy) ?=-1 Best 160um 1 {sigma} flux error (1)
real*4 a_F100B (nr__) ! (arcsec) ?=-1 Aperture radius used for F100B (1)
real*4 a_F160B (nr__) ! (arcsec) ?=-1 Aperture radius used for F160B (1)
integer*8 OBJID (nr__) ! ? SDSS DR7 objID (Cat. II/294)
real*8 RAsdeg (nr__) ! (deg) ?=-99 SDSS DR7 Right Ascension (J2000)
real*8 DEsdeg (nr__) ! (deg) ?=-99 SDSS DR7 Declination (J2000)
real*8 rmag (nr__) ! (mag) ?=-99 SDSS DR7 r-band model mag
real*8 e_rmag (nr__) ! (mag) ?=-99 SDSS DR7 r-band model mag error
real*8 IsoA (nr__) ! (arcsec) ?=-99 SDSS DR7 isophotal major axis, used for
* extended source SPIRE aperture fluxes
real*8 sepS (nr__) ! (arcsec) ?=-99 Separation between the 250um and
* SDSS DR7 positions
real*4 LR (nr__) ! ?=-99 Likelihood ratio
real*8 Rel (nr__) ! [0/1]?=-99 Reliability
integer*4 S_G (nr__) ! [0/1]?=-99 Star(1)/Galaxy(0) class
real*8 zph (nr__) ! ?=-99 Photometric redshifts from ANNz
* (Smith et al., 2011MNRAS.416..857S)
real*8 e_zph (nr__) ! ?=-99 Error on zph (2011MNRAS.416..857S)
character*9 n_zph (nr__) ! Bands used to derive zph
integer*4 f_zph (nr__) ! [0/511]?=-99 Integer representing n_zph
* (u=1, g=2, r=4, ... H=128, K=256)
integer*4 o_zph (nr__) ! ?=-99 Number of bands used in deriving photoz
real*8 zsp (nr__) ! ?=-99 The spectroscopic redshift
integer*4 q_zsp (nr__) ! [0/5]?=-99 quality flag of zsp (only use
* q_zsp>= 3, which have >=90% reliability)
integer*4 r_zsp (nr__) ! Flag indicating which the sources of zsp (2)
character*11 IRAS (nr__) ! Name of associated IRAS source in
* the IIFSCz, Cat. J/MNRAS/398/109
real*8 RAideg (nr__) ! (deg) ?=-99 IIFSCz Right ascension
real*8 DEideg (nr__) ! (deg) ?=-99 IIFSCz Declination
integer*4 fpos (nr__) ! [1/5]?=-99 Indicates the source of the
* position used in the IIFSCz (e.g. SDSS DR6,
* radio, IRAS, etc)
real*8 S12 (nr__) ! (Jy) ?=-99 IIFSCz 12um flux or upper limit
real*8 S25 (nr__) ! (Jy) ?=-99 IIFSCz 25um flux or upper limit
real*8 S60 (nr__) ! (Jy) ?=-99 IIFSCz 60um flux or upper limit
real*8 S100 (nr__) ! (Jy) ?=-99 IIFSCz 100um flux or upper limit
integer*4 q_S12 (nr__) ! ?=-99 Quality flag of IIFSCz S12 (3)
integer*4 q_S25 (nr__) ! ?=-99 Quality flag of IIFSCz S25 (3)
integer*4 q_S60 (nr__) ! ?=-99 Quality flag of IIFSCz S60 (3)
integer*4 q_S100 (nr__) ! ?=-99 Quality flag of IIFSCz S100 (3)
real*8 sepI (nr__) ! (arcsec) ?=-99 Separation between IIFSCz and 250um
* source positions
character*19 GAMA_IAU (nr__) ! GAMA IAU identifier, JHHMMSS.ss+DDMMSS.s
* Cat. J/MNRAS/413/971
integer*4 GAMA (nr__) ! ? GAMA catalogue index (GAMA_ID)
real*8 RAgdeg (nr__) ! (deg) ?=-99 GAMA Right ascension (J2000)
real*8 DEgdeg (nr__) ! (deg) ?=-99 GAMA Declination (J2000)
real*8 rpmag (nr__) ! (mag) ?=-99 Petrosian r-band mag from GAMA
real*8 rext (nr__) ! (mag) ?=-99 GAMA derived r-band extinction
real*8 rSmag (nr__) ! (mag) ?=-99 GAMA r-band sersic mag out to 10 Re
integer*4 FUV (nr__) ! (mag) [-99] Extinction-corrected Far UV GALEX mag
* (currently just a placeholder)
integer*4 NUV (nr__) ! (mag) [-99] Extinction-corrected Near UV GALEX mag
* (currently just a placeholder)
integer*4 e_FUV (nr__) ! (mag) [-99] Error on Far UV GALEX mag
* (currently just a placeholder)
integer*4 e_NUV (nr__) ! (mag) [-99] Error on Near UV GALEX mag
* (currently just a placeholder)
real*8 ukron (nr__) ! (mag) ?=-99 Extinction-corrected GAMA Kron u mag
* (either sd or rd depending on phsrc column)
real*8 gkron (nr__) ! (mag) ?=-99 Extinction-corrected GAMA Kron g mag
* (either sd or rd depending on phsrc column)
real*8 rkron (nr__) ! (mag) ?=-99 Extinction-corrected GAMA Kron r mag
* (either sd or rd depending on phsrc column)
real*8 ikron (nr__) ! (mag) ?=-99 Extinction-corrected GAMA Kron i mag
* (either sd or rd depending on phsrc column)
real*8 zkron (nr__) ! (mag) ?=-99 Extinction-corrected GAMA Kron z mag
* (either sd or rd depending on phsrc column)
real*8 ykron (nr__) ! (mag) ?=-99 Extinction-corrected GAMA Kron y mag
* (either sd or rd depending on phsrc column)
real*8 jkron (nr__) ! (mag) ?=-99 Extinction-corrected GAMA Kron j mag
* (either sd or rd depending on phsrc column)
real*8 hkron (nr__) ! (mag) ?=-99 Extinction-corrected GAMA Kron h mag
* (either sd or rd depending on phsrc column)
real*8 kkron (nr__) ! (mag) ?=-99 Extinction-corrected GAMA Kron k mag
* (either sd or rd depending on phsrc column)
real*8 e_ukron (nr__) ! (mag) ?=-99 GAMA Kron u mag error (4)
real*8 e_gkron (nr__) ! (mag) ?=-99 GAMA Kron g mag error (4)
real*8 e_rkron (nr__) ! (mag) ?=-99 GAMA Kron r mag error (4)
real*8 e_ikron (nr__) ! (mag) ?=-99 GAMA Kron i mag error (4)
real*8 e_zkron (nr__) ! (mag) ?=-99 GAMA Kron z mag error (4)
real*8 e_ykron (nr__) ! (mag) ?=-99 GAMA Kron y mag error (4)
real*8 e_jkron (nr__) ! (mag) ?=-99 GAMA Kron j mag error (4)
real*8 e_hkron (nr__) ! (mag) ?=-99 GAMA Kron h mag error (4)
real*8 e_kkron (nr__) ! (mag) ?=-99 GAMA Kron k mag error (4)
integer*8 LAS8 (nr__) ! ? Source ID for UKIDSS-LAS detections in
* YJHK (Cat. II/314)
character*2 phsrc (nr__) ! [rd sd] rd: defined, or sd: self-defined
*Note (1): A value of -99 means that a source wasn't detected.
* A value of -1 indicates a source that falls outside the PACS region.
*Note (2): additive flags with:
* 1 = SDSS
* 2 = 6dFGS
* 4 = 2SLAQ_QSO survey
* 8 = 2SLAQ_LRG survey
* 16 = GAMA
* (the GAMA redshift this will always take precedence over the other surveys)
*Note (3): Quality flag as:
* 1 = upper limit
* 2 = low quality
* 3 = good quality
*Note (4): corrected for factor of 4 in Hill et al., 2010MNRAS.404.1215H.
*Note (5): the "best" flux is identical to the point source flux for
* non-extended sources (its error is then -99); an aperture flux
* is computed for extended sources.
C=============================================================================
C Loading file 'hatlas.dat' ! Spire catalogue
C Format for file interpretation
1 format(
+ 7X,A16,1X,I5,1X,F9.5,1X,F9.6,1X,E22.15,1X,E22.15,1X,E22.15,1X,
+ E24.17,1X,E22.15,1X,E22.15,1X,E22.15,1X,E22.15,1X,E22.15,1X,
+ E22.15,1X,E22.15,1X,E22.15,1X,F10.6,1X,F10.6,1X,F5.1,1X,
+ E22.15,1X,E22.15,1X,E22.15,1X,E22.15,1X,E22.15,1X,E22.15,1X,
+ E22.15,1X,E22.15,1X,E22.15,1X,E22.15,1X,I18,1X,F10.5,1X,F9.5,
+ 1X,F8.4,1X,F10.6,1X,F10.6,1X,F10.6,1X,E12.6,1X,F10.6,1X,I3,1X,
+ F10.6,1X,F10.6,1X,A9,1X,I4,1X,I3,1X,F10.6,1X,I3,1X,I3,1X,A11,
+ 1X,F10.6,1X,F10.6,1X,I3,1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,I3,
+ 1X,I3,1X,I3,1X,I3,1X,F10.6,5X,A19,1X,I6,1X,F7.3,1X,F9.3,1X,
+ F10.6,1X,F10.6,1X,F10.6,1X,I3,1X,I3,1X,I3,1X,I3,1X,F10.6,1X,
+ F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,
+ 1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,
+ F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,I12,1X,A2)
C Effective file loading
open(unit=1,status='old',file=
+'hatlas.dat')
write(6,*) '....Loading file: hatlas.dat'
do i__=1,6876
read(1,'(A1164)')ar__
read(ar__,1)
+ HATLAS(i__),SDP(i__),RAdeg(i__),DEdeg(i__),F250P(i__),
+ F350P(i__),F500P(i__),e_F250P(i__),e_F350P(i__),e_F500P(i__),
+ F250B(i__),F350B(i__),F500B(i__),e_F250B(i__),e_F350B(i__),
+ e_F500B(i__),a_F250B(i__),a_F350B(i__),a_F500B(i__),
+ F100P(i__),F160P(i__),e_F100P(i__),e_F160P(i__),F100B(i__),
+ F160B(i__),e_F100B(i__),e_F160B(i__),a_F100B(i__),
+ a_F160B(i__),OBJID(i__),RAsdeg(i__),DEsdeg(i__),rmag(i__),
+ e_rmag(i__),IsoA(i__),sepS(i__),LR(i__),Rel(i__),S_G(i__),
+ zph(i__),e_zph(i__),n_zph(i__),f_zph(i__),o_zph(i__),zsp(i__),
+ q_zsp(i__),r_zsp(i__),IRAS(i__),RAideg(i__),DEideg(i__),
+ fpos(i__),S12(i__),S25(i__),S60(i__),S100(i__),q_S12(i__),
+ q_S25(i__),q_S60(i__),q_S100(i__),sepI(i__),GAMA_IAU(i__),
+ GAMA(i__),RAgdeg(i__),DEgdeg(i__),rpmag(i__),rext(i__),
+ rSmag(i__),FUV(i__),NUV(i__),e_FUV(i__),e_NUV(i__),ukron(i__),
+ gkron(i__),rkron(i__),ikron(i__),zkron(i__),ykron(i__),
+ jkron(i__),hkron(i__),kkron(i__),e_ukron(i__),e_gkron(i__),
+ e_rkron(i__),e_ikron(i__),e_zkron(i__),e_ykron(i__),
+ e_jkron(i__),e_hkron(i__),e_kkron(i__),LAS8(i__),phsrc(i__)
if(ar__(587:604) .EQ. '') OBJID(i__) = iNULL__
if(ar__(878:883) .EQ. '') GAMA(i__) = iNULL__
if(ar__(1150:1161) .EQ. '') LAS8(i__) = iNULL__
c ..............Just test output...........
write(6,1)
+ HATLAS(i__),SDP(i__),RAdeg(i__),DEdeg(i__),F250P(i__),
+ F350P(i__),F500P(i__),e_F250P(i__),e_F350P(i__),e_F500P(i__),
+ F250B(i__),F350B(i__),F500B(i__),e_F250B(i__),e_F350B(i__),
+ e_F500B(i__),a_F250B(i__),a_F350B(i__),a_F500B(i__),
+ F100P(i__),F160P(i__),e_F100P(i__),e_F160P(i__),F100B(i__),
+ F160B(i__),e_F100B(i__),e_F160B(i__),a_F100B(i__),
+ a_F160B(i__),OBJID(i__),RAsdeg(i__),DEsdeg(i__),rmag(i__),
+ e_rmag(i__),IsoA(i__),sepS(i__),LR(i__),Rel(i__),S_G(i__),
+ zph(i__),e_zph(i__),n_zph(i__),f_zph(i__),o_zph(i__),zsp(i__),
+ q_zsp(i__),r_zsp(i__),IRAS(i__),RAideg(i__),DEideg(i__),
+ fpos(i__),S12(i__),S25(i__),S60(i__),S100(i__),q_S12(i__),
+ q_S25(i__),q_S60(i__),q_S100(i__),sepI(i__),GAMA_IAU(i__),
+ GAMA(i__),RAgdeg(i__),DEgdeg(i__),rpmag(i__),rext(i__),
+ rSmag(i__),FUV(i__),NUV(i__),e_FUV(i__),e_NUV(i__),ukron(i__),
+ gkron(i__),rkron(i__),ikron(i__),zkron(i__),ykron(i__),
+ jkron(i__),hkron(i__),kkron(i__),e_ukron(i__),e_gkron(i__),
+ e_rkron(i__),e_ikron(i__),e_zkron(i__),e_ykron(i__),
+ e_jkron(i__),e_hkron(i__),e_kkron(i__),LAS8(i__),phsrc(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end