Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/818/116 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-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/818/116 Gould's Belt VLA survey. V. Perseus region (Pech+, 2016)
*================================================================================
*The Gould's Belt very large array survey.
*V. The Perseus region.
* Pech G., Loinard L., Dzib S.A., Mioduszewski A.J., Rodriguez L.F.,
* Ortiz-Leon G.N., Rivera J.L., Torres R.M., Boden A.F., Hartman L.,
* Kounkel M.A., Evans Ii N.J., Briceno C., Tobin J., Zapata L.A.
* <Astrophys. J., 818, 116 (2016)>
* =2016ApJ...818..116P (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! VLA Observations
integer*4 nr__
parameter (nr__=12) ! Number of records
character*50 ar__ ! Full-size record
character*8 Region (nr__) ! Region identifier
character*1 n_Region (nr__) ! [d] Note on Region (1)
character*1 Epoch (nr__) ! [123C] Epoch (2)
real*4 Bmaj4_5 (nr__) ! (arcsec) Synthesised beam major axis at 4.5GHz
real*4 Bmin4_5 (nr__) ! (arcsec) Synthesised beam minor axis at 4.5GHz
real*4 PA4_5 (nr__) ! (deg) [] Position angle of synthesized beam
* at 4.5GHz
real*4 Bmaj7_5 (nr__) ! (arcsec) Synthesised beam major axis at 7.5GHz
real*4 Bmin7_5 (nr__) ! (arcsec) Synthesised beam minor axis at 7.5GHz
real*4 PA7_5 (nr__) ! (deg) [] Position angle of synthesized beam
* at 7.5GHz
integer*4 rms4_5 (nr__) ! (uJy/beam) rms noise measured at the center of
* the Stokes I image, at 4.5GHz
integer*4 rms7_5 (nr__) ! (uJy/beam) rms noise measured at the center of
* the Stokes I image, at 7.5GHz
*Note (1): Note as follows:
* d = Average values in the seven individual fields.
*Note (2): C indicates parameters measured in the images after combining
* the epochs.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Radio Sources Detected in IC 348
integer*4 nr__1
parameter (nr__1=91) ! Number of records
character*71 ar__1 ! Full-size record
character*19 Name (nr__1) ! GBS-VLA Name (JHHMMSS.ss+DDMMSS.s)
character*1 New (nr__1) ! [YN] New source? (1)
real*4 F4_5 (nr__1) ! (mJy) Flux density at 4.5GHz
real*4 e_F4_5 (nr__1) ! (mJy) Uncertainty in F4.5
character*1 l_Var4_5 (nr__1) ! [>] Limit flag on Var4.5
integer*4 Var4_5 (nr__1) ! (%) [4/96]?=- Variability of the source at 4.5GHz
integer*4 e_Var4_5 (nr__1) ! (%) [1/68]? Uncertainty in Var4.5
character*1 n_Var4_5 (nr__1) ! [E] Note on Var4.5 (2)
character*1 l_F7_5 (nr__1) ! [<] Limit flag on F7.5
real*4 F7_5 (nr__1) ! (mJy) ?=- Flux density at 7.5GHz
real*4 e_F7_5 (nr__1) ! (mJy) ? Uncertainty in F7.5
character*1 l_Var7_5 (nr__1) ! [>] Limit flag on Var7.5
integer*4 Var7_5 (nr__1) ! (%) [2/88]?=- Variability of the source at 7.5GHz
integer*4 e_Var7_5 (nr__1) ! (%) [6/49]? Uncertainty in Var7.5
character*1 n_Var7_5 (nr__1) ! [AE] Note on Var7.5 (2)
character*1 l_alpha (nr__1) ! [<] Limit flag on alpha
real*4 alpha (nr__1) ! ?=- Spectral index
real*4 e_alpha (nr__1) ! ? Uncertainty in alpha
*Note (1): Source as follows:
* Y = Source without reported counterparts at any frequency;
* N = Source with known counterpart.
*Note (2): Note as follows:
* A = Source not detected at three times the noise level on individual epochs,
* but detected on the image of the concatenated epochs;
* E = Sources with extended emission features.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Radio Sources Detected in Single Fields
in Perseus
integer*4 nr__2
parameter (nr__2=41) ! Number of records
character*71 ar__2 ! Full-size record
character*19 Name_1 (nr__2) ! GBS-VLA Name (JHHMMSS.ss+DDMMSS.s)
character*1 New_1 (nr__2) ! [YN] New source? (1)
real*4 F4_5_1 (nr__2) ! (mJy) Flux density at 4.5GHz
real*4 e_F4_5_1 (nr__2) ! (mJy) Uncertainty in F4.5
character*1 l_Var4_5_1 (nr__2) ! [>] Limit flag on Var4.5
integer*4 Var4_5_1 (nr__2) ! (%) [4/96]?=- Variability of the source at 4.5GHz
integer*4 e_Var4_5_1 (nr__2) ! (%) [1/68]? Uncertainty in Var4.5
character*1 n_Var4_5_1 (nr__2) ! [E] Note on Var4.5 (2)
character*1 l_F7_5_1 (nr__2) ! [<] Limit flag on F7.5
real*4 F7_5_1 (nr__2) ! (mJy) ?=- Flux density at 7.5GHz
real*4 e_F7_5_1 (nr__2) ! (mJy) ? Uncertainty in F7.5
character*1 l_Var7_5_1 (nr__2) ! [>] Limit flag on Var7.5
integer*4 Var7_5_1 (nr__2) ! (%) [2/88]?=- Variability of the source at 7.5GHz
integer*4 e_Var7_5_1 (nr__2) ! (%) [6/49]? Uncertainty in Var7.5
character*1 n_Var7_5_1 (nr__2) ! [AE] Note on Var7.5 (2)
character*1 l_alpha_1 (nr__2) ! [<] Limit flag on alpha
real*4 alpha_1 (nr__2) ! ?=- Spectral index
real*4 e_alpha_1 (nr__2) ! ? Uncertainty in alpha
*Note (1): Source as follows:
* Y = Source without reported counterparts at any frequency;
* N = Source with known counterpart.
*Note (2): Note as follows:
* A = Source not detected at three times the noise level on individual epochs,
* but detected on the image of the concatenated epochs;
* E = Sources with extended emission features.
C=============================================================================
C Loading file 'table1.dat' ! VLA Observations
C Format for file interpretation
1 format(
+ A8,A1,1X,A1,1X,F4.2,1X,F4.2,1X,F5.1,1X,F4.2,1X,F4.2,1X,F6.2,
+ 1X,I2,1X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,12
read(1,'(A50)')ar__
read(ar__,1)
+ Region(i__),n_Region(i__),Epoch(i__),Bmaj4_5(i__),
+ Bmin4_5(i__),PA4_5(i__),Bmaj7_5(i__),Bmin7_5(i__),PA7_5(i__),
+ rms4_5(i__),rms7_5(i__)
c ..............Just test output...........
write(6,1)
+ Region(i__),n_Region(i__),Epoch(i__),Bmaj4_5(i__),
+ Bmin4_5(i__),PA4_5(i__),Bmaj7_5(i__),Bmin7_5(i__),PA7_5(i__),
+ rms4_5(i__),rms7_5(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Radio Sources Detected in IC 348
C Format for file interpretation
2 format(
+ A19,1X,A1,1X,F5.2,1X,F5.3,1X,A1,I2,1X,I2,A1,1X,A1,F5.2,1X,
+ F4.2,1X,A1,I2,1X,I2,A1,1X,A1,F4.1,1X,F3.1)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,91
read(1,'(A71)')ar__1
read(ar__1,2)
+ Name(i__),New(i__),F4_5(i__),e_F4_5(i__),l_Var4_5(i__),
+ Var4_5(i__),e_Var4_5(i__),n_Var4_5(i__),l_F7_5(i__),F7_5(i__),
+ e_F7_5(i__),l_Var7_5(i__),Var7_5(i__),e_Var7_5(i__),
+ n_Var7_5(i__),l_alpha(i__),alpha(i__),e_alpha(i__)
if (Var4_5(i__) .EQ. 45) Var4_5(i__) = iNULL__
if(ar__1(39:40) .EQ. '') e_Var4_5(i__) = iNULL__
if (idig(ar__1(44:48)).EQ.0) F7_5(i__) = rNULL__
if(ar__1(50:53) .EQ. '') e_F7_5(i__) = rNULL__
if (Var7_5(i__) .EQ. 45) Var7_5(i__) = iNULL__
if(ar__1(59:60) .EQ. '') e_Var7_5(i__) = iNULL__
if (idig(ar__1(64:67)).EQ.0) alpha(i__) = rNULL__
if(ar__1(69:71) .EQ. '') e_alpha(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Name(i__),New(i__),F4_5(i__),e_F4_5(i__),l_Var4_5(i__),
+ Var4_5(i__),e_Var4_5(i__),n_Var4_5(i__),l_F7_5(i__),F7_5(i__),
+ e_F7_5(i__),l_Var7_5(i__),Var7_5(i__),e_Var7_5(i__),
+ n_Var7_5(i__),l_alpha(i__),alpha(i__),e_alpha(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Radio Sources Detected in Single Fields
* in Perseus
C Format for file interpretation
3 format(
+ A19,1X,A1,1X,F5.2,1X,F5.3,1X,A1,I2,1X,I2,A1,1X,A1,F5.2,1X,
+ F4.2,1X,A1,I2,1X,I2,A1,1X,A1,F4.1,1X,F3.1)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,41
read(1,'(A71)')ar__2
read(ar__2,3)
+ Name_1(i__),New_1(i__),F4_5_1(i__),e_F4_5_1(i__),
+ l_Var4_5_1(i__),Var4_5_1(i__),e_Var4_5_1(i__),n_Var4_5_1(i__),
+ l_F7_5_1(i__),F7_5_1(i__),e_F7_5_1(i__),l_Var7_5_1(i__),
+ Var7_5_1(i__),e_Var7_5_1(i__),n_Var7_5_1(i__),l_alpha_1(i__),
+ alpha_1(i__),e_alpha_1(i__)
if (Var4_5_1(i__) .EQ. 45) Var4_5_1(i__) = iNULL__
if(ar__2(39:40) .EQ. '') e_Var4_5_1(i__) = iNULL__
if (idig(ar__2(44:48)).EQ.0) F7_5_1(i__) = rNULL__
if(ar__2(50:53) .EQ. '') e_F7_5_1(i__) = rNULL__
if (Var7_5_1(i__) .EQ. 45) Var7_5_1(i__) = iNULL__
if(ar__2(59:60) .EQ. '') e_Var7_5_1(i__) = iNULL__
if (idig(ar__2(64:67)).EQ.0) alpha_1(i__) = rNULL__
if(ar__2(69:71) .EQ. '') e_alpha_1(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ Name_1(i__),New_1(i__),F4_5_1(i__),e_F4_5_1(i__),
+ l_Var4_5_1(i__),Var4_5_1(i__),e_Var4_5_1(i__),n_Var4_5_1(i__),
+ l_F7_5_1(i__),F7_5_1(i__),e_F7_5_1(i__),l_Var7_5_1(i__),
+ Var7_5_1(i__),e_Var7_5_1(i__),n_Var7_5_1(i__),l_alpha_1(i__),
+ alpha_1(i__),e_alpha_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end
C Locate position of first digit in string; or return 0
integer function idig(c)
character*(*) c
character*1 c1
integer lc,i
lc=len(c)
idig=0
do i=1,lc
if(c(i:i).ne.' ') go to 1
end do
1 if(i.gt.lc) return
c1=c(i:i)
if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1
if(i.gt.lc) return
c1=c(i:i)
if(c1.ge.'0'.and.c1.le.'9') idig=i
return
end