Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/786/37 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-Apr-19
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 (nr__) ! [1/166] YSO ID number
character*16 Name (nr__) ! Name specifying J2000 coordinates
* (HHMMSSss+DDMMSSs)
character*3 Class (nr__) ! [I II III F] YSO SED Class
real*4 alpha (nr__) ! Infrared spectral slope
real*8 Flux3_6 (nr__) ! (mJy) ? Flux at IRAC1 3.6 micron band
real*4 e_Flux3_6 (nr__) ! (mJy) ? Error in Flux3.6
real*8 Flux4_5 (nr__) ! (mJy) ? Flux at IRAC2 4.5 micron band
real*4 e_Flux4_5 (nr__) ! (mJy) ? Error in Flux4.5
real*8 Flux5_8 (nr__) ! (mJy) ? Flux at IRAC3 5.8 micron band
real*4 e_Flux5_8 (nr__) ! (mJy) ? Error in Flux5.8
real*8 Flux8_0 (nr__) ! (mJy) ? Flux at IRAC4 8.0 micron band
real*4 e_Flux8_0 (nr__) ! (mJy) ? Error in Flux8.0
real*8 Flux3_4 (nr__) ! (mJy) ? Flux at WISE1 3.4 micron band
real*4 e_Flux3_4 (nr__) ! (mJy) ? Error in Flux3.4
real*8 Flux4_6 (nr__) ! (mJy) ? Flux at WISE2 4.6 micron band
real*4 e_Flux4_6 (nr__) ! (mJy) ? Error in Flux4.6
real*8 Flux12 (nr__) ! (mJy) ? Flux at WISE3 12 micron band
real*4 e_Flux12 (nr__) ! (mJy) ? Error in Flux12
character*1 l_Flux12 (nr__) ! [<] Limit flag on Flux12
real*4 Flux22 (nr__) ! (mJy) ? Flux at WISE4 22 micron band
real*4 e_Flux22 (nr__) ! (mJy) ? Error in Flux22
character*1 l_Flux24 (nr__) ! [<] Limit flag on Flux24
real*8 Flux24 (nr__) ! (mJy) ? Flux at MIPS1 24 micron band
real*8 e_Flux24 (nr__) ! (mJy) ? Error in Flux24
real*8 Flux70 (nr__) ! (mJy) ? Flux at MIPS2 70 micron band
real*8 e_Flux70 (nr__) ! (mJy) ? Error in Flux70
character*1 lowNH2 (nr__) ! [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 (nr__1) ! [1/166] YSO ID number
character*16 Name_1 (nr__1) ! Name specifying J2000 coordinates
* (HHMMSSss+DDMMSSs)
character*3 Class_1 (nr__1) ! [I II III F] YSO SED Class
real*4 alpha_1 (nr__1) ! Infrared spectral slope
real*8 Flux3_6_1 (nr__1) ! (mJy) ? Flux at IRAC1 3.6 micron band
real*4 e_Flux3_6_1(nr__1) ! (mJy) ? Error in Flux3.6
real*8 Flux4_5_1 (nr__1) ! (mJy) ? Flux at IRAC2 4.5 micron band
real*4 e_Flux4_5_1(nr__1) ! (mJy) ? Error in Flux4.5
real*8 Flux5_8_1 (nr__1) ! (mJy) ? Flux at IRAC3 5.8 micron band
real*4 e_Flux5_8_1(nr__1) ! (mJy) ? Error in Flux5.8
real*8 Flux8_0_1 (nr__1) ! (mJy) ? Flux at IRAC4 8.0 micron band
real*4 e_Flux8_0_1(nr__1) ! (mJy) ? Error in Flux8.0
real*8 Flux3_4_1 (nr__1) ! (mJy) ? Flux at WISE1 3.4 micron band
real*4 e_Flux3_4_1(nr__1) ! (mJy) ? Error in Flux3.4
real*8 Flux4_6_1 (nr__1) ! (mJy) ? Flux at WISE2 4.6 micron band
real*4 e_Flux4_6_1(nr__1) ! (mJy) ? Error in Flux4.6
real*8 Flux12_1 (nr__1) ! (mJy) ? Flux at WISE3 12 micron band
real*4 e_Flux12_1 (nr__1) ! (mJy) ? Error in Flux12
character*1 l_Flux12_1 (nr__1) ! [<] Limit flag on Flux12
real*4 Flux22_1 (nr__1) ! (mJy) ? Flux at WISE4 22 micron band
real*4 e_Flux22_1 (nr__1) ! (mJy) ? Error in Flux22
character*1 l_Flux24_1 (nr__1) ! [<] Limit flag on Flux24
real*8 Flux24_1 (nr__1) ! (mJy) ? Flux at MIPS1 24 micron band
real*8 e_Flux24_1 (nr__1) ! (mJy) ? Error in Flux24
real*8 Flux70_1 (nr__1) ! (mJy) ? Flux at MIPS2 70 micron band
real*8 e_Flux70_1 (nr__1) ! (mJy) ? Error in Flux70
character*1 lowNH2_1 (nr__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 (nr__2) ! [2/166] YSO ID number
character*2 Fitted (nr__2) ! [K7-A0] Fitted stellar spectrum to the SED
real*4 Av (nr__2) ! (mag) Extinction value fitted to the optical data
real*4 L_ (nr__2) ! (Lsun) Luminosity of the fitted stellar component
real*4 lambda_TO (nr__2) ! (um) Last wavelength before an excess >80%
* is measured (1)
real*4 alpha_ex (nr__2) ! ? Slope of SED across wavelengths longward
* of lambda-TO (2)
real*8 Ld_L_ (nr__2) ! 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 (nr__3) ! [2/166] YSO ID number
character*2 Fitted_1 (nr__3) ! [K7-A0] Fitted stellar spectrum to the SED
real*4 Av_1 (nr__3) ! (mag) Extinction value fitted to the optical data
real*4 L__1 (nr__3) ! (Lsun) Luminosity of the fitted stellar component
real*4 lambda_TO_1(nr__3) ! (um) Last wavelength before an excess >80%
* is measured (1)
real*4 alpha_ex_1 (nr__3) ! ? Slope of SED across wavelengths longward
* of lambda-TO (2)
real*8 Ld_L__1 (nr__3) ! 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(i__),Name(i__),Class(i__),alpha(i__),Flux3_6(i__),
+ e_Flux3_6(i__),Flux4_5(i__),e_Flux4_5(i__),Flux5_8(i__),
+ e_Flux5_8(i__),Flux8_0(i__),e_Flux8_0(i__),Flux3_4(i__),
+ e_Flux3_4(i__),Flux4_6(i__),e_Flux4_6(i__),Flux12(i__),
+ e_Flux12(i__),l_Flux12(i__),Flux22(i__),e_Flux22(i__),
+ l_Flux24(i__),Flux24(i__),e_Flux24(i__),Flux70(i__),
+ e_Flux70(i__),lowNH2(i__)
if(ar__(32:39) .EQ. '') Flux3_6(i__) = rNULL__
if(ar__(41:46) .EQ. '') e_Flux3_6(i__) = rNULL__
if(ar__(48:54) .EQ. '') Flux4_5(i__) = rNULL__
if(ar__(56:61) .EQ. '') e_Flux4_5(i__) = rNULL__
if(ar__(63:69) .EQ. '') Flux5_8(i__) = rNULL__
if(ar__(71:76) .EQ. '') e_Flux5_8(i__) = rNULL__
if(ar__(78:84) .EQ. '') Flux8_0(i__) = rNULL__
if(ar__(86:91) .EQ. '') e_Flux8_0(i__) = rNULL__
if(ar__(93:99) .EQ. '') Flux3_4(i__) = rNULL__
if(ar__(101:106) .EQ. '') e_Flux3_4(i__) = rNULL__
if(ar__(108:114) .EQ. '') Flux4_6(i__) = rNULL__
if(ar__(116:121) .EQ. '') e_Flux4_6(i__) = rNULL__
if(ar__(123:129) .EQ. '') Flux12(i__) = rNULL__
if(ar__(131:135) .EQ. '') e_Flux12(i__) = rNULL__
if(ar__(138:143) .EQ. '') Flux22(i__) = rNULL__
if(ar__(145:149) .EQ. '') e_Flux22(i__) = rNULL__
if(ar__(152:159) .EQ. '') Flux24(i__) = rNULL__
if(ar__(161:167) .EQ. '') e_Flux24(i__) = rNULL__
if(ar__(169:175) .EQ. '') Flux70(i__) = rNULL__
if(ar__(177:184) .EQ. '') e_Flux70(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ ID(i__),Name(i__),Class(i__),alpha(i__),Flux3_6(i__),
+ e_Flux3_6(i__),Flux4_5(i__),e_Flux4_5(i__),Flux5_8(i__),
+ e_Flux5_8(i__),Flux8_0(i__),e_Flux8_0(i__),Flux3_4(i__),
+ e_Flux3_4(i__),Flux4_6(i__),e_Flux4_6(i__),Flux12(i__),
+ e_Flux12(i__),l_Flux12(i__),Flux22(i__),e_Flux22(i__),
+ l_Flux24(i__),Flux24(i__),e_Flux24(i__),Flux70(i__),
+ e_Flux70(i__),lowNH2(i__)
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(i__),Name_1(i__),Class_1(i__),alpha_1(i__),
+ Flux3_6_1(i__),e_Flux3_6_1(i__),Flux4_5_1(i__),
+ e_Flux4_5_1(i__),Flux5_8_1(i__),e_Flux5_8_1(i__),
+ Flux8_0_1(i__),e_Flux8_0_1(i__),Flux3_4_1(i__),
+ e_Flux3_4_1(i__),Flux4_6_1(i__),e_Flux4_6_1(i__),
+ Flux12_1(i__),e_Flux12_1(i__),l_Flux12_1(i__),Flux22_1(i__),
+ e_Flux22_1(i__),l_Flux24_1(i__),Flux24_1(i__),e_Flux24_1(i__),
+ Flux70_1(i__),e_Flux70_1(i__),lowNH2_1(i__)
if(ar__1(32:39) .EQ. '') Flux3_6_1(i__) = rNULL__
if(ar__1(41:46) .EQ. '') e_Flux3_6_1(i__) = rNULL__
if(ar__1(48:54) .EQ. '') Flux4_5_1(i__) = rNULL__
if(ar__1(56:61) .EQ. '') e_Flux4_5_1(i__) = rNULL__
if(ar__1(63:69) .EQ. '') Flux5_8_1(i__) = rNULL__
if(ar__1(71:76) .EQ. '') e_Flux5_8_1(i__) = rNULL__
if(ar__1(78:84) .EQ. '') Flux8_0_1(i__) = rNULL__
if(ar__1(86:91) .EQ. '') e_Flux8_0_1(i__) = rNULL__
if(ar__1(93:99) .EQ. '') Flux3_4_1(i__) = rNULL__
if(ar__1(101:106) .EQ. '') e_Flux3_4_1(i__) = rNULL__
if(ar__1(108:114) .EQ. '') Flux4_6_1(i__) = rNULL__
if(ar__1(116:121) .EQ. '') e_Flux4_6_1(i__) = rNULL__
if(ar__1(123:129) .EQ. '') Flux12_1(i__) = rNULL__
if(ar__1(131:135) .EQ. '') e_Flux12_1(i__) = rNULL__
if(ar__1(138:143) .EQ. '') Flux22_1(i__) = rNULL__
if(ar__1(145:149) .EQ. '') e_Flux22_1(i__) = rNULL__
if(ar__1(152:159) .EQ. '') Flux24_1(i__) = rNULL__
if(ar__1(161:167) .EQ. '') e_Flux24_1(i__) = rNULL__
if(ar__1(169:175) .EQ. '') Flux70_1(i__) = rNULL__
if(ar__1(177:184) .EQ. '') e_Flux70_1(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ ID_1(i__),Name_1(i__),Class_1(i__),alpha_1(i__),
+ Flux3_6_1(i__),e_Flux3_6_1(i__),Flux4_5_1(i__),
+ e_Flux4_5_1(i__),Flux5_8_1(i__),e_Flux5_8_1(i__),
+ Flux8_0_1(i__),e_Flux8_0_1(i__),Flux3_4_1(i__),
+ e_Flux3_4_1(i__),Flux4_6_1(i__),e_Flux4_6_1(i__),
+ Flux12_1(i__),e_Flux12_1(i__),l_Flux12_1(i__),Flux22_1(i__),
+ e_Flux22_1(i__),l_Flux24_1(i__),Flux24_1(i__),e_Flux24_1(i__),
+ Flux70_1(i__),e_Flux70_1(i__),lowNH2_1(i__)
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(i__),Fitted(i__),Av(i__),L_(i__),lambda_TO(i__),
+ alpha_ex(i__),Ld_L_(i__)
if(ar__2(25:28) .EQ. '') alpha_ex(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ ID_2(i__),Fitted(i__),Av(i__),L_(i__),lambda_TO(i__),
+ alpha_ex(i__),Ld_L_(i__)
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(i__),Fitted_1(i__),Av_1(i__),L__1(i__),lambda_TO_1(i__),
+ alpha_ex_1(i__),Ld_L__1(i__)
if(ar__3(25:28) .EQ. '') alpha_ex_1(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ ID_3(i__),Fitted_1(i__),Av_1(i__),L__1(i__),lambda_TO_1(i__),
+ alpha_ex_1(i__),Ld_L__1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end