Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/925/76 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-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 (nr__) ! Object type
integer*4 Seq (nr__) ! [1/110] Running sequence number
character*19 LAMOST (nr__) ! LAMOST name (JHHMMSS.ss+DDMMSS.s)
character*10 Date (nr__) ! ("Y/M/D") Date of observation (UT)
real*8 RAdeg (nr__) ! (deg) Right ascension in decimal degrees (J2000)
real*8 DEdeg (nr__) ! (deg) Declination in decimal degrees (J2000)
character*7 r_OType (nr__) ! Classification reference or
* "new" (63 occurrences) (1)
character*4 n_Seq (nr__) ! 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 (nr__1) ! Object type
integer*4 Seq_1 (nr__1) ! [1/110] Running sequence number
real*4 cHb (nr__1) ! [0.001/0.73] Extinction coefficient, c(H{beta})
real*4 ne (nr__1) ! (cm-3) [1/1000]? Electron density
integer*4 TOIII (nr__1) ! (K) [7283/10568]? Electron temperature, [OIII] (3)
integer*4 TNII (nr__1) ! (K) [8237/11232]? Electron temperature, [NII] (3)
integer*4 TOII (nr__1) ! (K) [6986/6986]? Electron temperature, [OII] (3)
real*4 O_H_N2 (nr__1) ! [8/8.6]? Oxygen abundance derived with N2
* diagnostics
real*4 O_H_O3 (nr__1) ! [8.1/8.6]? Oxygen abundance derived with O3N2
* diagnostics
real*4 O_H (nr__1) ! [6.3/8.9] Oxygen abundance derived with direct
* method
real*4 N_H (nr__1) ! [6.18/8.1]? Nitrogen abundance derived with
* direct method
real*4 Ne_H (nr__1) ! [7.3/8.5]? Neon abundance derived with direct
* method
real*4 S_H (nr__1) ! [5.9/7]? Sulfur abundance derived with direct
* method
real*4 Cl_H (nr__1) ! [4.6/5.2]? Chlorine abundance derived with
* direct method
real*4 Ar_H (nr__1) ! [4.9/6.8]? Argon abundance derived with direct
* method
real*4 He_H (nr__1) ! [10/11.2]? Helium abundance derived with
* direct method from optical recombination lines
integer*4 RV (nr__1) ! (km/s) [-120/114] Radial velocity (4)
integer*4 e_RV (nr__1) ! (km/s) [0/90] RV standard deviation
real*4 R (nr__1) ! (kpc) [0.29/8.1] Distance to the center of M33 (5)
real*4 R_R25 (nr__1) ! [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 (nr__2) ! [1/110] Running sequence number
character*19 LAMOST_1 (nr__2) ! LAMOST name (JHHMMSS.ss+DDMMSS.s)
real*8 lambda (nr__2) ! (0.1nm) [3726.03/9068.6] Wavelength in {AA}
character*8 Line (nr__2) ! Ion identifier
real*8 Fline (nr__2) ! [0.086/1213] Emission line intensities on
* scale Hbeta=100
real*8 Iline (nr__2) ! [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(i__),Seq(i__),LAMOST(i__),Date(i__),RAdeg(i__),
+ DEdeg(i__),r_OType(i__),n_Seq(i__)
c ..............Just test output...........
write(6,1)
+ OType(i__),Seq(i__),LAMOST(i__),Date(i__),RAdeg(i__),
+ DEdeg(i__),r_OType(i__),n_Seq(i__)
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(i__),Seq_1(i__),cHb(i__),ne(i__),TOIII(i__),TNII(i__),
+ TOII(i__),O_H_N2(i__),O_H_O3(i__),O_H(i__),N_H(i__),Ne_H(i__),
+ S_H(i__),Cl_H(i__),Ar_H(i__),He_H(i__),RV(i__),e_RV(i__),
+ R(i__),R_R25(i__)
if(ar__1(15:20) .EQ. '') ne(i__) = rNULL__
if(ar__1(22:26) .EQ. '') TOIII(i__) = iNULL__
if(ar__1(28:32) .EQ. '') TNII(i__) = iNULL__
if(ar__1(34:37) .EQ. '') TOII(i__) = iNULL__
if(ar__1(39:42) .EQ. '') O_H_N2(i__) = rNULL__
if(ar__1(44:47) .EQ. '') O_H_O3(i__) = rNULL__
if(ar__1(54:57) .EQ. '') N_H(i__) = rNULL__
if(ar__1(59:62) .EQ. '') Ne_H(i__) = rNULL__
if(ar__1(64:67) .EQ. '') S_H(i__) = rNULL__
if(ar__1(69:72) .EQ. '') Cl_H(i__) = rNULL__
if(ar__1(74:77) .EQ. '') Ar_H(i__) = rNULL__
if(ar__1(79:83) .EQ. '') He_H(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ OType_1(i__),Seq_1(i__),cHb(i__),ne(i__),TOIII(i__),TNII(i__),
+ TOII(i__),O_H_N2(i__),O_H_O3(i__),O_H(i__),N_H(i__),Ne_H(i__),
+ S_H(i__),Cl_H(i__),Ar_H(i__),He_H(i__),RV(i__),e_RV(i__),
+ R(i__),R_R25(i__)
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(i__),LAMOST_1(i__),lambda(i__),Line(i__),Fline(i__),
+ Iline(i__)
c ..............Just test output...........
write(6,3)
+ Seq_2(i__),LAMOST_1(i__),lambda(i__),Line(i__),Fline(i__),
+ Iline(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end