CDS.headStuff2
FORTRAN Generation
Conversion of standardized ReadMe file for
catalog 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-26 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/ApJS/198/1 Photometry catalogs for the Lockman Hole (Fotopoulou+, 2012) *================================================================================ *Photometry and photometric redshift catalogs for the Lockman Hole deep field. * Fotopoulou S., Salvato M., Hasinger G., Rovilos E., Brusa M., Egami E., * Lutz D., Burwitz V., Henry J.P., Huang J.H., Rigopoulou D., Vaccari M. * <Astrophys. J. Suppl. Ser., 198, 1 (2012)> * =2012ApJS..198....1F C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table5.dat' ! Lockman Hole Photometric catalog integer*4 nr__ parameter (nr__=187611) ! Number of records character*576 ar__ ! Full-size record C J2000 position composed of: RAdeg DEdeg integer*4 ID (nr__) ! [1,206786] Catalog identifier * (<[FSH2012] NNNNNN> in Simbad) real*8 RAdeg (nr__) ! (deg) Optical Right Ascension (J2000) real*8 DEdeg (nr__) ! (deg) Optical Declination (J2000) real*8 RAGdeg (nr__) ! (deg) ?=-99 GALEX Right Ascension (J2000) real*8 DEGdeg (nr__) ! (deg) ?=-99 GALEX Declination (J2000) real*8 RAUdeg (nr__) ! (deg) ?=-99 UKIDSS Right Ascension (J2000) real*8 DEUdeg (nr__) ! (deg) ?=-99 UKIDSS Declination (J2000) real*8 RASdeg (nr__) ! (deg) ?=-99 SDSS Right Ascension (J2000) real*8 DESdeg (nr__) ! (deg) ?=-99 SDSS Declination (J2000) real*8 RAIdeg (nr__) ! (deg) ?=-99 IRAC Right Ascension (J2000) real*8 DEIdeg (nr__) ! (deg) ?=-99 IRAC Declination (J2000) real*8 FUVmag (nr__) ! (mag) ?=-99 GALEX FUV band AB magnitude (1) real*8 e_FUVmag (nr__) ! (mag) ?=-99 Error in FUVmag real*8 NUVmag (nr__) ! (mag) ?=-99 GALEX NUV band AB magnitude (1) real*8 e_NUVmag (nr__) ! (mag) ?=-99 Error in NUVmag real*8 Umag (nr__) ! (mag) ?=-99 LBT U band AB magnitude (1) real*8 e_Umag (nr__) ! (mag) ?=-99 Error in Umag real*8 Bmag (nr__) ! (mag) ?=-99 LBT B band AB magnitude (1) real*8 e_Bmag (nr__) ! (mag) ?=-99 Error in Bmag real*8 Vmag (nr__) ! (mag) ?=-99 LBT V band AB magnitude (1) real*8 e_Vmag (nr__) ! (mag) ?=-99 Error in Vmag real*8 z_Lmag (nr__) ! (mag) ?=-99 LBT z' band AB magnitude (1) real*8 e_z_Lmag (nr__) ! (mag) ?=-99 Error in z'LBTmag real*8 Ymag (nr__) ! (mag) ?=-99 LBT Y band AB magnitude (1) real*8 e_Ymag (nr__) ! (mag) ?=-99 Error in Ymag real*8 Rcmag (nr__) ! (mag) ?=-99 Subaru R_C_ band magnitude (1) real*8 e_Rcmag (nr__) ! (mag) ?=-99 Error in Rmag real*8 Icmag (nr__) ! (mag) ?=-99 Subaru I_C_ band magnitude (1) real*8 e_Icmag (nr__) ! (mag) ?=-99 Error in Imag real*8 z_Smag (nr__) ! (mag) ?=-99 Subaru z' band magnitude (1) real*8 e_z_Smag (nr__) ! (mag) ?=-99 Error in z'Submag real*8 u_mag (nr__) ! (mag) ?=-99 SDSS u' band magnitude (1) real*8 e_u_mag (nr__) ! (mag) ?=-99 Error in u'mag real*8 g_mag (nr__) ! (mag) ?=-99 SDSS g' band magnitude (1) real*8 e_g_mag (nr__) ! (mag) ?=-99 Error in g'mag real*8 r_mag (nr__) ! (mag) ?=-99 SDSS r' band magnitude (1) real*8 e_r_mag (nr__) ! (mag) ?=-99 Error in r'mag real*8 i_mag (nr__) ! (mag) ?=-99 SDSS i' band magnitude (1) real*8 e_i_mag (nr__) ! (mag) ?=-99 Error in i'mag real*8 z_mag (nr__) ! (mag) ?=-99 SDSS z' band magnitude (1) real*8 e_z_mag (nr__) ! (mag) ?=-99 Error in z'mag real*8 Jmag (nr__) ! (mag) ?=-99 UKIDSS J band magnitude (1) real*8 e_Jmag (nr__) ! (mag) ?=-99 Error in Jmag real*8 Kmag (nr__) ! (mag) ?=-99 UKIDSS K band magnitude (1) real*8 e_Kmag (nr__) ! (mag) ?=-99 Error in Kmag real*8 v_3_6_ (nr__) ! (mag) ?=-99 Spitzer/IRAC 3.6 micron band magnitude (1) real*8 e__3_6_ (nr__) ! (mag) ?=-99 Error in [3.6] real*8 v_4_5_ (nr__) ! (mag) ?=-99 Spitzer/IRAC 4.5 micron band magnitude (1) real*8 e__4_5_ (nr__) ! (mag) ?=-99 Error in [4.5] real*8 v_5_8_ (nr__) ! (mag) ?=-99 Spitzer/IRAC 5.8 micron band magnitude (1) real*8 e__5_8_ (nr__) ! (mag) ?=-99 Error in [5.8] real*8 v_8_0_ (nr__) ! (mag) ?=-99 Spitzer/IRAC 8.0 micron band magnitude (1) real*8 e__8_0_ (nr__) ! (mag) ?=-99 Error in [8.0] integer*4 Det (nr__) ! [1,6] Detection flag (2) integer*4 f_FUVmag (nr__) ! Photometry flag on FUVmag (3) integer*4 f_NUVmag (nr__) ! Photometry flag on NUVmag (3) integer*4 f_Umag (nr__) ! Photometry flag on Umag (3) integer*4 f_Bmag (nr__) ! Photometry flag on Bmag (3) integer*4 f_Vmag (nr__) ! Photometry flag on Vmag (3) integer*4 f_z_Lmag (nr__) ! Photometry flag on zLBTmag (3) integer*4 f_Ymag (nr__) ! Photometry flag on Ymag (3) integer*4 f_Rcmag (nr__) ! Photometry flag on Rmag (3) integer*4 f_Icmag (nr__) ! Photometry flag on Imag (3) integer*4 f_z_Smag (nr__) ! Photometry flag on z'Smag (3) integer*4 f__3_6_ (nr__) ! Photometry flag on [3.6] (3) integer*4 f__4_5_ (nr__) ! Photometry flag on [4.5] (3) integer*4 f__5_8_ (nr__) ! Photometry flag on [5.8] (3) integer*4 f__8_0_ (nr__) ! Photometry flag on [8.0] (3) integer*4 N (nr__) ! [0/1] Neighbor flag (1: neighbor within 1.5" or O) real*4 Class (nr__) ! [0,1] SExtractor star/galaxy classification * (1: star, or 0: galaxy) of the corresponding * detection image real*8 z (nr__) ! ?=-99 Spectroscopic redshift (4) integer*4 r_z (nr__) ! [1,27]? Reference for z (see refs.dat file) integer*4 Var (nr__) ! ?=-99 Variability flag (5) integer*4 Xray (nr__) ! X-ray detection flag (6) integer*4 XID (nr__) ! [5,2838]?=-99 X-ray identification (XMMNR number * for 44 sources) (7) integer*4 Morph (nr__) ! ?=-99 Morphology (8) *Note (1): LBT = Large Binocular Telescope. * Aperture photometry corrected to total (only for GALEX and IRAC * photometry) for point-like sources. *Note (2): Detection flag as follows: * 1 = R_c_ detection (5{sigma}); * 2 = z' detection (5{sigma}); * 3 = B detection (5{sigma}); * 4 = R_c_ detection (3{sigma}), * 5 = z' detection (3{sigma}); * 6 = B detection (3{sigma}). *Note (3): Complementary to the flag provided by SExtractor indicating * saturation, we masked problematic regions on the images which include * bad pixels and problematic areas close to stars, using 'Weight * Watcher' (Marmo & Bertin, 2008ASPC..394..619M). For more details, * refer to the associated description file of the catalog. * Photometry flag as follows: * -99 = 99 value produced by SExtractor either for the magnitude or the error, * or source outside of the field * -5 = Flag in the UKIDSS catalog marking noise in the JHK bands * -4 = Saturation or incomplete/corrupted data produced by SExtractor * -3 = Source inside a stripe, as marked by optical inspection of the images * -2 = Magnitude error was negative * -1 = Magnitude error was greater than 1 * 0 = Everything is OK * 1 = FWHM in the detection band is zero * 2 = Source inside wings of stars, potentially fake * 3 = Magnitude greater than the detection limit *Note (4): We use only the high quality spectroscopic redshifts from the * original catalogs according to Table 4. *Note (5): Sources with |Z_var,k_{+/-}{delta}Z_var,k_|>0.2 carry a flag 1, * while the rest carry a flag 0. *Note (6): Sources detected in the X-rays are flagged as 1, sources inside the * XMM area without X-ray detection are flagged as 0, while sources * outside the XMM area are flagged as -99. *Note (7): The sources identified as X-ray counterparts will have the * corresponding X-ray ID number from the XMM catalog (Brunner et al. * 2008, Cat. J/A+A/479/283, <[BCH2008] NNNN> for Simbad). *Note (8): As discussed in detail in Section 3.2.3, this column gives the * merged morphology classification information from the HST and ground * based images, provided only for the counterparts of the X-ray * sources. Morphology code as follows: * -2 = too faint; * -1 = photometry blended with nearby source; * 0 = extended source; * 1 = point-like source. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table10.dat' ! Lockman Hole Photometric Redshift catalog integer*4 nr__1 parameter (nr__1=187611) ! Number of records character*137 ar__1 ! Full-size record integer*4 ID_1 (nr__1) ! [1,206786] Catalog identifier * (<[FSH2012] NNNNNN> in Simbad) real*8 zph (nr__1) ! ?=-99 Best photometric redshift solution real*8 b_zph (nr__1) ! ?=-99 Lowest redshift at 68% significance real*8 B_zph_1 (nr__1) ! ?=-99 Highest redshift at 68% significance real*4 b_90z (nr__1) ! [0,6.8] Lowest redshift at 90% significance real*4 B_90z_1 (nr__1) ! [0.01,7] Highest redshift at 90% significance real*4 chi2z (nr__1) ! ?=-99 Lowest {chi}^2^ for best galaxy model real*8 pr_z (nr__1) ! (%) ? Probability zph is correct photometric redshift integer*4 mod_z (nr__1) ! [1,31]?=-99 Best model number (9) integer*4 Ext_z (nr__1) ! [0,4] Extinction-law used to compute zph * (unexplained) real*4 E_B_V_1 (nr__1) ! (mag) ?=-99 E(B-V) used to compute zph integer*4 Nband (nr__1) ! [1,21] Number of bands used in SED fit real*8 z2 (nr__1) ! ?=-99 Second best photometric redshift solution real*4 chi2z2 (nr__1) ! ?=-99 Lowest {chi}^2^ for second best galaxy * model real*4 pr_z2 (nr__1) ! (%) ? Probability z2 is correct photometric redshift integer*4 mod_z2 (nr__1) ! [1,31]?=-99 Second best model number (9) real*4 E_B_V_2 (nr__1) ! ?=-99 E(B-V) used to compute z2 real*4 chi2_ (nr__1) ! ?=-99 Lowest {chi}^2^ for best star model integer*4 cl (nr__1) ! [0/1] Star flag (1: star or 0: galaxy); * see section 5.3 *Note (9): Models from Ilbert et al. (2009ApJ...690.1236I) and * Salvato et al. (2009ApJ...690.1250S). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'refs.dat' ! Spectroscopic redshift references integer*4 nr__2 parameter (nr__2=27) ! Number of records character*86 ar__2 ! Full-size record integer*4 Ref (nr__2) ! Code reference character*19 BibCode (nr__2) ! bibcode character*17 Aut (nr__2) ! First author's name character*45 Comm (nr__2) ! Comment C============================================================================= C Loading file 'table5.dat' ! Lockman Hole Photometric catalog C Format for file interpretation 1 format( + I6,1X,F9.5,1X,F8.5,1X,F9.5,1X,F9.5,1X,F9.5,1X,F9.5,1X,F9.5,1X, + F9.5,1X,F9.5,1X,F9.5,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4, + 1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X, + F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4, + 1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X, + F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4, + 1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,I1, + 1X,I3,1X,I3,1X,I3,1X,I3,1X,I3,1X,I3,1X,I3,1X,I2,1X,I2,1X,I2, + 1X,I3,1X,I3,1X,I3,1X,I3,1X,I1,1X,F5.3,1X,F8.4,2X,I2,1X,I3,1X, + I3,1X,I4,1X,I3) C Effective file loading open(unit=1,status='old',file= +'table5.dat') write(6,*) '....Loading file: table5.dat' do i__=1,187611 read(1,'(A576)')ar__ read(ar__,1) + ID(i__),RAdeg(i__),DEdeg(i__),RAGdeg(i__),DEGdeg(i__), + RAUdeg(i__),DEUdeg(i__),RASdeg(i__),DESdeg(i__),RAIdeg(i__), + DEIdeg(i__),FUVmag(i__),e_FUVmag(i__),NUVmag(i__), + e_NUVmag(i__),Umag(i__),e_Umag(i__),Bmag(i__),e_Bmag(i__), + Vmag(i__),e_Vmag(i__),z_Lmag(i__),e_z_Lmag(i__),Ymag(i__), + e_Ymag(i__),Rcmag(i__),e_Rcmag(i__),Icmag(i__),e_Icmag(i__), + z_Smag(i__),e_z_Smag(i__),u_mag(i__),e_u_mag(i__),g_mag(i__), + e_g_mag(i__),r_mag(i__),e_r_mag(i__),i_mag(i__),e_i_mag(i__), + z_mag(i__),e_z_mag(i__),Jmag(i__),e_Jmag(i__),Kmag(i__), + e_Kmag(i__),v_3_6_(i__),e__3_6_(i__),v_4_5_(i__),e__4_5_(i__), + v_5_8_(i__),e__5_8_(i__),v_8_0_(i__),e__8_0_(i__),Det(i__), + f_FUVmag(i__),f_NUVmag(i__),f_Umag(i__),f_Bmag(i__), + f_Vmag(i__),f_z_Lmag(i__),f_Ymag(i__),f_Rcmag(i__), + f_Icmag(i__),f_z_Smag(i__),f__3_6_(i__),f__4_5_(i__), + f__5_8_(i__),f__8_0_(i__),N(i__),Class(i__),z(i__),r_z(i__), + Var(i__),Xray(i__),XID(i__),Morph(i__) if(ar__(558:559) .EQ. '') r_z(i__) = iNULL__ c ..............Just test output........... write(6,1) + ID(i__),RAdeg(i__),DEdeg(i__),RAGdeg(i__),DEGdeg(i__), + RAUdeg(i__),DEUdeg(i__),RASdeg(i__),DESdeg(i__),RAIdeg(i__), + DEIdeg(i__),FUVmag(i__),e_FUVmag(i__),NUVmag(i__), + e_NUVmag(i__),Umag(i__),e_Umag(i__),Bmag(i__),e_Bmag(i__), + Vmag(i__),e_Vmag(i__),z_Lmag(i__),e_z_Lmag(i__),Ymag(i__), + e_Ymag(i__),Rcmag(i__),e_Rcmag(i__),Icmag(i__),e_Icmag(i__), + z_Smag(i__),e_z_Smag(i__),u_mag(i__),e_u_mag(i__),g_mag(i__), + e_g_mag(i__),r_mag(i__),e_r_mag(i__),i_mag(i__),e_i_mag(i__), + z_mag(i__),e_z_mag(i__),Jmag(i__),e_Jmag(i__),Kmag(i__), + e_Kmag(i__),v_3_6_(i__),e__3_6_(i__),v_4_5_(i__),e__4_5_(i__), + v_5_8_(i__),e__5_8_(i__),v_8_0_(i__),e__8_0_(i__),Det(i__), + f_FUVmag(i__),f_NUVmag(i__),f_Umag(i__),f_Bmag(i__), + f_Vmag(i__),f_z_Lmag(i__),f_Ymag(i__),f_Rcmag(i__), + f_Icmag(i__),f_z_Smag(i__),f__3_6_(i__),f__4_5_(i__), + f__5_8_(i__),f__8_0_(i__),N(i__),Class(i__),z(i__),r_z(i__), + Var(i__),Xray(i__),XID(i__),Morph(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table10.dat' ! Lockman Hole Photometric Redshift catalog C Format for file interpretation 2 format( + I6,1X,F8.4,1X,F7.3,1X,F7.3,1X,F5.3,1X,F5.3,1X,E13.6,1X,F7.3, + 1X,I3,1X,I1,1X,F6.2,1X,I2,1X,F7.3,1X,E13.6,1X,F6.3,1X,I3,1X, + F6.2,1X,E13.6,1X,I1) C Effective file loading open(unit=1,status='old',file= +'table10.dat') write(6,*) '....Loading file: table10.dat' do i__=1,187611 read(1,'(A137)')ar__1 read(ar__1,2) + ID_1(i__),zph(i__),b_zph(i__),B_zph_1(i__),b_90z(i__), + B_90z_1(i__),chi2z(i__),pr_z(i__),mod_z(i__),Ext_z(i__), + E_B_V_1(i__),Nband(i__),z2(i__),chi2z2(i__),pr_z2(i__), + mod_z2(i__),E_B_V_2(i__),chi2_(i__),cl(i__) if(ar__1(59:65) .EQ. '') pr_z(i__) = rNULL__ if(ar__1(105:110) .EQ. '') pr_z2(i__) = rNULL__ c ..............Just test output........... write(6,2) + ID_1(i__),zph(i__),b_zph(i__),B_zph_1(i__),b_90z(i__), + B_90z_1(i__),chi2z(i__),pr_z(i__),mod_z(i__),Ext_z(i__), + E_B_V_1(i__),Nband(i__),z2(i__),chi2z2(i__),pr_z2(i__), + mod_z2(i__),E_B_V_2(i__),chi2_(i__),cl(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'refs.dat' ! Spectroscopic redshift references C Format for file interpretation 3 format(I2,1X,A19,1X,A17,1X,A45) C Effective file loading open(unit=1,status='old',file= +'refs.dat') write(6,*) '....Loading file: refs.dat' do i__=1,27 read(1,'(A86)')ar__2 read(ar__2,3)Ref(i__),BibCode(i__),Aut(i__),Comm(i__) c ..............Just test output........... write(6,3)Ref(i__),BibCode(i__),Aut(i__),Comm(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end