Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/786/37 into FORTRAN code for reading data files line by line.
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-Apr-18
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/786/37 The Auriga-California molecular cloud (Broekhoven-Fiene+, 2014)
*================================================================================
*The Spitzer survey of interstellar clouds in the Gould belt.
*VI. The Auriga-California molecular cloud observed with IRAC and MIPS.
* Broekhoven-Fiene H., Matthews B.C., Harvey P.M., Gutermuth R.A.,
* Huard T.L., Tothill N.F.H., Nutter D., Bourke T.L., DiFrancesco J.,
* Jorgensen J.K., Allen L.E., Chapman N.L., Cieza L.A., Dunham M.M.,
* Merin B., Miller J.F., Terebey S., Peterson D.E., Stapelfeldt K.R.
* <Astrophys. J., 786, 37 (2014)>
* =2014ApJ...786...37B (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! YSOs in the AMC Based on IRAC and MIPS
integer*4 nr__
parameter (nr__=138) ! Number of records
character*186 ar__ ! Full-size record
integer*4 ID ! [1/166] YSO ID number
character*16 Name ! Name specifying J2000 coordinates
* (HHMMSSss+DDMMSSs)
character*3 Class ! [I II III F] YSO SED Class
real*4 alpha ! Infrared spectral slope
real*8 Flux3_6 ! (mJy) ? Flux at IRAC1 3.6 micron band
real*4 e_Flux3_6 ! (mJy) ? Error in Flux3.6
real*8 Flux4_5 ! (mJy) ? Flux at IRAC2 4.5 micron band
real*4 e_Flux4_5 ! (mJy) ? Error in Flux4.5
real*8 Flux5_8 ! (mJy) ? Flux at IRAC3 5.8 micron band
real*4 e_Flux5_8 ! (mJy) ? Error in Flux5.8
real*8 Flux8_0 ! (mJy) ? Flux at IRAC4 8.0 micron band
real*4 e_Flux8_0 ! (mJy) ? Error in Flux8.0
real*8 Flux3_4 ! (mJy) ? Flux at WISE1 3.4 micron band
real*4 e_Flux3_4 ! (mJy) ? Error in Flux3.4
real*8 Flux4_6 ! (mJy) ? Flux at WISE2 4.6 micron band
real*4 e_Flux4_6 ! (mJy) ? Error in Flux4.6
real*8 Flux12 ! (mJy) ? Flux at WISE3 12 micron band
real*4 e_Flux12 ! (mJy) ? Error in Flux12
character*1 l_Flux12 ! [<] Limit flag on Flux12
real*4 Flux22 ! (mJy) ? Flux at WISE4 22 micron band
real*4 e_Flux22 ! (mJy) ? Error in Flux22
character*1 l_Flux24 ! [<] Limit flag on Flux24
real*8 Flux24 ! (mJy) ? Flux at MIPS1 24 micron band
real*8 e_Flux24 ! (mJy) ? Error in Flux24
real*8 Flux70 ! (mJy) ? Flux at MIPS2 70 micron band
real*8 e_Flux70 ! (mJy) ? Error in Flux70
character*1 lowNH2 ! [LN] Low column density flag (1)
*Note (1): Flag as follows:
* L = The YSO is in a region of low column density, N_H2_<5x10^21^cm^-2^, and
* so is a possible contaminant;
* N = The YSO lies beyond the N_H2_ column density map from Harvey et al.
* (2013, J/ApJ/764/133) and so N_H2_ at its position is unknown.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! YSO Candidates in the AMC Based on WISE and MIPS
integer*4 nr__1
parameter (nr__1=28) ! Number of records
character*186 ar__1 ! Full-size record
integer*4 ID_1 ! [1/166] YSO ID number
character*16 Name_1 ! Name specifying J2000 coordinates
* (HHMMSSss+DDMMSSs)
character*3 Class_1 ! [I II III F] YSO SED Class
real*4 alpha_1 ! Infrared spectral slope
real*8 Flux3_6_1 ! (mJy) ? Flux at IRAC1 3.6 micron band
real*4 e_Flux3_6_1 ! (mJy) ? Error in Flux3.6
real*8 Flux4_5_1 ! (mJy) ? Flux at IRAC2 4.5 micron band
real*4 e_Flux4_5_1 ! (mJy) ? Error in Flux4.5
real*8 Flux5_8_1 ! (mJy) ? Flux at IRAC3 5.8 micron band
real*4 e_Flux5_8_1 ! (mJy) ? Error in Flux5.8
real*8 Flux8_0_1 ! (mJy) ? Flux at IRAC4 8.0 micron band
real*4 e_Flux8_0_1 ! (mJy) ? Error in Flux8.0
real*8 Flux3_4_1 ! (mJy) ? Flux at WISE1 3.4 micron band
real*4 e_Flux3_4_1 ! (mJy) ? Error in Flux3.4
real*8 Flux4_6_1 ! (mJy) ? Flux at WISE2 4.6 micron band
real*4 e_Flux4_6_1 ! (mJy) ? Error in Flux4.6
real*8 Flux12_1 ! (mJy) ? Flux at WISE3 12 micron band
real*4 e_Flux12_1 ! (mJy) ? Error in Flux12
character*1 l_Flux12_1 ! [<] Limit flag on Flux12
real*4 Flux22_1 ! (mJy) ? Flux at WISE4 22 micron band
real*4 e_Flux22_1 ! (mJy) ? Error in Flux22
character*1 l_Flux24_1 ! [<] Limit flag on Flux24
real*8 Flux24_1 ! (mJy) ? Flux at MIPS1 24 micron band
real*8 e_Flux24_1 ! (mJy) ? Error in Flux24
real*8 Flux70_1 ! (mJy) ? Flux at MIPS2 70 micron band
real*8 e_Flux70_1 ! (mJy) ? Error in Flux70
character*1 lowNH2_1 ! [LN] Low column density flag (1)
*Note (1): Flag as follows:
* L = The YSO is in a region of low column density, N_H2_<5x10^21^cm^-2^, and
* so is a possible contaminant;
* N = The YSO lies beyond the N_H2_ column density map from Harvey et al.
* (2013, J/ApJ/764/133) and so N_H2_ at its position is unknown.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! SED modelling results in Class II sources
integer*4 nr__2
parameter (nr__2=91) ! Number of records
character*37 ar__2 ! Full-size record
integer*4 ID_2 ! [2/166] YSO ID number
character*2 Fitted ! [K7-A0] Fitted stellar spectrum to the SED
real*4 Av ! (mag) Extinction value fitted to the optical data
real*4 L_ ! (Lsun) Luminosity of the fitted stellar component
real*4 lambda_TO ! (um) Last wavelength before an excess >80%
* is measured (1)
real*4 alpha_ex ! ? Slope of SED across wavelengths longward
* of lambda-TO (2)
real*8 Ld_L_ ! Ratio of disk luminosity to luminosity of
* fitted stellar component
*Note (1): If no excess >80% is observed, then lambda-TO is set to 24 microns.
*Note (2): No alpha-ex value is listed for sources with lambda-TO at 24 micron
* because there are not enough data points to measure the slope of
* the excess.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table8.dat' ! SED modelling results in Class III sources
integer*4 nr__3
parameter (nr__3=17) ! Number of records
character*37 ar__3 ! Full-size record
integer*4 ID_3 ! [2/166] YSO ID number
character*2 Fitted_1 ! [K7-A0] Fitted stellar spectrum to the SED
real*4 Av_1 ! (mag) Extinction value fitted to the optical data
real*4 L__1 ! (Lsun) Luminosity of the fitted stellar component
real*4 lambda_TO_1 ! (um) Last wavelength before an excess >80%
* is measured (1)
real*4 alpha_ex_1 ! ? Slope of SED across wavelengths longward
* of lambda-TO (2)
real*8 Ld_L__1 ! Ratio of disk luminosity to luminosity of
* fitted stellar component
*Note (1): If no excess >80% is observed, then lambda-TO is set to 24 microns.
*Note (2): No alpha-ex value is listed for sources with lambda-TO at 24 micron
* because there are not enough data points to measure the slope of
* the excess.
C=============================================================================
C Loading file 'table4.dat' ! YSOs in the AMC Based on IRAC and MIPS
C Format for file interpretation
1 format(
+ I3,1X,A16,1X,A3,1X,F5.2,1X,F8.3,1X,F6.3,1X,F7.2,1X,F6.2,1X,
+ F7.2,1X,F6.2,1X,F7.2,1X,F6.2,1X,F7.2,1X,F6.2,1X,F7.2,1X,F6.2,
+ 1X,F7.2,1X,F5.2,1X,A1,F6.2,1X,F5.2,1X,A1,F8.2,1X,F7.2,1X,F7.1,
+ 1X,F8.3,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,138
read(1,'(A186)')ar__
read(ar__,1)
+ ID,Name,Class,alpha,Flux3_6,e_Flux3_6,Flux4_5,e_Flux4_5,
+ Flux5_8,e_Flux5_8,Flux8_0,e_Flux8_0,Flux3_4,e_Flux3_4,Flux4_6,
+ e_Flux4_6,Flux12,e_Flux12,l_Flux12,Flux22,e_Flux22,l_Flux24,
+ Flux24,e_Flux24,Flux70,e_Flux70,lowNH2
if(ar__(32:39) .EQ. '') Flux3_6 = rNULL__
if(ar__(41:46) .EQ. '') e_Flux3_6 = rNULL__
if(ar__(48:54) .EQ. '') Flux4_5 = rNULL__
if(ar__(56:61) .EQ. '') e_Flux4_5 = rNULL__
if(ar__(63:69) .EQ. '') Flux5_8 = rNULL__
if(ar__(71:76) .EQ. '') e_Flux5_8 = rNULL__
if(ar__(78:84) .EQ. '') Flux8_0 = rNULL__
if(ar__(86:91) .EQ. '') e_Flux8_0 = rNULL__
if(ar__(93:99) .EQ. '') Flux3_4 = rNULL__
if(ar__(101:106) .EQ. '') e_Flux3_4 = rNULL__
if(ar__(108:114) .EQ. '') Flux4_6 = rNULL__
if(ar__(116:121) .EQ. '') e_Flux4_6 = rNULL__
if(ar__(123:129) .EQ. '') Flux12 = rNULL__
if(ar__(131:135) .EQ. '') e_Flux12 = rNULL__
if(ar__(138:143) .EQ. '') Flux22 = rNULL__
if(ar__(145:149) .EQ. '') e_Flux22 = rNULL__
if(ar__(152:159) .EQ. '') Flux24 = rNULL__
if(ar__(161:167) .EQ. '') e_Flux24 = rNULL__
if(ar__(169:175) .EQ. '') Flux70 = rNULL__
if(ar__(177:184) .EQ. '') e_Flux70 = rNULL__
c ..............Just test output...........
write(6,1)
+ ID,Name,Class,alpha,Flux3_6,e_Flux3_6,Flux4_5,e_Flux4_5,
+ Flux5_8,e_Flux5_8,Flux8_0,e_Flux8_0,Flux3_4,e_Flux3_4,Flux4_6,
+ e_Flux4_6,Flux12,e_Flux12,l_Flux12,Flux22,e_Flux22,l_Flux24,
+ Flux24,e_Flux24,Flux70,e_Flux70,lowNH2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! YSO Candidates in the AMC Based on WISE and MIPS
C Format for file interpretation
2 format(
+ I3,1X,A16,1X,A3,1X,F5.2,1X,F8.3,1X,F6.3,1X,F7.2,1X,F6.2,1X,
+ F7.2,1X,F6.2,1X,F7.2,1X,F6.2,1X,F7.2,1X,F6.2,1X,F7.2,1X,F6.2,
+ 1X,F7.2,1X,F5.2,1X,A1,F6.2,1X,F5.2,1X,A1,F8.2,1X,F7.2,1X,F7.1,
+ 1X,F8.3,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,28
read(1,'(A186)')ar__1
read(ar__1,2)
+ ID_1,Name_1,Class_1,alpha_1,Flux3_6_1,e_Flux3_6_1,Flux4_5_1,
+ e_Flux4_5_1,Flux5_8_1,e_Flux5_8_1,Flux8_0_1,e_Flux8_0_1,
+ Flux3_4_1,e_Flux3_4_1,Flux4_6_1,e_Flux4_6_1,Flux12_1,
+ e_Flux12_1,l_Flux12_1,Flux22_1,e_Flux22_1,l_Flux24_1,Flux24_1,
+ e_Flux24_1,Flux70_1,e_Flux70_1,lowNH2_1
if(ar__1(32:39) .EQ. '') Flux3_6_1 = rNULL__
if(ar__1(41:46) .EQ. '') e_Flux3_6_1 = rNULL__
if(ar__1(48:54) .EQ. '') Flux4_5_1 = rNULL__
if(ar__1(56:61) .EQ. '') e_Flux4_5_1 = rNULL__
if(ar__1(63:69) .EQ. '') Flux5_8_1 = rNULL__
if(ar__1(71:76) .EQ. '') e_Flux5_8_1 = rNULL__
if(ar__1(78:84) .EQ. '') Flux8_0_1 = rNULL__
if(ar__1(86:91) .EQ. '') e_Flux8_0_1 = rNULL__
if(ar__1(93:99) .EQ. '') Flux3_4_1 = rNULL__
if(ar__1(101:106) .EQ. '') e_Flux3_4_1 = rNULL__
if(ar__1(108:114) .EQ. '') Flux4_6_1 = rNULL__
if(ar__1(116:121) .EQ. '') e_Flux4_6_1 = rNULL__
if(ar__1(123:129) .EQ. '') Flux12_1 = rNULL__
if(ar__1(131:135) .EQ. '') e_Flux12_1 = rNULL__
if(ar__1(138:143) .EQ. '') Flux22_1 = rNULL__
if(ar__1(145:149) .EQ. '') e_Flux22_1 = rNULL__
if(ar__1(152:159) .EQ. '') Flux24_1 = rNULL__
if(ar__1(161:167) .EQ. '') e_Flux24_1 = rNULL__
if(ar__1(169:175) .EQ. '') Flux70_1 = rNULL__
if(ar__1(177:184) .EQ. '') e_Flux70_1 = rNULL__
c ..............Just test output...........
write(6,2)
+ ID_1,Name_1,Class_1,alpha_1,Flux3_6_1,e_Flux3_6_1,Flux4_5_1,
+ e_Flux4_5_1,Flux5_8_1,e_Flux5_8_1,Flux8_0_1,e_Flux8_0_1,
+ Flux3_4_1,e_Flux3_4_1,Flux4_6_1,e_Flux4_6_1,Flux12_1,
+ e_Flux12_1,l_Flux12_1,Flux22_1,e_Flux22_1,l_Flux24_1,Flux24_1,
+ e_Flux24_1,Flux70_1,e_Flux70_1,lowNH2_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7.dat' ! SED modelling results in Class II sources
C Format for file interpretation
3 format(I3,1X,A2,1X,F4.1,1X,F6.2,1X,F4.1,1X,F4.1,1X,F8.6)
C Effective file loading
open(unit=1,status='old',file=
+'table7.dat')
write(6,*) '....Loading file: table7.dat'
do i__=1,91
read(1,'(A37)')ar__2
read(ar__2,3)ID_2,Fitted,Av,L_,lambda_TO,alpha_ex,Ld_L_
if(ar__2(25:28) .EQ. '') alpha_ex = rNULL__
c ..............Just test output...........
write(6,3)ID_2,Fitted,Av,L_,lambda_TO,alpha_ex,Ld_L_
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table8.dat' ! SED modelling results in Class III sources
C Format for file interpretation
4 format(I3,1X,A2,1X,F4.1,1X,F6.2,1X,F4.1,1X,F4.1,1X,F8.6)
C Effective file loading
open(unit=1,status='old',file=
+'table8.dat')
write(6,*) '....Loading file: table8.dat'
do i__=1,17
read(1,'(A37)')ar__3
read(ar__3,4)
+ ID_3,Fitted_1,Av_1,L__1,lambda_TO_1,alpha_ex_1,Ld_L__1
if(ar__3(25:28) .EQ. '') alpha_ex_1 = rNULL__
c ..............Just test output...........
write(6,4)
+ ID_3,Fitted_1,Av_1,L__1,lambda_TO_1,alpha_ex_1,Ld_L__1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end