Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/275/20 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-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 (nr__) ! Maser Line (G1)
integer*4 Epoch (nr__) ! [1/50] Epoch number
character*10 ObsDate (nr__) ! ("Y/M/D") Observation date (UT)
real*4 Phase (nr__) ! (d) [0.18/10.6] Phase
character*27 Speak (nr__) ! (Jy) Peak Flux density, S_{nu}_, list (1)
real*4 Snudnu (nr__) ! (Jy.km.s-1) [1.7/270.7]? Integral, S{nu}*d{nu}
real*4 rms (nr__) ! (Jy) [0.1/0.7]? Continuum root-mean-square noise
* in observation
character*24 Vpeak (nr__) ! (km/s) Peak V, list (1)
real*4 Vmean (nr__) ! (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 (nr__1) ! [1/50] Epoch number
character*10 ObsDate_1 (nr__1) ! ("Y/M/D") Observation date (UT)
real*4 Phase_1 (nr__1) ! (d) [0.18/10.6] Phase
real*4 Speak_B4 (nr__1) ! (Jy) [0.8/44]? Peak Flux density, Blue-shifted
* Line, 4
real*4 Speak_B3 (nr__1) ! (Jy) [0.7/4.5]? Peak Flux density,
* Blue-shifted Line, 3
real*4 Speak_B2 (nr__1) ! (Jy) [0.8/121]? Peak Flux density,
* Blue-shifted Line, 2
real*4 Speak_B1 (nr__1) ! (Jy) [0.7/26.4]? Peak Flux density,
* Blue-shifted Line, 1
real*4 Speak_R1 (nr__1) ! (Jy) [1/29.3]? Peak Flux density, Red-shifted
* Line, 1
real*4 Speak_R2 (nr__1) ! (Jy) [3/18.7]? Peak Flux density, Red-shifted
* Line, 2
real*4 Snudnu_1 (nr__1) ! (Jy.km.s-1) [4/224.6] Integral, S{nu}*d{nu}
real*4 rms_1 (nr__1) ! (Jy) [0.1/0.5] Continuum root-mean-square
* noise in observation
real*4 Vpeak_B4 (nr__1) ! (km.s-1) [-3/-1.5]? Peak V, Blue-shifted Line, 4
real*4 Vpeak_B3 (nr__1) ! (km.s-1) [-1.3/-0.3]? Peak V, Blue-shifted Line, 3
real*4 Vpeak_B2 (nr__1) ! (km.s-1) [0.4/2.2]? Peak V, Blue-shifted Line, 2
real*4 Vpeak_B1 (nr__1) ! (km.s-1) [2.5/4.7]? Peak V, Blue-shifted Line, 1
real*4 Vpeak_R1 (nr__1) ! (km.s-1) [6.0/7.9]? Peak V, Red-shifted Line,1
real*4 Vpeak_R2 (nr__1) ! (km.s-1) [8.3/9]? Peak V, Red-shifted Line,2
real*4 Vmean_1 (nr__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 (nr__2) ! Epoch number (1)
character*10 ObsDate_2 (nr__2) ! ("Y/M/D") Observation date (UT)
real*4 Phase_2 (nr__2) ! [0/1]? Optical Phase, normalized
real*4 Ratio1 (nr__2) ! [0.09/2.8] Maser Intensity Ratio,
* I.F.(H2O)/I.F.(SiO43) (G1)
real*4 Ratio2 (nr__2) ! [0.27/2.2] Maser Intensity Ratio,
* I.F.(SiO42)/I.F.(SiO43) (G1)
real*4 Ratio3 (nr__2) ! [0.18/1.51]? Maser Intensity Ratio,
* I.F.(SiO86)/I.F.(SiO43) (G1)
real*4 Ratio4 (nr__2) ! [0.01/0.31]? Maser Intensity Ratio,
* I.F.(SiO129)/I.F.(SiO43) (G1)
real*4 Ratio5 (nr__2) ! [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 (nr__3) ! [1/50] Epoch number
character*10 ObsDate_3 (nr__3) ! ("Y/M/D") Observation date (UT)
real*4 Phase_3 (nr__3) ! (d) [0.18/10.6] Phase
real*4 VB_SiO43 (nr__3) ! (km.s-1) [-3.1/2.6] Blue velocity edge, SiO43 (G1)
real*4 VR_SiO43 (nr__3) ! (km.s-1) [9.2/17.4] Red velocity edge, SiO43 (G1)
real*4 VB_SiO42 (nr__3) ! (km.s-1) [-4.0/3.5] Blue velocity edge, SiO42 (G1)
real*4 VR_SiO42 (nr__3) ! (km.s-1) [8.4/15.2] Red velocity edge, SiO42 (G1)
real*4 VB_SiO86 (nr__3) ! (km.s-1) [-3.1/3.4]? Blue velocity edge, SiO86 (G1)
real*4 VR_SiO86 (nr__3) ! (km.s-1) [10.9/15.6]? Red velocity edge, SiO86 (G1)
real*4 VB_SiO129 (nr__3) ! (km.s-1) [0.3/5.1]? Blue velocity edge, SiO129 (G1)
real*4 VR_SiO129 (nr__3) ! (km.s-1) [7.9/14.6]? Red velocity edge, SiO129 (G1)
real*4 VB_H2O (nr__3) ! (km.s-1) [-8.7/2.4] Blue velocity edge, H2O (G1)
real*4 VR_H2O (nr__3) ! (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(i__),Epoch(i__),ObsDate(i__),Phase(i__),Speak(i__),
+ Snudnu(i__),rms(i__),Vpeak(i__),Vmean(i__)
if(ar__(67:72) .EQ. '') Snudnu(i__) = rNULL__
if(ar__(74:78) .EQ. '') rms(i__) = rNULL__
if(ar__(105:107) .EQ. '') Vmean(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Maser(i__),Epoch(i__),ObsDate(i__),Phase(i__),Speak(i__),
+ Snudnu(i__),rms(i__),Vpeak(i__),Vmean(i__)
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(i__),ObsDate_1(i__),Phase_1(i__),Speak_B4(i__),
+ Speak_B3(i__),Speak_B2(i__),Speak_B1(i__),Speak_R1(i__),
+ Speak_R2(i__),Snudnu_1(i__),rms_1(i__),Vpeak_B4(i__),
+ Vpeak_B3(i__),Vpeak_B2(i__),Vpeak_B1(i__),Vpeak_R1(i__),
+ Vpeak_R2(i__),Vmean_1(i__)
if(ar__1(21:24) .EQ. '') Speak_B4(i__) = rNULL__
if(ar__1(26:28) .EQ. '') Speak_B3(i__) = rNULL__
if(ar__1(30:34) .EQ. '') Speak_B2(i__) = rNULL__
if(ar__1(36:39) .EQ. '') Speak_B1(i__) = rNULL__
if(ar__1(41:44) .EQ. '') Speak_R1(i__) = rNULL__
if(ar__1(46:49) .EQ. '') Speak_R2(i__) = rNULL__
if(ar__1(63:66) .EQ. '') Vpeak_B4(i__) = rNULL__
if(ar__1(68:71) .EQ. '') Vpeak_B3(i__) = rNULL__
if(ar__1(73:75) .EQ. '') Vpeak_B2(i__) = rNULL__
if(ar__1(77:79) .EQ. '') Vpeak_B1(i__) = rNULL__
if(ar__1(81:83) .EQ. '') Vpeak_R1(i__) = rNULL__
if(ar__1(85:87) .EQ. '') Vpeak_R2(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Epoch_1(i__),ObsDate_1(i__),Phase_1(i__),Speak_B4(i__),
+ Speak_B3(i__),Speak_B2(i__),Speak_B1(i__),Speak_R1(i__),
+ Speak_R2(i__),Snudnu_1(i__),rms_1(i__),Vpeak_B4(i__),
+ Vpeak_B3(i__),Vpeak_B2(i__),Vpeak_B1(i__),Vpeak_R1(i__),
+ Vpeak_R2(i__),Vmean_1(i__)
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(i__),ObsDate_2(i__),Phase_2(i__),Ratio1(i__),
+ Ratio2(i__),Ratio3(i__),Ratio4(i__),Ratio5(i__)
if(ar__2(17:20) .EQ. '') Phase_2(i__) = rNULL__
if(ar__2(32:35) .EQ. '') Ratio3(i__) = rNULL__
if(ar__2(37:40) .EQ. '') Ratio4(i__) = rNULL__
if(ar__2(42:45) .EQ. '') Ratio5(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ Epoch_2(i__),ObsDate_2(i__),Phase_2(i__),Ratio1(i__),
+ Ratio2(i__),Ratio3(i__),Ratio4(i__),Ratio5(i__)
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(i__),ObsDate_3(i__),Phase_3(i__),VB_SiO43(i__),
+ VR_SiO43(i__),VB_SiO42(i__),VR_SiO42(i__),VB_SiO86(i__),
+ VR_SiO86(i__),VB_SiO129(i__),VR_SiO129(i__),VB_H2O(i__),
+ VR_H2O(i__)
if(ar__3(41:44) .EQ. '') VB_SiO86(i__) = rNULL__
if(ar__3(46:49) .EQ. '') VR_SiO86(i__) = rNULL__
if(ar__3(51:53) .EQ. '') VB_SiO129(i__) = rNULL__
if(ar__3(55:58) .EQ. '') VR_SiO129(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ Epoch_3(i__),ObsDate_3(i__),Phase_3(i__),VB_SiO43(i__),
+ VR_SiO43(i__),VB_SiO42(i__),VR_SiO42(i__),VB_SiO86(i__),
+ VR_SiO86(i__),VB_SiO129(i__),VR_SiO129(i__),VB_H2O(i__),
+ VR_H2O(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end