Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/698/1989 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/698/1989 Spitzer survey of NGC 2451 (Balog+, 2009) *================================================================================ *Spitzer/IRAC-MIPS survey of NGC 2451A and B: debris disks at 50-80 million *years. * Balog Z., Kiss L.L., Vinko J., Rieke G.H., Muzerolle J., Gaspar A., * Young E.T., Gorlova N. * <Astrophys. J., 698, 1989-2013 (2009)> * =2009ApJ...698.1989B C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table3.dat' ! UBVRIJHK photometry of the members of NGC2451A integer*4 nr__ parameter (nr__=70) ! Number of records character*128 ar__ ! Full-size record C J2000 position composed of: RAdeg DEdeg character*1 Ref (nr__) ! [HP] Membership reference (1) real*8 RAdeg (nr__) ! (deg) Right Ascension in degrees (J2000) real*8 DEdeg (nr__) ! (deg) Declination in degrees (J2000) integer*4 e24 (nr__) ! [6,19]? The 24{mu}m excess source from Table 7 real*4 Umag (nr__) ! (mag) ? The Johnson U band magnitude real*4 e_Umag (nr__) ! (mag) ? The 1{sigma} error in Umag real*4 Bmag (nr__) ! (mag) ? The Johnson B band magnitude real*4 e_Bmag (nr__) ! (mag) ? The 1{sigma} error in Bmag real*4 Vmag (nr__) ! (mag) The Johnson V band magnitude real*4 e_Vmag (nr__) ! (mag) ? The 1{sigma} error in Vmag real*4 Rcmag (nr__) ! (mag) ? The Cousins R band magnitude real*4 e_Rcmag (nr__) ! (mag) ? The 1{sigma} error in Rcmag real*4 Icmag (nr__) ! (mag) ? The Cousins I band magnitude real*4 e_Icmag (nr__) ! (mag) ? The 1{sigma} error in Icmag real*4 Jmag (nr__) ! (mag) 2MASS Johnson J band magnitude real*4 e_Jmag (nr__) ! (mag) ? The 1{sigma} error in Jmag real*4 Hmag (nr__) ! (mag) 2MASS Johnson H band magnitude real*4 e_Hmag (nr__) ! (mag) ? The 1{sigma} error in Hmag real*4 Ksmag (nr__) ! (mag) 2MASS Johnson Ks band magnitude real*4 e_Ksmag (nr__) ! (mag) ? The 1{sigma} error in Ksmag *Note (1): Reference as follows: * P = Platais et al. (2001, Cat. J/AJ/122/1486; * <Cl* NGC 2451 PKB NNNN> in Simbad); * H = Hunsch et al. (2003, Cat. J/A+A/402/571; <[HWS2003] XNNN> in Simbad). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table4.dat' ! UBVRIJHK photometry of the members of NGC2451B integer*4 nr__1 parameter (nr__1=131) ! Number of records character*128 ar__1 ! Full-size record C J2000 position composed of: RAdeg DEdeg character*1 Ref_1 (nr__1) ! [HP] Membership reference (1) real*8 RAdeg_1 (nr__1) ! (deg) Right Ascension in degrees (J2000) real*8 DEdeg_1 (nr__1) ! (deg) Declination in degrees (J2000) integer*4 e24_1 (nr__1) ! [6,19]? The 24{mu}m excess source from Table 7 real*4 Umag_1 (nr__1) ! (mag) ? The Johnson U band magnitude real*4 e_Umag_1 (nr__1) ! (mag) ? The 1{sigma} error in Umag real*4 Bmag_1 (nr__1) ! (mag) ? The Johnson B band magnitude real*4 e_Bmag_1 (nr__1) ! (mag) ? The 1{sigma} error in Bmag real*4 Vmag_1 (nr__1) ! (mag) The Johnson V band magnitude real*4 e_Vmag_1 (nr__1) ! (mag) ? The 1{sigma} error in Vmag real*4 Rcmag_1 (nr__1) ! (mag) ? The Cousins R band magnitude real*4 e_Rcmag_1 (nr__1) ! (mag) ? The 1{sigma} error in Rcmag real*4 Icmag_1 (nr__1) ! (mag) ? The Cousins I band magnitude real*4 e_Icmag_1 (nr__1) ! (mag) ? The 1{sigma} error in Icmag real*4 Jmag_1 (nr__1) ! (mag) 2MASS Johnson J band magnitude real*4 e_Jmag_1 (nr__1) ! (mag) ? The 1{sigma} error in Jmag real*4 Hmag_1 (nr__1) ! (mag) 2MASS Johnson H band magnitude real*4 e_Hmag_1 (nr__1) ! (mag) ? The 1{sigma} error in Hmag real*4 Ksmag_1 (nr__1) ! (mag) 2MASS Johnson Ks band magnitude real*4 e_Ksmag_1 (nr__1) ! (mag) ? The 1{sigma} error in Ksmag *Note (1): Reference as follows: * P = Platais et al. (2001, Cat. J/AJ/122/1486; * <Cl* NGC 2451 PKB NNNN> in Simbad); * H = Hunsch et al. (2003, Cat. J/A+A/402/571; <[HWS2003] XNNN> in Simbad). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table5.dat' ! Spitzer 3.6-24um photometry of members of NGC2451A integer*4 nr__2 parameter (nr__2=70) ! Number of records character*89 ar__2 ! Full-size record C J2000 position composed of: RAdeg DEdeg character*1 Ref_2 (nr__2) ! [HP] Membership reference (1) real*8 RAdeg_2 (nr__2) ! (deg) Right Ascension in degrees (J2000) real*8 DEdeg_2 (nr__2) ! (deg) Declination in degrees (J2000) integer*4 e24_2 (nr__2) ! [6,19]? The 24{mu}m excess source from Table 7 real*4 v_3_6_ (nr__2) ! (mag) ? Spitzer/IRAC 3.6{mu}m band magnitude real*4 e__3_6_ (nr__2) ! (mag) ? The 1{sigma} error in 3.6{mu}m mag real*4 v_4_5_ (nr__2) ! (mag) ? Spitzer/IRAC 4.5{mu}m band magnitude real*4 e__4_5_ (nr__2) ! (mag) ? The 1{sigma} error in 4.5{mu}m mag real*4 v_5_8_ (nr__2) ! (mag) ? Spitzer/IRAC 5.8{mu}m band magnitude real*4 e__5_8_ (nr__2) ! (mag) ? The 1{sigma} error in 5.8{mu}m mag real*4 v_8_0_ (nr__2) ! (mag) ? Spitzer/IRAC 8.0{mu}m band magnitude real*4 e__8_0_ (nr__2) ! (mag) ? The 1{sigma} error in 8.0{mu}m mag real*4 v_24_ (nr__2) ! (mag) ? Spitzer/MIPS 24{mu}m band magnitude real*4 e__24_ (nr__2) ! (mag) ? The 1{sigma} error in 24{mu}m mag *Note (1): Reference as follows: * P = Platais et al. (2001, Cat. J/AJ/122/1486; <Cl* NGC 2451 PKB NNNN> in * Simbad); * H = Hunsch et al. (2003, Cat. J/A+A/402/571; <[HWS2003] XNNN> in Simbad). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table6.dat' ! Spitzer 3.6-24um photometry of members of NGC2451B integer*4 nr__3 parameter (nr__3=131) ! Number of records character*89 ar__3 ! Full-size record C J2000 position composed of: RAdeg DEdeg character*1 Ref_3 (nr__3) ! [HP] Membership reference (1) real*8 RAdeg_3 (nr__3) ! (deg) Right Ascension in degrees (J2000) real*8 DEdeg_3 (nr__3) ! (deg) Declination in degrees (J2000) integer*4 e24_3 (nr__3) ! [6,19]? The 24{mu}m excess source from Table 7 real*4 v_3_6__1 (nr__3) ! (mag) ? Spitzer/IRAC 3.6{mu}m band magnitude real*4 e__3_6__1 (nr__3) ! (mag) ? The 1{sigma} error in 3.6{mu}m mag real*4 v_4_5__1 (nr__3) ! (mag) ? Spitzer/IRAC 4.5{mu}m band magnitude real*4 e__4_5__1 (nr__3) ! (mag) ? The 1{sigma} error in 4.5{mu}m mag real*4 v_5_8__1 (nr__3) ! (mag) ? Spitzer/IRAC 5.8{mu}m band magnitude real*4 e__5_8__1 (nr__3) ! (mag) ? The 1{sigma} error in 5.8{mu}m mag real*4 v_8_0__1 (nr__3) ! (mag) ? Spitzer/IRAC 8.0{mu}m band magnitude real*4 e__8_0__1 (nr__3) ! (mag) ? The 1{sigma} error in 8.0{mu}m mag real*4 v_24__1 (nr__3) ! (mag) ? Spitzer/MIPS 24{mu}m band magnitude real*4 e__24__1 (nr__3) ! (mag) ? The 1{sigma} error in 24{mu}m mag *Note (1): Reference as follows: * P = Platais et al. (2001, Cat. J/AJ/122/1486; <Cl* NGC 2451 PKB NNNN> in * Simbad); * H = Hunsch et al. (2003, Cat. J/A+A/402/571; <[HWS2003] XNNN> in Simbad). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table7.dat' ! Measured and predicted photospheric (K-[24])_0_ colors of cluster members with 24um excesses integer*4 nr__4 parameter (nr__4=14) ! Number of records character*46 ar__4 ! Full-size record C J2000 position composed of: RAdeg DEdeg integer*4 e24_4 (nr__4) ! [6-19] Identification number * (<[BKV2009] NN> in Simbad) character*1 Comp (nr__4) ! [AB] NGC 2451 component real*8 RAdeg_4 (nr__4) ! (deg) Right ascension in decimal degrees (J2000) real*8 DEdeg_4 (nr__4) ! (deg) Declination in decimal degrees (J2000) real*4 K__24_ (nr__4) ! (mag) Measured (K-[24])_0_ color index real*4 e_K__24_ (nr__4) ! (mag) K-[24] uncertainty real*4 K__24_p (nr__4) ! (mag) Predicted photospheric K-[24] color index integer*4 Sigma (nr__4) ! {sigma} significance of the excess (1) *Note (1): The significance was calculated using the photometric errors * and the rms around this locus. C============================================================================= C Loading file 'table3.dat' ! UBVRIJHK photometry of the members of NGC2451A C Format for file interpretation 1 format( + A1,1X,F9.5,1X,F9.5,1X,I2,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X, + F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3, + 1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3) C Effective file loading open(unit=1,status='old',file= +'table3.dat') write(6,*) '....Loading file: table3.dat' do i__=1,70 read(1,'(A128)')ar__ read(ar__,1) + Ref(i__),RAdeg(i__),DEdeg(i__),e24(i__),Umag(i__),e_Umag(i__), + Bmag(i__),e_Bmag(i__),Vmag(i__),e_Vmag(i__),Rcmag(i__), + e_Rcmag(i__),Icmag(i__),e_Icmag(i__),Jmag(i__),e_Jmag(i__), + Hmag(i__),e_Hmag(i__),Ksmag(i__),e_Ksmag(i__) if(ar__(23:24) .EQ. '') e24(i__) = iNULL__ if(ar__(26:31) .EQ. '') Umag(i__) = rNULL__ if(ar__(33:37) .EQ. '') e_Umag(i__) = rNULL__ if(ar__(39:44) .EQ. '') Bmag(i__) = rNULL__ if(ar__(46:50) .EQ. '') e_Bmag(i__) = rNULL__ if(ar__(59:63) .EQ. '') e_Vmag(i__) = rNULL__ if(ar__(65:70) .EQ. '') Rcmag(i__) = rNULL__ if(ar__(72:76) .EQ. '') e_Rcmag(i__) = rNULL__ if(ar__(78:83) .EQ. '') Icmag(i__) = rNULL__ if(ar__(85:89) .EQ. '') e_Icmag(i__) = rNULL__ if(ar__(98:102) .EQ. '') e_Jmag(i__) = rNULL__ if(ar__(111:115) .EQ. '') e_Hmag(i__) = rNULL__ if(ar__(124:128) .EQ. '') e_Ksmag(i__) = rNULL__ c ..............Just test output........... write(6,1) + Ref(i__),RAdeg(i__),DEdeg(i__),e24(i__),Umag(i__),e_Umag(i__), + Bmag(i__),e_Bmag(i__),Vmag(i__),e_Vmag(i__),Rcmag(i__), + e_Rcmag(i__),Icmag(i__),e_Icmag(i__),Jmag(i__),e_Jmag(i__), + Hmag(i__),e_Hmag(i__),Ksmag(i__),e_Ksmag(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table4.dat' ! UBVRIJHK photometry of the members of NGC2451B C Format for file interpretation 2 format( + A1,1X,F9.5,1X,F9.5,1X,I2,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X, + F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3, + 1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3) C Effective file loading open(unit=1,status='old',file= +'table4.dat') write(6,*) '....Loading file: table4.dat' do i__=1,131 read(1,'(A128)')ar__1 read(ar__1,2) + Ref_1(i__),RAdeg_1(i__),DEdeg_1(i__),e24_1(i__),Umag_1(i__), + e_Umag_1(i__),Bmag_1(i__),e_Bmag_1(i__),Vmag_1(i__), + e_Vmag_1(i__),Rcmag_1(i__),e_Rcmag_1(i__),Icmag_1(i__), + e_Icmag_1(i__),Jmag_1(i__),e_Jmag_1(i__),Hmag_1(i__), + e_Hmag_1(i__),Ksmag_1(i__),e_Ksmag_1(i__) if(ar__1(23:24) .EQ. '') e24_1(i__) = iNULL__ if(ar__1(26:31) .EQ. '') Umag_1(i__) = rNULL__ if(ar__1(33:37) .EQ. '') e_Umag_1(i__) = rNULL__ if(ar__1(39:44) .EQ. '') Bmag_1(i__) = rNULL__ if(ar__1(46:50) .EQ. '') e_Bmag_1(i__) = rNULL__ if(ar__1(59:63) .EQ. '') e_Vmag_1(i__) = rNULL__ if(ar__1(65:70) .EQ. '') Rcmag_1(i__) = rNULL__ if(ar__1(72:76) .EQ. '') e_Rcmag_1(i__) = rNULL__ if(ar__1(78:83) .EQ. '') Icmag_1(i__) = rNULL__ if(ar__1(85:89) .EQ. '') e_Icmag_1(i__) = rNULL__ if(ar__1(98:102) .EQ. '') e_Jmag_1(i__) = rNULL__ if(ar__1(111:115) .EQ. '') e_Hmag_1(i__) = rNULL__ if(ar__1(124:128) .EQ. '') e_Ksmag_1(i__) = rNULL__ c ..............Just test output........... write(6,2) + Ref_1(i__),RAdeg_1(i__),DEdeg_1(i__),e24_1(i__),Umag_1(i__), + e_Umag_1(i__),Bmag_1(i__),e_Bmag_1(i__),Vmag_1(i__), + e_Vmag_1(i__),Rcmag_1(i__),e_Rcmag_1(i__),Icmag_1(i__), + e_Icmag_1(i__),Jmag_1(i__),e_Jmag_1(i__),Hmag_1(i__), + e_Hmag_1(i__),Ksmag_1(i__),e_Ksmag_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table5.dat' ! Spitzer 3.6-24um photometry of members of NGC2451A C Format for file interpretation 3 format( + A1,1X,F9.5,1X,F9.5,1X,I2,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X, + F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3) C Effective file loading open(unit=1,status='old',file= +'table5.dat') write(6,*) '....Loading file: table5.dat' do i__=1,70 read(1,'(A89)')ar__2 read(ar__2,3) + Ref_2(i__),RAdeg_2(i__),DEdeg_2(i__),e24_2(i__),v_3_6_(i__), + e__3_6_(i__),v_4_5_(i__),e__4_5_(i__),v_5_8_(i__), + e__5_8_(i__),v_8_0_(i__),e__8_0_(i__),v_24_(i__),e__24_(i__) if(ar__2(23:24) .EQ. '') e24_2(i__) = iNULL__ if(ar__2(26:31) .EQ. '') v_3_6_(i__) = rNULL__ if(ar__2(33:37) .EQ. '') e__3_6_(i__) = rNULL__ if(ar__2(39:44) .EQ. '') v_4_5_(i__) = rNULL__ if(ar__2(46:50) .EQ. '') e__4_5_(i__) = rNULL__ if(ar__2(52:57) .EQ. '') v_5_8_(i__) = rNULL__ if(ar__2(59:63) .EQ. '') e__5_8_(i__) = rNULL__ if(ar__2(65:70) .EQ. '') v_8_0_(i__) = rNULL__ if(ar__2(72:76) .EQ. '') e__8_0_(i__) = rNULL__ if(ar__2(78:83) .EQ. '') v_24_(i__) = rNULL__ if(ar__2(85:89) .EQ. '') e__24_(i__) = rNULL__ c ..............Just test output........... write(6,3) + Ref_2(i__),RAdeg_2(i__),DEdeg_2(i__),e24_2(i__),v_3_6_(i__), + e__3_6_(i__),v_4_5_(i__),e__4_5_(i__),v_5_8_(i__), + e__5_8_(i__),v_8_0_(i__),e__8_0_(i__),v_24_(i__),e__24_(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table6.dat' ! Spitzer 3.6-24um photometry of members of NGC2451B C Format for file interpretation 4 format( + A1,1X,F9.5,1X,F9.5,1X,I2,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X, + F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.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,131 read(1,'(A89)')ar__3 read(ar__3,4) + Ref_3(i__),RAdeg_3(i__),DEdeg_3(i__),e24_3(i__),v_3_6__1(i__), + e__3_6__1(i__),v_4_5__1(i__),e__4_5__1(i__),v_5_8__1(i__), + e__5_8__1(i__),v_8_0__1(i__),e__8_0__1(i__),v_24__1(i__), + e__24__1(i__) if(ar__3(23:24) .EQ. '') e24_3(i__) = iNULL__ if(ar__3(26:31) .EQ. '') v_3_6__1(i__) = rNULL__ if(ar__3(33:37) .EQ. '') e__3_6__1(i__) = rNULL__ if(ar__3(39:44) .EQ. '') v_4_5__1(i__) = rNULL__ if(ar__3(46:50) .EQ. '') e__4_5__1(i__) = rNULL__ if(ar__3(52:57) .EQ. '') v_5_8__1(i__) = rNULL__ if(ar__3(59:63) .EQ. '') e__5_8__1(i__) = rNULL__ if(ar__3(65:70) .EQ. '') v_8_0__1(i__) = rNULL__ if(ar__3(72:76) .EQ. '') e__8_0__1(i__) = rNULL__ if(ar__3(78:83) .EQ. '') v_24__1(i__) = rNULL__ if(ar__3(85:89) .EQ. '') e__24__1(i__) = rNULL__ c ..............Just test output........... write(6,4) + Ref_3(i__),RAdeg_3(i__),DEdeg_3(i__),e24_3(i__),v_3_6__1(i__), + e__3_6__1(i__),v_4_5__1(i__),e__4_5__1(i__),v_5_8__1(i__), + e__5_8__1(i__),v_8_0__1(i__),e__8_0__1(i__),v_24__1(i__), + e__24__1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table7.dat' ! Measured and predicted photospheric (K-[24])_0_ * colors of cluster members with 24um excesses C Format for file interpretation 5 format(I2,1X,A1,1X,F9.5,1X,F9.5,1X,F5.3,1X,F5.3,1X,F6.3,1X,I2) C Effective file loading open(unit=1,status='old',file= +'table7.dat') write(6,*) '....Loading file: table7.dat' do i__=1,14 read(1,'(A46)')ar__4 read(ar__4,5) + e24_4(i__),Comp(i__),RAdeg_4(i__),DEdeg_4(i__),K__24_(i__), + e_K__24_(i__),K__24_p(i__),Sigma(i__) c ..............Just test output........... write(6,5) + e24_4(i__),Comp(i__),RAdeg_4(i__),DEdeg_4(i__),K__24_(i__), + e_K__24_(i__),K__24_p(i__),Sigma(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end