Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/275/20 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 2026-Aug-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/ApJS/275/20 12yr KVN obs. of SiO and H2O masers toward WX Ser (Lim+, 2024)
*================================================================================
*Long-term simultaneous monitoring observations of SiO and H2O masers toward the
*Mira variable WX Serpentis.
* Lim J.-H., Kim J., Cho S.-H., Kim H., Yoon D.-H., Son S.-M., Suh K.-W.
* <Astrophys. J. Suppl. Ser., 275, 20 (2024)>
* =2024ApJS..275...20L
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Observational results of SiO transitions
integer*4 nr__
parameter (nr__=188) ! Number of records
character*107 ar__ ! Full-size record
character*17 Maser ! Maser Line (G1)
integer*4 Epoch ! [1/50] Epoch number
character*10 ObsDate ! ("Y/M/D") Observation date (UT)
real*4 Phase ! (d) [0.18/10.6] Phase
character*27 Speak ! (Jy) Peak Flux density, S_{nu}_, list (1)
real*4 Snudnu ! (Jy.km.s-1) [1.7/270.7]? Integral, S{nu}*d{nu}
real*4 rms ! (Jy) [0.1/0.7]? Continuum root-mean-square noise
* in observation
character*24 Vpeak ! (km/s) Peak V, list (1)
real*4 Vmean ! (km.s-1) [3.8/9.1]? Mean maser velocity
*Note (1): Speak and Vpeak values are given as lists of detected
* maser lines at that Epoch. The units of the Speak and Vpeak values are
* Jy and km.s-1, respectively.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Observational results of the H2O maser line
integer*4 nr__1
parameter (nr__1=50) ! Number of records
character*91 ar__1 ! Full-size record
integer*4 Epoch_1 ! [1/50] Epoch number
character*10 ObsDate_1 ! ("Y/M/D") Observation date (UT)
real*4 Phase_1 ! (d) [0.18/10.6] Phase
real*4 Speak_B4 ! (Jy) [0.8/44]? Peak Flux density, Blue-shifted
* Line, 4
real*4 Speak_B3 ! (Jy) [0.7/4.5]? Peak Flux density,
* Blue-shifted Line, 3
real*4 Speak_B2 ! (Jy) [0.8/121]? Peak Flux density,
* Blue-shifted Line, 2
real*4 Speak_B1 ! (Jy) [0.7/26.4]? Peak Flux density,
* Blue-shifted Line, 1
real*4 Speak_R1 ! (Jy) [1/29.3]? Peak Flux density, Red-shifted
* Line, 1
real*4 Speak_R2 ! (Jy) [3/18.7]? Peak Flux density, Red-shifted
* Line, 2
real*4 Snudnu_1 ! (Jy.km.s-1) [4/224.6] Integral, S{nu}*d{nu}
real*4 rms_1 ! (Jy) [0.1/0.5] Continuum root-mean-square
* noise in observation
real*4 Vpeak_B4 ! (km.s-1) [-3/-1.5]? Peak V, Blue-shifted Line, 4
real*4 Vpeak_B3 ! (km.s-1) [-1.3/-0.3]? Peak V, Blue-shifted Line, 3
real*4 Vpeak_B2 ! (km.s-1) [0.4/2.2]? Peak V, Blue-shifted Line, 2
real*4 Vpeak_B1 ! (km.s-1) [2.5/4.7]? Peak V, Blue-shifted Line, 1
real*4 Vpeak_R1 ! (km.s-1) [6.0/7.9]? Peak V, Red-shifted Line,1
real*4 Vpeak_R2 ! (km.s-1) [8.3/9]? Peak V, Red-shifted Line,2
real*4 Vmean_1 ! (km.s-1) [0.3/7.7] Mean maser velocity
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! Maser intensity ratios
integer*4 nr__2
parameter (nr__2=52) ! Number of records
character*45 ar__2 ! Full-size record
character*4 Epoch_2 ! Epoch number (1)
character*10 ObsDate_2 ! ("Y/M/D") Observation date (UT)
real*4 Phase_2 ! [0/1]? Optical Phase, normalized
real*4 Ratio1 ! [0.09/2.8] Maser Intensity Ratio,
* I.F.(H2O)/I.F.(SiO43) (G1)
real*4 Ratio2 ! [0.27/2.2] Maser Intensity Ratio,
* I.F.(SiO42)/I.F.(SiO43) (G1)
real*4 Ratio3 ! [0.18/1.51]? Maser Intensity Ratio,
* I.F.(SiO86)/I.F.(SiO43) (G1)
real*4 Ratio4 ! [0.01/0.31]? Maser Intensity Ratio,
* I.F.(SiO129)/I.F.(SiO43) (G1)
real*4 Ratio5 ! [0.02/0.5]? Maser Intensity Ratio,
* I.F.(SiO129)/I.F.(SiO86) (G1)
*Note (1): The last two rows of this table provide the mean and standard
* deviation of each column.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! Velocity extents of SiO and H2O maser lines
integer*4 nr__3
parameter (nr__3=50) ! Number of records
character*68 ar__3 ! Full-size record
integer*4 Epoch_3 ! [1/50] Epoch number
character*10 ObsDate_3 ! ("Y/M/D") Observation date (UT)
real*4 Phase_3 ! (d) [0.18/10.6] Phase
real*4 VB_SiO43 ! (km.s-1) [-3.1/2.6] Blue velocity edge, SiO43 (G1)
real*4 VR_SiO43 ! (km.s-1) [9.2/17.4] Red velocity edge, SiO43 (G1)
real*4 VB_SiO42 ! (km.s-1) [-4.0/3.5] Blue velocity edge, SiO42 (G1)
real*4 VR_SiO42 ! (km.s-1) [8.4/15.2] Red velocity edge, SiO42 (G1)
real*4 VB_SiO86 ! (km.s-1) [-3.1/3.4]? Blue velocity edge, SiO86 (G1)
real*4 VR_SiO86 ! (km.s-1) [10.9/15.6]? Red velocity edge, SiO86 (G1)
real*4 VB_SiO129 ! (km.s-1) [0.3/5.1]? Blue velocity edge, SiO129 (G1)
real*4 VR_SiO129 ! (km.s-1) [7.9/14.6]? Red velocity edge, SiO129 (G1)
real*4 VB_H2O ! (km.s-1) [-8.7/2.4] Blue velocity edge, H2O (G1)
real*4 VR_H2O ! (km.s-1) [9.2/16.6] Red velocity edge, H2O (G1)
C=============================================================================
C Loading file 'table3.dat' ! Observational results of SiO transitions
C Format for file interpretation
1 format(
+ A17,1X,I2,1X,A10,1X,F5.2,1X,A27,1X,F6.2,1X,F5.3,1X,A24,1X,F3.1)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,188
read(1,'(A107)')ar__
read(ar__,1)
+ Maser,Epoch,ObsDate,Phase,Speak,Snudnu,rms,Vpeak,Vmean
if(ar__(67:72) .EQ. '') Snudnu = rNULL__
if(ar__(74:78) .EQ. '') rms = rNULL__
if(ar__(105:107) .EQ. '') Vmean = rNULL__
c ..............Just test output...........
write(6,1)
+ Maser,Epoch,ObsDate,Phase,Speak,Snudnu,rms,Vpeak,Vmean
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Observational results of the H2O maser line
C Format for file interpretation
2 format(
+ I2,1X,A10,1X,F5.2,1X,F4.1,1X,F3.1,1X,F5.1,1X,F4.1,1X,F4.1,1X,
+ F4.1,1X,F5.1,1X,F5.3,1X,F4.1,1X,F4.1,1X,F3.1,1X,F3.1,1X,F3.1,
+ 1X,F3.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,50
read(1,'(A91)')ar__1
read(ar__1,2)
+ Epoch_1,ObsDate_1,Phase_1,Speak_B4,Speak_B3,Speak_B2,Speak_B1,
+ Speak_R1,Speak_R2,Snudnu_1,rms_1,Vpeak_B4,Vpeak_B3,Vpeak_B2,
+ Vpeak_B1,Vpeak_R1,Vpeak_R2,Vmean_1
if(ar__1(21:24) .EQ. '') Speak_B4 = rNULL__
if(ar__1(26:28) .EQ. '') Speak_B3 = rNULL__
if(ar__1(30:34) .EQ. '') Speak_B2 = rNULL__
if(ar__1(36:39) .EQ. '') Speak_B1 = rNULL__
if(ar__1(41:44) .EQ. '') Speak_R1 = rNULL__
if(ar__1(46:49) .EQ. '') Speak_R2 = rNULL__
if(ar__1(63:66) .EQ. '') Vpeak_B4 = rNULL__
if(ar__1(68:71) .EQ. '') Vpeak_B3 = rNULL__
if(ar__1(73:75) .EQ. '') Vpeak_B2 = rNULL__
if(ar__1(77:79) .EQ. '') Vpeak_B1 = rNULL__
if(ar__1(81:83) .EQ. '') Vpeak_R1 = rNULL__
if(ar__1(85:87) .EQ. '') Vpeak_R2 = rNULL__
c ..............Just test output...........
write(6,2)
+ Epoch_1,ObsDate_1,Phase_1,Speak_B4,Speak_B3,Speak_B2,Speak_B1,
+ Speak_R1,Speak_R2,Snudnu_1,rms_1,Vpeak_B4,Vpeak_B3,Vpeak_B2,
+ Vpeak_B1,Vpeak_R1,Vpeak_R2,Vmean_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! Maser intensity ratios
C Format for file interpretation
3 format(
+ A4,1X,A10,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,52
read(1,'(A45)')ar__2
read(ar__2,3)
+ Epoch_2,ObsDate_2,Phase_2,Ratio1,Ratio2,Ratio3,Ratio4,Ratio5
if(ar__2(17:20) .EQ. '') Phase_2 = rNULL__
if(ar__2(32:35) .EQ. '') Ratio3 = rNULL__
if(ar__2(37:40) .EQ. '') Ratio4 = rNULL__
if(ar__2(42:45) .EQ. '') Ratio5 = rNULL__
c ..............Just test output...........
write(6,3)
+ Epoch_2,ObsDate_2,Phase_2,Ratio1,Ratio2,Ratio3,Ratio4,Ratio5
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7.dat' ! Velocity extents of SiO and H2O maser lines
C Format for file interpretation
4 format(
+ I2,1X,A10,1X,F5.2,1X,F4.1,1X,F4.1,1X,F4.1,1X,F4.1,1X,F4.1,1X,
+ F4.1,1X,F3.1,1X,F4.1,1X,F4.1,1X,F4.1)
C Effective file loading
open(unit=1,status='old',file=
+'table7.dat')
write(6,*) '....Loading file: table7.dat'
do i__=1,50
read(1,'(A68)')ar__3
read(ar__3,4)
+ Epoch_3,ObsDate_3,Phase_3,VB_SiO43,VR_SiO43,VB_SiO42,VR_SiO42,
+ VB_SiO86,VR_SiO86,VB_SiO129,VR_SiO129,VB_H2O,VR_H2O
if(ar__3(41:44) .EQ. '') VB_SiO86 = rNULL__
if(ar__3(46:49) .EQ. '') VR_SiO86 = rNULL__
if(ar__3(51:53) .EQ. '') VB_SiO129 = rNULL__
if(ar__3(55:58) .EQ. '') VR_SiO129 = rNULL__
c ..............Just test output...........
write(6,4)
+ Epoch_3,ObsDate_3,Phase_3,VB_SiO43,VR_SiO43,VB_SiO42,VR_SiO42,
+ VB_SiO86,VR_SiO86,VB_SiO129,VR_SiO129,VB_H2O,VR_H2O
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end