Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/704/A276 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-11
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/704/A276 Nearby Evolved Stars Survey. III. (Wallstroem+, 2025)
*================================================================================
*The Nearby Evolved Stars Survey.
*III: First data release of JCMT CO-line observations.
* Wallstroem S.H.J., Scicluna P., Srinivasan S., Wouterloot J.G.A.,
* McDonald I., Decock L., Wijshoff M., Chen R., Torres D., Umans L.,
* Willebrords B., Kemper F., Rau G., Feng S., Jeste M., Kaminski T., Li D.,
* Liu F.C., Trejo-Cruz A., Chawner H., Goldman S., MacIsaac H., Tang J.,
* Zeegers S.T., Danilovich T., Matsuura M., Menten K.M., van Loon J.T.,
* Cami J., Clark C.J.R., Dharmawardena T.E., Greaves J., He J., Imai H.,
* Jones O.C., Kim H., Marshall J.P., Shinnaga H., Wesson R.,
* the NESS Collaboration
* <Astron. Astrophys. 704, A276 (2025)>
* =2025A&A...704A.276W (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablec1.dat' ! *CO observations and estimated parameters
integer*4 nr__
parameter (nr__=485) ! Number of records
character*319 ar__ ! Full-size record
character*11 IRASPSC ! IRAS PSC identifier (IRASPSC)
character*24 SIMBADID ! SIMBAD identifier (SIMBAD_ID)
character*1 Chemtype ! [O/C] Chemical type based on mid-IR
* spectra [O for O-rich, C for C-rich];
* see Section 2.4 (Chem_type)
character*12 Tier ! Grouping from Scicluna et al.
* (2022MNRAS.512.1091S,) based on
* location in distance-DPR space
* (very low, low, intermediate, high,
* or extreme) (Tier)
real*8 PeakCO_2_1 ! (K) ?=-999 Peak intensity of the CO(2-1)
* line (Peak_CO(2-1))
real*8 e_PeakCO_2_1 ! (K) ?=-999 Uncertainty in the peak
* intensity of the CO(2-1) line
* (Peak_CO(2-1)_error)
real*8 vinfCO_2_1 ! (km/s) ?=-999 Expansion velocity of the
* shell derived from the CO(2-1) line
* (uncertainty 1km/s) (v_inf_CO(2-1))
real*8 IntCO_2_1 ! (K.km/s) ?=-999 Integrated intensity of the
* CO(2-1) line (Int_CO(2-1))
real*8 e_IntCO_2_1 ! (K.km/s) []?=-999 Uncertainty in the integrated
* intensity of the CO(2-1) line
* (Int_CO(2-1)_error)
integer*4 NchanCO_2_1 ! ?=-999 Number of velocity channels in
* the CO(2-1) spectral band
* (nchan_CO(2-1))
real*8 CO_2_1_rms ! (K) ?=-999 RMS noise in the CO(2-1) line
* (CO(2-1)_rms)
real*8 Int13CO_2_1 ! (K.km/s) ?=-999 Integrated intensity of the
* 13CO(2-1) line (Int_13CO(2-1))
real*8 e_Int13CO_2_1 ! (K.km/s) ?=-999 Uncertainty in the integrated
* intensity of the 13CO(2-1) line
* (Int_13CO(2-1)_error)
real*8 PeakCO_3_2 ! (K) ?=-999 Peak intensity of the CO(3-2)
* line (Peak_CO(3-2))
real*8 e_PeakCO_3_2 ! (K) ?=-999 Uncertainty in the peak
* intensity of the CO(3-2) line
* (Peak_CO(3-2)_error)
real*8 vinfCO_3_2 ! (km/s) ?=-999 Expansion velocity of the shell
* derived from the CO(3-2) line
* (uncertainty 1km/s) (v_inf_CO(3-2))
real*8 IntCO_3_2 ! (K.km/s) ?=-999 Integrated intensity of the
* CO(3-2) line (Int_CO(3-2))
real*8 e_IntCO_3_2 ! (K.km/s) []?=-999 Uncertainty in the integrated
* intensity of the CO(3-2) line
* (Int_CO(3-2)_error)
integer*4 NchanCO_3_2 ! ?=-999 Number of velocity channels in
* the CO(3-2) spectral band
* (nchan_CO(3-2))
real*8 CO_3_2_rms ! ?=-999 RMS noise in the CO(3-2) line
* (CO(3-2)_rms)
real*8 Int13CO_3_2 ! (K.km/s) ?=-999 Integrated intensity of the
* 13CO(3-2) line (Int_13CO(3-2))
real*8 e_Int13CO_3_2 ! (K.km/s) ?=-999 Uncertainty in the integrated
* intensity of the 13CO(3-2) line
* (Int_13CO(3-2)_error)
real*4 MLRCO_2_1 ! (Msun/yr) ?=- Empirical mass-loss rate derived
* from the CO(2-1) line using the
* Ramstedt et al. (2008A&A...487..645R)
* formula (MLR_CO(2-1))
real*4 e_MLRCO_2_1 ! (Msun/yr) ?=- Uncertainty in the mass-loss
* rate derived from the CO(2-1) line
* (MLR_CO(2-1)_error)
integer*4 f_CO_2_1 ! [0/1] Flag designating whether the
* CO(2-1) line is detected (0) or is
* below 3 times the RMS (1)
* (CO(2-1)_upperlimit)
real*4 MLRCO_3_2 ! (Msun/yr) ?=- Empirical mass-loss rate derived
* from the CO(3-2) line using the
* Ramstedt et al. (2008A&A...487..645R)
* formula (MLR_CO(3-2))
real*4 e_MLRCO_3_2 ! (Msun/yr) ?=- Uncertainty in the mass-loss rate
* derived from the CO(3-2) line
* (MLR_CO(3-2)_error)
integer*4 f_CO_3_2 ! [0/1] Flag denoting whether the CO(3-2) line
* is detected (0) or is below 3 times
* the RMS (1) (CO(3-2)_upperlimit)
real*4 IsoRatio_2_1 ! (Msun/yr) Isotopic ratio 12CO/13CO derived from
* the (2-1) lines (IsoRatio(2-1))
real*4 IsoRatio_3_2 ! (Msun/yr) Isotopic ratio 12CO/13CO derived from
* the (3-2) lines (IsoRatio(3-2))
real*4 DPR ! Dust-production rate from SED fit
* (DPR)
real*4 e_DPR ! Uncertainty in the dust-production
* rate (DPR_error)
real*8 GasDustCO_2_1 ! ?=- Gas-to-dust ratio derived from the
* CO(2-1) line (GasDust_CO(2-1))
real*8 e_GasDustCO_2_1 ! ?=- Uncertainty in the gas-to-dust
* ratio derived from the CO(2-1) line
* (GasDust_CO(2-1)_error)
real*8 GasDustCO_3_2 ! ?=- Gas-to-dust ratio derived from the
* CO(3-2) line (GasDust_CO(3-2))
real*8 e_GasDustCO_3_2 ! ?=- Uncertainty in the gas-to-dust
* ratio derived from the CO(3-2) line
* (GasDust_CO(3-2)_error)
integer*4 SoftParaCO_2_1 ! ?=-999[0/1] Flag denoting whether the
* CO(2-1) line shape is a soft parabola
* (1) or not (0) (SoftParabola_CO(2-1))
integer*4 SoftParaCO_3_2 ! ?=-999[0/1] Flag denoting whether the
* CO(3-2) line shape is a soft parabola
* (1) or not (0) (SoftParabola_CO(3-2))
integer*4 OptThin12CO_2_1 ! [0/1] Flag denoting whether the
* CO(2-1) line is optically thin (1)
* or not (0) (OptThin_12CO(2-1))
integer*4 OptThin12CO_3_2 ! [0/1] Flag denoting whether the
* CO(3-2) line is optically thin (1)
* or not (0) (OptThin_12CO(3-2))
integer*4 FirstDet ! [0/1] Flag denoting whether the NESS
* observation is the first CO detection
* for the source (1) or not (0)
* (First_det)
C=============================================================================
C Loading file 'tablec1.dat' ! *CO observations and estimated parameters
C Format for file interpretation
1 format(
+ A11,1X,A24,1X,A1,1X,A12,1X,F7.2,1X,F7.2,1X,F7.2,1X,F7.2,1X,
+ F7.2,1X,I4,1X,F7.2,1X,F7.2,1X,F7.2,1X,F7.2,1X,F7.2,1X,F7.2,1X,
+ F7.2,1X,F7.2,1X,I4,1X,F7.2,1X,F7.2,1X,F7.2,1X,E10.5,1X,E10.5,
+ 1X,I1,1X,E10.5,1X,E10.5,1X,I1,1X,F5.2,1X,F5.2,1X,E10.5,E10.5,
+ F7.2,1X,F7.2,1X,F8.2,1X,F8.2,1X,I4,1X,I4,1X,I1,1X,I1,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'tablec1.dat')
write(6,*) '....Loading file: tablec1.dat'
do i__=1,485
read(1,'(A319)')ar__
read(ar__,1)
+ IRASPSC,SIMBADID,Chemtype,Tier,PeakCO_2_1,e_PeakCO_2_1,
+ vinfCO_2_1,IntCO_2_1,e_IntCO_2_1,NchanCO_2_1,CO_2_1_rms,
+ Int13CO_2_1,e_Int13CO_2_1,PeakCO_3_2,e_PeakCO_3_2,vinfCO_3_2,
+ IntCO_3_2,e_IntCO_3_2,NchanCO_3_2,CO_3_2_rms,Int13CO_3_2,
+ e_Int13CO_3_2,MLRCO_2_1,e_MLRCO_2_1,f_CO_2_1,MLRCO_3_2,
+ e_MLRCO_3_2,f_CO_3_2,IsoRatio_2_1,IsoRatio_3_2,DPR,e_DPR,
+ GasDustCO_2_1,e_GasDustCO_2_1,GasDustCO_3_2,e_GasDustCO_3_2,
+ SoftParaCO_2_1,SoftParaCO_3_2,OptThin12CO_2_1,OptThin12CO_3_2,
+ FirstDet
if (idig(ar__(191:200)).EQ.0) MLRCO_2_1 = rNULL__
if (idig(ar__(202:211)).EQ.0) e_MLRCO_2_1 = rNULL__
if (idig(ar__(215:224)).EQ.0) MLRCO_3_2 = rNULL__
if (idig(ar__(226:235)).EQ.0) e_MLRCO_3_2 = rNULL__
if (idig(ar__(271:277)).EQ.0) GasDustCO_2_1 = rNULL__
if (idig(ar__(279:285)).EQ.0) e_GasDustCO_2_1 = rNULL__
if (idig(ar__(287:294)).EQ.0) GasDustCO_3_2 = rNULL__
if (idig(ar__(296:303)).EQ.0) e_GasDustCO_3_2 = rNULL__
c ..............Just test output...........
write(6,1)
+ IRASPSC,SIMBADID,Chemtype,Tier,PeakCO_2_1,e_PeakCO_2_1,
+ vinfCO_2_1,IntCO_2_1,e_IntCO_2_1,NchanCO_2_1,CO_2_1_rms,
+ Int13CO_2_1,e_Int13CO_2_1,PeakCO_3_2,e_PeakCO_3_2,vinfCO_3_2,
+ IntCO_3_2,e_IntCO_3_2,NchanCO_3_2,CO_3_2_rms,Int13CO_3_2,
+ e_Int13CO_3_2,MLRCO_2_1,e_MLRCO_2_1,f_CO_2_1,MLRCO_3_2,
+ e_MLRCO_3_2,f_CO_3_2,IsoRatio_2_1,IsoRatio_3_2,DPR,e_DPR,
+ GasDustCO_2_1,e_GasDustCO_2_1,GasDustCO_3_2,e_GasDustCO_3_2,
+ SoftParaCO_2_1,SoftParaCO_3_2,OptThin12CO_2_1,OptThin12CO_3_2,
+ FirstDet
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end
C Locate position of first digit in string; or return 0
integer function idig(c)
character*(*) c
character*1 c1
integer lc,i
lc=len(c)
idig=0
do i=1,lc
if(c(i:i).ne.' ') go to 1
end do
1 if(i.gt.lc) return
c1=c(i:i)
if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1
if(i.gt.lc) return
c1=c(i:i)
if(c1.ge.'0'.and.c1.le.'9') idig=i
return
end