Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/675/491 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 2026-Jun-13
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/675/491 Lynds 988e JHKs and IRAC photometry (Allen+, 2008)
*================================================================================
*Spitzer, near-infrared, and submillimeter imaging of the relatively sparse
*young cluster, Lynds 988e.
* Allen T.S., Pipher J.L., Gutermuth R.A., Megeath S.T., Adams J.D.,
* Herter T.L., Williams J.P., Goetz-Bixby J.A., Allen L.E., Myers P.C.
* <Astrophys. J., 675, 491-506 (2008)>
* =2008ApJ...675..491A
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! *Objects with Mid-Infrared Excesses: Class I
integer*4 nr__
parameter (nr__=19) ! Number of records
character*153 ar__ ! Full-size record
C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg (nr__) ! (deg) Right Ascension J2000.0
real*8 DEdeg (nr__) ! (deg) Declination J2000.0
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 RAh (nr__) ! (h) Right ascension (J2000.0)
integer*4 RAm (nr__) ! (min) Right ascension (J2000.0)
real*4 RAs (nr__) ! (s) Right ascension (J2000.0)
real*4 e_RAs (nr__) ! (s) rms uncertainty on RA
character*1 DE_ (nr__) ! Declination sign (J2000.0)
integer*4 DEd (nr__) ! (deg) Declination (J2000.0)
integer*4 DEm (nr__) ! (arcmin) Declination (J2000.0)
real*4 DEs (nr__) ! (arcsec) Declination (J2000.0)
real*4 e_DEs (nr__) ! (arcsec) rms uncertainty on DE
real*4 Jmag (nr__) ! (mag) ? Dereddened J magnitude (1)
real*4 e_Jmag (nr__) ! (mag) ? rms uncertainty on Jmag
real*4 Hmag (nr__) ! (mag) ? Dereddened H magnitude (1)
real*4 e_Hmag (nr__) ! (mag) ? rms uncertainty on Hmag
real*4 Ksmag (nr__) ! (mag) ? Dereddened Ks magnitude (1)
real*4 e_Ksmag (nr__) ! (mag) ? rms uncertainty on Ksmag
real*4 v_3_6_ (nr__) ! (mag) ? Dereddened 3.6 magnitude (1)
real*4 e__3_6_ (nr__) ! (mag) ? rms uncertainty on [3.6]
real*4 v_4_5_ (nr__) ! (mag) ? Dereddened 4.5um magnitude (1)
real*4 e__4_5_ (nr__) ! (mag) ? rms uncertainty on [4.5]
real*4 v_5_8_ (nr__) ! (mag) ? Dereddened 5.8um magnitude (1)
real*4 e__5_8_ (nr__) ! (mag) ? rms uncertainty on [5.8]
real*4 v_8_0_ (nr__) ! (mag) ? Dereddened 8.0um magnitude (1)
real*4 e__8_0_ (nr__) ! (mag) ? rms uncertainty on [8.0]
real*4 v_24_ (nr__) ! (mag) ? Dereddened 24um magnitude (1)
real*4 e__24_ (nr__) ! (mag) ? rms uncertainty on [24]
real*4 AKs (nr__) ! (mag) ? Absoption in Ks band (2)
real*4 alpha (nr__) ! ? Slope {alpha} determined from IRAC photometry
character*1 n_alpha (nr__) ! [gf] how {alpha} is computed (3)
integer*4 IfAHA (nr__) ! ? Herbig & Dahm (2006, Cat. J/AJ/131/1530) ID
real*4 EWHa (nr__) ! (0.1nm) ? Halpha equivalent width (4)
*Note (1): Magnitudes are dereddened when extinction value AKs quoted.
*Note (2): AKs determined from method outlined in Sect.5.
*Note (3): Notes as follows:
* f = Slope alpha determined from [4.5] to [24] photometry.
* g = Slope alpha determined from [3.6], [5.8] and [24] photometry.
*Note (4): from Herbig & Dahm (2006, Cat. J/AJ/131/1530).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! *Objects with Mid-Infrared Excesses: Class II
integer*4 nr__1
parameter (nr__1=74) ! Number of records
character*153 ar__1 ! Full-size record
C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_1 (nr__1) ! (deg) Right Ascension J2000.0
real*8 DEdeg_1 (nr__1) ! (deg) Declination J2000.0
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 RAh_1 (nr__1) ! (h) Right ascension (J2000.0)
integer*4 RAm_1 (nr__1) ! (min) Right ascension (J2000.0)
real*4 RAs_1 (nr__1) ! (s) Right ascension (J2000.0)
real*4 e_RAs_1 (nr__1) ! (s) rms uncertainty on RA
character*1 DE__1 (nr__1) ! Declination sign (J2000.0)
integer*4 DEd_1 (nr__1) ! (deg) Declination (J2000.0)
integer*4 DEm_1 (nr__1) ! (arcmin) Declination (J2000.0)
real*4 DEs_1 (nr__1) ! (arcsec) Declination (J2000.0)
real*4 e_DEs_1 (nr__1) ! (arcsec) rms uncertainty on DE
real*4 Jmag_1 (nr__1) ! (mag) ? Dereddened J magnitude (1)
real*4 e_Jmag_1 (nr__1) ! (mag) ? rms uncertainty on Jmag
real*4 Hmag_1 (nr__1) ! (mag) ? Dereddened H magnitude (1)
real*4 e_Hmag_1 (nr__1) ! (mag) ? rms uncertainty on Hmag
real*4 Ksmag_1 (nr__1) ! (mag) ? Dereddened Ks magnitude (1)
real*4 e_Ksmag_1 (nr__1) ! (mag) ? rms uncertainty on Ksmag
real*4 v_3_6__1 (nr__1) ! (mag) ? Dereddened 3.6 magnitude (1)
real*4 e__3_6__1 (nr__1) ! (mag) ? rms uncertainty on [3.6]
real*4 v_4_5__1 (nr__1) ! (mag) ? Dereddened 4.5um magnitude (1)
real*4 e__4_5__1 (nr__1) ! (mag) ? rms uncertainty on [4.5]
real*4 v_5_8__1 (nr__1) ! (mag) ? Dereddened 5.8um magnitude (1)
real*4 e__5_8__1 (nr__1) ! (mag) ? rms uncertainty on [5.8]
real*4 v_8_0__1 (nr__1) ! (mag) ? Dereddened 8.0um magnitude (1)
real*4 e__8_0__1 (nr__1) ! (mag) ? rms uncertainty on [8.0]
real*4 v_24__1 (nr__1) ! (mag) ? Dereddened 24um magnitude (1)
real*4 e__24__1 (nr__1) ! (mag) ? rms uncertainty on [24]
real*4 AKs_1 (nr__1) ! (mag) ? Absoption in Ks band (2)
real*4 alpha_1 (nr__1) ! ? Slope {alpha} determined from IRAC photometry
character*1 n_alpha_1 (nr__1) ! [gf] how {alpha} is computed (3)
integer*4 IfAHA_1 (nr__1) ! ? Herbig & Dahm (2006, Cat. J/AJ/131/1530) ID
real*4 EWHa_1 (nr__1) ! (0.1nm) ? Halpha equivalent width (4)
*Note (1): Magnitudes are dereddened when extinction value AKs quoted.
*Note (2): AKs determined from method outlined in Sect.5.
*Note (3): Notes as follows:
* f = Slope alpha determined from [4.5] to [24] photometry.
* g = Slope alpha determined from [3.6], [5.8] and [24] photometry.
*Note (4): from Herbig & Dahm (2006, Cat. J/AJ/131/1530).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! *Objects with Mid-Infrared Excesses: Transition disks
integer*4 nr__2
parameter (nr__2=7) ! Number of records
character*153 ar__2 ! Full-size record
C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_2 (nr__2) ! (deg) Right Ascension J2000.0
real*8 DEdeg_2 (nr__2) ! (deg) Declination J2000.0
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 RAh_2 (nr__2) ! (h) Right ascension (J2000.0)
integer*4 RAm_2 (nr__2) ! (min) Right ascension (J2000.0)
real*4 RAs_2 (nr__2) ! (s) Right ascension (J2000.0)
real*4 e_RAs_2 (nr__2) ! (s) rms uncertainty on RA
character*1 DE__2 (nr__2) ! Declination sign (J2000.0)
integer*4 DEd_2 (nr__2) ! (deg) Declination (J2000.0)
integer*4 DEm_2 (nr__2) ! (arcmin) Declination (J2000.0)
real*4 DEs_2 (nr__2) ! (arcsec) Declination (J2000.0)
real*4 e_DEs_2 (nr__2) ! (arcsec) rms uncertainty on DE
real*4 Jmag_2 (nr__2) ! (mag) ? Dereddened J magnitude (1)
real*4 e_Jmag_2 (nr__2) ! (mag) ? rms uncertainty on Jmag
real*4 Hmag_2 (nr__2) ! (mag) ? Dereddened H magnitude (1)
real*4 e_Hmag_2 (nr__2) ! (mag) ? rms uncertainty on Hmag
real*4 Ksmag_2 (nr__2) ! (mag) ? Dereddened Ks magnitude (1)
real*4 e_Ksmag_2 (nr__2) ! (mag) ? rms uncertainty on Ksmag
real*4 v_3_6__2 (nr__2) ! (mag) ? Dereddened 3.6 magnitude (1)
real*4 e__3_6__2 (nr__2) ! (mag) ? rms uncertainty on [3.6]
real*4 v_4_5__2 (nr__2) ! (mag) ? Dereddened 4.5um magnitude (1)
real*4 e__4_5__2 (nr__2) ! (mag) ? rms uncertainty on [4.5]
real*4 v_5_8__2 (nr__2) ! (mag) ? Dereddened 5.8um magnitude (1)
real*4 e__5_8__2 (nr__2) ! (mag) ? rms uncertainty on [5.8]
real*4 v_8_0__2 (nr__2) ! (mag) ? Dereddened 8.0um magnitude (1)
real*4 e__8_0__2 (nr__2) ! (mag) ? rms uncertainty on [8.0]
real*4 v_24__2 (nr__2) ! (mag) ? Dereddened 24um magnitude (1)
real*4 e__24__2 (nr__2) ! (mag) ? rms uncertainty on [24]
real*4 AKs_2 (nr__2) ! (mag) ? Absoption in Ks band (2)
real*4 alpha_2 (nr__2) ! ? Slope {alpha} determined from IRAC photometry
character*1 n_alpha_2 (nr__2) ! [gf] how {alpha} is computed (3)
integer*4 IfAHA_2 (nr__2) ! ? Herbig & Dahm (2006, Cat. J/AJ/131/1530) ID
real*4 EWHa_2 (nr__2) ! (0.1nm) ? Halpha equivalent width (4)
*Note (1): Magnitudes are dereddened when extinction value AKs quoted.
*Note (2): AKs determined from method outlined in Sect.5.
*Note (3): Notes as follows:
* f = Slope alpha determined from [4.5] to [24] photometry.
* g = Slope alpha determined from [3.6], [5.8] and [24] photometry.
*Note (4): from Herbig & Dahm (2006, Cat. J/AJ/131/1530).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! *Objects without Mid-Infrared Excesses: Class III
integer*4 nr__3
parameter (nr__3=20) ! Number of records
character*153 ar__3 ! Full-size record
C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_3 (nr__3) ! (deg) Right Ascension J2000.0
real*8 DEdeg_3 (nr__3) ! (deg) Declination J2000.0
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 RAh_3 (nr__3) ! (h) Right ascension (J2000.0)
integer*4 RAm_3 (nr__3) ! (min) Right ascension (J2000.0)
real*4 RAs_3 (nr__3) ! (s) Right ascension (J2000.0)
real*4 e_RAs_3 (nr__3) ! (s) rms uncertainty on RA
character*1 DE__3 (nr__3) ! Declination sign (J2000.0)
integer*4 DEd_3 (nr__3) ! (deg) Declination (J2000.0)
integer*4 DEm_3 (nr__3) ! (arcmin) Declination (J2000.0)
real*4 DEs_3 (nr__3) ! (arcsec) Declination (J2000.0)
real*4 e_DEs_3 (nr__3) ! (arcsec) rms uncertainty on DE
real*4 Jmag_3 (nr__3) ! (mag) ? Dereddened J magnitude (1)
real*4 e_Jmag_3 (nr__3) ! (mag) ? rms uncertainty on Jmag
real*4 Hmag_3 (nr__3) ! (mag) ? Dereddened H magnitude (1)
real*4 e_Hmag_3 (nr__3) ! (mag) ? rms uncertainty on Hmag
real*4 Ksmag_3 (nr__3) ! (mag) ? Dereddened Ks magnitude (1)
real*4 e_Ksmag_3 (nr__3) ! (mag) ? rms uncertainty on Ksmag
real*4 v_3_6__3 (nr__3) ! (mag) ? Dereddened 3.6 magnitude (1)
real*4 e__3_6__3 (nr__3) ! (mag) ? rms uncertainty on [3.6]
real*4 v_4_5__3 (nr__3) ! (mag) ? Dereddened 4.5um magnitude (1)
real*4 e__4_5__3 (nr__3) ! (mag) ? rms uncertainty on [4.5]
real*4 v_5_8__3 (nr__3) ! (mag) ? Dereddened 5.8um magnitude (1)
real*4 e__5_8__3 (nr__3) ! (mag) ? rms uncertainty on [5.8]
real*4 v_8_0__3 (nr__3) ! (mag) ? Dereddened 8.0um magnitude (1)
real*4 e__8_0__3 (nr__3) ! (mag) ? rms uncertainty on [8.0]
real*4 v_24__3 (nr__3) ! (mag) ? Dereddened 24um magnitude (1)
real*4 e__24__3 (nr__3) ! (mag) ? rms uncertainty on [24]
real*4 AKs_3 (nr__3) ! (mag) ? Absoption in Ks band (2)
real*4 alpha_3 (nr__3) ! ? Slope {alpha} determined from IRAC photometry
character*1 n_alpha_3 (nr__3) ! [gf] how {alpha} is computed (3)
integer*4 IfAHA_3 (nr__3) ! ? Herbig & Dahm (2006, Cat. J/AJ/131/1530) ID
real*4 EWHa_3 (nr__3) ! (0.1nm) ? Halpha equivalent width (4)
*Note (1): Magnitudes are dereddened when extinction value AKs quoted.
*Note (2): AKs determined from method outlined in Sect.5.
*Note (3): Notes as follows:
* f = Slope alpha determined from [4.5] to [24] photometry.
* g = Slope alpha determined from [3.6], [5.8] and [24] photometry.
*Note (4): from Herbig & Dahm (2006, Cat. J/AJ/131/1530).
C=============================================================================
C Loading file 'table2.dat' ! *Objects with Mid-Infrared Excesses: Class I
C Format for file interpretation
1 format(
+ I2,1X,I2,1X,F5.2,1X,F4.2,1X,A1,I2,1X,I2,1X,F5.2,1X,F4.2,1X,
+ F5.2,1X,F6.3,1X,F5.2,1X,F5.3,1X,F5.2,1X,F5.3,1X,F5.2,1X,F5.3,
+ 1X,F5.2,1X,F5.3,1X,F5.2,1X,F5.3,1X,F5.2,1X,F5.3,1X,F4.2,1X,
+ F5.3,1X,F4.2,1X,F6.3,A1,1X,I3,1X,F5.1)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,19
read(1,'(A153)')ar__
read(ar__,1)
+ RAh(i__),RAm(i__),RAs(i__),e_RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),e_DEs(i__),Jmag(i__),e_Jmag(i__),Hmag(i__),
+ e_Hmag(i__),Ksmag(i__),e_Ksmag(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__),AKs(i__),alpha(i__),
+ n_alpha(i__),IfAHA(i__),EWHa(i__)
if(ar__(36:40) .EQ. '') Jmag(i__) = rNULL__
if(ar__(42:47) .EQ. '') e_Jmag(i__) = rNULL__
if(ar__(49:53) .EQ. '') Hmag(i__) = rNULL__
if(ar__(55:59) .EQ. '') e_Hmag(i__) = rNULL__
if(ar__(61:65) .EQ. '') Ksmag(i__) = rNULL__
if(ar__(67:71) .EQ. '') e_Ksmag(i__) = rNULL__
if(ar__(73:77) .EQ. '') v_3_6_(i__) = rNULL__
if(ar__(79:83) .EQ. '') e__3_6_(i__) = rNULL__
if(ar__(85:89) .EQ. '') v_4_5_(i__) = rNULL__
if(ar__(91:95) .EQ. '') e__4_5_(i__) = rNULL__
if(ar__(97:101) .EQ. '') v_5_8_(i__) = rNULL__
if(ar__(103:107) .EQ. '') e__5_8_(i__) = rNULL__
if(ar__(109:113) .EQ. '') v_8_0_(i__) = rNULL__
if(ar__(115:119) .EQ. '') e__8_0_(i__) = rNULL__
if(ar__(121:124) .EQ. '') v_24_(i__) = rNULL__
if(ar__(126:130) .EQ. '') e__24_(i__) = rNULL__
if(ar__(132:135) .EQ. '') AKs(i__) = rNULL__
if(ar__(137:142) .EQ. '') alpha(i__) = rNULL__
if(ar__(145:147) .EQ. '') IfAHA(i__) = iNULL__
if(ar__(149:153) .EQ. '') EWHa(i__) = rNULL__
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)
+ RAh(i__),RAm(i__),RAs(i__),e_RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),e_DEs(i__),Jmag(i__),e_Jmag(i__),Hmag(i__),
+ e_Hmag(i__),Ksmag(i__),e_Ksmag(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__),AKs(i__),alpha(i__),
+ n_alpha(i__),IfAHA(i__),EWHa(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' ! *Objects with Mid-Infrared Excesses: Class II
C Format for file interpretation
2 format(
+ I2,1X,I2,1X,F5.2,1X,F4.2,1X,A1,I2,1X,I2,1X,F5.2,1X,F4.2,1X,
+ F5.2,1X,F6.3,1X,F5.2,1X,F5.3,1X,F5.2,1X,F5.3,1X,F5.2,1X,F5.3,
+ 1X,F5.2,1X,F5.3,1X,F5.2,1X,F5.3,1X,F5.2,1X,F5.3,1X,F4.2,1X,
+ F5.3,1X,F4.2,1X,F6.3,A1,1X,I3,1X,F5.1)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,74
read(1,'(A153)')ar__1
read(ar__1,2)
+ RAh_1(i__),RAm_1(i__),RAs_1(i__),e_RAs_1(i__),DE__1(i__),
+ DEd_1(i__),DEm_1(i__),DEs_1(i__),e_DEs_1(i__),Jmag_1(i__),
+ e_Jmag_1(i__),Hmag_1(i__),e_Hmag_1(i__),Ksmag_1(i__),
+ e_Ksmag_1(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__),AKs_1(i__),
+ alpha_1(i__),n_alpha_1(i__),IfAHA_1(i__),EWHa_1(i__)
if(ar__1(36:40) .EQ. '') Jmag_1(i__) = rNULL__
if(ar__1(42:47) .EQ. '') e_Jmag_1(i__) = rNULL__
if(ar__1(49:53) .EQ. '') Hmag_1(i__) = rNULL__
if(ar__1(55:59) .EQ. '') e_Hmag_1(i__) = rNULL__
if(ar__1(61:65) .EQ. '') Ksmag_1(i__) = rNULL__
if(ar__1(67:71) .EQ. '') e_Ksmag_1(i__) = rNULL__
if(ar__1(73:77) .EQ. '') v_3_6__1(i__) = rNULL__
if(ar__1(79:83) .EQ. '') e__3_6__1(i__) = rNULL__
if(ar__1(85:89) .EQ. '') v_4_5__1(i__) = rNULL__
if(ar__1(91:95) .EQ. '') e__4_5__1(i__) = rNULL__
if(ar__1(97:101) .EQ. '') v_5_8__1(i__) = rNULL__
if(ar__1(103:107) .EQ. '') e__5_8__1(i__) = rNULL__
if(ar__1(109:113) .EQ. '') v_8_0__1(i__) = rNULL__
if(ar__1(115:119) .EQ. '') e__8_0__1(i__) = rNULL__
if(ar__1(121:124) .EQ. '') v_24__1(i__) = rNULL__
if(ar__1(126:130) .EQ. '') e__24__1(i__) = rNULL__
if(ar__1(132:135) .EQ. '') AKs_1(i__) = rNULL__
if(ar__1(137:142) .EQ. '') alpha_1(i__) = rNULL__
if(ar__1(145:147) .EQ. '') IfAHA_1(i__) = iNULL__
if(ar__1(149:153) .EQ. '') EWHa_1(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,2)
+ RAh_1(i__),RAm_1(i__),RAs_1(i__),e_RAs_1(i__),DE__1(i__),
+ DEd_1(i__),DEm_1(i__),DEs_1(i__),e_DEs_1(i__),Jmag_1(i__),
+ e_Jmag_1(i__),Hmag_1(i__),e_Hmag_1(i__),Ksmag_1(i__),
+ e_Ksmag_1(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__),AKs_1(i__),
+ alpha_1(i__),n_alpha_1(i__),IfAHA_1(i__),EWHa_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' ! *Objects with Mid-Infrared Excesses: Transition disks
C Format for file interpretation
3 format(
+ I2,1X,I2,1X,F5.2,1X,F4.2,1X,A1,I2,1X,I2,1X,F5.2,1X,F4.2,1X,
+ F5.2,1X,F6.3,1X,F5.2,1X,F5.3,1X,F5.2,1X,F5.3,1X,F5.2,1X,F5.3,
+ 1X,F5.2,1X,F5.3,1X,F5.2,1X,F5.3,1X,F5.2,1X,F5.3,1X,F4.2,1X,
+ F5.3,1X,F4.2,1X,F6.3,A1,1X,I3,1X,F5.1)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,7
read(1,'(A153)')ar__2
read(ar__2,3)
+ RAh_2(i__),RAm_2(i__),RAs_2(i__),e_RAs_2(i__),DE__2(i__),
+ DEd_2(i__),DEm_2(i__),DEs_2(i__),e_DEs_2(i__),Jmag_2(i__),
+ e_Jmag_2(i__),Hmag_2(i__),e_Hmag_2(i__),Ksmag_2(i__),
+ e_Ksmag_2(i__),v_3_6__2(i__),e__3_6__2(i__),v_4_5__2(i__),
+ e__4_5__2(i__),v_5_8__2(i__),e__5_8__2(i__),v_8_0__2(i__),
+ e__8_0__2(i__),v_24__2(i__),e__24__2(i__),AKs_2(i__),
+ alpha_2(i__),n_alpha_2(i__),IfAHA_2(i__),EWHa_2(i__)
if(ar__2(36:40) .EQ. '') Jmag_2(i__) = rNULL__
if(ar__2(42:47) .EQ. '') e_Jmag_2(i__) = rNULL__
if(ar__2(49:53) .EQ. '') Hmag_2(i__) = rNULL__
if(ar__2(55:59) .EQ. '') e_Hmag_2(i__) = rNULL__
if(ar__2(61:65) .EQ. '') Ksmag_2(i__) = rNULL__
if(ar__2(67:71) .EQ. '') e_Ksmag_2(i__) = rNULL__
if(ar__2(73:77) .EQ. '') v_3_6__2(i__) = rNULL__
if(ar__2(79:83) .EQ. '') e__3_6__2(i__) = rNULL__
if(ar__2(85:89) .EQ. '') v_4_5__2(i__) = rNULL__
if(ar__2(91:95) .EQ. '') e__4_5__2(i__) = rNULL__
if(ar__2(97:101) .EQ. '') v_5_8__2(i__) = rNULL__
if(ar__2(103:107) .EQ. '') e__5_8__2(i__) = rNULL__
if(ar__2(109:113) .EQ. '') v_8_0__2(i__) = rNULL__
if(ar__2(115:119) .EQ. '') e__8_0__2(i__) = rNULL__
if(ar__2(121:124) .EQ. '') v_24__2(i__) = rNULL__
if(ar__2(126:130) .EQ. '') e__24__2(i__) = rNULL__
if(ar__2(132:135) .EQ. '') AKs_2(i__) = rNULL__
if(ar__2(137:142) .EQ. '') alpha_2(i__) = rNULL__
if(ar__2(145:147) .EQ. '') IfAHA_2(i__) = iNULL__
if(ar__2(149:153) .EQ. '') EWHa_2(i__) = rNULL__
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)
+ RAh_2(i__),RAm_2(i__),RAs_2(i__),e_RAs_2(i__),DE__2(i__),
+ DEd_2(i__),DEm_2(i__),DEs_2(i__),e_DEs_2(i__),Jmag_2(i__),
+ e_Jmag_2(i__),Hmag_2(i__),e_Hmag_2(i__),Ksmag_2(i__),
+ e_Ksmag_2(i__),v_3_6__2(i__),e__3_6__2(i__),v_4_5__2(i__),
+ e__4_5__2(i__),v_5_8__2(i__),e__5_8__2(i__),v_8_0__2(i__),
+ e__8_0__2(i__),v_24__2(i__),e__24__2(i__),AKs_2(i__),
+ alpha_2(i__),n_alpha_2(i__),IfAHA_2(i__),EWHa_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' ! *Objects without Mid-Infrared Excesses: Class III
C Format for file interpretation
4 format(
+ I2,1X,I2,1X,F5.2,1X,F4.2,1X,A1,I2,1X,I2,1X,F5.2,1X,F4.2,1X,
+ F5.2,1X,F6.3,1X,F5.2,1X,F5.3,1X,F5.2,1X,F5.3,1X,F5.2,1X,F5.3,
+ 1X,F5.2,1X,F5.3,1X,F5.2,1X,F5.3,1X,F5.2,1X,F5.3,1X,F4.2,1X,
+ F5.3,1X,F4.2,1X,F6.3,A1,1X,I3,1X,F5.1)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,20
read(1,'(A153)')ar__3
read(ar__3,4)
+ RAh_3(i__),RAm_3(i__),RAs_3(i__),e_RAs_3(i__),DE__3(i__),
+ DEd_3(i__),DEm_3(i__),DEs_3(i__),e_DEs_3(i__),Jmag_3(i__),
+ e_Jmag_3(i__),Hmag_3(i__),e_Hmag_3(i__),Ksmag_3(i__),
+ e_Ksmag_3(i__),v_3_6__3(i__),e__3_6__3(i__),v_4_5__3(i__),
+ e__4_5__3(i__),v_5_8__3(i__),e__5_8__3(i__),v_8_0__3(i__),
+ e__8_0__3(i__),v_24__3(i__),e__24__3(i__),AKs_3(i__),
+ alpha_3(i__),n_alpha_3(i__),IfAHA_3(i__),EWHa_3(i__)
if(ar__3(36:40) .EQ. '') Jmag_3(i__) = rNULL__
if(ar__3(42:47) .EQ. '') e_Jmag_3(i__) = rNULL__
if(ar__3(49:53) .EQ. '') Hmag_3(i__) = rNULL__
if(ar__3(55:59) .EQ. '') e_Hmag_3(i__) = rNULL__
if(ar__3(61:65) .EQ. '') Ksmag_3(i__) = rNULL__
if(ar__3(67:71) .EQ. '') e_Ksmag_3(i__) = rNULL__
if(ar__3(73:77) .EQ. '') v_3_6__3(i__) = rNULL__
if(ar__3(79:83) .EQ. '') e__3_6__3(i__) = rNULL__
if(ar__3(85:89) .EQ. '') v_4_5__3(i__) = rNULL__
if(ar__3(91:95) .EQ. '') e__4_5__3(i__) = rNULL__
if(ar__3(97:101) .EQ. '') v_5_8__3(i__) = rNULL__
if(ar__3(103:107) .EQ. '') e__5_8__3(i__) = rNULL__
if(ar__3(109:113) .EQ. '') v_8_0__3(i__) = rNULL__
if(ar__3(115:119) .EQ. '') e__8_0__3(i__) = rNULL__
if(ar__3(121:124) .EQ. '') v_24__3(i__) = rNULL__
if(ar__3(126:130) .EQ. '') e__24__3(i__) = rNULL__
if(ar__3(132:135) .EQ. '') AKs_3(i__) = rNULL__
if(ar__3(137:142) .EQ. '') alpha_3(i__) = rNULL__
if(ar__3(145:147) .EQ. '') IfAHA_3(i__) = iNULL__
if(ar__3(149:153) .EQ. '') EWHa_3(i__) = rNULL__
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)
+ RAh_3(i__),RAm_3(i__),RAs_3(i__),e_RAs_3(i__),DE__3(i__),
+ DEd_3(i__),DEm_3(i__),DEs_3(i__),e_DEs_3(i__),Jmag_3(i__),
+ e_Jmag_3(i__),Hmag_3(i__),e_Hmag_3(i__),Ksmag_3(i__),
+ e_Ksmag_3(i__),v_3_6__3(i__),e__3_6__3(i__),v_4_5__3(i__),
+ e__4_5__3(i__),v_5_8__3(i__),e__5_8__3(i__),v_8_0__3(i__),
+ e__8_0__3(i__),v_24__3(i__),e__24__3(i__),AKs_3(i__),
+ alpha_3(i__),n_alpha_3(i__),IfAHA_3(i__),EWHa_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=============================================================================
stop
end