Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/175/147 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. II/295 Spitzer IRAC survey of the galactic center (Ramirez+, 2008) *================================================================================ *Point sources from a Spitzer IRAC survey of the galactic center. * Ramirez S.V., Arendt R.G., Sellgren K., Stolovy S.R., Cotera A., * Smith H.A., Yusef-Zadeh F. * <Astrophys. J. Suppl. Ser., 175, 147-164 (2008)> * =2008ApJS..175..147R C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table2.dat' ! Galactic Center IRAC 1 (3.6um) source lists integer*4 nr__ parameter (nr__=735020) ! Number of records character*106 ar__ ! Full-size record C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg (nr__) ! (deg) Right Ascension J2000 real*8 DEdeg (nr__) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) integer*4 Ch (nr__) ! [1/4] Channel number: * 1=3.6um, 2=4.5um, 3=5.8um, 4=8.0um integer*4 ID (nr__) ! Source identification number (1) integer*4 RAh (nr__) ! (h) Hour of Right Ascension (J2000) integer*4 RAm (nr__) ! (min) Minute of Right Ascension (J2000) real*4 RAs (nr__) ! (s) Second of Right Ascension (J2000) character*1 DE_ (nr__) ! Sign of the Declination (J2000) integer*4 DEd (nr__) ! (deg) Degree of Declination (J2000) integer*4 DEm (nr__) ! (arcmin) Arcminute of Declination (J2000) real*4 DEs (nr__) ! (arcsec) Arcsecond of Declination (J2000) real*8 GLON (nr__) ! (deg) Galactic longitude real*8 GLAT (nr__) ! (deg) Galactic latitude real*8 Flux (nr__) ! (mJy) Flux density in the channel wavelength real*8 e_Flux (nr__) ! (mJy) Uncertainty in Flux integer*4 N (nr__) ! Number of observations of source real*8 S_N (nr__) ! Signal-to-noise ratio character*3 Meth (nr__) ! Method used to compute the flux (2) *Note (1): Name in Simbad is [RAS2008] GC-IRACc-NNNNNN * where c is the channel number, and N the ID number *Note (2): Method as follows: * APC = Aperture Corrected Photometry. * PRF = Point Response Function Fitting. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table3.dat' ! Galactic Center IRAC 2 (4.5um) source lists integer*4 nr__1 parameter (nr__1=700923) ! Number of records character*106 ar__1 ! Full-size record C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg_1 (nr__1) ! (deg) Right Ascension J2000 real*8 DEdeg_1 (nr__1) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) integer*4 Ch_1 (nr__1) ! [1/4] Channel number: * 1=3.6um, 2=4.5um, 3=5.8um, 4=8.0um integer*4 ID_1 (nr__1) ! Source identification number (1) integer*4 RAh_1 (nr__1) ! (h) Hour of Right Ascension (J2000) integer*4 RAm_1 (nr__1) ! (min) Minute of Right Ascension (J2000) real*4 RAs_1 (nr__1) ! (s) Second of Right Ascension (J2000) character*1 DE__1 (nr__1) ! Sign of the Declination (J2000) integer*4 DEd_1 (nr__1) ! (deg) Degree of Declination (J2000) integer*4 DEm_1 (nr__1) ! (arcmin) Arcminute of Declination (J2000) real*4 DEs_1 (nr__1) ! (arcsec) Arcsecond of Declination (J2000) real*8 GLON_1 (nr__1) ! (deg) Galactic longitude real*8 GLAT_1 (nr__1) ! (deg) Galactic latitude real*8 Flux_1 (nr__1) ! (mJy) Flux density in the channel wavelength real*8 e_Flux_1 (nr__1) ! (mJy) Uncertainty in Flux integer*4 N_1 (nr__1) ! Number of observations of source real*8 S_N_1 (nr__1) ! Signal-to-noise ratio character*3 Meth_1 (nr__1) ! Method used to compute the flux (2) *Note (1): Name in Simbad is [RAS2008] GC-IRACc-NNNNNN * where c is the channel number, and N the ID number *Note (2): Method as follows: * APC = Aperture Corrected Photometry. * PRF = Point Response Function Fitting. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table4.dat' ! Galactic Center IRAC 3 (5.8um) source lists integer*4 nr__2 parameter (nr__2=493207) ! Number of records character*106 ar__2 ! Full-size record C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg_2 (nr__2) ! (deg) Right Ascension J2000 real*8 DEdeg_2 (nr__2) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) integer*4 Ch_2 (nr__2) ! [1/4] Channel number: * 1=3.6um, 2=4.5um, 3=5.8um, 4=8.0um integer*4 ID_2 (nr__2) ! Source identification number (1) integer*4 RAh_2 (nr__2) ! (h) Hour of Right Ascension (J2000) integer*4 RAm_2 (nr__2) ! (min) Minute of Right Ascension (J2000) real*4 RAs_2 (nr__2) ! (s) Second of Right Ascension (J2000) character*1 DE__2 (nr__2) ! Sign of the Declination (J2000) integer*4 DEd_2 (nr__2) ! (deg) Degree of Declination (J2000) integer*4 DEm_2 (nr__2) ! (arcmin) Arcminute of Declination (J2000) real*4 DEs_2 (nr__2) ! (arcsec) Arcsecond of Declination (J2000) real*8 GLON_2 (nr__2) ! (deg) Galactic longitude real*8 GLAT_2 (nr__2) ! (deg) Galactic latitude real*8 Flux_2 (nr__2) ! (mJy) Flux density in the channel wavelength real*8 e_Flux_2 (nr__2) ! (mJy) Uncertainty in Flux integer*4 N_2 (nr__2) ! Number of observations of source real*8 S_N_2 (nr__2) ! Signal-to-noise ratio character*3 Meth_2 (nr__2) ! Method used to compute the flux (2) *Note (1): Name in Simbad is [RAS2008] GC-IRACc-NNNNNN * where c is the channel number, and N the ID number *Note (2): Method as follows: * APC = Aperture Corrected Photometry. * PRF = Point Response Function Fitting. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table5.dat' ! Galactic Center IRAC 4 (8.0um) source lists integer*4 nr__3 parameter (nr__3=323512) ! Number of records character*106 ar__3 ! Full-size record C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg_3 (nr__3) ! (deg) Right Ascension J2000 real*8 DEdeg_3 (nr__3) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) integer*4 Ch_3 (nr__3) ! [1/4] Channel number: * 1=3.6um, 2=4.5um, 3=5.8um, 4=8.0um integer*4 ID_3 (nr__3) ! Source identification number (1) integer*4 RAh_3 (nr__3) ! (h) Hour of Right Ascension (J2000) integer*4 RAm_3 (nr__3) ! (min) Minute of Right Ascension (J2000) real*4 RAs_3 (nr__3) ! (s) Second of Right Ascension (J2000) character*1 DE__3 (nr__3) ! Sign of the Declination (J2000) integer*4 DEd_3 (nr__3) ! (deg) Degree of Declination (J2000) integer*4 DEm_3 (nr__3) ! (arcmin) Arcminute of Declination (J2000) real*4 DEs_3 (nr__3) ! (arcsec) Arcsecond of Declination (J2000) real*8 GLON_3 (nr__3) ! (deg) Galactic longitude real*8 GLAT_3 (nr__3) ! (deg) Galactic latitude real*8 Flux_3 (nr__3) ! (mJy) Flux density in the channel wavelength real*8 e_Flux_3 (nr__3) ! (mJy) Uncertainty in Flux integer*4 N_3 (nr__3) ! Number of observations of source real*8 S_N_3 (nr__3) ! Signal-to-noise ratio character*3 Meth_3 (nr__3) ! Method used to compute the flux (2) *Note (1): Name in Simbad is [RAS2008] GC-IRACc-NNNNNN * where c is the channel number, and N the ID number *Note (2): Method as follows: * APC = Aperture Corrected Photometry. * PRF = Point Response Function Fitting. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table6.dat' ! Galactic Center IRAC sub-array photometry integer*4 nr__4 parameter (nr__4=207) ! Number of records character*108 ar__4 ! Full-size record C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg_4 (nr__4) ! (deg) Right Ascension J2000 real*8 DEdeg_4 (nr__4) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) integer*4 ID_4 (nr__4) ! Source identification (NNN) character*5 Field (nr__4) ! IRAC sub-array field (1) integer*4 RAh_4 (nr__4) ! (h) Hour of Right Ascension (J2000) integer*4 RAm_4 (nr__4) ! (min) Minute of Right Ascension (J2000) real*4 RAs_4 (nr__4) ! (s) Second of Right Ascension (J2000) character*1 DE__4 (nr__4) ! Sign of the Declination (J2000) integer*4 DEd_4 (nr__4) ! (deg) Degree of Declination (J2000) integer*4 DEm_4 (nr__4) ! (arcmin) Arcminute of Declination (J2000) real*4 DEs_4 (nr__4) ! (arcsec) Arcsecond of Declination (J2000) real*8 GLON_4 (nr__4) ! (deg) Galactic longitude real*8 GLAT_4 (nr__4) ! (deg) Galactic latitude real*4 v3_6mag (nr__4) ! (mag) ? IRAC 3.6 micron band magnitude real*4 e_3_6mag (nr__4) ! (mag) ? Uncertainty in 3.6mag real*4 v4_5mag (nr__4) ! (mag) ? IRAC 4.5 micron band magnitude real*4 e_4_5mag (nr__4) ! (mag) ? Uncertainty in 4.5mag real*4 v5_8mag (nr__4) ! (mag) ? IRAC 5.8 micron band magnitude real*4 e_5_8mag (nr__4) ! (mag) ? Uncertainty in 5.8mag real*4 v8_0mag (nr__4) ! (mag) ? IRAC 8.0 micron band magnitude real*4 e_8_0mag (nr__4) ! (mag) ? Uncertainty in 8.0mag *Note (1): Note as follows: * sat = individual sub-array pointings. * sgra = sub-array field centered on SgrA* * quint = sub-array field centered on the Quintuplet cluster c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table7.dat' ! Galactic Center 2MASS/IRAC Catalog integer*4 nr__5 parameter (nr__5=1065565) ! Number of records character*340 ar__5 ! Full-size record C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg_5 (nr__5) ! (deg) Right Ascension J2000 real*8 DEdeg_5 (nr__5) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) integer*4 SSTGC (nr__5) ! Source identification integer*4 RAh_5 (nr__5) ! (h) Hour of Right Ascension (J2000) integer*4 RAm_5 (nr__5) ! (min) Minute of Right Ascension (J2000) real*4 RAs_5 (nr__5) ! (s) Second of Right Ascension (J2000) character*1 DE__5 (nr__5) ! Sign of the Declination (J2000) integer*4 DEd_5 (nr__5) ! (deg) Degree of Declination (J2000) integer*4 DEm_5 (nr__5) ! (arcmin) Arcminute of Declination (J2000) real*4 DEs_5 (nr__5) ! (arcsec) Arcsecond of Declination (J2000) real*8 GLON_5 (nr__5) ! (deg) Galactic longitude real*8 GLAT_5 (nr__5) ! (deg) Galactic latitude integer*4 v2_1 (nr__5) ! [0,1] satisfies the "2+1" criterion flag (1) integer*4 loc (nr__5) ! [0,3] Note about source location (2) character*17 v2MASS (nr__5) ! 2MASS identification (HHMMSSss+DDMMSSsA) real*4 Jmag (nr__5) ! (mag) ?=-9.999 2MASS J band magnitude real*4 e_Jmag (nr__5) ! (mag) ?=-9.999 Uncertainty in Jmag real*8 Hmag (nr__5) ! (mag) ?=-9.999 2MASS H band magnitude real*4 e_Hmag (nr__5) ! (mag) ?=-9.999 Uncertainty in Hmag real*8 Ksmag (nr__5) ! (mag) ?=-9.999 2MASS K_S_ band magnitude real*4 e_Ksmag (nr__5) ! (mag) ?=-9.999 Uncertainty in Ksmag character*3 q2M (nr__5) ! [ABCDEFUXZ] 2MASS quality flag in J H K (3) integer*4 ID1 (nr__5) ! ? Identification in Channel 1 (3.6um) real*4 v3_6mag_1 (nr__5) ! (mag) ?=-9.999 Channel 1 magnitude at 3.6um real*4 e_3_6mag_1 (nr__5) ! (mag) ?=-9.999 Uncertainty in 3.6mag real*8 SNR1 (nr__5) ! ?=-9.9 Channel 1 signal-to-noise integer*4 f1 (nr__5) ! [0,4] Channel 1 flag (4) integer*4 N1 (nr__5) ! Channel 1 Coverage (#frames at source position) real*4 M_N1 (nr__5) ! Channel 1 M/N ratio (5) integer*4 ID2 (nr__5) ! ? Identification in Channel 2 (4.5um) real*4 v4_5mag_1 (nr__5) ! (mag) ?=-9.999 Channel 2 magnitude at 4.5um real*4 e_4_5mag_1 (nr__5) ! (mag) ?=-9.999 Uncertainty in 4.5mag real*8 SNR2 (nr__5) ! ?=-9.9 Channel 2 signal-to-noise integer*4 f2 (nr__5) ! [0,4] Channel 2 flag (4) integer*4 N2 (nr__5) ! Channel 2 Coverage (#frames at source position) real*4 M_N2 (nr__5) ! Channel 2 M/N ratio (5) integer*4 ID3 (nr__5) ! ? Identification in Channel 3 (5.8um) real*4 v5_8mag_1 (nr__5) ! (mag) ?=-9.999 Channel 3 magnitude at 5.8um real*4 e_5_8mag_1 (nr__5) ! (mag) ?=-9.999 Uncertainty in 5.8mag real*8 SNR3 (nr__5) ! ?=-9.9 Channel 3 signal-to-noise integer*4 f3 (nr__5) ! [0,4] Channel 3 flag (4) integer*4 N3 (nr__5) ! Channel 3 Coverage (#frames at source position) real*4 M_N3 (nr__5) ! Channel 3 M/N ratio (5) integer*4 ID4 (nr__5) ! ? Identification in Channel 4 (8.0um) real*4 v8_0mag_1 (nr__5) ! (mag) ?=-9.999 Channel 4 magnitude at 8.0um real*4 e_8_0mag_1 (nr__5) ! (mag) ?=-9.999 Uncertainty in 8.0mag real*8 SNR4 (nr__5) ! ?=-9.9 Channel 4 signal-to-noise integer*4 f4 (nr__5) ! [0,4] Channel 4 flag (4) integer*4 N4 (nr__5) ! Channel 4 Coverage (#frames at source position) real*4 M_N4 (nr__5) ! Channel 4 M/N ratio (5) character*7 ID6 (nr__5) ! Name from table6 (sub-array photometry) (6) *Note (1): the value '1' indicates that the source satisfies the GLIMPSE * "2+1" criterion: M/N>=0.6 in one IRAC band, M/N>=0.4 in an adjacent * band where M is the number of detections and N, the number of possible * observations. *Note (2): Flag as follows: * 0 = source is located in areas of incomplete coverage. * 2 = source is located within 40 arcmin on the Central Cluster. * 3 = source is located within 42 arcmin of the Quintuplet Cluster. * 1 = none of the above *Note (3): JHK photometric quality flag from the 2MASS Catalog as follows: * A = SNR>10 * B = SNR>7 * C = SNR>5 * D = no SNR requirement * E = poor profile-fit photometry * F = detection without photometric uncertainty * U = detection with upper limit on magnitude * X = detection without brightness estimate * ZZZ = no 2MASS counterpart *Note (4): Flag as follows: * 0 = source is not detected at the corresponding IRAC band. * 1 = magnitude comes from PRF fitting algorithm. * 2 = magnitude comes from aperture corrected measurement. * 3 = magnitude is greater than corresponding saturation limit. * 4 = magnitude comes from sub-array photometry. *Note (5): Ratio between M number of detections over N number of possible * observations (coverage). *Note (6): for sources observed in the sub-array photometry, the name is * * satNNN = source in individual sub-array pointings * * sgraNNN = source in sub-array field centered on SgrA* * * qintNNN = source in sub-array field centered on Quintuplet cluster C============================================================================= C Loading file 'table2.dat' ! Galactic Center IRAC 1 (3.6um) source lists C Format for file interpretation 1 format( + 7X,I1,1X,I6,3X,I2,1X,I2,1X,F5.2,2X,A1,I2,1X,I2,1X,F4.1,2X, + F12.8,1X,F11.8,1X,F10.3,1X,F9.3,1X,I2,1X,F9.1,1X,A3) C Effective file loading open(unit=1,status='old',file= +'table2.dat') write(6,*) '....Loading file: table2.dat' do i__=1,735020 read(1,'(A106)')ar__ read(ar__,1) + Ch(i__),ID(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__), + DEm(i__),DEs(i__),GLON(i__),GLAT(i__),Flux(i__),e_Flux(i__), + N(i__),S_N(i__),Meth(i__) RAdeg(i__) = rNULL__ DEdeg(i__) = rNULL__ c Derive coordinates RAdeg and DEdeg from input data c (RAdeg and DEdeg are set to rNULL__ when unknown) if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15. if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4. if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240. if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__) if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60. if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600. if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__) c ..............Just test output........... write(6,1) + Ch(i__),ID(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__), + DEm(i__),DEs(i__),GLON(i__),GLAT(i__),Flux(i__),e_Flux(i__), + N(i__),S_N(i__),Meth(i__) write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table3.dat' ! Galactic Center IRAC 2 (4.5um) source lists C Format for file interpretation 2 format( + 7X,I1,1X,I6,3X,I2,1X,I2,1X,F5.2,2X,A1,I2,1X,I2,1X,F4.1,2X, + F12.8,1X,F11.8,1X,F10.3,1X,F9.3,1X,I2,1X,F9.1,1X,A3) C Effective file loading open(unit=1,status='old',file= +'table3.dat') write(6,*) '....Loading file: table3.dat' do i__=1,700923 read(1,'(A106)')ar__1 read(ar__1,2) + Ch_1(i__),ID_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__), + DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),GLON_1(i__), + GLAT_1(i__),Flux_1(i__),e_Flux_1(i__),N_1(i__),S_N_1(i__), + Meth_1(i__) RAdeg_1(i__) = rNULL__ DEdeg_1(i__) = rNULL__ c Derive coordinates RAdeg_1 and DEdeg_1 from input data c (RAdeg_1 and DEdeg_1 are set to rNULL__ when unknown) if(RAh_1(i__) .GT. -180) RAdeg_1(i__)=RAh_1(i__)*15. if(RAm_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAm_1(i__)/4. if(RAs_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAs_1(i__)/240. if(DEd_1(i__) .GE. 0) DEdeg_1(i__)=DEd_1(i__) if(DEm_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEm_1(i__)/60. if(DEs_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEs_1(i__)/3600. if(DE__1(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__) c ..............Just test output........... write(6,2) + Ch_1(i__),ID_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__), + DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),GLON_1(i__), + GLAT_1(i__),Flux_1(i__),e_Flux_1(i__),N_1(i__),S_N_1(i__), + Meth_1(i__) write(6,'(6H Pos: 2F8.4)') RAdeg_1(i__),DEdeg_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table4.dat' ! Galactic Center IRAC 3 (5.8um) source lists C Format for file interpretation 3 format( + 7X,I1,1X,I6,3X,I2,1X,I2,1X,F5.2,2X,A1,I2,1X,I2,1X,F4.1,2X, + F12.8,1X,F11.8,1X,F10.3,1X,F9.3,1X,I2,1X,F9.1,1X,A3) C Effective file loading open(unit=1,status='old',file= +'table4.dat') write(6,*) '....Loading file: table4.dat' do i__=1,493207 read(1,'(A106)')ar__2 read(ar__2,3) + Ch_2(i__),ID_2(i__),RAh_2(i__),RAm_2(i__),RAs_2(i__), + DE__2(i__),DEd_2(i__),DEm_2(i__),DEs_2(i__),GLON_2(i__), + GLAT_2(i__),Flux_2(i__),e_Flux_2(i__),N_2(i__),S_N_2(i__), + Meth_2(i__) RAdeg_2(i__) = rNULL__ DEdeg_2(i__) = rNULL__ c Derive coordinates RAdeg_2 and DEdeg_2 from input data c (RAdeg_2 and DEdeg_2 are set to rNULL__ when unknown) if(RAh_2(i__) .GT. -180) RAdeg_2(i__)=RAh_2(i__)*15. if(RAm_2(i__) .GT. -180) RAdeg_2(i__)=RAdeg_2(i__)+RAm_2(i__)/4. if(RAs_2(i__) .GT. -180) RAdeg_2(i__)=RAdeg_2(i__)+RAs_2(i__)/240. if(DEd_2(i__) .GE. 0) DEdeg_2(i__)=DEd_2(i__) if(DEm_2(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEm_2(i__)/60. if(DEs_2(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEs_2(i__)/3600. if(DE__2(i__).EQ.'-'.AND.DEdeg_2(i__).GE.0) DEdeg_2(i__)=-DEdeg_2(i__) c ..............Just test output........... write(6,3) + Ch_2(i__),ID_2(i__),RAh_2(i__),RAm_2(i__),RAs_2(i__), + DE__2(i__),DEd_2(i__),DEm_2(i__),DEs_2(i__),GLON_2(i__), + GLAT_2(i__),Flux_2(i__),e_Flux_2(i__),N_2(i__),S_N_2(i__), + Meth_2(i__) write(6,'(6H Pos: 2F8.4)') RAdeg_2(i__),DEdeg_2(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table5.dat' ! Galactic Center IRAC 4 (8.0um) source lists C Format for file interpretation 4 format( + 7X,I1,1X,I6,3X,I2,1X,I2,1X,F5.2,2X,A1,I2,1X,I2,1X,F4.1,2X, + F12.8,1X,F11.8,1X,F10.3,1X,F9.3,1X,I2,1X,F9.1,1X,A3) C Effective file loading open(unit=1,status='old',file= +'table5.dat') write(6,*) '....Loading file: table5.dat' do i__=1,323512 read(1,'(A106)')ar__3 read(ar__3,4) + Ch_3(i__),ID_3(i__),RAh_3(i__),RAm_3(i__),RAs_3(i__), + DE__3(i__),DEd_3(i__),DEm_3(i__),DEs_3(i__),GLON_3(i__), + GLAT_3(i__),Flux_3(i__),e_Flux_3(i__),N_3(i__),S_N_3(i__), + Meth_3(i__) RAdeg_3(i__) = rNULL__ DEdeg_3(i__) = rNULL__ c Derive coordinates RAdeg_3 and DEdeg_3 from input data c (RAdeg_3 and DEdeg_3 are set to rNULL__ when unknown) if(RAh_3(i__) .GT. -180) RAdeg_3(i__)=RAh_3(i__)*15. if(RAm_3(i__) .GT. -180) RAdeg_3(i__)=RAdeg_3(i__)+RAm_3(i__)/4. if(RAs_3(i__) .GT. -180) RAdeg_3(i__)=RAdeg_3(i__)+RAs_3(i__)/240. if(DEd_3(i__) .GE. 0) DEdeg_3(i__)=DEd_3(i__) if(DEm_3(i__) .GE. 0) DEdeg_3(i__)=DEdeg_3(i__)+DEm_3(i__)/60. if(DEs_3(i__) .GE. 0) DEdeg_3(i__)=DEdeg_3(i__)+DEs_3(i__)/3600. if(DE__3(i__).EQ.'-'.AND.DEdeg_3(i__).GE.0) DEdeg_3(i__)=-DEdeg_3(i__) c ..............Just test output........... write(6,4) + Ch_3(i__),ID_3(i__),RAh_3(i__),RAm_3(i__),RAs_3(i__), + DE__3(i__),DEd_3(i__),DEm_3(i__),DEs_3(i__),GLON_3(i__), + GLAT_3(i__),Flux_3(i__),e_Flux_3(i__),N_3(i__),S_N_3(i__), + Meth_3(i__) write(6,'(6H Pos: 2F8.4)') RAdeg_3(i__),DEdeg_3(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table6.dat' ! Galactic Center IRAC sub-array photometry C Format for file interpretation 5 format( + I3,1X,A5,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F12.8, + 1X,F11.8,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F5.3,1X,F5.3,1X, + F5.3,1X,F5.3) C Effective file loading open(unit=1,status='old',file= +'table6.dat') write(6,*) '....Loading file: table6.dat' do i__=1,207 read(1,'(A108)')ar__4 read(ar__4,5) + ID_4(i__),Field(i__),RAh_4(i__),RAm_4(i__),RAs_4(i__), + DE__4(i__),DEd_4(i__),DEm_4(i__),DEs_4(i__),GLON_4(i__), + GLAT_4(i__),v3_6mag(i__),e_3_6mag(i__),v4_5mag(i__), + e_4_5mag(i__),v5_8mag(i__),e_5_8mag(i__),v8_0mag(i__), + e_8_0mag(i__) if(ar__4(60:65) .EQ. '') v3_6mag(i__) = rNULL__ if(ar__4(67:71) .EQ. '') e_3_6mag(i__) = rNULL__ if(ar__4(73:78) .EQ. '') v4_5mag(i__) = rNULL__ if(ar__4(80:84) .EQ. '') e_4_5mag(i__) = rNULL__ if(ar__4(86:90) .EQ. '') v5_8mag(i__) = rNULL__ if(ar__4(92:96) .EQ. '') e_5_8mag(i__) = rNULL__ if(ar__4(98:102) .EQ. '') v8_0mag(i__) = rNULL__ if(ar__4(104:108) .EQ. '') e_8_0mag(i__) = rNULL__ RAdeg_4(i__) = rNULL__ DEdeg_4(i__) = rNULL__ c Derive coordinates RAdeg_4 and DEdeg_4 from input data c (RAdeg_4 and DEdeg_4 are set to rNULL__ when unknown) if(RAh_4(i__) .GT. -180) RAdeg_4(i__)=RAh_4(i__)*15. if(RAm_4(i__) .GT. -180) RAdeg_4(i__)=RAdeg_4(i__)+RAm_4(i__)/4. if(RAs_4(i__) .GT. -180) RAdeg_4(i__)=RAdeg_4(i__)+RAs_4(i__)/240. if(DEd_4(i__) .GE. 0) DEdeg_4(i__)=DEd_4(i__) if(DEm_4(i__) .GE. 0) DEdeg_4(i__)=DEdeg_4(i__)+DEm_4(i__)/60. if(DEs_4(i__) .GE. 0) DEdeg_4(i__)=DEdeg_4(i__)+DEs_4(i__)/3600. if(DE__4(i__).EQ.'-'.AND.DEdeg_4(i__).GE.0) DEdeg_4(i__)=-DEdeg_4(i__) c ..............Just test output........... write(6,5) + ID_4(i__),Field(i__),RAh_4(i__),RAm_4(i__),RAs_4(i__), + DE__4(i__),DEd_4(i__),DEm_4(i__),DEs_4(i__),GLON_4(i__), + GLAT_4(i__),v3_6mag(i__),e_3_6mag(i__),v4_5mag(i__), + e_4_5mag(i__),v5_8mag(i__),e_5_8mag(i__),v8_0mag(i__), + e_8_0mag(i__) write(6,'(6H Pos: 2F8.4)') RAdeg_4(i__),DEdeg_4(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table7.dat' ! Galactic Center 2MASS/IRAC Catalog C Format for file interpretation 6 format( + 6X,I7,1X,I2,1X,I2,1X,F5.2,2X,A1,I2,1X,I2,1X,F4.1,2X,F9.5,1X, + F8.5,1X,I1,1X,I1,1X,A17,1X,F6.3,1X,F6.3,1X,F8.3,1X,F6.3,1X, + F8.3,1X,F6.3,1X,A3,10X,I6,2X,F6.3,1X,F6.3,1X,F8.1,1X,I1,1X,I2, + 1X,F4.2,10X,I6,2X,F6.3,1X,F6.3,1X,F9.1,1X,I1,1X,I2,1X,F4.2, + 10X,I6,2X,F6.3,1X,F6.3,1X,F8.1,1X,I1,1X,I2,1X,F4.2,10X,I6,2X, + F6.3,1X,F6.3,1X,F8.1,1X,I1,1X,I2,1X,F4.2,2X,A7) C Effective file loading open(unit=1,status='old',file= +'table7.dat') write(6,*) '....Loading file: table7.dat' do i__=1,1065565 read(1,'(A340)')ar__5 read(ar__5,6) + SSTGC(i__),RAh_5(i__),RAm_5(i__),RAs_5(i__),DE__5(i__), + DEd_5(i__),DEm_5(i__),DEs_5(i__),GLON_5(i__),GLAT_5(i__), + v2_1(i__),loc(i__),v2MASS(i__),Jmag(i__),e_Jmag(i__), + Hmag(i__),e_Hmag(i__),Ksmag(i__),e_Ksmag(i__),q2M(i__), + ID1(i__),v3_6mag_1(i__),e_3_6mag_1(i__),SNR1(i__),f1(i__), + N1(i__),M_N1(i__),ID2(i__),v4_5mag_1(i__),e_4_5mag_1(i__), + SNR2(i__),f2(i__),N2(i__),M_N2(i__),ID3(i__),v5_8mag_1(i__), + e_5_8mag_1(i__),SNR3(i__),f3(i__),N3(i__),M_N3(i__),ID4(i__), + v8_0mag_1(i__),e_8_0mag_1(i__),SNR4(i__),f4(i__),N4(i__), + M_N4(i__),ID6(i__) if(ar__5(141:146) .EQ. '') ID1(i__) = iNULL__ if(ar__5(191:196) .EQ. '') ID2(i__) = iNULL__ if(ar__5(242:247) .EQ. '') ID3(i__) = iNULL__ if(ar__5(292:297) .EQ. '') ID4(i__) = iNULL__ RAdeg_5(i__) = rNULL__ DEdeg_5(i__) = rNULL__ c Derive coordinates RAdeg_5 and DEdeg_5 from input data c (RAdeg_5 and DEdeg_5 are set to rNULL__ when unknown) if(RAh_5(i__) .GT. -180) RAdeg_5(i__)=RAh_5(i__)*15. if(RAm_5(i__) .GT. -180) RAdeg_5(i__)=RAdeg_5(i__)+RAm_5(i__)/4. if(RAs_5(i__) .GT. -180) RAdeg_5(i__)=RAdeg_5(i__)+RAs_5(i__)/240. if(DEd_5(i__) .GE. 0) DEdeg_5(i__)=DEd_5(i__) if(DEm_5(i__) .GE. 0) DEdeg_5(i__)=DEdeg_5(i__)+DEm_5(i__)/60. if(DEs_5(i__) .GE. 0) DEdeg_5(i__)=DEdeg_5(i__)+DEs_5(i__)/3600. if(DE__5(i__).EQ.'-'.AND.DEdeg_5(i__).GE.0) DEdeg_5(i__)=-DEdeg_5(i__) c ..............Just test output........... write(6,6) + SSTGC(i__),RAh_5(i__),RAm_5(i__),RAs_5(i__),DE__5(i__), + DEd_5(i__),DEm_5(i__),DEs_5(i__),GLON_5(i__),GLAT_5(i__), + v2_1(i__),loc(i__),v2MASS(i__),Jmag(i__),e_Jmag(i__), + Hmag(i__),e_Hmag(i__),Ksmag(i__),e_Ksmag(i__),q2M(i__), + ID1(i__),v3_6mag_1(i__),e_3_6mag_1(i__),SNR1(i__),f1(i__), + N1(i__),M_N1(i__),ID2(i__),v4_5mag_1(i__),e_4_5mag_1(i__), + SNR2(i__),f2(i__),N2(i__),M_N2(i__),ID3(i__),v5_8mag_1(i__), + e_5_8mag_1(i__),SNR3(i__),f3(i__),N3(i__),M_N3(i__),ID4(i__), + v8_0mag_1(i__),e_8_0mag_1(i__),SNR4(i__),f4(i__),N4(i__), + M_N4(i__),ID6(i__) write(6,'(6H Pos: 2F8.4)') RAdeg_5(i__),DEdeg_5(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end