Conversion of standardized ReadMe file for
file /./ftp/cats/J/PASJ/62/391 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-20
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/PASJ/62/391 SiO and H2O masers in Stephenson 2 cluster (Deguchi+ 2010)
*================================================================================
*SiO and H2O maser observations of red supergiants in star clusters embedded in
*the Galactic disk.
* Deguchi S., Nakashima J., Zhang Y., Chong S.S.N., Koike K., Kwok S.
* <Publ. Astron. Soc. Jap., 62, 391-407 (2010)>
* =2010PASJ...62..391D
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Properties of observed objects in Stephenson
and embedded clusters
integer*4 nr__
parameter (nr__=36) ! Number of records
character*86 ar__ ! Full-size record
character*16 v2MASS ! 2MASS name (HHMMSSss+DDMMSSs)
character*6 Name ! Object name (McN-NN or St2-NN)
real*4 Kmag ! (mag) K magnitude
real*4 J_H ! (mag) J-H colour index
real*4 H_K ! (mag) H-K colour index
character*17 MSX ! MSX name (GLLL.llll+BB.bbbb)
real*4 Sep ! (arcsec) ? Separation 2MASS-MSX
real*4 FC ! (Jy) ? MSX C-band (12um) flux density
real*4 v_C_A_ ! ? Flux density ratio log(FC/FA),
* from MSX fluxes at 12um and 8um
real*4 v_E_C_ ! ? Flux density ratio log(FE/FC),
* from MSX fluxes at 21um and 12um
real*4 VSiO ! (km/s) ? SiO LSR velocity
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Observational results of the SiO maser search
integer*4 nr__1
parameter (nr__1=37) ! Number of records
character*85 ar__1 ! Full-size record
character*16 v2MASS_1 ! 2MASS name (HHMMSSss+DDMMSSs)
character*6 Name_1 ! Object name (McN-NN or St2-NN)
real*4 Ta1 ! (K) ?=- Antenna temperature for SiO
* (J=1-0 v=1, 43.122GHz)
real*4 Vlsr1 ! (km/s) ?=- LSR velocity for SiO (J=1-0 v=1)
real*4 F1 ! (K.km/s) ?=- Flux for SiO (J=1-0 v=1)
real*4 e_Ta1 ! (K) RMS value for SiO (J=1-0 v=1)
real*4 Ta2 ! (K) ?=- Antenna temperature for SiO
* (J=1-0 v=2, 42.821GHz)
real*4 Vlsr2 ! (km/s) ?=- LSR velocity for SiO (J=1-0 v=2)
real*4 F2 ! (K.km/s) ?=- Flux (SiO J=1-0 v=2)
real*4 e_Ta2 ! (K) RMS value for SiO (J=1-0 v=2)
character*12 ObsDate ! Observation date, in YYYY/MM/DD.d unit
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Observational results of the H2O maser search
integer*4 nr__2
parameter (nr__2=36) ! Number of records
character*86 ar__2 ! Full-size record
character*16 v2MASS_2 ! 2MASS name (HHMMSSss+DDMMSSs)
character*6 Name_2 ! Object name (McN-NN or St2-NN)
real*4 Ta1_1 ! (K) ?=- Antenna temperature for H2O
* (6_16_-5_23_, 22.235GHz)
character*1 n_Ta1 ! [a] contamination from Mc13-5
real*4 Vlsr1_1 ! (km/s) ?=- LSR velocity for H2O (6_16_-5_23_)
character*1 n_Vlsr1 ! [a] contamination from Mc13-5
real*4 F1_1 ! (K.km/s) ?=- Flux for H2O (6_16_-5_23_)
character*1 n_F1 ! [a] contamination from Mc13-5
real*4 e_Ta1_1 ! (K) RMS value for H2O (6_16_-5_23_)
real*4 Ta2_1 ! (K) ?=- Antenna temperature for H2O (second peak)
character*1 n_Ta2 ! [ab] contamination from Mc13-5(a) or Mc13-4(b)
real*4 Vlsr2_1 ! (km/s) ?=- LSR velocity for H2O (second peak)
character*1 n_Vlsr2 ! [ab] contamination from Mc13-5(a) or Mc13-4(b)
real*4 F2_1 ! (K.km/s) ?=- Flux for H2O (second peak)
character*1 n_F2 ! [ab] contamination from Mc13-5(a) or Mc13-4(b)
integer*4 e_Ta2_1 ! (K) ?=- RMS value for H2O (second peak)
character*12 ObsDate_1 ! Observation date, in YYYY/MM/DD.d unit
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Other names for observed objects
integer*4 nr__3
parameter (nr__3=36) ! Number of records
character*81 ar__3 ! Full-size record
character*16 v2MASS_3 ! 2MASS name (HHMMSSss+DDMMSSs)
character*6 Name_3 ! Object name (McN-NN or St2-NN)
integer*4 D07 ! ?=- Davis et al. (2007, Cat. J/ApJ/671/781)
* number, Cl* Stephenson 2 DFK NN in Simbad
character*17 MSX_1 ! MSX name (GLLL.llll+BB.bbbb)
character*31 ONames ! Other name(s)
character*5 Notes ! Notes (1)
*Note (1): 'a' for Ando et al. (2002ApJ...574..187A) sources.
C=============================================================================
C Loading file 'table1.dat' ! Properties of observed objects in Stephenson
* and embedded clusters
C Format for file interpretation
1 format(
+ A16,1X,A6,1X,F4.2,1X,F4.2,1X,F4.2,1X,A17,1X,F3.1,1X,F5.1,1X,
+ F6.3,1X,F6.3,1X,F5.1)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,36
read(1,'(A86)')ar__
read(ar__,1)
+ v2MASS,Name,Kmag,J_H,H_K,MSX,Sep,FC,v_C_A_,v_E_C_,VSiO
if(ar__(58:60) .EQ. '') Sep = rNULL__
if(ar__(62:66) .EQ. '') FC = rNULL__
if(ar__(68:73) .EQ. '') v_C_A_ = rNULL__
if(ar__(75:80) .EQ. '') v_E_C_ = rNULL__
if(ar__(82:86) .EQ. '') VSiO = rNULL__
c ..............Just test output...........
write(6,1)
+ v2MASS,Name,Kmag,J_H,H_K,MSX,Sep,FC,v_C_A_,v_E_C_,VSiO
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Observational results of the SiO maser search
C Format for file interpretation
2 format(
+ A16,1X,A6,1X,F5.3,1X,F5.1,1X,F6.3,1X,F5.3,1X,F5.3,1X,F5.1,1X,
+ F5.3,1X,F5.3,1X,A12)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,37
read(1,'(A85)')ar__1
read(ar__1,2)
+ v2MASS_1,Name_1,Ta1,Vlsr1,F1,e_Ta1,Ta2,Vlsr2,F2,e_Ta2,ObsDate
if (idig(ar__1(25:29)).EQ.0) Ta1 = rNULL__
if (idig(ar__1(31:35)).EQ.0) Vlsr1 = rNULL__
if (idig(ar__1(37:42)).EQ.0) F1 = rNULL__
if (idig(ar__1(50:54)).EQ.0) Ta2 = rNULL__
if (idig(ar__1(56:60)).EQ.0) Vlsr2 = rNULL__
if (idig(ar__1(62:66)).EQ.0) F2 = rNULL__
c ..............Just test output...........
write(6,2)
+ v2MASS_1,Name_1,Ta1,Vlsr1,F1,e_Ta1,Ta2,Vlsr2,F2,e_Ta2,ObsDate
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Observational results of the H2O maser search
C Format for file interpretation
3 format(
+ A16,1X,A6,1X,F5.3,A1,F5.1,A1,1X,F6.3,A1,1X,F5.3,1X,F5.3,A1,1X,
+ F4.1,A1,1X,F5.3,A1,1X,I2,1X,A12)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,36
read(1,'(A86)')ar__2
read(ar__2,3)
+ v2MASS_2,Name_2,Ta1_1,n_Ta1,Vlsr1_1,n_Vlsr1,F1_1,n_F1,e_Ta1_1,
+ Ta2_1,n_Ta2,Vlsr2_1,n_Vlsr2,F2_1,n_F2,e_Ta2_1,ObsDate_1
if (idig(ar__2(25:29)).EQ.0) Ta1_1 = rNULL__
if (idig(ar__2(31:35)).EQ.0) Vlsr1_1 = rNULL__
if (idig(ar__2(38:43)).EQ.0) F1_1 = rNULL__
if (idig(ar__2(52:56)).EQ.0) Ta2_1 = rNULL__
if (idig(ar__2(59:62)).EQ.0) Vlsr2_1 = rNULL__
if (idig(ar__2(65:69)).EQ.0) F2_1 = rNULL__
if (e_Ta2_1 .EQ. 45) e_Ta2_1 = iNULL__
c ..............Just test output...........
write(6,3)
+ v2MASS_2,Name_2,Ta1_1,n_Ta1,Vlsr1_1,n_Vlsr1,F1_1,n_F1,e_Ta1_1,
+ Ta2_1,n_Ta2,Vlsr2_1,n_Vlsr2,F2_1,n_F2,e_Ta2_1,ObsDate_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Other names for observed objects
C Format for file interpretation
4 format(A16,1X,A6,1X,I2,1X,A17,1X,A31,A5)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,36
read(1,'(A81)')ar__3
read(ar__3,4)v2MASS_3,Name_3,D07,MSX_1,ONames,Notes
if (D07 .EQ. 45) D07 = iNULL__
c ..............Just test output...........
write(6,4)v2MASS_3,Name_3,D07,MSX_1,ONames,Notes
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