Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/625/A19 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/625/A19 33 GHz continuum map of M51 (Querejeta+, 2019)
*================================================================================
*Dense gas is not enough: environmental variations in the star formation
*efficiency of dense molecular gas at 100pc scales in M51.
* Querejeta M., Schinnerer E., Schruba A., Murphy E., Meidt S., Usero A.,
* Leroy A. K., Pety J., Bigiel F., Chevance M., Faesi C. M., Gallagher M.,
* Garcia-Burillo S., Glover S. C. O., Hygate A. P. S., Jimenez-Donaire M. J.,
* Kruijssen J. M. D., Momjian E., Rosolowsky E., Utomo D.
* <Astron. Astrophys. 625, A19 (2019)>
* =2019A&A...625A..19Q (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea1.dat' ! Measurements for apertures on HCN peaks
integer*4 nr__
parameter (nr__=27) ! Number of records
character*122 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*8 Name ! Region in M51 (north, ring, south)
real*8 RAdeg ! (deg) Right ascension (J2000)
real*8 DEdeg ! (deg) Declination (J2000)
integer*4 SigStar ! (Msun/pc2) Stellar mass surface density
integer*4 e_SigStar ! (Msun/pc2) Uncertainty stellar mass surface density
real*4 TIR ! (10+6Lsun) Total Infrared Luminosity (TIR)
real*4 e_TIR ! (10+6Lsun) Uncertainty TIR
real*4 F33GHz ! (uJy) 33GHz Spectral Flux Density
real*4 e_F33GHz ! (uJy) Uncertainty 33GHz Flux Density
real*4 tfrac ! (%) Thermal fraction at 33GHz
real*4 e_tfrac ! (%) Uncertainty thermal fraction 33GHz
real*4 SigSFR ! (Msun/Myr/pc2) Star formation rate surface density
real*4 e_SigSFR ! (Msun/Myr/pc2) Uncertainty SFR surface density
real*4 IHCN ! (K.km/s) HCN(1-0) intensity
real*4 e_IHCN ! (K.km/s) Uncertainty HCN(1-0) intensity
real*4 ICO ! (K.km/s) ?=- CO(1-0) intensity
real*4 e_ICO ! (K.km/s) ?=- Uncertainty CO(1-0) intensity
real*4 sigHCN ! (km/s) ?=- HCN(1-0) velocity dispersion
real*4 e_sigHCN ! (km/s) ?=- Uncertainty HCN(1-0) velocity
* dispersion
real*4 sigCO ! (km/s) ?=- CO(1-0) velocity dispersion
real*4 e_sigCO ! (km/s) ?=- Uncertainty CO(1-0) velocity
* dispersion
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea2.dat' ! Measurements for Nyquist-sampled apertures
integer*4 nr__1
parameter (nr__1=4280) ! Number of records
character*122 ar__1 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*8 Name_1 ! Region in M51 (north, ring, south)
real*8 RAdeg_1 ! (deg) Right ascension (J2000)
real*8 DEdeg_1 ! (deg) Declination (J2000)
integer*4 SigStar_1 ! (Msun/pc2) Stellar mass surface density
integer*4 e_SigStar_1 ! (Msun/pc2) Uncertainty stellar mass surface density
real*4 TIR_1 ! (10+6Lsun) Total Infrared Luminosity (TIR)
real*4 e_TIR_1 ! (10+6Lsun) Uncertainty TIR
real*4 F33GHz_1 ! (uJy) 33GHz Spectral Flux Density
real*4 e_F33GHz_1 ! (uJy) Uncertainty 33GHz Flux Density
real*4 tfrac_1 ! (%) Thermal fraction at 33GHz
real*4 e_tfrac_1 ! (%) Uncertainty thermal fraction 33GHz
real*4 SigSFR_1 ! (Msun/Myr/pc2) Star formation rate surface density
real*4 e_SigSFR_1 ! (Msun/Myr/pc2) Uncertainty SFR surface density
real*4 IHCN_1 ! (K.km/s) HCN(1-0) intensity
real*4 e_IHCN_1 ! (K.km/s) Uncertainty HCN(1-0) intensity
real*4 ICO_1 ! (K.km/s) ?=- CO(1-0) intensity
real*4 e_ICO_1 ! (K.km/s) ?=- Uncertainty CO(1-0) intensity
real*4 sigHCN_1 ! (km/s) ?=- HCN(1-0) velocity dispersion
real*4 e_sigHCN_1 ! (km/s) ?=- Uncertainty HCN(1-0) velocity
* dispersion
real*4 sigCO_1 ! (km/s) ?=- CO(1-0) velocity dispersion
real*4 e_sigCO_1 ! (km/s) ?=- Uncertainty CO(1-0) velocity
* dispersion
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'list.dat' ! Information of fits datacube
integer*4 nr__2
parameter (nr__2=1) ! Number of records
character*88 ar__2 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
real*8 RAdeg_2 ! (deg) Right Ascension of center (J2000)
real*8 DEdeg_2 ! (deg) Declination of center (J2000)
integer*4 Nx ! Number of pixels along X-axis
integer*4 Ny ! Number of pixels along Y-axis
character*19 Obs_date ! ("datime") Observation date
integer*4 size ! (Kibyte) Size of FITS file
character*8 FileName ! Name of FITS file, in subdirectory fits
character*27 Title ! Title of the FITS file
C=============================================================================
C Loading file 'tablea1.dat' ! Measurements for apertures on HCN peaks
C Format for file interpretation
1 format(
+ A8,1X,F8.4,1X,F8.5,1X,I4,1X,I4,1X,F5.1,1X,F4.1,1X,F5.1,1X,
+ F3.1,1X,F6.1,1X,F4.1,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.1,
+ 1X,F4.1,1X,F4.1,1X,F3.1,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'tablea1.dat')
write(6,*) '....Loading file: tablea1.dat'
do i__=1,27
read(1,'(A122)')ar__
read(ar__,1)
+ Name,RAdeg,DEdeg,SigStar,e_SigStar,TIR,e_TIR,F33GHz,e_F33GHz,
+ tfrac,e_tfrac,SigSFR,e_SigSFR,IHCN,e_IHCN,ICO,e_ICO,sigHCN,
+ e_sigHCN,sigCO,e_sigCO
if (idig(ar__(93:97)).EQ.0) ICO = rNULL__
if (idig(ar__(99:102)).EQ.0) e_ICO = rNULL__
if (idig(ar__(104:107)).EQ.0) sigHCN = rNULL__
if (idig(ar__(109:111)).EQ.0) e_sigHCN = rNULL__
if (idig(ar__(113:117)).EQ.0) sigCO = rNULL__
if (idig(ar__(119:122)).EQ.0) e_sigCO = rNULL__
c ..............Just test output...........
write(6,1)
+ Name,RAdeg,DEdeg,SigStar,e_SigStar,TIR,e_TIR,F33GHz,e_F33GHz,
+ tfrac,e_tfrac,SigSFR,e_SigSFR,IHCN,e_IHCN,ICO,e_ICO,sigHCN,
+ e_sigHCN,sigCO,e_sigCO
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea2.dat' ! Measurements for Nyquist-sampled apertures
C Format for file interpretation
2 format(
+ A8,1X,F8.4,1X,F8.5,1X,I4,1X,I4,1X,F5.1,1X,F4.1,1X,F5.1,1X,
+ F3.1,1X,F6.1,1X,F4.1,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.1,
+ 1X,F4.1,1X,F4.1,1X,F3.1,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'tablea2.dat')
write(6,*) '....Loading file: tablea2.dat'
do i__=1,4280
read(1,'(A122)')ar__1
read(ar__1,2)
+ Name_1,RAdeg_1,DEdeg_1,SigStar_1,e_SigStar_1,TIR_1,e_TIR_1,
+ F33GHz_1,e_F33GHz_1,tfrac_1,e_tfrac_1,SigSFR_1,e_SigSFR_1,
+ IHCN_1,e_IHCN_1,ICO_1,e_ICO_1,sigHCN_1,e_sigHCN_1,sigCO_1,
+ e_sigCO_1
if (idig(ar__1(93:97)).EQ.0) ICO_1 = rNULL__
if (idig(ar__1(99:102)).EQ.0) e_ICO_1 = rNULL__
if (idig(ar__1(104:107)).EQ.0) sigHCN_1 = rNULL__
if (idig(ar__1(109:111)).EQ.0) e_sigHCN_1 = rNULL__
if (idig(ar__1(113:117)).EQ.0) sigCO_1 = rNULL__
if (idig(ar__1(119:122)).EQ.0) e_sigCO_1 = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1,RAdeg_1,DEdeg_1,SigStar_1,e_SigStar_1,TIR_1,e_TIR_1,
+ F33GHz_1,e_F33GHz_1,tfrac_1,e_tfrac_1,SigSFR_1,e_SigSFR_1,
+ IHCN_1,e_IHCN_1,ICO_1,e_ICO_1,sigHCN_1,e_sigHCN_1,sigCO_1,
+ e_sigCO_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'list.dat' ! Information of fits datacube
C Format for file interpretation
3 format(F9.5,F9.5,1X,I3,1X,I3,1X,A19,1X,I4,1X,A8,1X,A27)
C Effective file loading
open(unit=1,status='old',file=
+'list.dat')
write(6,*) '....Loading file: list.dat'
do i__=1,1
read(1,'(A88)')ar__2
read(ar__2,3)
+ RAdeg_2,DEdeg_2,Nx,Ny,Obs_date,size,FileName,Title
c ..............Just test output...........
write(6,3)
+ RAdeg_2,DEdeg_2,Nx,Ny,Obs_date,size,FileName,Title
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