Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/658/A17 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-Jun-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/A+A/658/A17 Post-common-envelope PN and mol. mass (Santander-Garcia+, 2022)
*================================================================================
*The ionised and molecular mass of post-common-envelope planetary nebulae.
*The missing mass problem.
* Santander-Garcia M., Jones D., Alcolea J., Bujarrabal V., Wesson R.
* <Astron. Astrophys. 658, A17 (2022)>
* =2022A&A...658A..17S (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Computed ionised and molecular masses of
post-common-envelope PNe
integer*4 nr__
parameter (nr__=21) ! Number of records
character*99 ar__ ! Full-size record
character*11 PNG (nr__) ! Galactic PN designation
character*14 Name (nr__) ! Common PN name
real*4 Dist (nr__) ! (kpc) Distance to the PN
real*4 E_Dist (nr__) ! (kpc) Positive (upper) error on distance
real*4 e_Dist_1 (nr__) ! (kpc) Negative (lower) error on distance
real*4 Mion (nr__) ! (Msun) Ionised mass
real*4 E_Mion (nr__) ! (Msun) Positive (upper) error on ionised mass
real*4 e_Mion_1 (nr__) ! (Msun) Negative (lower) error on ionised mass
character*1 l_Mmol (nr__) ! [<] Upper limit flag on Mmol
real*4 Mmol (nr__) ! (Msun) Molecular mass (or upper limit, see flag)
real*4 e_Mmol (nr__) ! (Msun) ? Error on molecular mass, only if
* Mmol is not an upper limit
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea1.dat' ! Computed ionised and molecular masses of
the regular comparison sample
integer*4 nr__1
parameter (nr__1=97) ! Number of records
character*99 ar__1 ! Full-size record
character*11 PNG_1 (nr__1) ! Galactic PN designation
character*14 Name_1 (nr__1) ! Common PN name
real*4 Dist_1 (nr__1) ! (kpc) Distance to the PN
real*4 E_Dist_2 (nr__1) ! (kpc) Positive (upper) error on distance
real*4 e_Dist_3 (nr__1) ! (kpc) Negative (lower) error on distance
real*4 Mion_1 (nr__1) ! (Msun) Ionised mass
real*4 E_Mion_2 (nr__1) ! (Msun) Positive (upper) error on ionised mass
real*4 e_Mion_3 (nr__1) ! (Msun) Negative (lower) error on ionised mass
character*1 l_Mmol_1 (nr__1) ! [<] Upper limit flag on Mmol
real*4 Mmol_1 (nr__1) ! (Msun) Molecular mass (or upper limit, see flag)
real*4 e_Mmol_1 (nr__1) ! (Msun) ? Error on molecular mass, only if
* Mmol is not an upper limit
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tableb1.dat' ! Parameters used to estimate the masses of the
whole sample
integer*4 nr__2
parameter (nr__2=118) ! Number of records
character*248 ar__2 ! Full-size record
character*14 Name_2 (nr__2) ! Common PN name
integer*4 n (nr__2) ! (cm-3) Electronic density
integer*4 E_n (nr__2) ! (cm-3) ? Positive (upper) error on density
integer*4 e_n_1 (nr__2) ! (cm-3) ? Negative (lower) error on density
real*4 Te (nr__2) ! (kK) Electronic temperature (1)
real*4 E_Te (nr__2) ! (kK) ? Positive (upper) error on temperature
real*4 e_Te_1 (nr__2) ! (kK) ? Negative (lower) error on temperature
real*4 Dist_2 (nr__2) ! (kpc) Distance to the PN
real*4 E_Dist_4 (nr__2) ! (kpc) Positive (upper) error on distance
real*4 e_Dist_5 (nr__2) ! (kpc) Negative (lower) error on distance
real*4 Diama (nr__2) ! (arcsec) PN diameter along major axis
real*4 Diamb (nr__2) ! (arcsec) PN diameter along minor axis
real*4 vexp (nr__2) ! (km/s) ? PN characteristic expansion velocity
real*4 logS_Ha (nr__2) ! ([mW/m2/sr]) Logarithm of H-alpha dereddened flux
real*4 e_logS_Ha (nr__2) ! ([mW/m2/sr]) Error on the logarithmic H-alpha flux
character*8 Tel (nr__2) ! Telescope used for CO observations
character*3 Trans (nr__2) ! CO rotational (J) transition used
real*4 rms (nr__2) ! (mK) ? Sensitivity on CO non-detection
real*4 Int (nr__2) ! (K.km/s) ? CO intensity on single-pointing
real*4 Flux (nr__2) ! (K.km/s.arcsec2) ? CO total flux
character*78 Ref (nr__2) ! References (in refs.dat file)
*Note (1): Te assumed to be 10kK when no data available in the literature
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'refs.dat' ! References
integer*4 nr__3
parameter (nr__3=74) ! Number of records
character*114 ar__3 ! Full-size record
integer*4 Ref_1 (nr__3) ! Reference code
character*19 BibCode (nr__3) ! BibCode
character*28 Aut (nr__3) ! Authors' name
character*63 Com (nr__3) ! Comments
C=============================================================================
C Loading file 'table4.dat' ! Computed ionised and molecular masses of
* post-common-envelope PNe
C Format for file interpretation
1 format(
+ A11,3X,A14,3X,F6.3,3X,F5.3,3X,F5.3,3X,F5.3,3X,F6.3,3X,F6.3,3X,
+ A1,3X,F5.3,3X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,21
read(1,'(A99)')ar__
read(ar__,1)
+ PNG(i__),Name(i__),Dist(i__),E_Dist(i__),e_Dist_1(i__),
+ Mion(i__),E_Mion(i__),e_Mion_1(i__),l_Mmol(i__),Mmol(i__),
+ e_Mmol(i__)
if(ar__(95:99) .EQ. '') e_Mmol(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ PNG(i__),Name(i__),Dist(i__),E_Dist(i__),e_Dist_1(i__),
+ Mion(i__),E_Mion(i__),e_Mion_1(i__),l_Mmol(i__),Mmol(i__),
+ e_Mmol(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea1.dat' ! Computed ionised and molecular masses of
* the regular comparison sample
C Format for file interpretation
2 format(
+ A11,3X,A14,3X,F6.3,3X,F5.3,3X,F5.3,3X,F5.3,3X,F6.3,3X,F6.3,3X,
+ A1,3X,F5.3,3X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'tablea1.dat')
write(6,*) '....Loading file: tablea1.dat'
do i__=1,97
read(1,'(A99)')ar__1
read(ar__1,2)
+ PNG_1(i__),Name_1(i__),Dist_1(i__),E_Dist_2(i__),
+ e_Dist_3(i__),Mion_1(i__),E_Mion_2(i__),e_Mion_3(i__),
+ l_Mmol_1(i__),Mmol_1(i__),e_Mmol_1(i__)
if(ar__1(95:99) .EQ. '') e_Mmol_1(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ PNG_1(i__),Name_1(i__),Dist_1(i__),E_Dist_2(i__),
+ e_Dist_3(i__),Mion_1(i__),E_Mion_2(i__),e_Mion_3(i__),
+ l_Mmol_1(i__),Mmol_1(i__),e_Mmol_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tableb1.dat' ! Parameters used to estimate the masses of the
* whole sample
C Format for file interpretation
3 format(
+ A14,3X,I5,3X,I5,3X,I5,3X,F5.2,3X,F4.2,3X,F4.2,3X,F6.3,3X,F5.3,
+ 3X,F5.3,3X,F5.1,3X,F5.1,3X,F5.1,3X,F5.2,3X,F4.2,3X,A8,3X,A3,
+ 3X,F5.1,3X,F5.1,3X,E7.1,3X,A78)
C Effective file loading
open(unit=1,status='old',file=
+'tableb1.dat')
write(6,*) '....Loading file: tableb1.dat'
do i__=1,118
read(1,'(A248)')ar__2
read(ar__2,3)
+ Name_2(i__),n(i__),E_n(i__),e_n_1(i__),Te(i__),E_Te(i__),
+ e_Te_1(i__),Dist_2(i__),E_Dist_4(i__),e_Dist_5(i__),
+ Diama(i__),Diamb(i__),vexp(i__),logS_Ha(i__),e_logS_Ha(i__),
+ Tel(i__),Trans(i__),rms(i__),Int(i__),Flux(i__),Ref(i__)
if(ar__2(26:30) .EQ. '') E_n(i__) = iNULL__
if(ar__2(34:38) .EQ. '') e_n_1(i__) = iNULL__
if(ar__2(50:53) .EQ. '') E_Te(i__) = rNULL__
if(ar__2(57:60) .EQ. '') e_Te_1(i__) = rNULL__
if(ar__2(105:109) .EQ. '') vexp(i__) = rNULL__
if(ar__2(145:149) .EQ. '') rms(i__) = rNULL__
if(ar__2(153:157) .EQ. '') Int(i__) = rNULL__
if(ar__2(161:167) .EQ. '') Flux(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ Name_2(i__),n(i__),E_n(i__),e_n_1(i__),Te(i__),E_Te(i__),
+ e_Te_1(i__),Dist_2(i__),E_Dist_4(i__),e_Dist_5(i__),
+ Diama(i__),Diamb(i__),vexp(i__),logS_Ha(i__),e_logS_Ha(i__),
+ Tel(i__),Trans(i__),rms(i__),Int(i__),Flux(i__),Ref(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'refs.dat' ! References
C Format for file interpretation
4 format(I2,1X,A19,1X,A28,A63)
C Effective file loading
open(unit=1,status='old',file=
+'refs.dat')
write(6,*) '....Loading file: refs.dat'
do i__=1,74
read(1,'(A114)')ar__3
read(ar__3,4)Ref_1(i__),BibCode(i__),Aut(i__),Com(i__)
c ..............Just test output...........
write(6,4)Ref_1(i__),BibCode(i__),Aut(i__),Com(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end