Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/702/178 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/ApJ/702/178 YSOs in the central 400pc of the Galaxy (Yusef-Zadeh+, 2009) *================================================================================ *Star formation in the central 400 pc of the Milky Way: *evidence for a population of massive young stellar objects. * Yusef-Zadeh F., Hewitt J.W., Arendt R.G., Whitney B., Rieke G., Wardle M., * Hinz J.L., Stolovy S., Lang C.C., Burton M.G., Ramirez S. * <Astrophys. J., 702, 178-225 (2009)> * =2009ApJ...702..178Y C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table1.dat' ! Galactic coordinates and flux density of sources in G359.43+0.02 integer*4 nr__ parameter (nr__=18) ! Number of records character*61 ar__ ! Full-size record C Position composed of: GLON GLAT integer*4 v_YHA2009_ (nr__) ! Source number character*1 m__YHA2009_(nr__) ! [e] e: Extended component character*1 n__YHA2009_(nr__) ! [abc] Individual notes (1) real*8 GLON (nr__) ! (deg) Galactic longitude real*8 GLAT (nr__) ! (deg) Galactic latitude real*4 S6cm (nr__) ! ? VLA 6cm flux density real*4 e_S6cm (nr__) ! ? S6cm uncertainty real*4 S20cm (nr__) ! ? VLA 20cm flux density real*4 e_S20cm (nr__) ! ? S20cm uncertainty real*4 amaj (nr__) ! (arcsec) ? Major axis real*4 amin (nr__) ! (arcsec) ? Minor axis integer*4 Surf (nr__) ! (arcsec2) ? Surface when circular *Note (1): Flag as follows: * a = Infrared counterpart 300758. * b = Infrared counterpart 284291. * c = Source 11 has a central compact source and an extended component. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table3.dat' ! parameters of the SED fits to candidate YSOs in l<0{deg} integer*4 nr__1 parameter (nr__1=223) ! Number of records character*88 ar__1 ! Full-size record C Position composed of: GLON GLAT integer*4 SSTGC (nr__1) ! SST (Spitzer Space Telescope) Galactic * center (GC) number (1) character*3 n_SSTGC (nr__1) ! [abc, ] Notes for individual sources (2) real*8 GLON_1 (nr__1) ! (deg) Galactic longitude real*8 GLAT_1 (nr__1) ! (deg) Galactic latitude real*4 Chi2 (nr__1) ! {chi}^2^ value integer*4 Nfits (nr__1) ! Number of acceptably fitting models real*4 Av (nr__1) ! (mag) Visual extinction A_v_ real*4 Mass (nr__1) ! (Msun) Average mass of the protostar real*4 e_Mass (nr__1) ! (Msun) Mass uncertainty real*4 Lum (nr__1) ! (Lsun) Average luminosity of the protostar real*4 e_Lum (nr__1) ! (Lsun) Lum uncertainty real*4 Menv (nr__1) ! (Msun/yr) Average mass accretion rate real*4 e_Menv (nr__1) ! (Msun/yr) Menv 1{sigma} standard deviation character*3 Stage (nr__1) ! Evolutionary stage of the fitted source * (I, II or III) *Note (1): Based on the IRAC catalogs (Ramirez et al. 2008, Cat. II/295; * <SSTGC NNNNNNN> in Simbad) *Note (2): Flag as follows: * a = radio counterpart 6 * b = radio counterpart 5 * c = Candidate YSO in the G359.43+0.02 cluster (table 2 of the paper; note * that "GLAT" for SSTGC 307175 in table 2 is not 0.02637 but, indeed, * -0.00026 as in table 3). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table4.dat' ! parameters of the SED fits to candidate YSOs in l>0{deg} integer*4 nr__2 parameter (nr__2=137) ! Number of records character*88 ar__2 ! Full-size record C Position composed of: GLON GLAT integer*4 SSTGC_1 (nr__2) ! SST (Spitzer Space Telescope) Galactic * center (GC) number (1) character*3 n_SSTGC_1 (nr__2) ! [abc, ] Notes for individual sources (2) real*8 GLON_2 (nr__2) ! (deg) Galactic longitude real*8 GLAT_2 (nr__2) ! (deg) Galactic latitude real*4 Chi2_1 (nr__2) ! {chi}^2^ value integer*4 Nfits_1 (nr__2) ! Number of acceptably fitting models real*4 Av_1 (nr__2) ! (mag) Visual extinction A_v_ real*4 Mass_1 (nr__2) ! (Msun) Average mass of the protostar real*4 e_Mass_1 (nr__2) ! (Msun) Mass uncertainty real*4 Lum_1 (nr__2) ! (Lsun) Average luminosity of the protostar real*4 e_Lum_1 (nr__2) ! (Lsun) Lum uncertainty real*4 Menv_1 (nr__2) ! (Msun/yr) Average mass accretion rate real*4 e_Menv_1 (nr__2) ! (Msun/yr) Menv 1{sigma} standard deviation character*3 Stage_1 (nr__2) ! Evolutionary stage of the fitted source * (I, II or III) *Note (1): Based on the IRAC catalogs (Ramirez et al. 2008, Cat. II/295; * <SSTGC NNNNNNN> in Simbad) *Note (2): Flag as follows: * a = radio counterpart 6 * b = radio counterpart 5 * c = Candidate YSO in the G359.43+0.02 cluster (table 2 of the paper; note * that "GLAT" for SSTGC 307175 in table 2 is not 0.02637 but, indeed, * -0.00026 as in table 3). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table5.dat' ! Table of sources with 4.5um Excess Emission integer*4 nr__3 parameter (nr__3=33) ! Number of records character*69 ar__3 ! Full-size record C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg (nr__3) ! (deg) Right Ascension J2000 real*8 DEdeg (nr__3) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) integer*4 ID (nr__3) ! Green source identification * (<[YHA2009] gNN> in Simbad) (1) real*8 GLON_3 (nr__3) ! (deg) Galactic longitude real*4 GLAT_3 (nr__3) ! (deg) Galactic latitude integer*4 RAh (nr__3) ! (h) Hour of right ascension (J2000) integer*4 RAm (nr__3) ! (min) Minute of right ascension (J2000) real*4 RAs (nr__3) ! (s) Second of right ascension (J2000) character*1 DE_ (nr__3) ! Sign of declination (J2000) integer*4 DEd (nr__3) ! (deg) Degree of declination (J2000) integer*4 DEm (nr__3) ! (arcmin) Arcminute of declination (J2000) real*4 DEs (nr__3) ! (arcsec) Arcsecond of declination (J2000) character*27 Comm (nr__3) ! Description *Note (1): Visually identified from IRAC data. Many point sources may be * multiple or extended sources at higher resolution. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table6.dat' ! Table of 70um Sources integer*4 nr__4 parameter (nr__4=9) ! Number of records character*58 ar__4 ! Full-size record C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg_1 (nr__4) ! (deg) Right Ascension J2000 real*8 DEdeg_1 (nr__4) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) integer*4 ID_1 (nr__4) ! Source identification from table 5, or from * table8 when n_ID is 'a' (1) character*1 n_ID (nr__4) ! [a] Methanol class II maser source (2) real*8 GLON_4 (nr__4) ! (deg) Galactic longitude real*8 GLAT_4 (nr__4) ! (deg) Galactic latitude integer*4 RAh_1 (nr__4) ! (h) Hour of right ascension (J2000) integer*4 RAm_1 (nr__4) ! (min) Minute of right ascension (J2000) real*4 RAs_1 (nr__4) ! (s) Second of right ascension (J2000) character*1 DE__1 (nr__4) ! Sign of declination (J2000) integer*4 DEd_1 (nr__4) ! (deg) Degree of declination (J2000) integer*4 DEm_1 (nr__4) ! (arcmin) Arcminute of declination (J2000) real*4 DEs_1 (nr__4) ! (arcsec) Arcsecond of declination (J2000) integer*4 Flux (nr__4) ! (mJy) ? Flux real*4 S_N (nr__4) ! ? Signal to noise ratio *Note (1): Subset of 4.5um excess sources that were observed with MIPS at 70um. *Note (2): Source No. given is for table of methanol class II maser sources * (table 8). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table7.dat' ! Table of Methanol class I maser sources integer*4 nr__5 parameter (nr__5=21) ! Number of records character*55 ar__5 ! Full-size record C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg_2 (nr__5) ! (deg) Right Ascension J2000 real*8 DEdeg_2 (nr__5) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) integer*4 ID_2 (nr__5) ! Source identification real*8 GLON_5 (nr__5) ! (deg) Galactic longitude real*4 GLAT_5 (nr__5) ! (deg) Galactic latitude integer*4 RAh_2 (nr__5) ! (h) Hour of right ascension (J2000) integer*4 RAm_2 (nr__5) ! (min) Minute of right ascension (J2000) real*4 RAs_2 (nr__5) ! (s) Second of right ascension (J2000) character*1 DE__2 (nr__5) ! Sign of declination (J2000) integer*4 DEd_2 (nr__5) ! (deg) Degree of declination (J2000) integer*4 DEm_2 (nr__5) ! (arcmin) Arcminute of declination (J2000) real*4 DEs_2 (nr__5) ! (arcsec) Arcsecond of declination (J2000) real*4 Vlsr (nr__5) ! (km/s) LSR velocity character*7 Ref (nr__5) ! Reference(s) (1) *Note (1): References as follows: * H90 = Haschick et al. 1990ApJ...354..556H (44GHz); * S94 = Slysh et al. 1994MNRAS.268..464S (44GHz); * M97 = Mehringer & Menten 1997ApJ...474..346M (44GHz, multiple sources in * Sgr B2); * V00 = Val'tts et al. 2000MNRAS.317..315V (95GHz). * C96 = Caswell 1996MNRAS.283..606C; * W98 = Walsh et al. 1998, Cat. J/MNRAS/301/640 c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table8.dat' ! Table of Methanol class II maser sources integer*4 nr__6 parameter (nr__6=23) ! Number of records character*55 ar__6 ! Full-size record C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg_3 (nr__6) ! (deg) Right Ascension J2000 real*8 DEdeg_3 (nr__6) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) integer*4 ID_3 (nr__6) ! Source identification real*8 GLON_6 (nr__6) ! (deg) Galactic longitude real*4 GLAT_6 (nr__6) ! (deg) Galactic latitude integer*4 RAh_3 (nr__6) ! (h) Hour of right ascension (J2000) integer*4 RAm_3 (nr__6) ! (min) Minute of right ascension (J2000) real*4 RAs_3 (nr__6) ! (s) Second of right ascension (J2000) character*1 DE__3 (nr__6) ! Sign of declination (J2000) integer*4 DEd_3 (nr__6) ! (deg) Degree of declination (J2000) integer*4 DEm_3 (nr__6) ! (arcmin) Arcminute of declination (J2000) real*4 DEs_3 (nr__6) ! (arcsec) Arcsecond of declination (J2000) real*4 Vlsr_1 (nr__6) ! (km/s) LSR velocity character*7 Ref_1 (nr__6) ! Reference(s) (1) *Note (1): References as follows: * H90 = Haschick et al. 1990ApJ...354..556H (44GHz); * S94 = Slysh et al. 1994MNRAS.268..464S (44GHz); * M97 = Mehringer & Menten 1997ApJ...474..346M (44GHz, multiple sources in * Sgr B2); * V00 = Val'tts et al. 2000MNRAS.317..315V (95GHz). * C96 = Caswell 1996MNRAS.283..606C; * W98 = Walsh et al. 1998, Cat. J/MNRAS/301/640 c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table9.dat' ! Parameters of SED fits to 4.5um excess sources near the Galactic Center integer*4 nr__7 parameter (nr__7=14) ! Number of records character*90 ar__7 ! Full-size record C Position composed of: GLON GLAT character*4 v_YHA2009__1(nr__7) ! Green number (gNNA) integer*4 SSTGC_2 (nr__7) ! SSTGC number real*8 GLON_7 (nr__7) ! (deg) Galactic longitude real*8 GLAT_7 (nr__7) ! (deg) Galactic latitude real*4 Chi2_2 (nr__7) ! {chi}^2^ value integer*4 Nfits_2 (nr__7) ! Number of acceptably fitting models real*4 Av_2 (nr__7) ! (mag) Visual extinction A_v_ real*4 Mass_2 (nr__7) ! (Msun) Average mass of the protostar real*4 e_Mass_2 (nr__7) ! (Msun) Mass uncertainty real*4 Lum_2 (nr__7) ! (Lsun) Average luminosity of the protostar real*4 e_Lum_2 (nr__7) ! (Lsun) Lum uncertainty real*4 Menv_2 (nr__7) ! (Msun/yr) Average mass accretion rate real*4 e_Menv_2 (nr__7) ! (Msun/yr) Menv 1{sigma} standard deviation character*3 Stage_2 (nr__7) ! Evolutionary stage of the fitted source c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table10.dat' ! parameters of SED fits to 4.5um excess sources foreground to the Galactic Center integer*4 nr__8 parameter (nr__8=19) ! Number of records character*90 ar__8 ! Full-size record C Position composed of: GLON GLAT character*4 v_YHA2009__2(nr__8) ! Green number (gNNA) integer*4 SSTGC_3 (nr__8) ! SSTGC number real*8 GLON_8 (nr__8) ! (deg) Galactic longitude real*8 GLAT_8 (nr__8) ! (deg) Galactic latitude real*4 Chi2_3 (nr__8) ! {chi}^2^ value integer*4 Nfits_3 (nr__8) ! Number of acceptably fitting models real*4 Av_3 (nr__8) ! (mag) Visual extinction A_v_ real*4 Mass_3 (nr__8) ! (Msun) Average mass of the protostar real*4 e_Mass_3 (nr__8) ! (Msun) Mass uncertainty real*4 Lum_3 (nr__8) ! (Lsun) Average luminosity of the protostar real*4 e_Lum_3 (nr__8) ! (Lsun) Lum uncertainty real*4 Menv_3 (nr__8) ! (Msun/yr) Average mass accretion rate real*4 e_Menv_3 (nr__8) ! (Msun/yr) Menv 1{sigma} standard deviation character*3 Stage_3 (nr__8) ! Evolutionary stage of the fitted source C============================================================================= C Loading file 'table1.dat' ! Galactic coordinates and flux density of * sources in G359.43+0.02 C Format for file interpretation 1 format( + I2,A1,1X,A1,1X,F10.6,1X,F9.6,1X,F4.1,1X,F4.1,1X,F4.1,1X,F4.1, + 1X,F4.1,1X,F4.1,1X,I4) C Effective file loading open(unit=1,status='old',file= +'table1.dat') write(6,*) '....Loading file: table1.dat' do i__=1,18 read(1,'(A61)')ar__ read(ar__,1) + v_YHA2009_(i__),m__YHA2009_(i__),n__YHA2009_(i__),GLON(i__), + GLAT(i__),S6cm(i__),e_S6cm(i__),S20cm(i__),e_S20cm(i__), + amaj(i__),amin(i__),Surf(i__) if(ar__(28:31) .EQ. '') S6cm(i__) = rNULL__ if(ar__(33:36) .EQ. '') e_S6cm(i__) = rNULL__ if(ar__(38:41) .EQ. '') S20cm(i__) = rNULL__ if(ar__(43:46) .EQ. '') e_S20cm(i__) = rNULL__ if(ar__(48:51) .EQ. '') amaj(i__) = rNULL__ if(ar__(53:56) .EQ. '') amin(i__) = rNULL__ if(ar__(58:61) .EQ. '') Surf(i__) = iNULL__ c ..............Just test output........... write(6,1) + v_YHA2009_(i__),m__YHA2009_(i__),n__YHA2009_(i__),GLON(i__), + GLAT(i__),S6cm(i__),e_S6cm(i__),S20cm(i__),e_S20cm(i__), + amaj(i__),amin(i__),Surf(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table3.dat' ! parameters of the SED fits to candidate YSOs * in l<0{deg} C Format for file interpretation 2 format( + I6,A3,1X,F8.4,1X,F8.5,1X,F5.2,1X,I4,1X,F4.1,1X,F4.1,1X,F3.1, + 1X,E7.2,1X,E7.2,1X,E7.2,1X,E7.2,1X,A3) C Effective file loading open(unit=1,status='old',file= +'table3.dat') write(6,*) '....Loading file: table3.dat' do i__=1,223 read(1,'(A88)')ar__1 read(ar__1,2) + SSTGC(i__),n_SSTGC(i__),GLON_1(i__),GLAT_1(i__),Chi2(i__), + Nfits(i__),Av(i__),Mass(i__),e_Mass(i__),Lum(i__),e_Lum(i__), + Menv(i__),e_Menv(i__),Stage(i__) c ..............Just test output........... write(6,2) + SSTGC(i__),n_SSTGC(i__),GLON_1(i__),GLAT_1(i__),Chi2(i__), + Nfits(i__),Av(i__),Mass(i__),e_Mass(i__),Lum(i__),e_Lum(i__), + Menv(i__),e_Menv(i__),Stage(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table4.dat' ! parameters of the SED fits to candidate YSOs * in l>0{deg} C Format for file interpretation 3 format( + I6,A3,1X,F8.4,1X,F8.5,1X,F5.2,1X,I4,1X,F4.1,1X,F4.1,1X,F3.1, + 1X,E7.2,1X,E7.2,1X,E7.2,1X,E7.2,1X,A3) C Effective file loading open(unit=1,status='old',file= +'table4.dat') write(6,*) '....Loading file: table4.dat' do i__=1,137 read(1,'(A88)')ar__2 read(ar__2,3) + SSTGC_1(i__),n_SSTGC_1(i__),GLON_2(i__),GLAT_2(i__), + Chi2_1(i__),Nfits_1(i__),Av_1(i__),Mass_1(i__),e_Mass_1(i__), + Lum_1(i__),e_Lum_1(i__),Menv_1(i__),e_Menv_1(i__),Stage_1(i__) c ..............Just test output........... write(6,3) + SSTGC_1(i__),n_SSTGC_1(i__),GLON_2(i__),GLAT_2(i__), + Chi2_1(i__),Nfits_1(i__),Av_1(i__),Mass_1(i__),e_Mass_1(i__), + Lum_1(i__),e_Lum_1(i__),Menv_1(i__),e_Menv_1(i__),Stage_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table5.dat' ! Table of sources with 4.5um Excess Emission C Format for file interpretation 4 format( + I2,1X,F7.3,1X,F6.3,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1, + 1X,A27) C Effective file loading open(unit=1,status='old',file= +'table5.dat') write(6,*) '....Loading file: table5.dat' do i__=1,33 read(1,'(A69)')ar__3 read(ar__3,4) + ID(i__),GLON_3(i__),GLAT_3(i__),RAh(i__),RAm(i__),RAs(i__), + DE_(i__),DEd(i__),DEm(i__),DEs(i__),Comm(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,4) + ID(i__),GLON_3(i__),GLAT_3(i__),RAh(i__),RAm(i__),RAs(i__), + DE_(i__),DEd(i__),DEm(i__),DEs(i__),Comm(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 'table6.dat' ! Table of 70um Sources C Format for file interpretation 5 format( + I2,A1,1X,F9.5,1X,F8.5,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X, + F4.1,1X,I5,1X,F5.2) C Effective file loading open(unit=1,status='old',file= +'table6.dat') write(6,*) '....Loading file: table6.dat' do i__=1,9 read(1,'(A58)')ar__4 read(ar__4,5) + ID_1(i__),n_ID(i__),GLON_4(i__),GLAT_4(i__),RAh_1(i__), + RAm_1(i__),RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__), + DEs_1(i__),Flux(i__),S_N(i__) if(ar__4(48:52) .EQ. '') Flux(i__) = iNULL__ if(ar__4(54:58) .EQ. '') S_N(i__) = rNULL__ 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,5) + ID_1(i__),n_ID(i__),GLON_4(i__),GLAT_4(i__),RAh_1(i__), + RAm_1(i__),RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__), + DEs_1(i__),Flux(i__),S_N(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 'table7.dat' ! Table of Methanol class I maser sources C Format for file interpretation 6 format( + I2,1X,F7.3,1X,F6.3,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1, + 1X,F5.1,1X,A7) C Effective file loading open(unit=1,status='old',file= +'table7.dat') write(6,*) '....Loading file: table7.dat' do i__=1,21 read(1,'(A55)')ar__5 read(ar__5,6) + ID_2(i__),GLON_5(i__),GLAT_5(i__),RAh_2(i__),RAm_2(i__), + RAs_2(i__),DE__2(i__),DEd_2(i__),DEm_2(i__),DEs_2(i__), + Vlsr(i__),Ref(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,6) + ID_2(i__),GLON_5(i__),GLAT_5(i__),RAh_2(i__),RAm_2(i__), + RAs_2(i__),DE__2(i__),DEd_2(i__),DEm_2(i__),DEs_2(i__), + Vlsr(i__),Ref(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 'table8.dat' ! Table of Methanol class II maser sources C Format for file interpretation 7 format( + I2,1X,F7.3,1X,F6.3,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1, + 1X,F5.1,1X,A7) C Effective file loading open(unit=1,status='old',file= +'table8.dat') write(6,*) '....Loading file: table8.dat' do i__=1,23 read(1,'(A55)')ar__6 read(ar__6,7) + ID_3(i__),GLON_6(i__),GLAT_6(i__),RAh_3(i__),RAm_3(i__), + RAs_3(i__),DE__3(i__),DEd_3(i__),DEm_3(i__),DEs_3(i__), + Vlsr_1(i__),Ref_1(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,7) + ID_3(i__),GLON_6(i__),GLAT_6(i__),RAh_3(i__),RAm_3(i__), + RAs_3(i__),DE__3(i__),DEd_3(i__),DEm_3(i__),DEs_3(i__), + Vlsr_1(i__),Ref_1(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 'table9.dat' ! Parameters of SED fits to 4.5um excess sources * near the Galactic Center C Format for file interpretation 8 format( + A4,1X,I7,1X,F8.4,1X,F7.4,1X,F5.2,1X,I4,1X,F4.1,1X,F4.1,1X, + F3.1,1X,E7.2,1X,E7.2,1X,E7.2,1X,E7.2,1X,A3) C Effective file loading open(unit=1,status='old',file= +'table9.dat') write(6,*) '....Loading file: table9.dat' do i__=1,14 read(1,'(A90)')ar__7 read(ar__7,8) + v_YHA2009__1(i__),SSTGC_2(i__),GLON_7(i__),GLAT_7(i__), + Chi2_2(i__),Nfits_2(i__),Av_2(i__),Mass_2(i__),e_Mass_2(i__), + Lum_2(i__),e_Lum_2(i__),Menv_2(i__),e_Menv_2(i__),Stage_2(i__) c ..............Just test output........... write(6,8) + v_YHA2009__1(i__),SSTGC_2(i__),GLON_7(i__),GLAT_7(i__), + Chi2_2(i__),Nfits_2(i__),Av_2(i__),Mass_2(i__),e_Mass_2(i__), + Lum_2(i__),e_Lum_2(i__),Menv_2(i__),e_Menv_2(i__),Stage_2(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table10.dat' ! parameters of SED fits to 4.5um excess sources * foreground to the Galactic Center C Format for file interpretation 9 format( + A4,1X,I7,1X,F8.4,1X,F7.4,1X,F5.2,1X,I4,1X,F4.1,1X,F4.1,1X, + F3.1,1X,E7.2,1X,E7.2,1X,E7.2,1X,E7.2,1X,A3) C Effective file loading open(unit=1,status='old',file= +'table10.dat') write(6,*) '....Loading file: table10.dat' do i__=1,19 read(1,'(A90)')ar__8 read(ar__8,9) + v_YHA2009__2(i__),SSTGC_3(i__),GLON_8(i__),GLAT_8(i__), + Chi2_3(i__),Nfits_3(i__),Av_3(i__),Mass_3(i__),e_Mass_3(i__), + Lum_3(i__),e_Lum_3(i__),Menv_3(i__),e_Menv_3(i__),Stage_3(i__) c ..............Just test output........... write(6,9) + v_YHA2009__2(i__),SSTGC_3(i__),GLON_8(i__),GLAT_8(i__), + Chi2_3(i__),Nfits_3(i__),Av_3(i__),Mass_3(i__),e_Mass_3(i__), + Lum_3(i__),e_Lum_3(i__),Menv_3(i__),e_Menv_3(i__),Stage_3(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end