Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/831/155 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-May-15
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/831/155 VLA survey for faint compact radio sources in ONC (Sheehan+, 2016)
*================================================================================
*A VLA survey for faint compact radio sources in the Orion Nebula Cluster.
* Sheehan P.D., Eisner J.A., Mann R.K., Williams J.P.
* <Astrophys. J., 831, 155-155 (2016)>
* =2016ApJ...831..155S (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Source detections, identifications, and fluxes
integer*4 nr__
parameter (nr__=909) ! Number of records
character*151 ar__ ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg (nr__) ! (deg) Right Ascension J2000
real*8 DEdeg (nr__) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 Seq (nr__) ! Running sequence number
character*8 v_RRS2008_ (nr__) ! Proplyd name (from Ricci+, 2008, J/AJ/136/2136)
integer*4 v_HC2000_ (nr__) ! ? Hillenbrand & Carpenter 2000, J/ApJ/540/236
* identifier
character*2 GMR (nr__) ! Garay+ 1987ApJ...314..535G identifier
integer*4 v_ZRK2004_ (nr__) ! ? Zapata+ 2004AJ....127.2252Z identifier
character*7 OName (nr__) ! Other Names
integer*4 RAh (nr__) ! (h) [5] Hour of Right Ascension (J2000)
integer*4 RAm (nr__) ! (min) [34/35] Minute of Right Ascension (J2000)
real*4 RAs (nr__) ! (s) Second of Right Ascension (J2000)
character*1 DE_ (nr__) ! [-] Sign of the Declination (J2000)
integer*4 DEd (nr__) ! (deg) [5] Degree of Declination (J2000)
integer*4 DEm (nr__) ! (arcmin) [18/28] Arcminute of Declination (J2000)
real*4 DEs (nr__) ! (arcsec) Arcsecond of Declination (J2000)
real*4 F6 (nr__) ! (mJy) [-0.4/36] The 6cm flux density
real*4 e_F6 (nr__) ! (mJy) [0.04/4] Uncertainty in F6
real*4 F3_6 (nr__) ! (mJy) [-1/37]? The 3.6cm flux density
real*4 e_F3_6 (nr__) ! (mJy) [0.03/4]? Uncertainty in F3.6
real*4 F1_3a (nr__) ! (mJy) [-1.5/29]? Nov. 10, 2013 1.3cm flux density
real*4 e_F1_3a (nr__) ! (mJy) [0.03/38]? Uncertainty in F1.3,1
real*4 F1_3b (nr__) ! (mJy) [-1/41]? Mar. 3, 2014 1.3cm flux density
real*4 e_F1_3b (nr__) ! (mJy) [0.03/25]? Uncertainty in F1.3,2
real*4 F1_3c (nr__) ! (mJy) [-1.2/30]? Mar. 7, 2014 1.3cm flux density
real*4 e_F1_3c (nr__) ! (mJy) [0.03/70]? Uncertainty in F1.3,3
real*4 F1_3d (nr__) ! (mJy) [-1/23]? May 3, 2014 1.3cm flux density
real*4 e_F1_3d (nr__) ! (mJy) [0.03/3]? Uncertainty in F1.3,4
real*4 v_F1_3_ (nr__) ! (mJy) [-0.9/19]? Mean 1.3cm flux density
real*8 e__F1_3_ (nr__) ! (mJy) [0.01/10335]? Uncertainty in <F1.3>
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Free-free emission model parameters and ALMA
band fluxes for detected sources
integer*4 nr__1
parameter (nr__1=175) ! Number of records
character*107 ar__1 ! Full-size record
integer*4 Seq_1 (nr__1) ! [11/904] Index identifier
character*1 l_nuturn (nr__1) ! Limit flag on nuturn
real*8 nuturn (nr__1) ! (GHz) [5/300]? Frequency where the wind becomes
* partially optically thick
real*4 e_nuturn (nr__1) ! (GHz) [0.9/60]? Uncertainty in {nu}turn [GHz]
real*4 Fturn (nr__1) ! (mJy) [0.05/24.5] Flux density at nuturn
real*4 e_Fturn (nr__1) ! (mJy) [0.01/1] Uncertainty in Fturn
real*8 Fdust (nr__1) ! (mJy) [0/274]? Dust flux density
real*4 e_Fdust (nr__1) ! (mJy) [0.06/18]? Uncertainty in Fdust
real*4 S35 (nr__1) ! (mJy) [0.05/20.4] Flux density at 35GHz; Band 1
real*4 S80 (nr__1) ! (mJy) [0.05/18.7] Flux density at 80GHz; Band 2
real*4 S100 (nr__1) ! (mJy) [0.04/18.3] Flux density at 100GHz; Band 3
real*4 S150 (nr__1) ! (mJy) [0.04/17.6] Flux density at 150GHz; Band 4
real*4 S200 (nr__1) ! (mJy) [0.04/17.1] Flux density at 200GHz; Band 5
real*4 S230 (nr__1) ! (mJy) [0.04/16.9] Flux density at 230GHz; Band 6
real*4 S345 (nr__1) ! (mJy) [0.04/16.2] Flux density at 345GHz; Band 7
real*4 S450 (nr__1) ! (mJy) [0.04/15.8] Flux density at 450GHz; Band 8
real*4 S650 (nr__1) ! (mJy) [0.04/15.2] Flux density at 650GHz; Band 9
real*4 S850 (nr__1) ! (mJy) [0.04/14.8] Flux density at 850GHz; Band 10
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Variability of ONC sources
integer*4 nr__2
parameter (nr__2=147) ! Number of records
character*31 ar__2 ! Full-size record
integer*4 Seq_2 (nr__2) ! [11/897] Index number
character*1 v1_3Var (nr__2) ! 1.3cm variable?
integer*4 DelF1_3 (nr__2) ! (%) [19/905]? {Delta}F/F at 1.3cm
integer*4 e_DelF1_3 (nr__2) ! (%) [3/574]? DelF1.3 uncertainty
character*1 v3_6Var (nr__2) ! 3.6cm variable?
integer*4 DelF3_6 (nr__2) ! (%) [20/215]? {Delta}F/F at 3.6cm
integer*4 e_DelF3_6 (nr__2) ! (%) [3/46]? DelF3.6 uncertainty
character*1 v6Var (nr__2) ! 6cm variable?
integer*4 DelF6 (nr__2) ! (%) [47/102]? {Delta}F/F at 6cm
integer*4 e_DelF6 (nr__2) ! (%) [3/11]? DelF6 uncertainty
C=============================================================================
C Loading file 'table2.dat' ! Source detections, identifications, and fluxes
C Format for file interpretation
1 format(
+ I3,1X,A8,1X,I3,1X,A2,1X,I2,1X,A7,1X,I1,1X,I2,1X,F5.2,1X,A1,I1,
+ 1X,I2,1X,F5.2,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F6.3,
+ 1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,
+ F9.3)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,909
read(1,'(A151)')ar__
read(ar__,1)
+ Seq(i__),v_RRS2008_(i__),v_HC2000_(i__),GMR(i__),
+ v_ZRK2004_(i__),OName(i__),RAh(i__),RAm(i__),RAs(i__),
+ DE_(i__),DEd(i__),DEm(i__),DEs(i__),F6(i__),e_F6(i__),
+ F3_6(i__),e_F3_6(i__),F1_3a(i__),e_F1_3a(i__),F1_3b(i__),
+ e_F1_3b(i__),F1_3c(i__),e_F1_3c(i__),F1_3d(i__),e_F1_3d(i__),
+ v_F1_3_(i__),e__F1_3_(i__)
if(ar__(14:16) .EQ. '') v_HC2000_(i__) = iNULL__
if(ar__(21:22) .EQ. '') v_ZRK2004_(i__) = iNULL__
if(ar__(68:73) .EQ. '') F3_6(i__) = rNULL__
if(ar__(75:79) .EQ. '') e_F3_6(i__) = rNULL__
if(ar__(81:86) .EQ. '') F1_3a(i__) = rNULL__
if(ar__(88:93) .EQ. '') e_F1_3a(i__) = rNULL__
if(ar__(95:100) .EQ. '') F1_3b(i__) = rNULL__
if(ar__(102:107) .EQ. '') e_F1_3b(i__) = rNULL__
if(ar__(109:114) .EQ. '') F1_3c(i__) = rNULL__
if(ar__(116:121) .EQ. '') e_F1_3c(i__) = rNULL__
if(ar__(123:128) .EQ. '') F1_3d(i__) = rNULL__
if(ar__(130:134) .EQ. '') e_F1_3d(i__) = rNULL__
if(ar__(136:141) .EQ. '') v_F1_3_(i__) = rNULL__
if(ar__(143:151) .EQ. '') e__F1_3_(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)
+ Seq(i__),v_RRS2008_(i__),v_HC2000_(i__),GMR(i__),
+ v_ZRK2004_(i__),OName(i__),RAh(i__),RAm(i__),RAs(i__),
+ DE_(i__),DEd(i__),DEm(i__),DEs(i__),F6(i__),e_F6(i__),
+ F3_6(i__),e_F3_6(i__),F1_3a(i__),e_F1_3a(i__),F1_3b(i__),
+ e_F1_3b(i__),F1_3c(i__),e_F1_3c(i__),F1_3d(i__),e_F1_3d(i__),
+ v_F1_3_(i__),e__F1_3_(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' ! Free-free emission model parameters and ALMA
* band fluxes for detected sources
C Format for file interpretation
2 format(
+ I3,1X,A1,F7.3,1X,F6.3,1X,F6.3,1X,F5.3,1X,F7.3,1X,F6.3,1X,F5.2,
+ 1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,
+ F5.2,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,175
read(1,'(A107)')ar__1
read(ar__1,2)
+ Seq_1(i__),l_nuturn(i__),nuturn(i__),e_nuturn(i__),Fturn(i__),
+ e_Fturn(i__),Fdust(i__),e_Fdust(i__),S35(i__),S80(i__),
+ S100(i__),S150(i__),S200(i__),S230(i__),S345(i__),S450(i__),
+ S650(i__),S850(i__)
if(ar__1(6:12) .EQ. '') nuturn(i__) = rNULL__
if(ar__1(14:19) .EQ. '') e_nuturn(i__) = rNULL__
if(ar__1(34:40) .EQ. '') Fdust(i__) = rNULL__
if(ar__1(42:47) .EQ. '') e_Fdust(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Seq_1(i__),l_nuturn(i__),nuturn(i__),e_nuturn(i__),Fturn(i__),
+ e_Fturn(i__),Fdust(i__),e_Fdust(i__),S35(i__),S80(i__),
+ S100(i__),S150(i__),S200(i__),S230(i__),S345(i__),S450(i__),
+ S650(i__),S850(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Variability of ONC sources
C Format for file interpretation
3 format(I3,1X,A1,1X,I3,1X,I3,1X,A1,1X,I3,1X,I2,1X,A1,1X,I3,1X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,147
read(1,'(A31)')ar__2
read(ar__2,3)
+ Seq_2(i__),v1_3Var(i__),DelF1_3(i__),e_DelF1_3(i__),
+ v3_6Var(i__),DelF3_6(i__),e_DelF3_6(i__),v6Var(i__),
+ DelF6(i__),e_DelF6(i__)
if(ar__2(7:9) .EQ. '') DelF1_3(i__) = iNULL__
if(ar__2(11:13) .EQ. '') e_DelF1_3(i__) = iNULL__
if(ar__2(17:19) .EQ. '') DelF3_6(i__) = iNULL__
if(ar__2(21:22) .EQ. '') e_DelF3_6(i__) = iNULL__
if(ar__2(26:28) .EQ. '') DelF6(i__) = iNULL__
if(ar__2(30:31) .EQ. '') e_DelF6(i__) = iNULL__
c ..............Just test output...........
write(6,3)
+ Seq_2(i__),v1_3Var(i__),DelF1_3(i__),e_DelF1_3(i__),
+ v3_6Var(i__),DelF3_6(i__),e_DelF3_6(i__),v6Var(i__),
+ DelF6(i__),e_DelF6(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end