Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/459/342 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 2024-Apr-20
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/MNRAS/459/342 Dense cores in Taurus L1495 cloud (Marsh+, 2016)
*================================================================================
*A census of dense cores in the Taurus L1495 cloud from the
*Herschel Gould Belt Survey.
* Marsh K.A., Kirk J.M., Andre P., Griffin M.J., Konyves V., Palmeirim P.,
* Men'shchikov A., Ward-thompson D., Benedettini M., Bresnahan D.W.,
* Di Francesco J., Elia D., Motte F., Peretto N., Pezzuto S., Roy A.,
* Sadavoy S., Schneider N., Spinoglio L., White G.J.
* <Mon. Not. R. Astron. Soc., 459, 342-356 (2016)>
* =2016MNRAS.459..342M (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tableb1.dat' ! TAURUS-L1495 observed core catalog
integer*4 nr__
parameter (nr__=536) ! Number of records
character*467 ar__ ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg ! (deg) Right Ascension J2000
real*8 DEdeg ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 Seq ! [1/536] Source number (Source_number)
character*15 Name ! Source name (Source_name)
integer*4 RAh ! (h) Right ascension (J2000) (RA_2000_h)
integer*4 RAm ! (min) Right ascension (J2000) (RA_2000_m)
real*4 RAs ! (s) Right ascension (J2000) (RA_2000_s)
character*1 DE_ ! Declination sign (J2000)
integer*4 DEd ! (deg) Declination (J2000) (DEC_2000_d)
integer*4 DEm ! (arcmin) Declination (J2000) (DEC_2000_m)
real*4 DEs ! (arcsec) Declination (J2000) (DEC_2000_m)
real*4 sig70 ! Detection significance at 70um (Sig_70)
real*4 I70peak ! (Jy/beam) Estimate of the peak intensity at 70um
* ((I_70)peak)
real*4 e_I70peak ! (Jy/beam) Error in the estimate of the peak intensity
* (e_(I_70)peak)
real*4 C70 ! Contrast over the local background at 70um
* ((I_70)peak/Ibg)
real*4 I70conv ! (Jy/beam) Peak intensity at 70um after smoothing to a
* resolution corresponding to the 500um
* ((I_70)conv500)
real*4 S70 ! (Jy) Estimate of the total flux at 70um (S_70)
real*4 e_S70 ! (Jy) Error in estimate of the total flux at 70um
* (e_S_70)
integer*4 a70 ! (arcsec) Estimate of size along major axis at 70um
* (a_70)
integer*4 b70 ! (arcsec) Estimate of size along minor axis at 70um
* (b_70)
integer*4 PA70 ! (deg) Position angle of source major axis at 70um
* (E to N) (PA_70)
real*4 sig160 ! Detection significance at 160um (Sig_160)
real*4 I160peak ! (Jy/beam) Estimate of the peak intensity at 160um
* ((I_160)peak)
real*4 e_I160peak ! (Jy/beam) Error in the estimate of the peak intensity
* at 160um (e_(I_160)peak)
real*4 C160 ! Contrast of the peak intensity at 160um over
* the local background ((I_160)peak/Ibg)
real*4 I160conv ! (Jy/beam) Peak intensity at 160um after smoothing to
* a resolution corresponding to the 500um
* ((I_160)conv500)
real*4 S160 ! (Jy) Estimate of the total flux at 160um (S_160)
real*4 e_S160 ! (Jy) Error in estimate of the total flux at
* 160um (e_S_160)
integer*4 a160 ! (arcsec) Estimate of size along major axis at 160um
* (a_160)
integer*4 b160 ! (arcsec) Estimate of size along minor axis at 160um
* (b_160)
integer*4 PA160 ! (deg) Position angle of source major axis at 160um
* (E to N) (PA_160)
real*4 sig250 ! Detection significance at 250um (Sig_250)
real*8 I250peak ! (Jy/beam) Estimate of the peak intensity at 250um
* ((I_250)peak)
real*4 e_I250peak ! (Jy/beam) Error in the estimate of the peak intensity
* at 250um (e_(I_250)peak)
real*4 C250 ! Contrast over the local background at 250um
* ((I_250)peak/Ibg)
real*4 I250conv ! (Jy/beam) Peak intensity at 250um after smoothing to a
* resolution corresponding to the 500um
* ((I_250)conv500)
real*8 S250 ! (Jy) Estimate of the total flux at 250um (S_250)
real*4 e_S250 ! (Jy) Error in estimate of the total flux at
* 250um (e_S_250)
integer*4 a250 ! (arcsec) Estimate of size along major axis at 250um
* (a_250)
integer*4 b250 ! (arcsec) Estimate of size along minor axis at 250um
* (b_250)
integer*4 PA250 ! (deg) Position angle of source major axis at 250um
* (E to N) (PA_250)
real*4 sig350 ! Detection significance at 350um (Sig_350)
real*4 I350peak ! (Jy/beam) Estimate of the peak intensity at 350um
* ((I_350)peak)
real*4 e_I350peak ! (Jy/beam) Error in the estimate of the peak intensity
* at 350um (e_(I_350)peak)
real*4 C350 ! Contrast over the local background at 350um
* ((I_350)peak/Ibg)
real*4 I350conv ! (Jy/beam) Peak intensity at 350um after smoothing to
* a resolution corresponding to the 500um
* ((I_350)conv500)
real*8 S350 ! (Jy) Estimate of the total flux at 350um (S_350)
real*4 e_S350 ! (Jy) Error in estimate of the total flux at
* 350um (e_S_350)
integer*4 a350 ! (arcsec) Estimate of size along major axis at 350um
* (a_350)
integer*4 b350 ! (arcsec) Estimate of size along minor axis at 350um
* (b_350)
integer*4 PA350 ! (deg) Position angle of source major axis at 350um
* (E to N) (PA_350)
real*4 sig500 ! Detection significance at 500um (Sig_500)
real*8 I500peak ! (Jy/beam) Estimate of the peak intensity at 500um
* ((I_500)peak)
real*4 e_I500peak ! (Jy/beam) Error in the estimate of the peak intensity
* at 500um (e_(I_5000)peak)
real*4 C500 ! Contrast over the local background at 500um
* ((I_500)peak/Ibg)
real*8 S500 ! (Jy) Estimate of the total flux at 500um (S_500)
real*4 e_S500 ! (Jy) Error in estimate of the total flux at
* 500um (e_S_500)
integer*4 a500 ! (arcsec) Estimate of size along major axis at 500um
* (a_500)
integer*4 b500 ! (arcsec) Estimate of size along minor axis at 500um
* (b_500)
integer*4 PA500 ! (deg) Position angle of source major axis at 500um
* (E to N) (PA_500)
real*4 sigNH2 ! Detection significance from single scales on
* high-res column density map (Sig(N_H2))
real*4 NH2peak ! (10+21cm-2) Peak column density at 18.2" resolution,
* estimated by getsources ((N_H2)peak)
real*4 CNH2 ! Contrast of peak column density over the
* local background ((N_H2)peak/(N_H2)bg)
real*4 NH2conv ! (10+21cm-2) Peak column density after smoothing to a
* resolution corresponding to the 500um
* ((N_H2)conv500)
real*4 NH2bg ! (10+21cm-2) Column density of local background
* ((N_H2)bg)
integer*4 aNH2 ! (arcsec) Estimate of source size along major axis in
* column density map (a_(N_H2))
integer*4 bNH2 ! (arcsec) Estimate of source size along minor axis in
* column density map (b_(N_H2))
integer*4 PANH2 ! (deg) Position angle of source major axis in
* column density map (PA_(N_H2))
integer*4 Nsed ! Number of bands in which the source is
* significant (Sig_lambda > 5) (N_SED)
integer*4 Cflag ! [0/2] CSAR flag (CSAR_flag) (1)
integer*4 CType ! [1/4] Core type (Core_type) (G1)
character*24 Simbad ! Simbad ID of nearest match if within 1
* arcmin of Herschel position (SIMBAD_assoc)
character*37 Com ! Comments
*Note on (1): CSAR flag as follows:
* 2 = if the getsources core has a counterpart detected by CSAR
* (Kirk et al. 2013MNRAS.432.1424K) within 6 arcsec of its peak position
* 1 = if no close CSAR counterpart was found but the peak position of a CSAR
* source lies within the FWHM contour of the getsources core in the
* high-resolution column density map
* 0 = otherwise
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tableb2.dat' ! TAURUS-L1495 derived core catalog
integer*4 nr__1
parameter (nr__1=536) ! Number of records
character*156 ar__1 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_1 ! (deg) Right Ascension J2000
real*8 DEdeg_1 ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 Seq_1 ! [1/536] Source number (Source_number)
character*15 Name_1 ! Source name (Source_name)
integer*4 RAh_1 ! (h) Right ascension (J2000) (RA_2000_h)
integer*4 RAm_1 ! (min) Right ascension (J2000) (RA_2000_m)
real*4 RAs_1 ! (s) Right ascension (J2000) (RA_2000_s)
character*1 DE__1 ! Declination sign (J2000)
integer*4 DEd_1 ! (deg) Declination (J2000) (DEC_2000_d)
integer*4 DEm_1 ! (arcmin) Declination (J2000) (DEC_2000_m)
real*4 DEs_1 ! (arcsec) Declination (J2000) (DEC_2000_m)
real*4 Rad ! (pc) Core radius as observed (radius)
real*4 RadBeam ! Core radius with 18.2 arcsec beam
* (radius_beam)
real*4 Mcore ! (Msun) Core mass (Core_mass)
real*4 e_Mcore ! (Msun) Uncertainty in core mass (e_Core_mass)
real*4 Tdust ! (K) Dust temperature (Dust_T)
real*4 e_Tdust ! (K) Uncertainty in dust temperature (e_Dust_T)
real*4 NH2peak_1 ! (10+21cm-2) Peak H2 column density at 500um ((N_H2)peak)
real*4 v_NH2_o ! (10+21cm-2) Average column density observed
* (Ave_column_a)
real*4 v_NH2_d ! (10+21cm-2) Average column density deconvolved
* (Ave_column_b)
real*4 nH2peak_2 ! (10+4cm-3) Beam-averaged peak volume density
* ((n_H2)peak)
real*4 v_nH2_o_1 ! (10+4cm-3) Average volume density observed
* (Ave_volume_a)
real*4 v_nH2_d_1 ! (10+4cm-3) Average volume density deconvolved
* (Ave_volume_b)
real*4 Mratio ! Bonnor-Ebert mass ratio (Bonnor-Ebert)
integer*4 CType_1 ! [1/4] Core type (Core_type) (G1)
character*34 Com_1 ! Comments
C=============================================================================
C Loading file 'tableb1.dat' ! TAURUS-L1495 observed core catalog
C Format for file interpretation
1 format(
+ I3,1X,A15,1X,I1,1X,I2,1X,F5.2,A1,I2,1X,I2,1X,F5.2,1X,F6.1,1X,
+ E8.5,1X,E7.4,1X,F6.2,1X,E8.5,1X,E9.5,1X,E7.4,1X,I3,1X,I3,1X,
+ I3,1X,F5.1,1X,E8.5,1X,F6.4,1X,F5.2,1X,E8.5,1X,E8.5,1X,F6.4,1X,
+ I3,1X,I3,1X,I3,1X,F5.1,1X,F7.4,1X,E7.4,1X,F5.2,1X,F6.3,1X,
+ F7.4,1X,E7.4,1X,I3,1X,I3,1X,I3,1X,F5.1,1X,F6.4,1X,E7.4,1X,
+ F6.2,1X,F6.4,1X,F7.4,1X,E7.4,1X,I3,1X,I3,1X,I3,1X,F5.1,1X,
+ F7.4,1X,E7.4,1X,F5.2,1X,F8.5,1X,E7.4,1X,I3,1X,I3,1X,I3,1X,
+ F5.1,1X,F4.1,1X,F4.2,1X,F4.1,1X,F4.1,1X,I3,1X,I3,1X,I3,1X,I1,
+ 1X,I1,1X,I1,1X,A24,1X,A37)
C Effective file loading
open(unit=1,status='old',file=
+'tableb1.dat')
write(6,*) '....Loading file: tableb1.dat'
do i__=1,536
read(1,'(A467)')ar__
read(ar__,1)
+ Seq,Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,sig70,I70peak,e_I70peak,
+ C70,I70conv,S70,e_S70,a70,b70,PA70,sig160,I160peak,e_I160peak,
+ C160,I160conv,S160,e_S160,a160,b160,PA160,sig250,I250peak,
+ e_I250peak,C250,I250conv,S250,e_S250,a250,b250,PA250,sig350,
+ I350peak,e_I350peak,C350,I350conv,S350,e_S350,a350,b350,PA350,
+ sig500,I500peak,e_I500peak,C500,S500,e_S500,a500,b500,PA500,
+ sigNH2,NH2peak,CNH2,NH2conv,NH2bg,aNH2,bNH2,PANH2,Nsed,Cflag,
+ CType,Simbad,Com
RAdeg = rNULL__
DEdeg = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh .GT. -180) RAdeg=RAh*15.
if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
if(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
if(DEd .GE. 0) DEdeg=DEd
if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c ..............Just test output...........
write(6,1)
+ Seq,Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,sig70,I70peak,e_I70peak,
+ C70,I70conv,S70,e_S70,a70,b70,PA70,sig160,I160peak,e_I160peak,
+ C160,I160conv,S160,e_S160,a160,b160,PA160,sig250,I250peak,
+ e_I250peak,C250,I250conv,S250,e_S250,a250,b250,PA250,sig350,
+ I350peak,e_I350peak,C350,I350conv,S350,e_S350,a350,b350,PA350,
+ sig500,I500peak,e_I500peak,C500,S500,e_S500,a500,b500,PA500,
+ sigNH2,NH2peak,CNH2,NH2conv,NH2bg,aNH2,bNH2,PANH2,Nsed,Cflag,
+ CType,Simbad,Com
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tableb2.dat' ! TAURUS-L1495 derived core catalog
C Format for file interpretation
2 format(
+ I3,1X,A15,1X,I1,1X,I2,1X,F5.2,A1,I2,1X,I2,1X,F5.2,1X,F5.3,1X,
+ F5.3,1X,F5.3,1X,F5.3,1X,F4.1,1X,F4.1,1X,F4.1,1X,F5.2,1X,F5.2,
+ 1X,F5.1,1X,F5.2,1X,F6.2,1X,F6.1,1X,I1,1X,A34)
C Effective file loading
open(unit=1,status='old',file=
+'tableb2.dat')
write(6,*) '....Loading file: tableb2.dat'
do i__=1,536
read(1,'(A156)')ar__1
read(ar__1,2)
+ Seq_1,Name_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,Rad,
+ RadBeam,Mcore,e_Mcore,Tdust,e_Tdust,NH2peak_1,v_NH2_o,v_NH2_d,
+ nH2peak_2,v_nH2_o_1,v_nH2_d_1,Mratio,CType_1,Com_1
RAdeg_1 = rNULL__
DEdeg_1 = rNULL__
c Derive coordinates RAdeg_1 and DEdeg_1 from input data
c (RAdeg_1 and DEdeg_1 are set to rNULL__ when unknown)
if(RAh_1 .GT. -180) RAdeg_1=RAh_1*15.
if(RAm_1 .GT. -180) RAdeg_1=RAdeg_1+RAm_1/4.
if(RAs_1 .GT. -180) RAdeg_1=RAdeg_1+RAs_1/240.
if(DEd_1 .GE. 0) DEdeg_1=DEd_1
if(DEm_1 .GE. 0) DEdeg_1=DEdeg_1+DEm_1/60.
if(DEs_1 .GE. 0) DEdeg_1=DEdeg_1+DEs_1/3600.
if(DE__1.EQ.'-'.AND.DEdeg_1.GE.0) DEdeg_1=-DEdeg_1
c ..............Just test output...........
write(6,2)
+ Seq_1,Name_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,Rad,
+ RadBeam,Mcore,e_Mcore,Tdust,e_Tdust,NH2peak_1,v_NH2_o,v_NH2_d,
+ nH2peak_2,v_nH2_o_1,v_nH2_d_1,Mratio,CType_1,Com_1
write(6,'(6H Pos: 2F8.4)') RAdeg_1,DEdeg_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end