Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/925/76 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-Oct-13
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/925/76 M33 with LAMOST survey sp. I. PNe & HII regions (Alexeeva+, 2022)
*================================================================================
*Spectroscopic study of M33 with the LAMOST survey.
*I. Chemical gradients from nebulae.
* Alexeeva S., Zhao G.
* <Astrophys. J., 925, 76 (2022)>
* =2022ApJ...925...76A
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Sample of PNe and HII regions in M33 from
LAMOST survey
integer*4 nr__
parameter (nr__=110) ! Number of records
character*71 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*3 OType ! Object type
integer*4 Seq ! [1/110] Running sequence number
character*19 LAMOST ! LAMOST name (JHHMMSS.ss+DDMMSS.s)
character*10 Date ! ("Y/M/D") Date of observation (UT)
real*8 RAdeg ! (deg) Right ascension in decimal degrees (J2000)
real*8 DEdeg ! (deg) Declination in decimal degrees (J2000)
character*7 r_OType ! Classification reference or
* "new" (63 occurrences) (1)
character*4 n_Seq ! Object type uncertainty for #84 (2)
*Note (1): We call regions "new" if there is no previous classification as
* an HII region or PN in the literature for objects within the radius
* of 2" (the site seeing of the telescope). For regions that were
* classified in previous works, we provide references as follows:
* ZCH2020 = Zhang et al. (2020RAA....20...97Z)
* CDL2004 = Ciardullo et al. (2004ApJ...614..167C)
* MCM2000 = Magrini et al. (2000A&A...355..713M)
* VGHC86 = Viallefond et al. (1986A&AS...64..237V)
* DMS93 = Drissen et al. (1993AJ....105.1400D)
* BCL74 = Boulesteix et al. (1974A&A....37...33B)
* RVP2013 = Relano et al. (2013, J/A+A/552/A140)
* LHK2017 = Lin et al. (2017, J/ApJ/842/97)
* B2011 = Bresolin (2011ApJ...730..129B)
* MVM2007 = Magrini et al. (2007A&A...470..865M)
* CPS87 = Courtes et al. (1987A&A...174...28C)
* LJD2011 = Lagrois et al. (2011MNRAS.413..705L)
* CS82 = Christian & Schommer (1982ApJS...49..405C)
* MOH2006 = Massey et al. (2006, J/AJ/131/2478)
*Note (2): The nebula J013340.36+302001.6 (#84) was classified as an HII region
* in our study. However, this nebula is located very close to
* J013340.37+302001.6, which was classified as a new SNR candidate in
* Zhang+ (2020RAA....20...97Z), who employed an additional
* classification based on Frew & Parker (2010PASA...27..129F).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! The chemical abundances and radial velocities
of nebulae
integer*4 nr__1
parameter (nr__1=110) ! Number of records
character*104 ar__1 ! Full-size record
character*3 OType_1 ! Object type
integer*4 Seq_1 ! [1/110] Running sequence number
real*4 cHb ! [0.001/0.73] Extinction coefficient, c(H{beta})
real*4 ne ! (cm-3) [1/1000]? Electron density
integer*4 TOIII ! (K) [7283/10568]? Electron temperature, [OIII] (3)
integer*4 TNII ! (K) [8237/11232]? Electron temperature, [NII] (3)
integer*4 TOII ! (K) [6986/6986]? Electron temperature, [OII] (3)
real*4 O_H_N2 ! [8/8.6]? Oxygen abundance derived with N2
* diagnostics
real*4 O_H_O3 ! [8.1/8.6]? Oxygen abundance derived with O3N2
* diagnostics
real*4 O_H ! [6.3/8.9] Oxygen abundance derived with direct
* method
real*4 N_H ! [6.18/8.1]? Nitrogen abundance derived with
* direct method
real*4 Ne_H ! [7.3/8.5]? Neon abundance derived with direct
* method
real*4 S_H ! [5.9/7]? Sulfur abundance derived with direct
* method
real*4 Cl_H ! [4.6/5.2]? Chlorine abundance derived with
* direct method
real*4 Ar_H ! [4.9/6.8]? Argon abundance derived with direct
* method
real*4 He_H ! [10/11.2]? Helium abundance derived with
* direct method from optical recombination lines
integer*4 RV ! (km/s) [-120/114] Radial velocity (4)
integer*4 e_RV ! (km/s) [0/90] RV standard deviation
real*4 R ! (kpc) [0.29/8.1] Distance to the center of M33 (5)
real*4 R_R25 ! [0.04/1.2] R to R25 ratio (6)
*Note (3): A value of Te=10000K was assumed to calculate the ionic and
* total abundances, when Te is not available.
*Note (4): Radial velocity was calculated from Halfa, Hbeta, and doublet
* of [SII]6716/6731{AA}. Systemic velocity of M33, -179km/s, has been
* subtracted.
*Note (5): With M33 at RA(J2000)=01:33:50.9,DE(J2000)=+30:39:36.8
* (from Plucinsky+ 2008, J/ApJS/174/366);
* Position angle (major axis)=23{deg} and Inclination=56{deg}
* (from Zaritsky+ 1989AJ.....97...97Z); see Table 2.
*Note (6): Assumed an R25 of 28.2 arcmin (de Vaucouleurs 1959ApJ...130..728D).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'flux.dat' ! Measured fluxes, fluxes corrected for extinction
for particular ions
integer*4 nr__2
parameter (nr__2=2052) ! Number of records
character*60 ar__2 ! Full-size record
integer*4 Seq_2 ! [1/110] Running sequence number
character*19 LAMOST_1 ! LAMOST name (JHHMMSS.ss+DDMMSS.s)
real*8 lambda ! (0.1nm) [3726.03/9068.6] Wavelength in {AA}
character*8 Line ! Ion identifier
real*8 Fline ! [0.086/1213] Emission line intensities on
* scale Hbeta=100
real*8 Iline ! [0.076/1213] Emission line intensities on
* scale Hbeta=100 (corrected for extinction)
C=============================================================================
C Loading file 'table1.dat' ! Sample of PNe and HII regions in M33 from
* LAMOST survey
C Format for file interpretation
1 format(A3,1X,I3,1X,A19,1X,A10,1X,F9.6,1X,F9.6,1X,A7,1X,A4)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,110
read(1,'(A71)')ar__
read(ar__,1)OType,Seq,LAMOST,Date,RAdeg,DEdeg,r_OType,n_Seq
c ..............Just test output...........
write(6,1)OType,Seq,LAMOST,Date,RAdeg,DEdeg,r_OType,n_Seq
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! The chemical abundances and radial velocities
* of nebulae
C Format for file interpretation
2 format(
+ A3,1X,I3,1X,F5.3,1X,F6.1,1X,I5,1X,I5,1X,I4,1X,F4.2,1X,F4.2,1X,
+ F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,I4,1X,
+ I2,1X,F5.3,1X,F6.4)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,110
read(1,'(A104)')ar__1
read(ar__1,2)
+ OType_1,Seq_1,cHb,ne,TOIII,TNII,TOII,O_H_N2,O_H_O3,O_H,N_H,
+ Ne_H,S_H,Cl_H,Ar_H,He_H,RV,e_RV,R,R_R25
if(ar__1(15:20) .EQ. '') ne = rNULL__
if(ar__1(22:26) .EQ. '') TOIII = iNULL__
if(ar__1(28:32) .EQ. '') TNII = iNULL__
if(ar__1(34:37) .EQ. '') TOII = iNULL__
if(ar__1(39:42) .EQ. '') O_H_N2 = rNULL__
if(ar__1(44:47) .EQ. '') O_H_O3 = rNULL__
if(ar__1(54:57) .EQ. '') N_H = rNULL__
if(ar__1(59:62) .EQ. '') Ne_H = rNULL__
if(ar__1(64:67) .EQ. '') S_H = rNULL__
if(ar__1(69:72) .EQ. '') Cl_H = rNULL__
if(ar__1(74:77) .EQ. '') Ar_H = rNULL__
if(ar__1(79:83) .EQ. '') He_H = rNULL__
c ..............Just test output...........
write(6,2)
+ OType_1,Seq_1,cHb,ne,TOIII,TNII,TOII,O_H_N2,O_H_O3,O_H,N_H,
+ Ne_H,S_H,Cl_H,Ar_H,He_H,RV,e_RV,R,R_R25
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'flux.dat' ! Measured fluxes, fluxes corrected for extinction
* for particular ions
C Format for file interpretation
3 format(I3,1X,A19,1X,F7.2,1X,A8,3X,F8.3,1X,F8.3)
C Effective file loading
open(unit=1,status='old',file=
+'flux.dat')
write(6,*) '....Loading file: flux.dat'
do i__=1,2052
read(1,'(A60)')ar__2
read(ar__2,3)Seq_2,LAMOST_1,lambda,Line,Fline,Iline
c ..............Just test output...........
write(6,3)Seq_2,LAMOST_1,lambda,Line,Fline,Iline
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end