Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/517/A44 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-Mar-28 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/517/A44 Spitzer sources in the Galactic bulge (Uttenthaler+, 2010) *================================================================================ *Galactic bulge giants: probing stellar and galactic evolution. *I. Catalogue of Spitzer IRAC and MIPS sources. * Uttenthaler S., Stute M., Sahai R., Blommaert J.A.D.L., Schultheis M., * Kraemer K.E., Groenewegen M.A.T., Price S.D. * <Astron. Astrophys. 517, A44 (2010)> * =2010A&A...517A..44U C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'fits.dat' ! Summary of FITS images integer*4 nr__ parameter (nr__=31) ! Number of records character*46 ar__ ! Full-size record C J2000 position composed of: RAdeg DEdeg character*11 FileName (nr__) ! Name of FITS file, in "fig" subdirectory character*1 Ch (nr__) ! [1324M] Channel integer*4 Xpix (nr__) ! (pix) X image size integer*4 Ypix (nr__) ! (pix) Y image size real*8 RAdeg (nr__) ! (deg) Right Ascension (J2000) of image centre real*8 DEdeg (nr__) ! (deg) Declination (J2000) of image centre c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'b1.dat' ! Catalogue of point sources in the Bulge 1 field integer*4 nr__1 parameter (nr__1=69047) ! Number of records character*466 ar__1 ! Full-size record C J2000 position composed of: RAdeg DEdeg character*18 USB (nr__1) ! Source identification (from position) real*8 RAdeg_1 (nr__1) ! (deg) Right ascension in degrees (J2000) real*8 DEdeg_1 (nr__1) ! (deg) Declination in degrees (J2000) real*4 IRAC1 (nr__1) ! (uJy) ?=-9.99 IRAC 1 (3.6um) flux real*4 IRAC2 (nr__1) ! (uJy) ?=-9.99 IRAC 2 (4.5um) flux real*4 IRAC3 (nr__1) ! (uJy) ?=-9.99 IRAC 3 (5.8um) flux real*4 IRAC4 (nr__1) ! (uJy) ?=-9.99 IRAC 4 (8.0um) flux real*4 MIPS (nr__1) ! (uJy) ?=-9.99 MIPS 24 (24um) flux character*1 q_IRAC1 (nr__1) ! [ABCD-] Quality flag of IRAC 1 (Sect.2.2.2) character*1 q_IRAC2 (nr__1) ! [ABCD-] Quality flag of IRAC 2 (Sect.2.2.2) character*1 q_IRAC3 (nr__1) ! [ABCD-] Quality flag of IRAC 3 (Sect.2.2.2) character*1 q_IRAC4 (nr__1) ! [ABCD-] Quality flag of IRAC 4 (Sect.2.2.2) character*1 q_MIPS (nr__1) ! [ABCD-] Quality flag of MIPS 24 (Sect.2.2.2) real*4 e_IRAC1 (nr__1) ! (uJy) ?=-9.99 Uncertainty IRAC 1 flux real*4 e_IRAC2 (nr__1) ! (uJy) ?=-9.99 Uncertainty IRAC 2 flux real*4 e_IRAC3 (nr__1) ! (uJy) ?=-9.99 Uncertainty IRAC 3 flux real*4 e_IRAC4 (nr__1) ! (uJy) ?=-9.99 Uncertainty IRAC 4 flux real*4 e_MIPS (nr__1) ! (uJy) ?=-9.99 Uncertainty MIPS 24 flux integer*4 f_IRAC1 (nr__1) ! [0/1] Outside/Inside IRAC 1 field of view integer*4 f_IRAC2 (nr__1) ! [0/1] Outside/Inside IRAC 2 field of view integer*4 f_IRAC3 (nr__1) ! [0/1] Outside/Inside IRAC 3 field of view integer*4 f_IRAC4 (nr__1) ! [0/1] Outside/Inside IRAC 4 field of view integer*4 f_MIPS (nr__1) ! [0/1] Outside/Inside MIPS 24 field of view real*4 dDEN (nr__1) ! (arcsec) ? Distance to closest DENIS source real*4 Imag (nr__1) ! (mag) ? I-band magnitude of DENIS counterpart real*4 Jmag (nr__1) ! (mag) ? J-band magnitude of DENIS counterpart real*4 Kmag (nr__1) ! (mag) ? K-band magnitude of DENIS counterpart real*4 d2M (nr__1) ! (arcsec) ? Distance to closest 2MASS source real*4 J2mag (nr__1) ! (mag) ?=- (?=-9.999) J-band magnitude of * 2MASS counterpart real*4 H2mag (nr__1) ! (mag) ?=- (?=-9.999) H-band magnitude of * 2MASS counterpart real*4 K2mag (nr__1) ! (mag) ?=- (?=-9.999) K-band magnitude of * 2MASS counterpart character*1 q_J2mag (nr__1) ! [ABCDEFU-] 2MASS J-band quality flag character*1 q_H2mag (nr__1) ! [ABCDEFU-] 2MASS H-band quality flag character*1 q_K2mag (nr__1) ! [ABCDEFU-] 2MASS K-band quality flag real*4 dIRAS (nr__1) ! (arcsec) ? Distance to closest IRAS source real*4 F12 (nr__1) ! (Jy) ?=- (?=-9.9) IRAS 12 micron flux integer*4 q_F12 (nr__1) ! [1/3]?=- IRAS 12 micron quality flag real*4 F25 (nr__1) ! (Jy) ?=- (?=-9.9) IRAS 25 micron flux integer*4 q_F25 (nr__1) ! [1/3]?=- IRAS 25 micron quality flag real*4 F60 (nr__1) ! (Jy) ?=- (?=-9.9) IRAS 60 micron flux integer*4 q_F60 (nr__1) ! [1/3]?=- IRAS 60 micron quality flag real*4 F100 (nr__1) ! (Jy) ?=- (?=-9.9) IRAS 100 micron flux integer*4 q_F100 (nr__1) ! [1/3]?=- IRAS 100 micron quality flag real*4 dMSX (nr__1) ! (arcsec) ? Distance to closest MSX source real*4 B1 (nr__1) ! (Jy) ?=- (?=-9.99) MSX B1 (4.29um) band flux integer*4 q_B1 (nr__1) ! [0/4]?=- MSX B1 band quality flag real*4 B2 (nr__1) ! (Jy) ?=- (?=-9.99) MSX B2 (4.25um) band flux integer*4 q_B2 (nr__1) ! [0/4]?=- MSX B2 band quality flag real*4 A (nr__1) ! (Jy) ?=- (?=-9.99) MSX A (8.28um) band flux integer*4 q_A (nr__1) ! [0/4]?=- MSX A band quality flag real*4 C (nr__1) ! (Jy) ?=- (?=-9.99) MSX C (12.13um) band flux integer*4 q_C (nr__1) ! [0/4]?=- MSX C band quality flag real*4 D (nr__1) ! (Jy) ?=- (?=-9.99) MSX D (14.65um) band flux integer*4 q_D (nr__1) ! [0/4]?=- MSX D band quality flag real*4 E (nr__1) ! (Jy) ?=- (?=-9.99) MSX E (21.34um) band flux integer*4 q_E (nr__1) ! [0/4]?=- MSX E band quality flag c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'b2.dat' ! Catalogue of point sources in the Bulge 2 field integer*4 nr__2 parameter (nr__2=60262) ! Number of records character*466 ar__2 ! Full-size record C J2000 position composed of: RAdeg DEdeg character*18 USB_1 (nr__2) ! Source identification (from position) real*8 RAdeg_2 (nr__2) ! (deg) Right ascension in degrees (J2000) real*8 DEdeg_2 (nr__2) ! (deg) Declination in degrees (J2000) real*4 IRAC1_1 (nr__2) ! (uJy) ?=-9.99 IRAC 1 (3.6um) flux real*4 IRAC2_1 (nr__2) ! (uJy) ?=-9.99 IRAC 2 (4.5um) flux real*4 IRAC3_1 (nr__2) ! (uJy) ?=-9.99 IRAC 3 (5.8um) flux real*4 IRAC4_1 (nr__2) ! (uJy) ?=-9.99 IRAC 4 (8.0um) flux real*4 MIPS_1 (nr__2) ! (uJy) ?=-9.99 MIPS 24 (24um) flux character*1 q_IRAC1_1 (nr__2) ! [ABCD-] Quality flag of IRAC 1 (Sect.2.2.2) character*1 q_IRAC2_1 (nr__2) ! [ABCD-] Quality flag of IRAC 2 (Sect.2.2.2) character*1 q_IRAC3_1 (nr__2) ! [ABCD-] Quality flag of IRAC 3 (Sect.2.2.2) character*1 q_IRAC4_1 (nr__2) ! [ABCD-] Quality flag of IRAC 4 (Sect.2.2.2) character*1 q_MIPS_1 (nr__2) ! [ABCD-] Quality flag of MIPS 24 (Sect.2.2.2) real*4 e_IRAC1_1 (nr__2) ! (uJy) ?=-9.99 Uncertainty IRAC 1 flux real*4 e_IRAC2_1 (nr__2) ! (uJy) ?=-9.99 Uncertainty IRAC 2 flux real*4 e_IRAC3_1 (nr__2) ! (uJy) ?=-9.99 Uncertainty IRAC 3 flux real*4 e_IRAC4_1 (nr__2) ! (uJy) ?=-9.99 Uncertainty IRAC 4 flux real*4 e_MIPS_1 (nr__2) ! (uJy) ?=-9.99 Uncertainty MIPS 24 flux integer*4 f_IRAC1_1 (nr__2) ! [0/1] Outside/Inside IRAC 1 field of view integer*4 f_IRAC2_1 (nr__2) ! [0/1] Outside/Inside IRAC 2 field of view integer*4 f_IRAC3_1 (nr__2) ! [0/1] Outside/Inside IRAC 3 field of view integer*4 f_IRAC4_1 (nr__2) ! [0/1] Outside/Inside IRAC 4 field of view integer*4 f_MIPS_1 (nr__2) ! [0/1] Outside/Inside MIPS 24 field of view real*4 dDEN_1 (nr__2) ! (arcsec) ? Distance to closest DENIS source real*4 Imag_1 (nr__2) ! (mag) ? I-band magnitude of DENIS counterpart real*4 Jmag_1 (nr__2) ! (mag) ? J-band magnitude of DENIS counterpart real*4 Kmag_1 (nr__2) ! (mag) ? K-band magnitude of DENIS counterpart real*4 d2M_1 (nr__2) ! (arcsec) ? Distance to closest 2MASS source real*4 J2mag_1 (nr__2) ! (mag) ?=- (?=-9.999) J-band magnitude of * 2MASS counterpart real*4 H2mag_1 (nr__2) ! (mag) ?=- (?=-9.999) H-band magnitude of * 2MASS counterpart real*4 K2mag_1 (nr__2) ! (mag) ?=- (?=-9.999) K-band magnitude of * 2MASS counterpart character*1 q_J2mag_1 (nr__2) ! [ABCDEFU-] 2MASS J-band quality flag character*1 q_H2mag_1 (nr__2) ! [ABCDEFU-] 2MASS H-band quality flag character*1 q_K2mag_1 (nr__2) ! [ABCDEFU-] 2MASS K-band quality flag real*4 dIRAS_1 (nr__2) ! (arcsec) ? Distance to closest IRAS source real*4 F12_1 (nr__2) ! (Jy) ?=- (?=-9.9) IRAS 12 micron flux integer*4 q_F12_1 (nr__2) ! [1/3]?=- IRAS 12 micron quality flag real*4 F25_1 (nr__2) ! (Jy) ?=- (?=-9.9) IRAS 25 micron flux integer*4 q_F25_1 (nr__2) ! [1/3]?=- IRAS 25 micron quality flag real*4 F60_1 (nr__2) ! (Jy) ?=- (?=-9.9) IRAS 60 micron flux integer*4 q_F60_1 (nr__2) ! [1/3]?=- IRAS 60 micron quality flag real*4 F100_1 (nr__2) ! (Jy) ?=- (?=-9.9) IRAS 100 micron flux integer*4 q_F100_1 (nr__2) ! [1/3]?=- IRAS 100 micron quality flag real*4 dMSX_1 (nr__2) ! (arcsec) ? Distance to closest MSX source real*4 B1_1 (nr__2) ! (Jy) ?=- (?=-9.99) MSX B1 (4.29um) band flux integer*4 q_B1_1 (nr__2) ! [0/4]?=- MSX B1 band quality flag real*4 B2_1 (nr__2) ! (Jy) ?=- (?=-9.99) MSX B2 (4.25um) band flux integer*4 q_B2_1 (nr__2) ! [0/4]?=- MSX B2 band quality flag real*4 A_1 (nr__2) ! (Jy) ?=- (?=-9.99) MSX A (8.28um) band flux integer*4 q_A_1 (nr__2) ! [0/4]?=- MSX A band quality flag real*4 C_1 (nr__2) ! (Jy) ?=- (?=-9.99) MSX C (12.13um) band flux integer*4 q_C_1 (nr__2) ! [0/4]?=- MSX C band quality flag real*4 D_1 (nr__2) ! (Jy) ?=- (?=-9.99) MSX D (14.65um) band flux integer*4 q_D_1 (nr__2) ! [0/4]?=- MSX D band quality flag real*4 E_1 (nr__2) ! (Jy) ?=- (?=-9.99) MSX E (21.34um) band flux integer*4 q_E_1 (nr__2) ! [0/4]?=- MSX E band quality flag c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'b3.dat' ! Catalogue of point sources in the Bulge 3 field integer*4 nr__3 parameter (nr__3=52576) ! Number of records character*466 ar__3 ! Full-size record C J2000 position composed of: RAdeg DEdeg character*18 USB_2 (nr__3) ! Source identification (from position) real*8 RAdeg_3 (nr__3) ! (deg) Right ascension in degrees (J2000) real*8 DEdeg_3 (nr__3) ! (deg) Declination in degrees (J2000) real*4 IRAC1_2 (nr__3) ! (uJy) ?=-9.99 IRAC 1 (3.6um) flux real*4 IRAC2_2 (nr__3) ! (uJy) ?=-9.99 IRAC 2 (4.5um) flux real*4 IRAC3_2 (nr__3) ! (uJy) ?=-9.99 IRAC 3 (5.8um) flux real*4 IRAC4_2 (nr__3) ! (uJy) ?=-9.99 IRAC 4 (8.0um) flux real*4 MIPS_2 (nr__3) ! (uJy) ?=-9.99 MIPS 24 (24um) flux character*1 q_IRAC1_2 (nr__3) ! [ABCD-] Quality flag of IRAC 1 (Sect.2.2.2) character*1 q_IRAC2_2 (nr__3) ! [ABCD-] Quality flag of IRAC 2 (Sect.2.2.2) character*1 q_IRAC3_2 (nr__3) ! [ABCD-] Quality flag of IRAC 3 (Sect.2.2.2) character*1 q_IRAC4_2 (nr__3) ! [ABCD-] Quality flag of IRAC 4 (Sect.2.2.2) character*1 q_MIPS_2 (nr__3) ! [ABCD-] Quality flag of MIPS 24 (Sect.2.2.2) real*4 e_IRAC1_2 (nr__3) ! (uJy) ?=-9.99 Uncertainty IRAC 1 flux real*4 e_IRAC2_2 (nr__3) ! (uJy) ?=-9.99 Uncertainty IRAC 2 flux real*4 e_IRAC3_2 (nr__3) ! (uJy) ?=-9.99 Uncertainty IRAC 3 flux real*4 e_IRAC4_2 (nr__3) ! (uJy) ?=-9.99 Uncertainty IRAC 4 flux real*4 e_MIPS_2 (nr__3) ! (uJy) ?=-9.99 Uncertainty MIPS 24 flux integer*4 f_IRAC1_2 (nr__3) ! [0/1] Outside/Inside IRAC 1 field of view integer*4 f_IRAC2_2 (nr__3) ! [0/1] Outside/Inside IRAC 2 field of view integer*4 f_IRAC3_2 (nr__3) ! [0/1] Outside/Inside IRAC 3 field of view integer*4 f_IRAC4_2 (nr__3) ! [0/1] Outside/Inside IRAC 4 field of view integer*4 f_MIPS_2 (nr__3) ! [0/1] Outside/Inside MIPS 24 field of view real*4 dDEN_2 (nr__3) ! (arcsec) ? Distance to closest DENIS source real*4 Imag_2 (nr__3) ! (mag) ? I-band magnitude of DENIS counterpart real*4 Jmag_2 (nr__3) ! (mag) ? J-band magnitude of DENIS counterpart real*4 Kmag_2 (nr__3) ! (mag) ? K-band magnitude of DENIS counterpart real*4 d2M_2 (nr__3) ! (arcsec) ? Distance to closest 2MASS source real*4 J2mag_2 (nr__3) ! (mag) ?=- (?=-9.999) J-band magnitude of * 2MASS counterpart real*4 H2mag_2 (nr__3) ! (mag) ?=- (?=-9.999) H-band magnitude of * 2MASS counterpart real*4 K2mag_2 (nr__3) ! (mag) ?=- (?=-9.999) K-band magnitude of * 2MASS counterpart character*1 q_J2mag_2 (nr__3) ! [ABCDEFU-] 2MASS J-band quality flag character*1 q_H2mag_2 (nr__3) ! [ABCDEFU-] 2MASS H-band quality flag character*1 q_K2mag_2 (nr__3) ! [ABCDEFU-] 2MASS K-band quality flag real*4 dIRAS_2 (nr__3) ! (arcsec) ? Distance to closest IRAS source real*4 F12_2 (nr__3) ! (Jy) ?=- (?=-9.9) IRAS 12 micron flux integer*4 q_F12_2 (nr__3) ! [1/3]?=- IRAS 12 micron quality flag real*4 F25_2 (nr__3) ! (Jy) ?=- (?=-9.9) IRAS 25 micron flux integer*4 q_F25_2 (nr__3) ! [1/3]?=- IRAS 25 micron quality flag real*4 F60_2 (nr__3) ! (Jy) ?=- (?=-9.9) IRAS 60 micron flux integer*4 q_F60_2 (nr__3) ! [1/3]?=- IRAS 60 micron quality flag real*4 F100_2 (nr__3) ! (Jy) ?=- (?=-9.9) IRAS 100 micron flux integer*4 q_F100_2 (nr__3) ! [1/3]?=- IRAS 100 micron quality flag real*4 dMSX_2 (nr__3) ! (arcsec) ? Distance to closest MSX source real*4 B1_2 (nr__3) ! (Jy) ?=- (?=-9.99) MSX B1 (4.29um) band flux integer*4 q_B1_2 (nr__3) ! [0/4]?=- MSX B1 band quality flag real*4 B2_2 (nr__3) ! (Jy) ?=- (?=-9.99) MSX B2 (4.25um) band flux integer*4 q_B2_2 (nr__3) ! [0/4]?=- MSX B2 band quality flag real*4 A_2 (nr__3) ! (Jy) ?=- (?=-9.99) MSX A (8.28um) band flux integer*4 q_A_2 (nr__3) ! [0/4]?=- MSX A band quality flag real*4 C_2 (nr__3) ! (Jy) ?=- (?=-9.99) MSX C (12.13um) band flux integer*4 q_C_2 (nr__3) ! [0/4]?=- MSX C band quality flag real*4 D_2 (nr__3) ! (Jy) ?=- (?=-9.99) MSX D (14.65um) band flux integer*4 q_D_2 (nr__3) ! [0/4]?=- MSX D band quality flag real*4 E_2 (nr__3) ! (Jy) ?=- (?=-9.99) MSX E (21.34um) band flux integer*4 q_E_2 (nr__3) ! [0/4]?=- MSX E band quality flag c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'b4.dat' ! Catalogue of point sources in the Bulge 4 field integer*4 nr__4 parameter (nr__4=58258) ! Number of records character*466 ar__4 ! Full-size record C J2000 position composed of: RAdeg DEdeg character*18 USB_3 (nr__4) ! Source identification (from position) real*8 RAdeg_4 (nr__4) ! (deg) Right ascension in degrees (J2000) real*8 DEdeg_4 (nr__4) ! (deg) Declination in degrees (J2000) real*4 IRAC1_3 (nr__4) ! (uJy) ?=-9.99 IRAC 1 (3.6um) flux real*4 IRAC2_3 (nr__4) ! (uJy) ?=-9.99 IRAC 2 (4.5um) flux real*4 IRAC3_3 (nr__4) ! (uJy) ?=-9.99 IRAC 3 (5.8um) flux real*4 IRAC4_3 (nr__4) ! (uJy) ?=-9.99 IRAC 4 (8.0um) flux real*4 MIPS_3 (nr__4) ! (uJy) ?=-9.99 MIPS 24 (24um) flux character*1 q_IRAC1_3 (nr__4) ! [ABCD-] Quality flag of IRAC 1 (Sect.2.2.2) character*1 q_IRAC2_3 (nr__4) ! [ABCD-] Quality flag of IRAC 2 (Sect.2.2.2) character*1 q_IRAC3_3 (nr__4) ! [ABCD-] Quality flag of IRAC 3 (Sect.2.2.2) character*1 q_IRAC4_3 (nr__4) ! [ABCD-] Quality flag of IRAC 4 (Sect.2.2.2) character*1 q_MIPS_3 (nr__4) ! [ABCD-] Quality flag of MIPS 24 (Sect.2.2.2) real*4 e_IRAC1_3 (nr__4) ! (uJy) ?=-9.99 Uncertainty IRAC 1 flux real*4 e_IRAC2_3 (nr__4) ! (uJy) ?=-9.99 Uncertainty IRAC 2 flux real*4 e_IRAC3_3 (nr__4) ! (uJy) ?=-9.99 Uncertainty IRAC 3 flux real*4 e_IRAC4_3 (nr__4) ! (uJy) ?=-9.99 Uncertainty IRAC 4 flux real*4 e_MIPS_3 (nr__4) ! (uJy) ?=-9.99 Uncertainty MIPS 24 flux integer*4 f_IRAC1_3 (nr__4) ! [0/1] Outside/Inside IRAC 1 field of view integer*4 f_IRAC2_3 (nr__4) ! [0/1] Outside/Inside IRAC 2 field of view integer*4 f_IRAC3_3 (nr__4) ! [0/1] Outside/Inside IRAC 3 field of view integer*4 f_IRAC4_3 (nr__4) ! [0/1] Outside/Inside IRAC 4 field of view integer*4 f_MIPS_3 (nr__4) ! [0/1] Outside/Inside MIPS 24 field of view real*4 dDEN_3 (nr__4) ! (arcsec) ? Distance to closest DENIS source real*4 Imag_3 (nr__4) ! (mag) ? I-band magnitude of DENIS counterpart real*4 Jmag_3 (nr__4) ! (mag) ? J-band magnitude of DENIS counterpart real*4 Kmag_3 (nr__4) ! (mag) ? K-band magnitude of DENIS counterpart real*4 d2M_3 (nr__4) ! (arcsec) ? Distance to closest 2MASS source real*4 J2mag_3 (nr__4) ! (mag) ?=- (?=-9.999) J-band magnitude of * 2MASS counterpart real*4 H2mag_3 (nr__4) ! (mag) ?=- (?=-9.999) H-band magnitude of * 2MASS counterpart real*4 K2mag_3 (nr__4) ! (mag) ?=- (?=-9.999) K-band magnitude of * 2MASS counterpart character*1 q_J2mag_3 (nr__4) ! [ABCDEFU-] 2MASS J-band quality flag character*1 q_H2mag_3 (nr__4) ! [ABCDEFU-] 2MASS H-band quality flag character*1 q_K2mag_3 (nr__4) ! [ABCDEFU-] 2MASS K-band quality flag real*4 dIRAS_3 (nr__4) ! (arcsec) ? Distance to closest IRAS source real*4 F12_3 (nr__4) ! (Jy) ?=- (?=-9.9) IRAS 12 micron flux integer*4 q_F12_3 (nr__4) ! [1/3]?=- IRAS 12 micron quality flag real*4 F25_3 (nr__4) ! (Jy) ?=- (?=-9.9) IRAS 25 micron flux integer*4 q_F25_3 (nr__4) ! [1/3]?=- IRAS 25 micron quality flag real*4 F60_3 (nr__4) ! (Jy) ?=- (?=-9.9) IRAS 60 micron flux integer*4 q_F60_3 (nr__4) ! [1/3]?=- IRAS 60 micron quality flag real*4 F100_3 (nr__4) ! (Jy) ?=- (?=-9.9) IRAS 100 micron flux integer*4 q_F100_3 (nr__4) ! [1/3]?=- IRAS 100 micron quality flag real*4 dMSX_3 (nr__4) ! (arcsec) ? Distance to closest MSX source real*4 B1_3 (nr__4) ! (Jy) ?=- (?=-9.99) MSX B1 (4.29um) band flux integer*4 q_B1_3 (nr__4) ! [0/4]?=- MSX B1 band quality flag real*4 B2_3 (nr__4) ! (Jy) ?=- (?=-9.99) MSX B2 (4.25um) band flux integer*4 q_B2_3 (nr__4) ! [0/4]?=- MSX B2 band quality flag real*4 A_3 (nr__4) ! (Jy) ?=- (?=-9.99) MSX A (8.28um) band flux integer*4 q_A_3 (nr__4) ! [0/4]?=- MSX A band quality flag real*4 C_3 (nr__4) ! (Jy) ?=- (?=-9.99) MSX C (12.13um) band flux integer*4 q_C_3 (nr__4) ! [0/4]?=- MSX C band quality flag real*4 D_3 (nr__4) ! (Jy) ?=- (?=-9.99) MSX D (14.65um) band flux integer*4 q_D_3 (nr__4) ! [0/4]?=- MSX D band quality flag real*4 E_3 (nr__4) ! (Jy) ?=- (?=-9.99) MSX E (21.34um) band flux integer*4 q_E_3 (nr__4) ! [0/4]?=- MSX E band quality flag c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'bn1.dat' ! Catalogue of point sources in the Bulge N1 field integer*4 nr__5 parameter (nr__5=64798) ! Number of records character*466 ar__5 ! Full-size record C J2000 position composed of: RAdeg DEdeg character*18 USB_4 (nr__5) ! Source identification (from position) real*8 RAdeg_5 (nr__5) ! (deg) Right ascension in degrees (J2000) real*8 DEdeg_5 (nr__5) ! (deg) Declination in degrees (J2000) real*4 IRAC1_4 (nr__5) ! (uJy) ?=-9.99 IRAC 1 (3.6um) flux real*4 IRAC2_4 (nr__5) ! (uJy) ?=-9.99 IRAC 2 (4.5um) flux real*4 IRAC3_4 (nr__5) ! (uJy) ?=-9.99 IRAC 3 (5.8um) flux real*4 IRAC4_4 (nr__5) ! (uJy) ?=-9.99 IRAC 4 (8.0um) flux real*4 MIPS_4 (nr__5) ! (uJy) ?=-9.99 MIPS 24 (24um) flux character*1 q_IRAC1_4 (nr__5) ! [ABCD-] Quality flag of IRAC 1 (Sect.2.2.2) character*1 q_IRAC2_4 (nr__5) ! [ABCD-] Quality flag of IRAC 2 (Sect.2.2.2) character*1 q_IRAC3_4 (nr__5) ! [ABCD-] Quality flag of IRAC 3 (Sect.2.2.2) character*1 q_IRAC4_4 (nr__5) ! [ABCD-] Quality flag of IRAC 4 (Sect.2.2.2) character*1 q_MIPS_4 (nr__5) ! [ABCD-] Quality flag of MIPS 24 (Sect.2.2.2) real*4 e_IRAC1_4 (nr__5) ! (uJy) ?=-9.99 Uncertainty IRAC 1 flux real*4 e_IRAC2_4 (nr__5) ! (uJy) ?=-9.99 Uncertainty IRAC 2 flux real*4 e_IRAC3_4 (nr__5) ! (uJy) ?=-9.99 Uncertainty IRAC 3 flux real*4 e_IRAC4_4 (nr__5) ! (uJy) ?=-9.99 Uncertainty IRAC 4 flux real*4 e_MIPS_4 (nr__5) ! (uJy) ?=-9.99 Uncertainty MIPS 24 flux integer*4 f_IRAC1_4 (nr__5) ! [0/1] Outside/Inside IRAC 1 field of view integer*4 f_IRAC2_4 (nr__5) ! [0/1] Outside/Inside IRAC 2 field of view integer*4 f_IRAC3_4 (nr__5) ! [0/1] Outside/Inside IRAC 3 field of view integer*4 f_IRAC4_4 (nr__5) ! [0/1] Outside/Inside IRAC 4 field of view integer*4 f_MIPS_4 (nr__5) ! [0/1] Outside/Inside MIPS 24 field of view real*4 dDEN_4 (nr__5) ! (arcsec) ? Distance to closest DENIS source real*4 Imag_4 (nr__5) ! (mag) ? I-band magnitude of DENIS counterpart real*4 Jmag_4 (nr__5) ! (mag) ? J-band magnitude of DENIS counterpart real*4 Kmag_4 (nr__5) ! (mag) ? K-band magnitude of DENIS counterpart real*4 d2M_4 (nr__5) ! (arcsec) ? Distance to closest 2MASS source real*4 J2mag_4 (nr__5) ! (mag) ?=- (?=-9.999) J-band magnitude of * 2MASS counterpart real*4 H2mag_4 (nr__5) ! (mag) ?=- (?=-9.999) H-band magnitude of * 2MASS counterpart real*4 K2mag_4 (nr__5) ! (mag) ?=- (?=-9.999) K-band magnitude of * 2MASS counterpart character*1 q_J2mag_4 (nr__5) ! [ABCDEFU-] 2MASS J-band quality flag character*1 q_H2mag_4 (nr__5) ! [ABCDEFU-] 2MASS H-band quality flag character*1 q_K2mag_4 (nr__5) ! [ABCDEFU-] 2MASS K-band quality flag real*4 dIRAS_4 (nr__5) ! (arcsec) ? Distance to closest IRAS source real*4 F12_4 (nr__5) ! (Jy) ?=- (?=-9.9) IRAS 12 micron flux integer*4 q_F12_4 (nr__5) ! [1/3]?=- IRAS 12 micron quality flag real*4 F25_4 (nr__5) ! (Jy) ?=- (?=-9.9) IRAS 25 micron flux integer*4 q_F25_4 (nr__5) ! [1/3]?=- IRAS 25 micron quality flag real*4 F60_4 (nr__5) ! (Jy) ?=- (?=-9.9) IRAS 60 micron flux integer*4 q_F60_4 (nr__5) ! [1/3]?=- IRAS 60 micron quality flag real*4 F100_4 (nr__5) ! (Jy) ?=- (?=-9.9) IRAS 100 micron flux integer*4 q_F100_4 (nr__5) ! [1/3]?=- IRAS 100 micron quality flag real*4 dMSX_4 (nr__5) ! (arcsec) ? Distance to closest MSX source real*4 B1_4 (nr__5) ! (Jy) ?=- (?=-9.99) MSX B1 (4.29um) band flux integer*4 q_B1_4 (nr__5) ! [0/4]?=- MSX B1 band quality flag real*4 B2_4 (nr__5) ! (Jy) ?=- (?=-9.99) MSX B2 (4.25um) band flux integer*4 q_B2_4 (nr__5) ! [0/4]?=- MSX B2 band quality flag real*4 A_4 (nr__5) ! (Jy) ?=- (?=-9.99) MSX A (8.28um) band flux integer*4 q_A_4 (nr__5) ! [0/4]?=- MSX A band quality flag real*4 C_4 (nr__5) ! (Jy) ?=- (?=-9.99) MSX C (12.13um) band flux integer*4 q_C_4 (nr__5) ! [0/4]?=- MSX C band quality flag real*4 D_4 (nr__5) ! (Jy) ?=- (?=-9.99) MSX D (14.65um) band flux integer*4 q_D_4 (nr__5) ! [0/4]?=- MSX D band quality flag real*4 E_4 (nr__5) ! (Jy) ?=- (?=-9.99) MSX E (21.34um) band flux integer*4 q_E_4 (nr__5) ! [0/4]?=- MSX E band quality flag c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'bn2.dat' ! Catalogue of point sources in the Bulge N2 field integer*4 nr__6 parameter (nr__6=428) ! Number of records character*466 ar__6 ! Full-size record C J2000 position composed of: RAdeg DEdeg character*18 USB_5 (nr__6) ! Source identification (from position) real*8 RAdeg_6 (nr__6) ! (deg) Right ascension in degrees (J2000) real*8 DEdeg_6 (nr__6) ! (deg) Declination in degrees (J2000) real*4 IRAC1_5 (nr__6) ! (uJy) ?=-9.99 IRAC 1 (3.6um) flux real*4 IRAC2_5 (nr__6) ! (uJy) ?=-9.99 IRAC 2 (4.5um) flux real*4 IRAC3_5 (nr__6) ! (uJy) ?=-9.99 IRAC 3 (5.8um) flux real*4 IRAC4_5 (nr__6) ! (uJy) ?=-9.99 IRAC 4 (8.0um) flux real*4 MIPS_5 (nr__6) ! (uJy) ?=-9.99 MIPS 24 (24um) flux character*1 q_IRAC1_5 (nr__6) ! [ABCD-] Quality flag of IRAC 1 (Sect.2.2.2) character*1 q_IRAC2_5 (nr__6) ! [ABCD-] Quality flag of IRAC 2 (Sect.2.2.2) character*1 q_IRAC3_5 (nr__6) ! [ABCD-] Quality flag of IRAC 3 (Sect.2.2.2) character*1 q_IRAC4_5 (nr__6) ! [ABCD-] Quality flag of IRAC 4 (Sect.2.2.2) character*1 q_MIPS_5 (nr__6) ! [ABCD-] Quality flag of MIPS 24 (Sect.2.2.2) real*4 e_IRAC1_5 (nr__6) ! (uJy) ?=-9.99 Uncertainty IRAC 1 flux real*4 e_IRAC2_5 (nr__6) ! (uJy) ?=-9.99 Uncertainty IRAC 2 flux real*4 e_IRAC3_5 (nr__6) ! (uJy) ?=-9.99 Uncertainty IRAC 3 flux real*4 e_IRAC4_5 (nr__6) ! (uJy) ?=-9.99 Uncertainty IRAC 4 flux real*4 e_MIPS_5 (nr__6) ! (uJy) ?=-9.99 Uncertainty MIPS 24 flux integer*4 f_IRAC1_5 (nr__6) ! [0/1] Outside/Inside IRAC 1 field of view integer*4 f_IRAC2_5 (nr__6) ! [0/1] Outside/Inside IRAC 2 field of view integer*4 f_IRAC3_5 (nr__6) ! [0/1] Outside/Inside IRAC 3 field of view integer*4 f_IRAC4_5 (nr__6) ! [0/1] Outside/Inside IRAC 4 field of view integer*4 f_MIPS_5 (nr__6) ! [0/1] Outside/Inside MIPS 24 field of view real*4 dDEN_5 (nr__6) ! (arcsec) ? Distance to closest DENIS source real*4 Imag_5 (nr__6) ! (mag) ? I-band magnitude of DENIS counterpart real*4 Jmag_5 (nr__6) ! (mag) ? J-band magnitude of DENIS counterpart real*4 Kmag_5 (nr__6) ! (mag) ? K-band magnitude of DENIS counterpart real*4 d2M_5 (nr__6) ! (arcsec) ? Distance to closest 2MASS source real*4 J2mag_5 (nr__6) ! (mag) ?=- (?=-9.999) J-band magnitude of * 2MASS counterpart real*4 H2mag_5 (nr__6) ! (mag) ?=- (?=-9.999) H-band magnitude of * 2MASS counterpart real*4 K2mag_5 (nr__6) ! (mag) ?=- (?=-9.999) K-band magnitude of * 2MASS counterpart character*1 q_J2mag_5 (nr__6) ! [ABCDEFU-] 2MASS J-band quality flag character*1 q_H2mag_5 (nr__6) ! [ABCDEFU-] 2MASS H-band quality flag character*1 q_K2mag_5 (nr__6) ! [ABCDEFU-] 2MASS K-band quality flag real*4 dIRAS_5 (nr__6) ! (arcsec) ? Distance to closest IRAS source real*4 F12_5 (nr__6) ! (Jy) ?=- (?=-9.9) IRAS 12 micron flux integer*4 q_F12_5 (nr__6) ! [1/3]?=- IRAS 12 micron quality flag real*4 F25_5 (nr__6) ! (Jy) ?=- (?=-9.9) IRAS 25 micron flux integer*4 q_F25_5 (nr__6) ! [1/3]?=- IRAS 25 micron quality flag real*4 F60_5 (nr__6) ! (Jy) ?=- (?=-9.9) IRAS 60 micron flux integer*4 q_F60_5 (nr__6) ! [1/3]?=- IRAS 60 micron quality flag real*4 F100_5 (nr__6) ! (Jy) ?=- (?=-9.9) IRAS 100 micron flux integer*4 q_F100_5 (nr__6) ! [1/3]?=- IRAS 100 micron quality flag real*4 dMSX_5 (nr__6) ! (arcsec) ? Distance to closest MSX source real*4 B1_5 (nr__6) ! (Jy) ?=- (?=-9.99) MSX B1 (4.29um) band flux integer*4 q_B1_5 (nr__6) ! [0/4]?=- MSX B1 band quality flag real*4 B2_5 (nr__6) ! (Jy) ?=- (?=-9.99) MSX B2 (4.25um) band flux integer*4 q_B2_5 (nr__6) ! [0/4]?=- MSX B2 band quality flag real*4 A_5 (nr__6) ! (Jy) ?=- (?=-9.99) MSX A (8.28um) band flux integer*4 q_A_5 (nr__6) ! [0/4]?=- MSX A band quality flag real*4 C_5 (nr__6) ! (Jy) ?=- (?=-9.99) MSX C (12.13um) band flux integer*4 q_C_5 (nr__6) ! [0/4]?=- MSX C band quality flag real*4 D_5 (nr__6) ! (Jy) ?=- (?=-9.99) MSX D (14.65um) band flux integer*4 q_D_5 (nr__6) ! [0/4]?=- MSX D band quality flag real*4 E_5 (nr__6) ! (Jy) ?=- (?=-9.99) MSX E (21.34um) band flux integer*4 q_E_5 (nr__6) ! [0/4]?=- MSX E band quality flag c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'bngc6522.dat' ! Catalogue of point sources in the NGC 6522 field integer*4 nr__7 parameter (nr__7=47232) ! Number of records character*466 ar__7 ! Full-size record C J2000 position composed of: RAdeg DEdeg character*18 USB_6 (nr__7) ! Source identification (from position) real*8 RAdeg_7 (nr__7) ! (deg) Right ascension in degrees (J2000) real*8 DEdeg_7 (nr__7) ! (deg) Declination in degrees (J2000) real*4 IRAC1_6 (nr__7) ! (uJy) ?=-9.99 IRAC 1 (3.6um) flux real*4 IRAC2_6 (nr__7) ! (uJy) ?=-9.99 IRAC 2 (4.5um) flux real*4 IRAC3_6 (nr__7) ! (uJy) ?=-9.99 IRAC 3 (5.8um) flux real*4 IRAC4_6 (nr__7) ! (uJy) ?=-9.99 IRAC 4 (8.0um) flux real*4 MIPS_6 (nr__7) ! (uJy) ?=-9.99 MIPS 24 (24um) flux character*1 q_IRAC1_6 (nr__7) ! [ABCD-] Quality flag of IRAC 1 (Sect.2.2.2) character*1 q_IRAC2_6 (nr__7) ! [ABCD-] Quality flag of IRAC 2 (Sect.2.2.2) character*1 q_IRAC3_6 (nr__7) ! [ABCD-] Quality flag of IRAC 3 (Sect.2.2.2) character*1 q_IRAC4_6 (nr__7) ! [ABCD-] Quality flag of IRAC 4 (Sect.2.2.2) character*1 q_MIPS_6 (nr__7) ! [ABCD-] Quality flag of MIPS 24 (Sect.2.2.2) real*4 e_IRAC1_6 (nr__7) ! (uJy) ?=-9.99 Uncertainty IRAC 1 flux real*4 e_IRAC2_6 (nr__7) ! (uJy) ?=-9.99 Uncertainty IRAC 2 flux real*4 e_IRAC3_6 (nr__7) ! (uJy) ?=-9.99 Uncertainty IRAC 3 flux real*4 e_IRAC4_6 (nr__7) ! (uJy) ?=-9.99 Uncertainty IRAC 4 flux real*4 e_MIPS_6 (nr__7) ! (uJy) ?=-9.99 Uncertainty MIPS 24 flux integer*4 f_IRAC1_6 (nr__7) ! [0/1] Outside/Inside IRAC 1 field of view integer*4 f_IRAC2_6 (nr__7) ! [0/1] Outside/Inside IRAC 2 field of view integer*4 f_IRAC3_6 (nr__7) ! [0/1] Outside/Inside IRAC 3 field of view integer*4 f_IRAC4_6 (nr__7) ! [0/1] Outside/Inside IRAC 4 field of view integer*4 f_MIPS_6 (nr__7) ! [0/1] Outside/Inside MIPS 24 field of view real*4 dDEN_6 (nr__7) ! (arcsec) ? Distance to closest DENIS source real*4 Imag_6 (nr__7) ! (mag) ? I-band magnitude of DENIS counterpart real*4 Jmag_6 (nr__7) ! (mag) ? J-band magnitude of DENIS counterpart real*4 Kmag_6 (nr__7) ! (mag) ? K-band magnitude of DENIS counterpart real*4 d2M_6 (nr__7) ! (arcsec) ? Distance to closest 2MASS source real*4 J2mag_6 (nr__7) ! (mag) ?=- (?=-9.999) J-band magnitude of * 2MASS counterpart real*4 H2mag_6 (nr__7) ! (mag) ?=- (?=-9.999) H-band magnitude of * 2MASS counterpart real*4 K2mag_6 (nr__7) ! (mag) ?=- (?=-9.999) K-band magnitude of * 2MASS counterpart character*1 q_J2mag_6 (nr__7) ! [ABCDEFU-] 2MASS J-band quality flag character*1 q_H2mag_6 (nr__7) ! [ABCDEFU-] 2MASS H-band quality flag character*1 q_K2mag_6 (nr__7) ! [ABCDEFU-] 2MASS K-band quality flag real*4 dIRAS_6 (nr__7) ! (arcsec) ? Distance to closest IRAS source real*4 F12_6 (nr__7) ! (Jy) ?=- (?=-9.9) IRAS 12 micron flux integer*4 q_F12_6 (nr__7) ! [1/3]?=- IRAS 12 micron quality flag real*4 F25_6 (nr__7) ! (Jy) ?=- (?=-9.9) IRAS 25 micron flux integer*4 q_F25_6 (nr__7) ! [1/3]?=- IRAS 25 micron quality flag real*4 F60_6 (nr__7) ! (Jy) ?=- (?=-9.9) IRAS 60 micron flux integer*4 q_F60_6 (nr__7) ! [1/3]?=- IRAS 60 micron quality flag real*4 F100_6 (nr__7) ! (Jy) ?=- (?=-9.9) IRAS 100 micron flux integer*4 q_F100_6 (nr__7) ! [1/3]?=- IRAS 100 micron quality flag real*4 dMSX_6 (nr__7) ! (arcsec) ? Distance to closest MSX source real*4 B1_6 (nr__7) ! (Jy) ?=- (?=-9.99) MSX B1 (4.29um) band flux integer*4 q_B1_6 (nr__7) ! [0/4]?=- MSX B1 band quality flag real*4 B2_6 (nr__7) ! (Jy) ?=- (?=-9.99) MSX B2 (4.25um) band flux integer*4 q_B2_6 (nr__7) ! [0/4]?=- MSX B2 band quality flag real*4 A_6 (nr__7) ! (Jy) ?=- (?=-9.99) MSX A (8.28um) band flux integer*4 q_A_6 (nr__7) ! [0/4]?=- MSX A band quality flag real*4 C_6 (nr__7) ! (Jy) ?=- (?=-9.99) MSX C (12.13um) band flux integer*4 q_C_6 (nr__7) ! [0/4]?=- MSX C band quality flag real*4 D_6 (nr__7) ! (Jy) ?=- (?=-9.99) MSX D (14.65um) band flux integer*4 q_D_6 (nr__7) ! [0/4]?=- MSX D band quality flag real*4 E_6 (nr__7) ! (Jy) ?=- (?=-9.99) MSX E (21.34um) band flux integer*4 q_E_6 (nr__7) ! [0/4]?=- MSX E band quality flag C============================================================================= C Loading file 'fits.dat' ! Summary of FITS images C Format for file interpretation 1 format(A11,3X,A1,2X,I4,1X,I4,2X,F9.5,F9.5) C Effective file loading open(unit=1,status='old',file= +'fits.dat') write(6,*) '....Loading file: fits.dat' do i__=1,31 read(1,'(A46)')ar__ read(ar__,1) + FileName(i__),Ch(i__),Xpix(i__),Ypix(i__),RAdeg(i__), + DEdeg(i__) c ..............Just test output........... write(6,1) + FileName(i__),Ch(i__),Xpix(i__),Ypix(i__),RAdeg(i__), + DEdeg(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'b1.dat' ! Catalogue of point sources in the Bulge 1 field C Format for file interpretation 2 format( + 3X,A18,3X,F10.6,3X,F10.6,3X,E12.6,3X,E12.6,3X,E12.6,3X,E12.6, + 3X,E12.6,3X,A1,A1,A1,A1,A1,3X,E12.6,3X,E12.6,3X,E12.6,3X, + E12.6,3X,E12.6,3X,I1,I1,I1,I1,I1,3X,F5.3,3X,F6.3,3X,F6.3,3X, + F6.3,3X,F5.3,3X,F6.3,3X,F6.3,3X,F6.3,3X,A1,A1,A1,3X,F5.1,3X, + E8.2,3X,I1,3X,E8.2,3X,I1,3X,E8.2,3X,I1,3X,E8.2,3X,I1,3X,F4.1, + 3X,E10.3,3X,I1,3X,E10.3,3X,I1,3X,E10.3,3X,I1,3X,E10.3,3X,I1, + 3X,E10.3,3X,I1,3X,E10.3,3X,I1) C Effective file loading open(unit=1,status='old',file= +'b1.dat') write(6,*) '....Loading file: b1.dat' do i__=1,69047 read(1,'(A466)')ar__1 read(ar__1,2) + USB(i__),RAdeg_1(i__),DEdeg_1(i__),IRAC1(i__),IRAC2(i__), + IRAC3(i__),IRAC4(i__),MIPS(i__),q_IRAC1(i__),q_IRAC2(i__), + q_IRAC3(i__),q_IRAC4(i__),q_MIPS(i__),e_IRAC1(i__), + e_IRAC2(i__),e_IRAC3(i__),e_IRAC4(i__),e_MIPS(i__), + f_IRAC1(i__),f_IRAC2(i__),f_IRAC3(i__),f_IRAC4(i__), + f_MIPS(i__),dDEN(i__),Imag(i__),Jmag(i__),Kmag(i__),d2M(i__), + J2mag(i__),H2mag(i__),K2mag(i__),q_J2mag(i__),q_H2mag(i__), + q_K2mag(i__),dIRAS(i__),F12(i__),q_F12(i__),F25(i__), + q_F25(i__),F60(i__),q_F60(i__),F100(i__),q_F100(i__), + dMSX(i__),B1(i__),q_B1(i__),B2(i__),q_B2(i__),A(i__),q_A(i__), + C(i__),q_C(i__),D(i__),q_D(i__),E(i__),q_E(i__) if(ar__1(217:221) .EQ. '') dDEN(i__) = rNULL__ if(ar__1(225:230) .EQ. '') Imag(i__) = rNULL__ if(ar__1(234:239) .EQ. '') Jmag(i__) = rNULL__ if(ar__1(243:248) .EQ. '') Kmag(i__) = rNULL__ if(ar__1(252:256) .EQ. '') d2M(i__) = rNULL__ if (idig(ar__1(260:265)).EQ.0) J2mag(i__) = rNULL__ if (idig(ar__1(269:274)).EQ.0) H2mag(i__) = rNULL__ if (idig(ar__1(278:283)).EQ.0) K2mag(i__) = rNULL__ if(ar__1(293:297) .EQ. '') dIRAS(i__) = rNULL__ if (idig(ar__1(301:308)).EQ.0) F12(i__) = rNULL__ if (q_F12(i__) .EQ. 45) q_F12(i__) = iNULL__ if (idig(ar__1(316:323)).EQ.0) F25(i__) = rNULL__ if (q_F25(i__) .EQ. 45) q_F25(i__) = iNULL__ if (idig(ar__1(331:338)).EQ.0) F60(i__) = rNULL__ if (q_F60(i__) .EQ. 45) q_F60(i__) = iNULL__ if (idig(ar__1(346:353)).EQ.0) F100(i__) = rNULL__ if (q_F100(i__) .EQ. 45) q_F100(i__) = iNULL__ if(ar__1(361:364) .EQ. '') dMSX(i__) = rNULL__ if (idig(ar__1(368:377)).EQ.0) B1(i__) = rNULL__ if (q_B1(i__) .EQ. 45) q_B1(i__) = iNULL__ if (idig(ar__1(385:394)).EQ.0) B2(i__) = rNULL__ if (q_B2(i__) .EQ. 45) q_B2(i__) = iNULL__ if (idig(ar__1(402:411)).EQ.0) A(i__) = rNULL__ if (q_A(i__) .EQ. 45) q_A(i__) = iNULL__ if (idig(ar__1(419:428)).EQ.0) C(i__) = rNULL__ if (q_C(i__) .EQ. 45) q_C(i__) = iNULL__ if (idig(ar__1(436:445)).EQ.0) D(i__) = rNULL__ if (q_D(i__) .EQ. 45) q_D(i__) = iNULL__ if (idig(ar__1(453:462)).EQ.0) E(i__) = rNULL__ if (q_E(i__) .EQ. 45) q_E(i__) = iNULL__ c ..............Just test output........... write(6,2) + USB(i__),RAdeg_1(i__),DEdeg_1(i__),IRAC1(i__),IRAC2(i__), + IRAC3(i__),IRAC4(i__),MIPS(i__),q_IRAC1(i__),q_IRAC2(i__), + q_IRAC3(i__),q_IRAC4(i__),q_MIPS(i__),e_IRAC1(i__), + e_IRAC2(i__),e_IRAC3(i__),e_IRAC4(i__),e_MIPS(i__), + f_IRAC1(i__),f_IRAC2(i__),f_IRAC3(i__),f_IRAC4(i__), + f_MIPS(i__),dDEN(i__),Imag(i__),Jmag(i__),Kmag(i__),d2M(i__), + J2mag(i__),H2mag(i__),K2mag(i__),q_J2mag(i__),q_H2mag(i__), + q_K2mag(i__),dIRAS(i__),F12(i__),q_F12(i__),F25(i__), + q_F25(i__),F60(i__),q_F60(i__),F100(i__),q_F100(i__), + dMSX(i__),B1(i__),q_B1(i__),B2(i__),q_B2(i__),A(i__),q_A(i__), + C(i__),q_C(i__),D(i__),q_D(i__),E(i__),q_E(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'b2.dat' ! Catalogue of point sources in the Bulge 2 field C Format for file interpretation 3 format( + 3X,A18,3X,F10.6,3X,F10.6,3X,E12.6,3X,E12.6,3X,E12.6,3X,E12.6, + 3X,E12.6,3X,A1,A1,A1,A1,A1,3X,E12.6,3X,E12.6,3X,E12.6,3X, + E12.6,3X,E12.6,3X,I1,I1,I1,I1,I1,3X,F5.3,3X,F6.3,3X,F6.3,3X, + F6.3,3X,F5.3,3X,F6.3,3X,F6.3,3X,F6.3,3X,A1,A1,A1,3X,F5.1,3X, + E8.2,3X,I1,3X,E8.2,3X,I1,3X,E8.2,3X,I1,3X,E8.2,3X,I1,3X,F4.1, + 3X,E10.3,3X,I1,3X,E10.3,3X,I1,3X,E10.3,3X,I1,3X,E10.3,3X,I1, + 3X,E10.3,3X,I1,3X,E10.3,3X,I1) C Effective file loading open(unit=1,status='old',file= +'b2.dat') write(6,*) '....Loading file: b2.dat' do i__=1,60262 read(1,'(A466)')ar__2 read(ar__2,3) + USB_1(i__),RAdeg_2(i__),DEdeg_2(i__),IRAC1_1(i__), + IRAC2_1(i__),IRAC3_1(i__),IRAC4_1(i__),MIPS_1(i__), + q_IRAC1_1(i__),q_IRAC2_1(i__),q_IRAC3_1(i__),q_IRAC4_1(i__), + q_MIPS_1(i__),e_IRAC1_1(i__),e_IRAC2_1(i__),e_IRAC3_1(i__), + e_IRAC4_1(i__),e_MIPS_1(i__),f_IRAC1_1(i__),f_IRAC2_1(i__), + f_IRAC3_1(i__),f_IRAC4_1(i__),f_MIPS_1(i__),dDEN_1(i__), + Imag_1(i__),Jmag_1(i__),Kmag_1(i__),d2M_1(i__),J2mag_1(i__), + H2mag_1(i__),K2mag_1(i__),q_J2mag_1(i__),q_H2mag_1(i__), + q_K2mag_1(i__),dIRAS_1(i__),F12_1(i__),q_F12_1(i__), + F25_1(i__),q_F25_1(i__),F60_1(i__),q_F60_1(i__),F100_1(i__), + q_F100_1(i__),dMSX_1(i__),B1_1(i__),q_B1_1(i__),B2_1(i__), + q_B2_1(i__),A_1(i__),q_A_1(i__),C_1(i__),q_C_1(i__),D_1(i__), + q_D_1(i__),E_1(i__),q_E_1(i__) if(ar__2(217:221) .EQ. '') dDEN_1(i__) = rNULL__ if(ar__2(225:230) .EQ. '') Imag_1(i__) = rNULL__ if(ar__2(234:239) .EQ. '') Jmag_1(i__) = rNULL__ if(ar__2(243:248) .EQ. '') Kmag_1(i__) = rNULL__ if(ar__2(252:256) .EQ. '') d2M_1(i__) = rNULL__ if (idig(ar__2(260:265)).EQ.0) J2mag_1(i__) = rNULL__ if (idig(ar__2(269:274)).EQ.0) H2mag_1(i__) = rNULL__ if (idig(ar__2(278:283)).EQ.0) K2mag_1(i__) = rNULL__ if(ar__2(293:297) .EQ. '') dIRAS_1(i__) = rNULL__ if (idig(ar__2(301:308)).EQ.0) F12_1(i__) = rNULL__ if (q_F12_1(i__) .EQ. 45) q_F12_1(i__) = iNULL__ if (idig(ar__2(316:323)).EQ.0) F25_1(i__) = rNULL__ if (q_F25_1(i__) .EQ. 45) q_F25_1(i__) = iNULL__ if (idig(ar__2(331:338)).EQ.0) F60_1(i__) = rNULL__ if (q_F60_1(i__) .EQ. 45) q_F60_1(i__) = iNULL__ if (idig(ar__2(346:353)).EQ.0) F100_1(i__) = rNULL__ if (q_F100_1(i__) .EQ. 45) q_F100_1(i__) = iNULL__ if(ar__2(361:364) .EQ. '') dMSX_1(i__) = rNULL__ if (idig(ar__2(368:377)).EQ.0) B1_1(i__) = rNULL__ if (q_B1_1(i__) .EQ. 45) q_B1_1(i__) = iNULL__ if (idig(ar__2(385:394)).EQ.0) B2_1(i__) = rNULL__ if (q_B2_1(i__) .EQ. 45) q_B2_1(i__) = iNULL__ if (idig(ar__2(402:411)).EQ.0) A_1(i__) = rNULL__ if (q_A_1(i__) .EQ. 45) q_A_1(i__) = iNULL__ if (idig(ar__2(419:428)).EQ.0) C_1(i__) = rNULL__ if (q_C_1(i__) .EQ. 45) q_C_1(i__) = iNULL__ if (idig(ar__2(436:445)).EQ.0) D_1(i__) = rNULL__ if (q_D_1(i__) .EQ. 45) q_D_1(i__) = iNULL__ if (idig(ar__2(453:462)).EQ.0) E_1(i__) = rNULL__ if (q_E_1(i__) .EQ. 45) q_E_1(i__) = iNULL__ c ..............Just test output........... write(6,3) + USB_1(i__),RAdeg_2(i__),DEdeg_2(i__),IRAC1_1(i__), + IRAC2_1(i__),IRAC3_1(i__),IRAC4_1(i__),MIPS_1(i__), + q_IRAC1_1(i__),q_IRAC2_1(i__),q_IRAC3_1(i__),q_IRAC4_1(i__), + q_MIPS_1(i__),e_IRAC1_1(i__),e_IRAC2_1(i__),e_IRAC3_1(i__), + e_IRAC4_1(i__),e_MIPS_1(i__),f_IRAC1_1(i__),f_IRAC2_1(i__), + f_IRAC3_1(i__),f_IRAC4_1(i__),f_MIPS_1(i__),dDEN_1(i__), + Imag_1(i__),Jmag_1(i__),Kmag_1(i__),d2M_1(i__),J2mag_1(i__), + H2mag_1(i__),K2mag_1(i__),q_J2mag_1(i__),q_H2mag_1(i__), + q_K2mag_1(i__),dIRAS_1(i__),F12_1(i__),q_F12_1(i__), + F25_1(i__),q_F25_1(i__),F60_1(i__),q_F60_1(i__),F100_1(i__), + q_F100_1(i__),dMSX_1(i__),B1_1(i__),q_B1_1(i__),B2_1(i__), + q_B2_1(i__),A_1(i__),q_A_1(i__),C_1(i__),q_C_1(i__),D_1(i__), + q_D_1(i__),E_1(i__),q_E_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'b3.dat' ! Catalogue of point sources in the Bulge 3 field C Format for file interpretation 4 format( + 3X,A18,3X,F10.6,3X,F10.6,3X,E12.6,3X,E12.6,3X,E12.6,3X,E12.6, + 3X,E12.6,3X,A1,A1,A1,A1,A1,3X,E12.6,3X,E12.6,3X,E12.6,3X, + E12.6,3X,E12.6,3X,I1,I1,I1,I1,I1,3X,F5.3,3X,F6.3,3X,F6.3,3X, + F6.3,3X,F5.3,3X,F6.3,3X,F6.3,3X,F6.3,3X,A1,A1,A1,3X,F5.1,3X, + E8.2,3X,I1,3X,E8.2,3X,I1,3X,E8.2,3X,I1,3X,E8.2,3X,I1,3X,F4.1, + 3X,E10.3,3X,I1,3X,E10.3,3X,I1,3X,E10.3,3X,I1,3X,E10.3,3X,I1, + 3X,E10.3,3X,I1,3X,E10.3,3X,I1) C Effective file loading open(unit=1,status='old',file= +'b3.dat') write(6,*) '....Loading file: b3.dat' do i__=1,52576 read(1,'(A466)')ar__3 read(ar__3,4) + USB_2(i__),RAdeg_3(i__),DEdeg_3(i__),IRAC1_2(i__), + IRAC2_2(i__),IRAC3_2(i__),IRAC4_2(i__),MIPS_2(i__), + q_IRAC1_2(i__),q_IRAC2_2(i__),q_IRAC3_2(i__),q_IRAC4_2(i__), + q_MIPS_2(i__),e_IRAC1_2(i__),e_IRAC2_2(i__),e_IRAC3_2(i__), + e_IRAC4_2(i__),e_MIPS_2(i__),f_IRAC1_2(i__),f_IRAC2_2(i__), + f_IRAC3_2(i__),f_IRAC4_2(i__),f_MIPS_2(i__),dDEN_2(i__), + Imag_2(i__),Jmag_2(i__),Kmag_2(i__),d2M_2(i__),J2mag_2(i__), + H2mag_2(i__),K2mag_2(i__),q_J2mag_2(i__),q_H2mag_2(i__), + q_K2mag_2(i__),dIRAS_2(i__),F12_2(i__),q_F12_2(i__), + F25_2(i__),q_F25_2(i__),F60_2(i__),q_F60_2(i__),F100_2(i__), + q_F100_2(i__),dMSX_2(i__),B1_2(i__),q_B1_2(i__),B2_2(i__), + q_B2_2(i__),A_2(i__),q_A_2(i__),C_2(i__),q_C_2(i__),D_2(i__), + q_D_2(i__),E_2(i__),q_E_2(i__) if(ar__3(217:221) .EQ. '') dDEN_2(i__) = rNULL__ if(ar__3(225:230) .EQ. '') Imag_2(i__) = rNULL__ if(ar__3(234:239) .EQ. '') Jmag_2(i__) = rNULL__ if(ar__3(243:248) .EQ. '') Kmag_2(i__) = rNULL__ if(ar__3(252:256) .EQ. '') d2M_2(i__) = rNULL__ if (idig(ar__3(260:265)).EQ.0) J2mag_2(i__) = rNULL__ if (idig(ar__3(269:274)).EQ.0) H2mag_2(i__) = rNULL__ if (idig(ar__3(278:283)).EQ.0) K2mag_2(i__) = rNULL__ if(ar__3(293:297) .EQ. '') dIRAS_2(i__) = rNULL__ if (idig(ar__3(301:308)).EQ.0) F12_2(i__) = rNULL__ if (q_F12_2(i__) .EQ. 45) q_F12_2(i__) = iNULL__ if (idig(ar__3(316:323)).EQ.0) F25_2(i__) = rNULL__ if (q_F25_2(i__) .EQ. 45) q_F25_2(i__) = iNULL__ if (idig(ar__3(331:338)).EQ.0) F60_2(i__) = rNULL__ if (q_F60_2(i__) .EQ. 45) q_F60_2(i__) = iNULL__ if (idig(ar__3(346:353)).EQ.0) F100_2(i__) = rNULL__ if (q_F100_2(i__) .EQ. 45) q_F100_2(i__) = iNULL__ if(ar__3(361:364) .EQ. '') dMSX_2(i__) = rNULL__ if (idig(ar__3(368:377)).EQ.0) B1_2(i__) = rNULL__ if (q_B1_2(i__) .EQ. 45) q_B1_2(i__) = iNULL__ if (idig(ar__3(385:394)).EQ.0) B2_2(i__) = rNULL__ if (q_B2_2(i__) .EQ. 45) q_B2_2(i__) = iNULL__ if (idig(ar__3(402:411)).EQ.0) A_2(i__) = rNULL__ if (q_A_2(i__) .EQ. 45) q_A_2(i__) = iNULL__ if (idig(ar__3(419:428)).EQ.0) C_2(i__) = rNULL__ if (q_C_2(i__) .EQ. 45) q_C_2(i__) = iNULL__ if (idig(ar__3(436:445)).EQ.0) D_2(i__) = rNULL__ if (q_D_2(i__) .EQ. 45) q_D_2(i__) = iNULL__ if (idig(ar__3(453:462)).EQ.0) E_2(i__) = rNULL__ if (q_E_2(i__) .EQ. 45) q_E_2(i__) = iNULL__ c ..............Just test output........... write(6,4) + USB_2(i__),RAdeg_3(i__),DEdeg_3(i__),IRAC1_2(i__), + IRAC2_2(i__),IRAC3_2(i__),IRAC4_2(i__),MIPS_2(i__), + q_IRAC1_2(i__),q_IRAC2_2(i__),q_IRAC3_2(i__),q_IRAC4_2(i__), + q_MIPS_2(i__),e_IRAC1_2(i__),e_IRAC2_2(i__),e_IRAC3_2(i__), + e_IRAC4_2(i__),e_MIPS_2(i__),f_IRAC1_2(i__),f_IRAC2_2(i__), + f_IRAC3_2(i__),f_IRAC4_2(i__),f_MIPS_2(i__),dDEN_2(i__), + Imag_2(i__),Jmag_2(i__),Kmag_2(i__),d2M_2(i__),J2mag_2(i__), + H2mag_2(i__),K2mag_2(i__),q_J2mag_2(i__),q_H2mag_2(i__), + q_K2mag_2(i__),dIRAS_2(i__),F12_2(i__),q_F12_2(i__), + F25_2(i__),q_F25_2(i__),F60_2(i__),q_F60_2(i__),F100_2(i__), + q_F100_2(i__),dMSX_2(i__),B1_2(i__),q_B1_2(i__),B2_2(i__), + q_B2_2(i__),A_2(i__),q_A_2(i__),C_2(i__),q_C_2(i__),D_2(i__), + q_D_2(i__),E_2(i__),q_E_2(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'b4.dat' ! Catalogue of point sources in the Bulge 4 field C Format for file interpretation 5 format( + 3X,A18,3X,F10.6,3X,F10.6,3X,E12.6,3X,E12.6,3X,E12.6,3X,E12.6, + 3X,E12.6,3X,A1,A1,A1,A1,A1,3X,E12.6,3X,E12.6,3X,E12.6,3X, + E12.6,3X,E12.6,3X,I1,I1,I1,I1,I1,3X,F5.3,3X,F6.3,3X,F6.3,3X, + F6.3,3X,F5.3,3X,F6.3,3X,F6.3,3X,F6.3,3X,A1,A1,A1,3X,F5.1,3X, + E8.2,3X,I1,3X,E8.2,3X,I1,3X,E8.2,3X,I1,3X,E8.2,3X,I1,3X,F4.1, + 3X,E10.3,3X,I1,3X,E10.3,3X,I1,3X,E10.3,3X,I1,3X,E10.3,3X,I1, + 3X,E10.3,3X,I1,3X,E10.3,3X,I1) C Effective file loading open(unit=1,status='old',file= +'b4.dat') write(6,*) '....Loading file: b4.dat' do i__=1,58258 read(1,'(A466)')ar__4 read(ar__4,5) + USB_3(i__),RAdeg_4(i__),DEdeg_4(i__),IRAC1_3(i__), + IRAC2_3(i__),IRAC3_3(i__),IRAC4_3(i__),MIPS_3(i__), + q_IRAC1_3(i__),q_IRAC2_3(i__),q_IRAC3_3(i__),q_IRAC4_3(i__), + q_MIPS_3(i__),e_IRAC1_3(i__),e_IRAC2_3(i__),e_IRAC3_3(i__), + e_IRAC4_3(i__),e_MIPS_3(i__),f_IRAC1_3(i__),f_IRAC2_3(i__), + f_IRAC3_3(i__),f_IRAC4_3(i__),f_MIPS_3(i__),dDEN_3(i__), + Imag_3(i__),Jmag_3(i__),Kmag_3(i__),d2M_3(i__),J2mag_3(i__), + H2mag_3(i__),K2mag_3(i__),q_J2mag_3(i__),q_H2mag_3(i__), + q_K2mag_3(i__),dIRAS_3(i__),F12_3(i__),q_F12_3(i__), + F25_3(i__),q_F25_3(i__),F60_3(i__),q_F60_3(i__),F100_3(i__), + q_F100_3(i__),dMSX_3(i__),B1_3(i__),q_B1_3(i__),B2_3(i__), + q_B2_3(i__),A_3(i__),q_A_3(i__),C_3(i__),q_C_3(i__),D_3(i__), + q_D_3(i__),E_3(i__),q_E_3(i__) if(ar__4(217:221) .EQ. '') dDEN_3(i__) = rNULL__ if(ar__4(225:230) .EQ. '') Imag_3(i__) = rNULL__ if(ar__4(234:239) .EQ. '') Jmag_3(i__) = rNULL__ if(ar__4(243:248) .EQ. '') Kmag_3(i__) = rNULL__ if(ar__4(252:256) .EQ. '') d2M_3(i__) = rNULL__ if (idig(ar__4(260:265)).EQ.0) J2mag_3(i__) = rNULL__ if (idig(ar__4(269:274)).EQ.0) H2mag_3(i__) = rNULL__ if (idig(ar__4(278:283)).EQ.0) K2mag_3(i__) = rNULL__ if(ar__4(293:297) .EQ. '') dIRAS_3(i__) = rNULL__ if (idig(ar__4(301:308)).EQ.0) F12_3(i__) = rNULL__ if (q_F12_3(i__) .EQ. 45) q_F12_3(i__) = iNULL__ if (idig(ar__4(316:323)).EQ.0) F25_3(i__) = rNULL__ if (q_F25_3(i__) .EQ. 45) q_F25_3(i__) = iNULL__ if (idig(ar__4(331:338)).EQ.0) F60_3(i__) = rNULL__ if (q_F60_3(i__) .EQ. 45) q_F60_3(i__) = iNULL__ if (idig(ar__4(346:353)).EQ.0) F100_3(i__) = rNULL__ if (q_F100_3(i__) .EQ. 45) q_F100_3(i__) = iNULL__ if(ar__4(361:364) .EQ. '') dMSX_3(i__) = rNULL__ if (idig(ar__4(368:377)).EQ.0) B1_3(i__) = rNULL__ if (q_B1_3(i__) .EQ. 45) q_B1_3(i__) = iNULL__ if (idig(ar__4(385:394)).EQ.0) B2_3(i__) = rNULL__ if (q_B2_3(i__) .EQ. 45) q_B2_3(i__) = iNULL__ if (idig(ar__4(402:411)).EQ.0) A_3(i__) = rNULL__ if (q_A_3(i__) .EQ. 45) q_A_3(i__) = iNULL__ if (idig(ar__4(419:428)).EQ.0) C_3(i__) = rNULL__ if (q_C_3(i__) .EQ. 45) q_C_3(i__) = iNULL__ if (idig(ar__4(436:445)).EQ.0) D_3(i__) = rNULL__ if (q_D_3(i__) .EQ. 45) q_D_3(i__) = iNULL__ if (idig(ar__4(453:462)).EQ.0) E_3(i__) = rNULL__ if (q_E_3(i__) .EQ. 45) q_E_3(i__) = iNULL__ c ..............Just test output........... write(6,5) + USB_3(i__),RAdeg_4(i__),DEdeg_4(i__),IRAC1_3(i__), + IRAC2_3(i__),IRAC3_3(i__),IRAC4_3(i__),MIPS_3(i__), + q_IRAC1_3(i__),q_IRAC2_3(i__),q_IRAC3_3(i__),q_IRAC4_3(i__), + q_MIPS_3(i__),e_IRAC1_3(i__),e_IRAC2_3(i__),e_IRAC3_3(i__), + e_IRAC4_3(i__),e_MIPS_3(i__),f_IRAC1_3(i__),f_IRAC2_3(i__), + f_IRAC3_3(i__),f_IRAC4_3(i__),f_MIPS_3(i__),dDEN_3(i__), + Imag_3(i__),Jmag_3(i__),Kmag_3(i__),d2M_3(i__),J2mag_3(i__), + H2mag_3(i__),K2mag_3(i__),q_J2mag_3(i__),q_H2mag_3(i__), + q_K2mag_3(i__),dIRAS_3(i__),F12_3(i__),q_F12_3(i__), + F25_3(i__),q_F25_3(i__),F60_3(i__),q_F60_3(i__),F100_3(i__), + q_F100_3(i__),dMSX_3(i__),B1_3(i__),q_B1_3(i__),B2_3(i__), + q_B2_3(i__),A_3(i__),q_A_3(i__),C_3(i__),q_C_3(i__),D_3(i__), + q_D_3(i__),E_3(i__),q_E_3(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'bn1.dat' ! Catalogue of point sources in the Bulge N1 field C Format for file interpretation 6 format( + 3X,A18,3X,F10.6,3X,F10.6,3X,E12.6,3X,E12.6,3X,E12.6,3X,E12.6, + 3X,E12.6,3X,A1,A1,A1,A1,A1,3X,E12.6,3X,E12.6,3X,E12.6,3X, + E12.6,3X,E12.6,3X,I1,I1,I1,I1,I1,3X,F5.3,3X,F6.3,3X,F6.3,3X, + F6.3,3X,F5.3,3X,F6.3,3X,F6.3,3X,F6.3,3X,A1,A1,A1,3X,F5.1,3X, + E8.2,3X,I1,3X,E8.2,3X,I1,3X,E8.2,3X,I1,3X,E8.2,3X,I1,3X,F4.1, + 3X,E10.3,3X,I1,3X,E10.3,3X,I1,3X,E10.3,3X,I1,3X,E10.3,3X,I1, + 3X,E10.3,3X,I1,3X,E10.3,3X,I1) C Effective file loading open(unit=1,status='old',file= +'bn1.dat') write(6,*) '....Loading file: bn1.dat' do i__=1,64798 read(1,'(A466)')ar__5 read(ar__5,6) + USB_4(i__),RAdeg_5(i__),DEdeg_5(i__),IRAC1_4(i__), + IRAC2_4(i__),IRAC3_4(i__),IRAC4_4(i__),MIPS_4(i__), + q_IRAC1_4(i__),q_IRAC2_4(i__),q_IRAC3_4(i__),q_IRAC4_4(i__), + q_MIPS_4(i__),e_IRAC1_4(i__),e_IRAC2_4(i__),e_IRAC3_4(i__), + e_IRAC4_4(i__),e_MIPS_4(i__),f_IRAC1_4(i__),f_IRAC2_4(i__), + f_IRAC3_4(i__),f_IRAC4_4(i__),f_MIPS_4(i__),dDEN_4(i__), + Imag_4(i__),Jmag_4(i__),Kmag_4(i__),d2M_4(i__),J2mag_4(i__), + H2mag_4(i__),K2mag_4(i__),q_J2mag_4(i__),q_H2mag_4(i__), + q_K2mag_4(i__),dIRAS_4(i__),F12_4(i__),q_F12_4(i__), + F25_4(i__),q_F25_4(i__),F60_4(i__),q_F60_4(i__),F100_4(i__), + q_F100_4(i__),dMSX_4(i__),B1_4(i__),q_B1_4(i__),B2_4(i__), + q_B2_4(i__),A_4(i__),q_A_4(i__),C_4(i__),q_C_4(i__),D_4(i__), + q_D_4(i__),E_4(i__),q_E_4(i__) if(ar__5(217:221) .EQ. '') dDEN_4(i__) = rNULL__ if(ar__5(225:230) .EQ. '') Imag_4(i__) = rNULL__ if(ar__5(234:239) .EQ. '') Jmag_4(i__) = rNULL__ if(ar__5(243:248) .EQ. '') Kmag_4(i__) = rNULL__ if(ar__5(252:256) .EQ. '') d2M_4(i__) = rNULL__ if (idig(ar__5(260:265)).EQ.0) J2mag_4(i__) = rNULL__ if (idig(ar__5(269:274)).EQ.0) H2mag_4(i__) = rNULL__ if (idig(ar__5(278:283)).EQ.0) K2mag_4(i__) = rNULL__ if(ar__5(293:297) .EQ. '') dIRAS_4(i__) = rNULL__ if (idig(ar__5(301:308)).EQ.0) F12_4(i__) = rNULL__ if (q_F12_4(i__) .EQ. 45) q_F12_4(i__) = iNULL__ if (idig(ar__5(316:323)).EQ.0) F25_4(i__) = rNULL__ if (q_F25_4(i__) .EQ. 45) q_F25_4(i__) = iNULL__ if (idig(ar__5(331:338)).EQ.0) F60_4(i__) = rNULL__ if (q_F60_4(i__) .EQ. 45) q_F60_4(i__) = iNULL__ if (idig(ar__5(346:353)).EQ.0) F100_4(i__) = rNULL__ if (q_F100_4(i__) .EQ. 45) q_F100_4(i__) = iNULL__ if(ar__5(361:364) .EQ. '') dMSX_4(i__) = rNULL__ if (idig(ar__5(368:377)).EQ.0) B1_4(i__) = rNULL__ if (q_B1_4(i__) .EQ. 45) q_B1_4(i__) = iNULL__ if (idig(ar__5(385:394)).EQ.0) B2_4(i__) = rNULL__ if (q_B2_4(i__) .EQ. 45) q_B2_4(i__) = iNULL__ if (idig(ar__5(402:411)).EQ.0) A_4(i__) = rNULL__ if (q_A_4(i__) .EQ. 45) q_A_4(i__) = iNULL__ if (idig(ar__5(419:428)).EQ.0) C_4(i__) = rNULL__ if (q_C_4(i__) .EQ. 45) q_C_4(i__) = iNULL__ if (idig(ar__5(436:445)).EQ.0) D_4(i__) = rNULL__ if (q_D_4(i__) .EQ. 45) q_D_4(i__) = iNULL__ if (idig(ar__5(453:462)).EQ.0) E_4(i__) = rNULL__ if (q_E_4(i__) .EQ. 45) q_E_4(i__) = iNULL__ c ..............Just test output........... write(6,6) + USB_4(i__),RAdeg_5(i__),DEdeg_5(i__),IRAC1_4(i__), + IRAC2_4(i__),IRAC3_4(i__),IRAC4_4(i__),MIPS_4(i__), + q_IRAC1_4(i__),q_IRAC2_4(i__),q_IRAC3_4(i__),q_IRAC4_4(i__), + q_MIPS_4(i__),e_IRAC1_4(i__),e_IRAC2_4(i__),e_IRAC3_4(i__), + e_IRAC4_4(i__),e_MIPS_4(i__),f_IRAC1_4(i__),f_IRAC2_4(i__), + f_IRAC3_4(i__),f_IRAC4_4(i__),f_MIPS_4(i__),dDEN_4(i__), + Imag_4(i__),Jmag_4(i__),Kmag_4(i__),d2M_4(i__),J2mag_4(i__), + H2mag_4(i__),K2mag_4(i__),q_J2mag_4(i__),q_H2mag_4(i__), + q_K2mag_4(i__),dIRAS_4(i__),F12_4(i__),q_F12_4(i__), + F25_4(i__),q_F25_4(i__),F60_4(i__),q_F60_4(i__),F100_4(i__), + q_F100_4(i__),dMSX_4(i__),B1_4(i__),q_B1_4(i__),B2_4(i__), + q_B2_4(i__),A_4(i__),q_A_4(i__),C_4(i__),q_C_4(i__),D_4(i__), + q_D_4(i__),E_4(i__),q_E_4(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'bn2.dat' ! Catalogue of point sources in the Bulge N2 field C Format for file interpretation 7 format( + 3X,A18,3X,F10.6,3X,F10.6,3X,E12.6,3X,E12.6,3X,E12.6,3X,E12.6, + 3X,E12.6,3X,A1,A1,A1,A1,A1,3X,E12.6,3X,E12.6,3X,E12.6,3X, + E12.6,3X,E12.6,3X,I1,I1,I1,I1,I1,3X,F5.3,3X,F6.3,3X,F6.3,3X, + F6.3,3X,F5.3,3X,F6.3,3X,F6.3,3X,F6.3,3X,A1,A1,A1,3X,F5.1,3X, + E8.2,3X,I1,3X,E8.2,3X,I1,3X,E8.2,3X,I1,3X,E8.2,3X,I1,3X,F4.1, + 3X,E10.3,3X,I1,3X,E10.3,3X,I1,3X,E10.3,3X,I1,3X,E10.3,3X,I1, + 3X,E10.3,3X,I1,3X,E10.3,3X,I1) C Effective file loading open(unit=1,status='old',file= +'bn2.dat') write(6,*) '....Loading file: bn2.dat' do i__=1,428 read(1,'(A466)')ar__6 read(ar__6,7) + USB_5(i__),RAdeg_6(i__),DEdeg_6(i__),IRAC1_5(i__), + IRAC2_5(i__),IRAC3_5(i__),IRAC4_5(i__),MIPS_5(i__), + q_IRAC1_5(i__),q_IRAC2_5(i__),q_IRAC3_5(i__),q_IRAC4_5(i__), + q_MIPS_5(i__),e_IRAC1_5(i__),e_IRAC2_5(i__),e_IRAC3_5(i__), + e_IRAC4_5(i__),e_MIPS_5(i__),f_IRAC1_5(i__),f_IRAC2_5(i__), + f_IRAC3_5(i__),f_IRAC4_5(i__),f_MIPS_5(i__),dDEN_5(i__), + Imag_5(i__),Jmag_5(i__),Kmag_5(i__),d2M_5(i__),J2mag_5(i__), + H2mag_5(i__),K2mag_5(i__),q_J2mag_5(i__),q_H2mag_5(i__), + q_K2mag_5(i__),dIRAS_5(i__),F12_5(i__),q_F12_5(i__), + F25_5(i__),q_F25_5(i__),F60_5(i__),q_F60_5(i__),F100_5(i__), + q_F100_5(i__),dMSX_5(i__),B1_5(i__),q_B1_5(i__),B2_5(i__), + q_B2_5(i__),A_5(i__),q_A_5(i__),C_5(i__),q_C_5(i__),D_5(i__), + q_D_5(i__),E_5(i__),q_E_5(i__) if(ar__6(217:221) .EQ. '') dDEN_5(i__) = rNULL__ if(ar__6(225:230) .EQ. '') Imag_5(i__) = rNULL__ if(ar__6(234:239) .EQ. '') Jmag_5(i__) = rNULL__ if(ar__6(243:248) .EQ. '') Kmag_5(i__) = rNULL__ if(ar__6(252:256) .EQ. '') d2M_5(i__) = rNULL__ if (idig(ar__6(260:265)).EQ.0) J2mag_5(i__) = rNULL__ if (idig(ar__6(269:274)).EQ.0) H2mag_5(i__) = rNULL__ if (idig(ar__6(278:283)).EQ.0) K2mag_5(i__) = rNULL__ if(ar__6(293:297) .EQ. '') dIRAS_5(i__) = rNULL__ if (idig(ar__6(301:308)).EQ.0) F12_5(i__) = rNULL__ if (q_F12_5(i__) .EQ. 45) q_F12_5(i__) = iNULL__ if (idig(ar__6(316:323)).EQ.0) F25_5(i__) = rNULL__ if (q_F25_5(i__) .EQ. 45) q_F25_5(i__) = iNULL__ if (idig(ar__6(331:338)).EQ.0) F60_5(i__) = rNULL__ if (q_F60_5(i__) .EQ. 45) q_F60_5(i__) = iNULL__ if (idig(ar__6(346:353)).EQ.0) F100_5(i__) = rNULL__ if (q_F100_5(i__) .EQ. 45) q_F100_5(i__) = iNULL__ if(ar__6(361:364) .EQ. '') dMSX_5(i__) = rNULL__ if (idig(ar__6(368:377)).EQ.0) B1_5(i__) = rNULL__ if (q_B1_5(i__) .EQ. 45) q_B1_5(i__) = iNULL__ if (idig(ar__6(385:394)).EQ.0) B2_5(i__) = rNULL__ if (q_B2_5(i__) .EQ. 45) q_B2_5(i__) = iNULL__ if (idig(ar__6(402:411)).EQ.0) A_5(i__) = rNULL__ if (q_A_5(i__) .EQ. 45) q_A_5(i__) = iNULL__ if (idig(ar__6(419:428)).EQ.0) C_5(i__) = rNULL__ if (q_C_5(i__) .EQ. 45) q_C_5(i__) = iNULL__ if (idig(ar__6(436:445)).EQ.0) D_5(i__) = rNULL__ if (q_D_5(i__) .EQ. 45) q_D_5(i__) = iNULL__ if (idig(ar__6(453:462)).EQ.0) E_5(i__) = rNULL__ if (q_E_5(i__) .EQ. 45) q_E_5(i__) = iNULL__ c ..............Just test output........... write(6,7) + USB_5(i__),RAdeg_6(i__),DEdeg_6(i__),IRAC1_5(i__), + IRAC2_5(i__),IRAC3_5(i__),IRAC4_5(i__),MIPS_5(i__), + q_IRAC1_5(i__),q_IRAC2_5(i__),q_IRAC3_5(i__),q_IRAC4_5(i__), + q_MIPS_5(i__),e_IRAC1_5(i__),e_IRAC2_5(i__),e_IRAC3_5(i__), + e_IRAC4_5(i__),e_MIPS_5(i__),f_IRAC1_5(i__),f_IRAC2_5(i__), + f_IRAC3_5(i__),f_IRAC4_5(i__),f_MIPS_5(i__),dDEN_5(i__), + Imag_5(i__),Jmag_5(i__),Kmag_5(i__),d2M_5(i__),J2mag_5(i__), + H2mag_5(i__),K2mag_5(i__),q_J2mag_5(i__),q_H2mag_5(i__), + q_K2mag_5(i__),dIRAS_5(i__),F12_5(i__),q_F12_5(i__), + F25_5(i__),q_F25_5(i__),F60_5(i__),q_F60_5(i__),F100_5(i__), + q_F100_5(i__),dMSX_5(i__),B1_5(i__),q_B1_5(i__),B2_5(i__), + q_B2_5(i__),A_5(i__),q_A_5(i__),C_5(i__),q_C_5(i__),D_5(i__), + q_D_5(i__),E_5(i__),q_E_5(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'bngc6522.dat' ! Catalogue of point sources in the NGC 6522 field C Format for file interpretation 8 format( + 3X,A18,3X,F10.6,3X,F10.6,3X,E12.6,3X,E12.6,3X,E12.6,3X,E12.6, + 3X,E12.6,3X,A1,A1,A1,A1,A1,3X,E12.6,3X,E12.6,3X,E12.6,3X, + E12.6,3X,E12.6,3X,I1,I1,I1,I1,I1,3X,F5.3,3X,F6.3,3X,F6.3,3X, + F6.3,3X,F5.3,3X,F6.3,3X,F6.3,3X,F6.3,3X,A1,A1,A1,3X,F5.1,3X, + E8.2,3X,I1,3X,E8.2,3X,I1,3X,E8.2,3X,I1,3X,E8.2,3X,I1,3X,F4.1, + 3X,E10.3,3X,I1,3X,E10.3,3X,I1,3X,E10.3,3X,I1,3X,E10.3,3X,I1, + 3X,E10.3,3X,I1,3X,E10.3,3X,I1) C Effective file loading open(unit=1,status='old',file= +'bngc6522.dat') write(6,*) '....Loading file: bngc6522.dat' do i__=1,47232 read(1,'(A466)')ar__7 read(ar__7,8) + USB_6(i__),RAdeg_7(i__),DEdeg_7(i__),IRAC1_6(i__), + IRAC2_6(i__),IRAC3_6(i__),IRAC4_6(i__),MIPS_6(i__), + q_IRAC1_6(i__),q_IRAC2_6(i__),q_IRAC3_6(i__),q_IRAC4_6(i__), + q_MIPS_6(i__),e_IRAC1_6(i__),e_IRAC2_6(i__),e_IRAC3_6(i__), + e_IRAC4_6(i__),e_MIPS_6(i__),f_IRAC1_6(i__),f_IRAC2_6(i__), + f_IRAC3_6(i__),f_IRAC4_6(i__),f_MIPS_6(i__),dDEN_6(i__), + Imag_6(i__),Jmag_6(i__),Kmag_6(i__),d2M_6(i__),J2mag_6(i__), + H2mag_6(i__),K2mag_6(i__),q_J2mag_6(i__),q_H2mag_6(i__), + q_K2mag_6(i__),dIRAS_6(i__),F12_6(i__),q_F12_6(i__), + F25_6(i__),q_F25_6(i__),F60_6(i__),q_F60_6(i__),F100_6(i__), + q_F100_6(i__),dMSX_6(i__),B1_6(i__),q_B1_6(i__),B2_6(i__), + q_B2_6(i__),A_6(i__),q_A_6(i__),C_6(i__),q_C_6(i__),D_6(i__), + q_D_6(i__),E_6(i__),q_E_6(i__) if(ar__7(217:221) .EQ. '') dDEN_6(i__) = rNULL__ if(ar__7(225:230) .EQ. '') Imag_6(i__) = rNULL__ if(ar__7(234:239) .EQ. '') Jmag_6(i__) = rNULL__ if(ar__7(243:248) .EQ. '') Kmag_6(i__) = rNULL__ if(ar__7(252:256) .EQ. '') d2M_6(i__) = rNULL__ if (idig(ar__7(260:265)).EQ.0) J2mag_6(i__) = rNULL__ if (idig(ar__7(269:274)).EQ.0) H2mag_6(i__) = rNULL__ if (idig(ar__7(278:283)).EQ.0) K2mag_6(i__) = rNULL__ if(ar__7(293:297) .EQ. '') dIRAS_6(i__) = rNULL__ if (idig(ar__7(301:308)).EQ.0) F12_6(i__) = rNULL__ if (q_F12_6(i__) .EQ. 45) q_F12_6(i__) = iNULL__ if (idig(ar__7(316:323)).EQ.0) F25_6(i__) = rNULL__ if (q_F25_6(i__) .EQ. 45) q_F25_6(i__) = iNULL__ if (idig(ar__7(331:338)).EQ.0) F60_6(i__) = rNULL__ if (q_F60_6(i__) .EQ. 45) q_F60_6(i__) = iNULL__ if (idig(ar__7(346:353)).EQ.0) F100_6(i__) = rNULL__ if (q_F100_6(i__) .EQ. 45) q_F100_6(i__) = iNULL__ if(ar__7(361:364) .EQ. '') dMSX_6(i__) = rNULL__ if (idig(ar__7(368:377)).EQ.0) B1_6(i__) = rNULL__ if (q_B1_6(i__) .EQ. 45) q_B1_6(i__) = iNULL__ if (idig(ar__7(385:394)).EQ.0) B2_6(i__) = rNULL__ if (q_B2_6(i__) .EQ. 45) q_B2_6(i__) = iNULL__ if (idig(ar__7(402:411)).EQ.0) A_6(i__) = rNULL__ if (q_A_6(i__) .EQ. 45) q_A_6(i__) = iNULL__ if (idig(ar__7(419:428)).EQ.0) C_6(i__) = rNULL__ if (q_C_6(i__) .EQ. 45) q_C_6(i__) = iNULL__ if (idig(ar__7(436:445)).EQ.0) D_6(i__) = rNULL__ if (q_D_6(i__) .EQ. 45) q_D_6(i__) = iNULL__ if (idig(ar__7(453:462)).EQ.0) E_6(i__) = rNULL__ if (q_E_6(i__) .EQ. 45) q_E_6(i__) = iNULL__ c ..............Just test output........... write(6,8) + USB_6(i__),RAdeg_7(i__),DEdeg_7(i__),IRAC1_6(i__), + IRAC2_6(i__),IRAC3_6(i__),IRAC4_6(i__),MIPS_6(i__), + q_IRAC1_6(i__),q_IRAC2_6(i__),q_IRAC3_6(i__),q_IRAC4_6(i__), + q_MIPS_6(i__),e_IRAC1_6(i__),e_IRAC2_6(i__),e_IRAC3_6(i__), + e_IRAC4_6(i__),e_MIPS_6(i__),f_IRAC1_6(i__),f_IRAC2_6(i__), + f_IRAC3_6(i__),f_IRAC4_6(i__),f_MIPS_6(i__),dDEN_6(i__), + Imag_6(i__),Jmag_6(i__),Kmag_6(i__),d2M_6(i__),J2mag_6(i__), + H2mag_6(i__),K2mag_6(i__),q_J2mag_6(i__),q_H2mag_6(i__), + q_K2mag_6(i__),dIRAS_6(i__),F12_6(i__),q_F12_6(i__), + F25_6(i__),q_F25_6(i__),F60_6(i__),q_F60_6(i__),F100_6(i__), + q_F100_6(i__),dMSX_6(i__),B1_6(i__),q_B1_6(i__),B2_6(i__), + q_B2_6(i__),A_6(i__),q_A_6(i__),C_6(i__),q_C_6(i__),D_6(i__), + q_D_6(i__),E_6(i__),q_E_6(i__) 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